Changeset 19 for palm/trunk/SOURCE/diffusion_e.f90
 Timestamp:
 Feb 23, 2007 4:53:48 AM (17 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/diffusion_e.f90
r4 r19 4 4 ! Actual revisions: 5 5 !  6 ! 6 ! Calculation extended for gridpoint nzt 7 7 ! 8 8 ! Former revisions: … … 56 56 REAL, DIMENSION(:,:), POINTER :: rif 57 57 REAL, DIMENSION(:,:,:), POINTER :: e, km, theta 58 REAL, DIMENSION(nzb+1:nzt 1,nys:nyn) :: dissipation, l, ll58 REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: dissipation, l, ll 59 59 60 60 … … 72 72 ENDIF 73 73 74 DO k = nzb_s_inner(j,i)+1, nzt 174 DO k = nzb_s_inner(j,i)+1, nzt 75 75 ! 76 76 ! Calculate the mixing length (for dissipation) … … 102 102 ! Calculate the tendency terms 103 103 DO j = nys, nyn 104 DO k = nzb_s_inner(j,i)+1, nzt 1104 DO k = nzb_s_inner(j,i)+1, nzt 105 105 106 106 dissipation(k,j) = ( 0.19 + 0.74 * l(k,j) / ll(k,j) ) * & … … 130 130 IF ( use_sgs_for_particles ) THEN 131 131 DO j = nys, nyn 132 DO k = nzb_s_inner(j,i)+1, nzt 1132 DO k = nzb_s_inner(j,i)+1, nzt 133 133 diss(k,j,i) = dissipation(k,j) 134 134 ENDDO … … 171 171 REAL, DIMENSION(:,:), POINTER :: rif 172 172 REAL, DIMENSION(:,:,:), POINTER :: e, km, theta 173 REAL, DIMENSION(nzb+1:nzt 1) :: dissipation, l, ll173 REAL, DIMENSION(nzb+1:nzt) :: dissipation, l, ll 174 174 175 175 … … 187 187 ! 188 188 ! Calculate the mixing length (for dissipation) 189 DO k = nzb_s_inner(j,i)+1, nzt 1189 DO k = nzb_s_inner(j,i)+1, nzt 190 190 dpt_dz = ( theta(k+1,j,i)  theta(k1,j,i) ) * dd2zu(k) 191 191 IF ( dpt_dz > 0.0 ) THEN … … 234 234 ! Store dissipation if needed for calculating the sgs particle velocities 235 235 IF ( use_sgs_for_particles ) THEN 236 DO k = nzb_s_inner(j,i)+1, nzt 1236 DO k = nzb_s_inner(j,i)+1, nzt 237 237 diss(k,j,i) = dissipation(k) 238 238 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.