Radar-derived microphysics temperature tendencies similar to operatio…

…nal HRRR (#986) This PR implements a feature of the operational HRRR, radar-derived microphysics temperature tendencies applied in the first N minutes of the forecast to improve clouds in the first few hours. Two regression tests are added.
ufs-community · Jan 11, 2022 · 1eef87e · 1eef87e
1 parent afa30f3
commit 1eef87e
Show file tree

Hide file tree

Showing 19 changed files with 4,714 additions and 3,396 deletions.
diff --git a/FV3 b/FV3
diff --git a/tests/RegressionTests_cheyenne.gnu.log b/tests/RegressionTests_cheyenne.gnu.log
diff --git a/tests/RegressionTests_cheyenne.intel.log b/tests/RegressionTests_cheyenne.intel.log
diff --git a/tests/RegressionTests_gaea.intel.log b/tests/RegressionTests_gaea.intel.log
diff --git a/tests/RegressionTests_hera.gnu.log b/tests/RegressionTests_hera.gnu.log
diff --git a/tests/RegressionTests_hera.intel.log b/tests/RegressionTests_hera.intel.log
diff --git a/tests/RegressionTests_jet.intel.log b/tests/RegressionTests_jet.intel.log
diff --git a/tests/RegressionTests_orion.intel.log b/tests/RegressionTests_orion.intel.log
diff --git a/tests/RegressionTests_wcoss_cray.log b/tests/RegressionTests_wcoss_cray.log
diff --git a/tests/RegressionTests_wcoss_dell_p3.log b/tests/RegressionTests_wcoss_dell_p3.log
diff --git a/tests/default_vars.sh b/tests/default_vars.sh
@@ -935,6 +935,7 @@ export LNDP_PRT_LIST=-999
 #IAU
 export IAU_INC_FILES="''"
 
+export FH_DFI_RADAR='-2e10'
 
 #Cellular automata
 export DO_CA=.false.

diff --git a/tests/fv3_conf/rrfs_warm_run.IN b/tests/fv3_conf/rrfs_warm_run.IN
@@ -0,0 +1,29 @@
+rm -fr INPUT RESTART
+mkdir INPUT RESTART
+
+OPNREQ_TEST=${OPNREQ_TEST:-false}
+SUFFIX=${RT_SUFFIX}
+cp -r @[INPUTDATA_ROOT]/FV3_input_data_conus13km/INPUT/* INPUT/
+
+for x in global_glacier.2x2.grb global_h2oprdlos.f77 global_maxice.2x2.grb \
+         global_o3prdlos.f77 global_snoclim.1.875.grb global_zorclim.1x1.grb \
+         CFSR.SEAICE.1982.2012.monthly.clim.grb RTGSST.1982.2012.monthly.clim.grb ; do
+  cp -f @[INPUTDATA_ROOT]/FV3_fix/"$x" .
+done
+
+cp    @[INPUTDATA_ROOT]/FV3_input_data/global_soilmgldas.t126.384.190.grb .
+cp    @[INPUTDATA_ROOT]/FV3_fix/seaice_newland.grb .
+
+touch data_table
+
+# Thompson MP lookup tables - copy standard and SIONlib tables of precomputed tables
+if [ $IMP_PHYSICS = 8 ]; then
+  cp    @[INPUTDATA_ROOT]/FV3_input_data_gsd/qr_acr_qsV2.dat .
+  cp    @[INPUTDATA_ROOT]/FV3_input_data_gsd/qr_acr_qgV2.dat .
+  cp    @[INPUTDATA_ROOT]/FV3_input_data_gsd/freezeH2O.dat .
+  cp    @[INPUTDATA_ROOT]/FV3_input_data_gsd/CCN_ACTIVATE.BIN .
+fi
+
+if [ $DO_RRTMGP = .true. ]; then
+  cp @[INPUTDATA_ROOT]/FV3_input_data_RRTMGP/* .
+fi
diff --git a/tests/parm/diag_table/diag_table_hrrr b/tests/parm/diag_table/diag_table_hrrr
diff --git a/tests/parm/model_configure_rrfs_conus13km.IN b/tests/parm/model_configure_rrfs_conus13km.IN
@@ -0,0 +1,105 @@
+total_member:            1
+start_year:              @[SYEAR]
+start_month:             @[SMONTH]
+start_day:               @[SDAY]
+start_hour:              @[SHOUR]
+start_minute:            0
+start_second:            0
+nhours_fcst:             @[FHMAX]
+
+dt_atmos:                @[DT_ATMOS]
+calendar:                'julian'
+restart_interval:        @[RESTART_INTERVAL]
+output_1st_tstep_rst:    .false.
+
+memuse_verbose:          .false.
+print_esmf:              .true.
+RUN_CONTINUE:            .false.
+ENS_SPS:                 .false.
+cpl:                     .false.
+
+quilting:                @[QUILTING]
+write_groups:            @[WRITE_GROUP]
+write_tasks_per_group:   @[WRTTASK_PER_GROUP]
+output_history:          @[OUTPUT_HISTORY]
+write_dopost:            @[WRITE_DOPOST]
+num_files:               @[NUM_FILES]
+filename_base:           @[FILENAME_BASE]
+output_file:             @[OUTPUT_FILE]
+ideflate:                @[IDEFLATE]
+nbits:                   @[NBITS]
+ichunk3d:                @[ICHUNK3D]
+jchunk3d:                @[JCHUNK3D]
+kchunk3d:                @[KCHUNK3D]
+imo:                     @[IMO]
+jmo:                     @[JMO]
+
+#
+# Write-component output frequency parameter definitions:
+#
+# nfhout:
+# Output frequency in hours after forecast hour "nfhmax_hf".
+#
+# nfhmax_hf:
+# Number of forecast hours until output frequency "nfhout" takes affect.
+#
+# nfhout_hf:
+# Output frequency in hours until forecast hour "nfhmax_hf".
+#
+# nsout:
+# Output frequency in time steps (positive values override "nfhout" and 
+# "nfhout_hf").
+#
+nfhout:     1
+nfhmax_hf:  60
+nfhout_hf:  1
+nsout:      -1
+
+#
+# Coordinate system used by the output grid.
+#
+output_grid:  @[OUTPUT_GRID]
+
+#
+# Parameter definitions for an output grid of type "lambert_conformal":
+#
+# cen_lon:
+# Longitude of center of grid (degrees).
+#
+# cen_lat:
+# Latitude of center of grid (degrees).
+#
+# stdlat1:
+# Latitude of first standard parallel (degrees).
+#
+# stdlat2:
+# Latitude of second standard parallel (degrees).
+#
+# nx:
+# Number of grid cells along x-axis in Lambert conformal (x,y) plane.
+#
+# ny:
+# Number of grid cells along y-axis in Lambert conformal (x,y) plane.
+#
+# lon1:
+# Longitude of center of grid cell at bottom-left corner of grid (degrees).
+#
+# lat1:
+# Latitude of center of grid cell at bottom-left corner of grid (degrees).
+#
+# dx:
+# Grid cell size in x direction (meters).
+#
+# dy:
+# Grid cell size in y direction (meters).
+#
+cen_lon:  -97.5
+cen_lat:  38.5
+stdlat1:  38.5
+stdlat2:  38.5
+nx:       393
+ny:       225
+lon1:     -121.70231097
+lat1:     22.57417972
+dx:       13000.0
+dy:       13000.0
+133 −7		ccpp/data/GFS_typedefs.F90
+66 −0		ccpp/data/GFS_typedefs.meta
+20 −1		ccpp/driver/GFS_diagnostics.F90
+42 −1		ccpp/driver/GFS_restart.F90
+1 −1		ccpp/physics