Changeset 1320 for palm/trunk/SOURCE/advec_v_pw.f90
- Timestamp:
- Mar 20, 2014 8:40:49 AM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_v_pw.f90
r1310 r1320 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! ONLY-attribute added to USE-statements, 23 ! kind-parameters added to all INTEGER and REAL declaration statements, 24 ! kinds are defined in new module kinds, 25 ! revision history before 2012 removed, 26 ! comment fields (!:) to be used for variable explanations added to 27 ! all variable declaration statements 23 28 ! 24 29 ! Former revisions: … … 28 33 ! 1036 2012-10-22 13:43:42Z raasch 29 34 ! code put under GPL (PALM 3.9) 30 !31 ! 106 2007-08-16 14:30:26Z raasch32 ! j loop is starting from nysv (needed for non-cyclic boundary conditions)33 !34 ! 75 2007-03-22 09:54:05Z raasch35 ! vynp eliminated36 !37 ! RCS Log replace by Id keyword, revision history cleaned up38 !39 ! Revision 1.15 2006/02/23 09:46:08 raasch40 ! nzb_2d replaced by nzb_v_inner41 35 ! 42 36 ! Revision 1.1 1997/08/11 06:09:57 raasch … … 68 62 SUBROUTINE advec_v_pw 69 63 70 USE arrays_3d 71 USE control_parameters 72 USE grid_variables 73 USE indices 64 USE arrays_3d, & 65 ONLY: ddzw, tend, u, v, w 66 67 USE control_parameters, & 68 ONLY: u_gtrans, v_gtrans 69 70 USE grid_variables, & 71 ONLY: ddx, ddy 72 73 USE indices, & 74 ONLY: nxl, nxr, nyn, nysv, nzb_v_inner, nzt 75 76 USE kinds 77 74 78 75 79 IMPLICIT NONE 76 80 77 INTEGER :: i, j, k 78 REAL :: gu, gv 81 INTEGER(iwp) :: i !: 82 INTEGER(iwp) :: j !: 83 INTEGER(iwp) :: k !: 84 85 REAL(wp) :: gu !: 86 REAL(wp) :: gv !: 79 87 80 88 … … 105 113 SUBROUTINE advec_v_pw_ij( i, j ) 106 114 107 USE arrays_3d 108 USE control_parameters 109 USE grid_variables 110 USE indices 115 USE arrays_3d, & 116 ONLY: ddzw, tend, u, v, w 117 118 USE control_parameters, & 119 ONLY: u_gtrans, v_gtrans 120 121 USE grid_variables, & 122 ONLY: ddx, ddy 123 124 USE indices, & 125 ONLY: nzb_v_inner, nzt 126 127 USE kinds 128 111 129 112 130 IMPLICIT NONE 113 131 114 INTEGER :: i, j, k 115 REAL :: gu, gv 132 INTEGER(iwp) :: i !: 133 INTEGER(iwp) :: j !: 134 INTEGER(iwp) :: k !: 135 136 REAL(wp) :: gu !: 137 REAL(wp) :: gv !: 116 138 117 139 … … 119 141 gv = 2.0 * v_gtrans 120 142 DO k = nzb_v_inner(j,i)+1, nzt 121 tend(k,j,i) = tend(k,j,i) - 0.25 * ( &143 tend(k,j,i) = tend(k,j,i) - 0.25 * ( & 122 144 ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu ) & 123 145 - v(k,j,i-1) * ( u(k,j-1,i) + u(k,j,i) - gu ) ) * ddx &
Note: See TracChangeset
for help on using the changeset viewer.