Merge HWRF CCPP physics into HAFS (#57)

This Pull Request enables using the HWRF CCPP physics in the HAFS system.

*Initial workflow level changes came from the feature/hwrf_physics branch of @mzhangw
*Import hwrf_physics_fix: DETAMPNEW_DATA.expanded_rain_LE is needed for the F-A microphysics scheme and the following files are needed for the HWRF version of NOAH LSM (note: other versions of NOAH LSM may need different versions of these files): GENPARM.TBL, SOILPARM.TBL, VEGPARM.TBL
*Add two hwrf physics suite related workflow level regression tests
*Workflow level (parm/script/build) clean ups and unification to support both HWRF ccpp suite and the existing HAFS suite.

All workflow level regression test passed. And this PR addresses issue #58.

Note: Enabling using the HWRF CCPP Suite in HAFS is a collaborative group effort among DTC, GSL, EMC, including: Man Zhang, Grant Firl, Mrinal Biswas, Dominikus Heinzeller, Chunxi Zhang, Eric Aligo, Weiguo Wang, Qingfu Liu, Bin Liu, etc.

parm/forecast/globnest/diag_table.tmp

@@ -303,7 +303,7 @@
 "gfs_phys",    "nddsf_ave",    "nddsf_ave",     "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_phys",    "trans_ave",    "trans_ave",     "fv3_history2d",  "all",  .false.,  "none",  2
 
-# Reflectivity from GFDL microphysics
+# Reflectivity from microphysics
 "gfs_phys", "refl_10cm", "refl_10cm", "fv3_history2d", "all", .false., "none", 2
 
 #Max/Min must be kept in separate files.  Time is controlled by diag_table

parm/forecast/globnest/input.nml.tmp

@@ -177,7 +177,8 @@
        redrag         = .true.
        dspheat        = .true.
        hybedmf        = .true.
-       moninq_fac     = -1.0
+       hurr_pbl       = .T.        ! HWRF moninedmf
+       moninq_fac     = -1.0       ! HWRF_moninedmf
        satmedmf       = .false.
        random_clds    = .false.
        trans_trac     = .true.

parm/forecast/regional_hafsv0p1a/diag_table.tmp → parm/forecast/globnest_hwrf/diag_table.tmp

@@ -119,9 +119,9 @@
 "gfs_dyn",     "rainwat",     "rwmr",      "fv3_history",    "all",  .false.,  "none",  2
 "gfs_dyn",     "snowwat",     "snmr",      "fv3_history",    "all",  .false.,  "none",  2
 "gfs_dyn",     "graupel",     "grle",      "fv3_history",    "all",  .false.,  "none",  2
+"gfs_dyn",     "q_rimef",     "qrim",      "fv3_history",    "all",  .false.,  "none",  2
 "gfs_dyn",     "ps",          "pressfc",   "fv3_history",    "all",  .false.,  "none",  2
 "gfs_dyn",     "hs",          "hgtsfc",    "fv3_history",    "all",  .false.,  "none",  2
-"gfs_dyn",     "cld_amt",     "cld_amt",   "fv3_history",    "all",  .false.,  "none",  2
 #"gfs_dyn",     "ice_nc",      "nicp",      "fv3_history",  "all",  .false.,  "none",  2
 #"gfs_dyn",     "rain_nc",     "ntrnc",     "fv3_history",  "all",  .false.,  "none",  2
 
@@ -201,6 +201,9 @@
 "gfs_phys",    "u10m",        "ugrd10m",   "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_phys",    "v10m",        "vgrd10m",   "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_phys",    "crain",       "crain",     "fv3_history2d",  "all",  .false.,  "none",  2
+#"gfs_phys",  "f_ice",         "f_ice",    "fv3_history2d",         "all",  .false.,  "none",  2
+#"gfs_phys",  "f_rain",        "f_rain",   "fv3_history2d",         "all",  .false.,  "none",  2
+#"gfs_phys",  "f_rimef",       "f_rimef",  "fv3_history2d",         "all",  .false.,  "none",  2
 "gfs_sfc",     "tprcp",       "tprcp",     "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_sfc",     "hgtsfc",      "orog",      "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_sfc",     "weasd",       "weasd",     "fv3_history2d",  "all",  .false.,  "none",  2
@@ -302,8 +305,12 @@
 "gfs_phys",    "nbdsf_ave",    "nbdsf_ave",     "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_phys",    "nddsf_ave",    "nddsf_ave",     "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_phys",    "trans_ave",    "trans_ave",     "fv3_history2d",  "all",  .false.,  "none",  2
-
-# Reflectivity from GFDL microphysics
+#HWRF
+#"gfs_phys",    "cleffr",       "cleffr",        "fv3_history2d",  "all",  .false.,  "none",  2 
+#"gfs_phys",    "cieffr",       "cieffr",        "fv3_history2d",  "all",  .false.,  "none",  2
+#"gfs_phys",    "cseffr",       "cseffr",        "fv3_history2d",  "all",  .false.,  "none",  2
+#
+# Reflectivity from microphysics
 "gfs_phys", "refl_10cm", "refl_10cm", "fv3_history2d", "all", .false., "none", 2
 
 #Max/Min must be kept in separate files.  Time is controlled by diag_table

parm/forecast/regional_hafsv0p1a/field_table → parm/forecast/globnest_hwrf/field_table

@@ -9,20 +9,18 @@
            "longname",     "cloud water mixing ratio"
            "units",        "kg/kg"
        "profile_type", "fixed", "surface_value=1.e30" /
- "TRACER", "atmos_mod", "rainwat"
-           "longname",     "rain mixing ratio"
-           "units",        "kg/kg"
-       "profile_type", "fixed", "surface_value=1.e30" /
  "TRACER", "atmos_mod", "ice_wat"
            "longname",     "cloud ice mixing ratio"
            "units",        "kg/kg"
        "profile_type", "fixed", "surface_value=1.e30" /
- "TRACER", "atmos_mod", "snowwat"
-           "longname",     "snow mixing ratio"
+# prognostic rain water mixing ratio
+ "TRACER", "atmos_mod", "rainwat"
+           "longname",     "rain mixing ratio"
            "units",        "kg/kg"
        "profile_type", "fixed", "surface_value=1.e30" /
- "TRACER", "atmos_mod", "graupel"
-           "longname",     "graupel mixing ratio"
+# prognostic mass weighted rime factor
+ "TRACER", "atmos_mod", "q_rimef"
+           "longname",     "mass weighted rime factor"
            "units",        "kg/kg"
        "profile_type", "fixed", "surface_value=1.e30" /
 # prognostic ozone mixing ratio tracer