Changeset 94 for palm/trunk/SOURCE/init_grid.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/init_grid.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Grid definition for ocean version
 !
 ! Former revisions:
@@ -70 +70 @@
        CALL local_stop
     ENDIF
-!
-!-- Since the w-level lies on the surface, the first u-level (staggered!) lies
-!-- below the surface (used for "mirror" boundary condition).
-!-- The first u-level above the surface corresponds to the top of the 
-!-- Prandtl-layer.
-    zu(0) = - dz * 0.5
-    zu(1) =   dz * 0.5
-
-    dz_stretch_level_index = nzt+1
-    dz_stretched = dz
-    DO  k = 2, nzt+1
-       IF ( dz_stretch_level <= zu(k-1)  .AND.  dz_stretched < dz_max )  THEN
-          dz_stretched = dz_stretched * dz_stretch_factor
-          dz_stretched = MIN( dz_stretched, dz_max )
-          IF ( dz_stretch_level_index == nzt+1 )  dz_stretch_level_index = k-1
-       ENDIF
-       zu(k) = zu(k-1) + dz_stretched
-    ENDDO
-
-!
-!-- Compute the u-levels. They are always staggered half-way between the 
-!-- corresponding w-levels. The top w-level is extrapolated linearly.
-    zw(0) = 0.0
-    DO  k = 1, nzt
-       zw(k) = ( zu(k) + zu(k+1) ) * 0.5
-    ENDDO
-    zw(nzt+1) = zw(nzt) + 2.0 * ( zu(nzt+1) - zw(nzt) )
+
+!
+!-- Define the vertical grid levels
+    IF ( .NOT. ocean )  THEN
+!
+!--    Grid for atmosphere with surface at z=0 (k=0, w-grid).
+!--    Since the w-level lies on the surface, the first u-level (staggered!)
+!--    lies below the surface (used for "mirror" boundary condition).
+!--    The first u-level above the surface corresponds to the top of the 
+!--    Prandtl-layer.
+       zu(0) = - dz * 0.5
+       zu(1) =   dz * 0.5
+
+       dz_stretch_level_index = nzt+1
+       dz_stretched = dz
+       DO  k = 2, nzt+1
+          IF ( dz_stretch_level <= zu(k-1)  .AND.  dz_stretched < dz_max )  THEN
+             dz_stretched = dz_stretched * dz_stretch_factor
+             dz_stretched = MIN( dz_stretched, dz_max )
+             IF ( dz_stretch_level_index == nzt+1 ) dz_stretch_level_index = k-1
+          ENDIF
+          zu(k) = zu(k-1) + dz_stretched
+       ENDDO
+
+!
+!--    Compute the w-levels. They are always staggered half-way between the 
+!--    corresponding u-levels. The top w-level is extrapolated linearly.
+       zw(0) = 0.0
+       DO  k = 1, nzt
+          zw(k) = ( zu(k) + zu(k+1) ) * 0.5
+       ENDDO
+       zw(nzt+1) = zw(nzt) + 2.0 * ( zu(nzt+1) - zw(nzt) )
+
+    ELSE
+!
+!--    Grid for ocean with solid surface at z=0 (k=0, w-grid). The free water
+!--    surface is at k=nzt (w-grid).
+!--    Since the w-level lies always on the surface, the first/last u-level
+!--    (staggered!) lies below the bottom surface / above the free surface.
+!--    It is used for "mirror" boundary condition.
+!--    The first u-level above the bottom surface corresponds to the top of the 
+!--    Prandtl-layer.
+       zu(nzt+1) =   dz * 0.5
+       zu(nzt)   = - dz * 0.5
+
+       dz_stretch_level_index = 0
+       dz_stretched = dz
+       DO  k = nzt-1, 0, -1
+          IF ( dz_stretch_level <= ABS( zu(k+1) )  .AND.  &
+               dz_stretched < dz_max )  THEN
+             dz_stretched = dz_stretched * dz_stretch_factor
+             dz_stretched = MIN( dz_stretched, dz_max )
+             IF ( dz_stretch_level_index == 0 ) dz_stretch_level_index = k+1
+          ENDIF
+          zu(k) = zu(k+1) - dz_stretched
+       ENDDO
+
+!
+!--    Compute the w-levels. They are always staggered half-way between the 
+!--    corresponding u-levels.
+!--    The top w-level (nzt+1) is not used but set for consistency, since
+!--    w and all scalar variables are defined up tp nzt+1.
+       zw(nzt+1) = dz
+       zw(nzt)   = 0.0
+       DO  k = 0, nzt
+          zw(k) = ( zu(k) + zu(k+1) ) * 0.5
+       ENDDO
+
+    ENDIF
 
 !