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

Timestamp:

Jun 6, 2014 1:05:08 PM (11 years ago)

Author:

fricke

Message:

Bugfixes concerning grid stretching for the ocean and calculation of the salinity flux in routine surface_coupler

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: Change if-condition for stretched grid in the ocean, with the old
+!          condition and a negative value for dz_stretch_level the condition
+!          was always true for the whole model domain
 ! 
 ! Former revisions:
@@ -272 +274 @@
        dz_stretched = dz
        DO  k = nzt-1, 0, -1
-          IF ( dz_stretch_level <= ABS( zu(k+1) )  .AND.  &
-               dz_stretched < dz_max )  THEN
+!
+!--       The default value of dz_stretch_level is positive, thus the first 
+!--       condition is always true. Hence, the second condition is necessary.
+          IF ( dz_stretch_level >= zu(k+1)  .AND.  dz_stretch_level <= 0.0  &
+               .AND.  dz_stretched < dz_max )  THEN
              dz_stretched = dz_stretched * dz_stretch_factor
              dz_stretched = MIN( dz_stretched, dz_max )