Delete obsolete module_initialize_*.F files and modify z_scale (#526)

TYPE: enhancement

KEYWORDS: ideal initialization z_scale calculation added

SOURCE: internal

DESCRIPTION OF CHANGES:
The parameter z_scale was originally hardwired in each of the ideal initialization modules. This switch was used to stretch the vertical grid levels approximately with pressure when stretch_grid=.true. (internal switch).

New code now replaces the hardwired value with a calculation so that stretching stays reasonable when ztop is changed. The calculation for zscale:

z_scale = 8000/ztop 
where 8000 is a pressure scale height.

To document these changes, modifications are made to both run/README.namelist and test/em_real/examples.namelist.

This modification will alter results for some of the idealized cases (hill2d, grav2d, convrad) that use the original hardwired values associated with stretched_grid=.true.
Old z_scale values are left in and commented out.

LIST OF MODIFIED FILES:
D dyn_em/module_initialize_b_wave.F
D dyn_em/module_initialize_convrad.F
D dyn_em/module_initialize_grav2d_x.F
D dyn_em/module_initialize_hill2d_x.F
M dyn_em/module_initialize_ideal.F
D dyn_em/module_initialize_les.F
D dyn_em/module_initialize_quarter_ss.F
D dyn_em/module_initialize_seabreeze2d_x.F
D dyn_em/module_initialize_squall2d_x.F
D dyn_em/module_initialize_squall2d_y.F
M run/README.namelist
M test/em_real/examples.namelist
M test/em_seabreeze2d_x/README.seabreeze

TESTS CONDUCTED:
WTF OK
Ideal cases not run with new levels.

Author: dudhia

Diff for: dyn_em/module_initialize_b_wave.F

@@ -176,6 +176,8 @@ SUBROUTINE init_domain_rk ( grid &
 
    END SELECT
 
+   CALL model_to_grid_config_rec ( grid%id , model_config_rec , config_flags )
+
   ideal_constants: SELECT CASE ( model_config_rec%ideal_case )
   CASE ( hill2d_x )
    hm = 100.
@@ -193,8 +195,8 @@ SUBROUTINE init_domain_rk ( grid &
 
    stretch_grid = .true.
    delt = 0.
-!   z_scale = .50
-   z_scale = .40
+!   z_scale = .40
+   z_scale = 8000./config_flags%ztop
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/2
@@ -204,8 +206,8 @@ SUBROUTINE init_domain_rk ( grid &
 
    stretch_grid = .true.
    delt = 3.
-!   z_scale = .50
-   z_scale = .40
+!   z_scale = .40
+   z_scale = 8000./config_flags%ztop
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/2
@@ -214,8 +216,8 @@ SUBROUTINE init_domain_rk ( grid &
   CASE (grav2d_x)
 
    stretch_grid = .true.
-!   z_scale = .50
-   z_scale = 1.675
+!    z_scale = 1.675
+   z_scale = 8000./config_flags%ztop
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/2
@@ -224,8 +226,8 @@ SUBROUTINE init_domain_rk ( grid &
   CASE (convrad)
    stretch_grid = .true.
    delt = 1.
-!   z_scale = .50
-   z_scale = .32
+!   z_scale = .32
+   z_scale = 8000./config_flags%ztop
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/2
@@ -245,7 +247,8 @@ SUBROUTINE init_domain_rk ( grid &
 
    stretch_grid = .true.
    delt = 0.
-   z_scale = .50
+!   z_scale = .50
+   z_scale = 8000./config_flags%ztop
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/4
@@ -255,8 +258,8 @@ SUBROUTINE init_domain_rk ( grid &
 
    stretch_grid = .true.
    delt = 6.
-!   z_scale = .50
-   z_scale = .40
+!   z_scale = .40
+   z_scale = 8000./config_flags%ztop
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/2
@@ -270,8 +273,8 @@ SUBROUTINE init_domain_rk ( grid &
   !  FOR LES, set stretch to false
    stretch_grid = .false.
    delt = 3.
-!   z_scale = .50
-   z_scale = .40
+!   z_scale = .40
+   z_scale = 8000./config_flags%ztop  ! not used if stretch_grid = .false.
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
    nxc = (ide-ids)/2
@@ -287,7 +290,6 @@ SUBROUTINE init_domain_rk ( grid &
 
 
 
-   CALL model_to_grid_config_rec ( grid%id , model_config_rec , config_flags )
 
 ! here we check to see if the boundary conditions are set properly