Changeset 57 for palm/trunk/SOURCE/diffusivities.f90

Timestamp:

Mar 9, 2007 12:05:41 PM (17 years ago)

Author:

raasch

Message:

preliminary update of further changes, advec_particles is not running!

File:

• palm/trunk/SOURCE/diffusivities.f90 (modified) (2 diffs)

palm/trunk/SOURCE/diffusivities.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Reference temperature pt_reference can be used in buoyancy term
 ! 
 ! Former revisions:
@@ -90 +90 @@
              dpt_dz = ( theta(k+1,j,i) - theta(k-1,j,i) ) * dd2zu(k)
              IF ( dpt_dz > 0.0 ) THEN
-                l_stable = 0.76 * sqrt_e(k) / &
-                                  SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
+                IF ( use_pt_reference )  THEN
+                   l_stable = 0.76 * sqrt_e(k) / &
+                                     SQRT( g / pt_reference * dpt_dz ) + 1E-5
+                ELSE
+                   l_stable = 0.76 * sqrt_e(k) / &
+                                     SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
+                ENDIF
              ELSE
                 l_stable = l_grid(k)