Changeset 843 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Feb 29, 2012 3:16:21 PM (12 years ago)

Author:

gryschka

Message:

Redefiend first u-level (k=0) in ocean in case of bc_uv_b='dirichlet'

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (6 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! gryschka
+! In case of ocean and dirichlet bc for u and v at the bottom
+! the first u-level ist defined at same height as the first w-level
 !
 ! Former revisions:
@@ -135 +137 @@
 !
 !--    Grid for atmosphere with surface at z=0 (k=0, w-grid).
-!--    The first u-level above the surface corresponds to the top of the 
+!--    The second u-level (k=1) corresponds to the top of the 
 !--    Prandtl-layer.
 
@@ -159 +161 @@
 !
 !--    Compute the w-levels. They are always staggered half-way between the 
-!--    corresponding u-levels. The top w-level is extrapolated linearly.
+!--    corresponding u-levels. In case of dirichlet bc for u and v at the 
+!--    ground the first u- and w-level (k=0) are defined at same height (z=0). 
+!--    The top w-level is extrapolated linearly.
        zw(0) = 0.0
        DO  k = 1, nzt
@@ -168 +172 @@
     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.
+!--    Grid for ocean with free water surface is at k=nzt (w-grid).
+!--    In case of neumann bc at the ground the first first u-level (k=0) lies
+!--    below the first w-level (k=0). In case of dirichlet bc the first u- and
+!--    w-level are defined at same height, but staggered from the second level.
+!--    The second u-level (k=1) corresponds to the top of the Prandtl-layer.
        zu(nzt+1) =   dz * 0.5
        zu(nzt)   = - dz * 0.5
@@ -192 +194 @@
 !
 !--    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.
+!--    corresponding u-levels, except in case of dirichlet bc for u and v
+!--    at the ground. In this case the first u- and w-level are defined at 
+!--    same height. 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
@@ -200 +203 @@
           zw(k) = ( zu(k) + zu(k+1) ) * 0.5
        ENDDO
+
+!
+!--    In case of dirichlet bc for u and v the first u- and w-level are defined
+!--    at same height.
+       IF ( ibc_uv_b == 0 ) THEN
+          zu(0) = zw(0)
+       ENDIF
 
     ENDIF