First attempt at adding the massden parameter.

DavidHuber-NOAA · Sep 5, 2023 · c95800b · c95800b
1 parent 4656e9f
commit c95800b
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 5 deletions.
diff --git a/parm/varmap_tables/GFSphys_var_map.txt b/parm/varmap_tables/GFSphys_var_map.txt
@@ -21,3 +21,4 @@ cnwat               cnwat     set_to_fill   0.0      S
 hice                icetk     set_to_fill   1.5      S
 weasd               weasd     set_to_fill   0.0      S
 snod                snod      set_to_fill   0.0      S       	
+massden             smoke     set_to_fill   1E-12    T
diff --git a/sorc/chgres_cube.fd/atm_input_data.F90 b/sorc/chgres_cube.fd/atm_input_data.F90
@@ -1978,7 +1978,7 @@ subroutine read_input_atm_grib2_file(localpet)
 
  integer, intent(in)                   :: localpet
 
- integer, parameter                    :: ntrac_max=14
+ integer, parameter                    :: ntrac_max=15
  integer, parameter                    :: max_levs=1000
 
  character(len=300)                    :: the_file
@@ -2030,18 +2030,18 @@ subroutine read_input_atm_grib2_file(localpet)
 
  tracers(:) = "NULL"
 
- trac_names_oct10 = (/1,  1,  14,  1,  1,  1,  1, 6,  6,   1,  6,  13,  13, 2 /)
- trac_names_oct11 = (/0, 22, 192, 23, 24, 25, 32, 1, 29, 100, 28, 193, 192, 2 /)
+ trac_names_oct10 = (/1,  1,  14,  1,  1,  1,  1, 6,  6,   1,  6,  13,  13, 2, 20 /)
+ trac_names_oct11 = (/0, 22, 192, 23, 24, 25, 32, 1, 29, 100, 28, 193, 192, 2, 0 /)
 
  trac_names_vmap = (/"sphum   ", "liq_wat ", "o3mr    ", "ice_wat ", &
                      "rainwat ", "snowwat ", "graupel ", "cld_amt ", "ice_nc  ", &
                      "rain_nc ", "water_nc", "liq_aero", "ice_aero", &
-                     "sgs_tke "/)
+                     "sgs_tke ", "massden "/)
 
  tracers_default = (/"sphum   ", "liq_wat ", "o3mr    ", "ice_wat ", &
                      "rainwat ", "snowwat ", "graupel ", "cld_amt ", "ice_nc  ", &
                      "rain_nc ", "water_nc", "liq_aero", "ice_aero", &
-                     "sgs_tke "/)
+                     "sgs_tke ", "massden "/)
 
  the_file = trim(data_dir_input_grid) // "/" // trim(grib2_file_input_grid)