Changeset 94 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Jun 1, 2007 3:25:22 PM (17 years ago)

Author:

raasch

Message:

preliminary uncomplete changes for ocean version

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Calculation of Deardorff velocity scale adjusted to be used with the ocean
+! version (HAS STILL TO BE COMPLETED!!!)
 !
 ! Former revisions:
@@ -784 +785 @@
 !
 !--    Determine the boundary layer height using two different schemes.
-!--    First scheme: Starting from the Earth's surface, look for the first 
-!--    relative minimum of the total heat flux. The corresponding height is
-!--    accepted as the boundary layer height, if it is less than 1.5 times the
-!--    height where the heat flux becomes negative for the first time.
+!--    First scheme: Starting from the Earth's (Ocean's) surface, look for the
+!--    first relative minimum (maximum) of the total heat flux.
+!--    The corresponding height is assumed as the boundary layer height, if it
+!--    is less than 1.5 times the height where the heat flux becomes negative
+!--    (positive) for the first time.
 !--    NOTE: This criterion is still capable of improving!
        z_i(1) = 0.0
        first = .TRUE.
-       DO  k = nzb, nzt-1
-          IF ( first .AND. hom(k,1,18,sr) < 0.0 )  THEN
-             first = .FALSE.
-             height = zw(k)
-          ENDIF
-          IF ( hom(k,1,18,sr) < 0.0  .AND. &
-               hom(k+1,1,18,sr) > hom(k,1,18,sr) )  THEN
-             IF ( zw(k) < 1.5 * height )  THEN
-                z_i(1) = zw(k)
-             ELSE
-                z_i(1) = height
+!       IF ( .NOT. ocean )  THEN
+          DO  k = nzb, nzt-1
+             IF ( first .AND. hom(k,1,18,sr) < 0.0 )  THEN
+                first = .FALSE.
+                height = zw(k)
              ENDIF
-             EXIT
-          ENDIF
-       ENDDO
+             IF ( hom(k,1,18,sr) < 0.0  .AND. &
+                  hom(k+1,1,18,sr) > hom(k,1,18,sr) )  THEN
+                IF ( zw(k) < 1.5 * height )  THEN
+                   z_i(1) = zw(k)
+                ELSE
+                   z_i(1) = height
+                ENDIF
+                EXIT
+             ENDIF
+          ENDDO
+!       ELSE
+!       ENDIF
 
 !
@@ -829 +834 @@
        IF ( hom(nzb,1,18,sr) > 0.0  .AND.  z_i(1) /= 0.0 )  THEN
           hom(nzb+8,1,pr_palm,sr)  = ( g / hom(nzb+1,1,4,sr) * &
-                                       hom(nzb,1,18,sr) * z_i(1) )**0.333333333
+                                       hom(nzb,1,18,sr) *      &
+                                       ABS( z_i(1) ) )**0.333333333
 !--       so far this only works if Prandtl layer is used
           hom(nzb+11,1,pr_palm,sr) = hom(nzb,1,16,sr) / hom(nzb+8,1,pr_palm,sr)