Changeset 57 for palm/trunk/SOURCE/diffusion_e.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/diffusion_e.f90 (modified) (3 diffs)

palm/trunk/SOURCE/diffusion_e.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Reference temperature pt_reference can be used in buoyancy term
 !
 ! Former revisions:
@@ -82 +82 @@
                 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,j,i) ) / &
-                                     SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
+                   IF ( use_pt_reference )  THEN
+                      l_stable = 0.76 * SQRT( e(k,j,i) ) / &
+                                        SQRT( g / pt_reference * dpt_dz ) + 1E-5
+                   ELSE
+                      l_stable = 0.76 * SQRT( e(k,j,i) ) / &
+                                        SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
+                   ENDIF
                 ELSE
                    l_stable = l_grid(k)
@@ -195 +200 @@
           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,j,i) ) / &
-                               SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
+             IF ( use_pt_reference )  THEN
+                l_stable = 0.76 * SQRT( e(k,j,i) ) / &
+                                  SQRT( g / pt_reference * dpt_dz ) + 1E-5
+             ELSE
+                l_stable = 0.76 * SQRT( e(k,j,i) ) / &
+                                  SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
+             ENDIF
           ELSE
              l_stable = l_grid(k)