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

palm/trunk/SOURCE/modules.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Variables concerning stretching introduced or revised
+! 
+! 3045 2018-05-28 07:55:41Z Giersch
 ! z_max_do2d removed
 ! 
@@ -1148 +1151 @@
     INTEGER(iwp) ::  dp_level_ind_b = 0                !< lowest grid index for external pressure gradient forcing
     INTEGER(iwp) ::  dvrp_filecount = 0                !< parameter for dvr visualization software
-    INTEGER(iwp) ::  dz_stretch_level_index            !< vertical grid level index above which the vertical grid spacing is stretched
     INTEGER(iwp) ::  ensemble_member_nr = 0            !< namelist parameter
     INTEGER(iwp) ::  gamma_mg                          !< switch for steering the multigrid cycle: 1: v-cycle, 2: w-cycle
@@ -1190 +1192 @@
     INTEGER(iwp) ::  num_var_fl                        !< number of sampling/output variables in virtual flight measurements
     INTEGER(iwp) ::  num_var_fl_user=0                 !< number of user-defined sampling/output variables in virtual flight measurements
+    INTEGER(iwp) ::  number_stretch_level_start        !< number of user-specified start levels for stretching
+    INTEGER(iwp) ::  number_stretch_level_end          !< number of user-specified end levels for stretching
     INTEGER(iwp) ::  nz_do3d = -9999                   !< namelist parameter
     INTEGER(iwp) ::  prt_time_count = 0                !< number of output intervals for particle data output
@@ -1199 +1203 @@
     INTEGER(iwp) ::  timestep_count = 0                !< number of timesteps carried out since the beginning of the initial run
     INTEGER(iwp) ::  y_shift = 0                       !< namelist parameter
+    
     INTEGER(iwp) ::  dist_nxl(0:1)                               !< left boundary of disturbance region
     INTEGER(iwp) ::  dist_nxr(0:1)                               !< right boundary of disturbance region
@@ -1211 +1216 @@
     INTEGER(iwp) ::  domask_no(max_masks,0:1) = 0                !< number of masked output quantities
     INTEGER(iwp) ::  domask_time_count(max_masks,0:1)            !< number of output intervals for masked data
+    INTEGER(iwp) ::  dz_stretch_level_end_index(9)               !< vertical grid level index until which the vertical grid spacing is stretched
+    INTEGER(iwp) ::  dz_stretch_level_start_index(9)             !< vertical grid level index above which the vertical grid spacing is stretched
     INTEGER(iwp) ::  mask_size(max_masks,3) = -1                 !< size of mask array per mask and dimension (for netcdf output)
     INTEGER(iwp) ::  mask_size_l(max_masks,3) = -1               !< subdomain size of mask array per mask and dimension (for netcdf output)
@@ -1416 +1423 @@
     REAL(wp) ::  dt_spinup = 60.0_wp                           !< namelist parameter
     REAL(wp) ::  dt_3d = 1.0_wp                                !< time step
-    REAL(wp) ::  dz = -1.0_wp                                  !< namelist parameter
-    REAL(wp) ::  dz_max = 9999999.9_wp                         !< namelist parameter
+    REAL(wp) ::  dz_max = 1000.0_wp                            !< namelist parameter
     REAL(wp) ::  dz_stretch_factor = 1.08_wp                   !< namelist parameter
-    REAL(wp) ::  dz_stretch_level = 100000.0_wp                !< namelist parameter
+    REAL(wp) ::  dz_stretch_level = -9999999.9_wp              !< namelist parameter
     REAL(wp) ::  e_init = 0.0_wp                               !< namelist parameter
     REAL(wp) ::  e_min = 0.0_wp                                !< namelist parameter
@@ -1529 +1535 @@
     REAL(wp) ::  dpdxy(1:2) = 0.0_wp                               !< namelist parameter
     REAL(wp) ::  dt_domask(max_masks) = 9999999.9_wp               !< namelist parameter
+    REAL(wp) ::  dz(10) = -1.0_wp                                  !< namelist parameter
+    REAL(wp) ::  dz_stretch_level_start(9) = -9999999.9_wp         !< namelist parameter
+    REAL(wp) ::  dz_stretch_level_end(9) = 9999999.9_wp            !< namelist parameter
+    REAL(wp) ::  dz_stretch_factor_array(9) = 1.08_wp              !< namelist parameter
     REAL(wp) ::  mask_scale(3)                                     !< collective array for mask_scale_x/y/z
     REAL(wp) ::  pt_vertical_gradient(10) = 0.0_wp                 !< namelist parameter