Changeset 3065 for palm/trunk/SOURCE/init_masks.f90

Timestamp:

Jun 12, 2018 7:03:02 AM (6 years ago)

Author:

Giersch

Message:

New vertical stretching procedure has been introduced

File:

• palm/trunk/SOURCE/init_masks.f90 (modified) (5 diffs)

palm/trunk/SOURCE/init_masks.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! dz_stretch_level was replaced by dz_stretch_level_start
+! 
+! 3049 2018-05-29 13:52:36Z Giersch
 ! Error messages revised
 ! 
@@ -133 +136 @@
         ONLY:  constant_diffusion, cloud_droplets, cloud_physics,              &
                data_output_masks, data_output_masks_user,                      &
-               doav, doav_n, domask, domask_no, dz, dz_stretch_level, humidity,&
-               mask, masks, mask_scale, mask_i,                                &
+               doav, doav_n, domask, domask_no, dz, dz_stretch_level_start,    &
+               humidity, mask, masks, mask_scale, mask_i,                      &
                mask_i_global, mask_j, mask_j_global, mask_k, mask_k_global,    &
                mask_loop, mask_size, mask_size_l, mask_start_l, mask_x,        &
@@ -141 +144 @@
                microphysics_morrison, microphysics_seifert, passive_scalar,    &
                ocean, varnamelength
+               
 
     USE grid_variables,                                                        &
@@ -490 +494 @@
        CALL set_mask_locations( 1, dx, 'dx', nx, 'nx', nxl, nxr )
        CALL set_mask_locations( 2, dy, 'dy', ny, 'ny', nys, nyn )
-       CALL set_mask_locations( 3, dz, 'dz', nz, 'nz', nzb, nzt )
+       CALL set_mask_locations( 3, dz(1), 'dz', nz, 'nz', nzb, nzt )
 !
 !--    Set global masks along all three dimensions (required by
@@ -727 +731 @@
 !--          The following line assumes a constant vertical grid spacing within
 !--          the vertical mask range; it fails for vertical grid stretching.
-!--          Maybe revise later. Issue warning but continue execution.
+!--          Maybe revise later. Issue warning but continue execution. ABS(...)
+!--          within the IF statement is necessary because the default value of
+!--          dz_stretch_level_start is -9999999.9_wp.
              loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim) * ddxyz )
 
-             IF ( mask_loop(mid,dim,2) * mask_scale(dim) > dz_stretch_level )  &
-                  THEN
+             IF ( mask_loop(mid,dim,2) * mask_scale(dim) >                     &
+                  ABS( dz_stretch_level_start(1) ) )  THEN
                 WRITE ( message_string, '(A,I3,A,I1,A,F9.3,A,F8.2,3A)' )       &
                      'mask_loop(',mid,',',dim,',2)=', mask_loop(mid,dim,2),    &
-                     ' exceeds dz_stretch_level=',dz_stretch_level,            &
+                     ' exceeds dz_stretch_level=',dz_stretch_level_start(1),   &
                      '.&Vertical mask locations will not ',                    &
                      'match the desired heights within the stretching ',       &