Changeset 57 for palm/trunk/SOURCE/diffusion_e.f90
 Timestamp:
 Mar 9, 2007 12:05:41 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/diffusion_e.f90
r39 r57 4 4 ! Actual revisions: 5 5 !  6 ! 6 ! Reference temperature pt_reference can be used in buoyancy term 7 7 ! 8 8 ! Former revisions: … … 82 82 dpt_dz = ( theta(k+1,j,i)  theta(k1,j,i) ) * dd2zu(k) 83 83 IF ( dpt_dz > 0.0 ) THEN 84 l_stable = 0.76 * SQRT( e(k,j,i) ) / & 85 SQRT( g / theta(k,j,i) * dpt_dz ) + 1E5 84 IF ( use_pt_reference ) THEN 85 l_stable = 0.76 * SQRT( e(k,j,i) ) / & 86 SQRT( g / pt_reference * dpt_dz ) + 1E5 87 ELSE 88 l_stable = 0.76 * SQRT( e(k,j,i) ) / & 89 SQRT( g / theta(k,j,i) * dpt_dz ) + 1E5 90 ENDIF 86 91 ELSE 87 92 l_stable = l_grid(k) … … 195 200 dpt_dz = ( theta(k+1,j,i)  theta(k1,j,i) ) * dd2zu(k) 196 201 IF ( dpt_dz > 0.0 ) THEN 197 l_stable = 0.76 * SQRT( e(k,j,i) ) / & 198 SQRT( g / theta(k,j,i) * dpt_dz ) + 1E5 202 IF ( use_pt_reference ) THEN 203 l_stable = 0.76 * SQRT( e(k,j,i) ) / & 204 SQRT( g / pt_reference * dpt_dz ) + 1E5 205 ELSE 206 l_stable = 0.76 * SQRT( e(k,j,i) ) / & 207 SQRT( g / theta(k,j,i) * dpt_dz ) + 1E5 208 ENDIF 199 209 ELSE 200 210 l_stable = l_grid(k)
Note: See TracChangeset
for help on using the changeset viewer.