Updates to P3 scheme (#1648)

TYPE: enhancement and bug fix 

KEYWORDS: P3, microphysics

SOURCE: Hugh Morrison (NCAR MMM)

DESCRIPTION OF CHANGES:
SUMMARY: Updates to improve the P3 scheme (re-factor 3-moment lookup table and minor bug fixes)

1. Reduce 3-moment lookup table file size by an order of magnitude. Re-factoring the lookup table 
dimensions for 3-moment greatly reduces file size with little loss of accuracy. Re-factoring of the 
3-moment lookup table so the dimension is in the space of mu (gamma size distribution shape parameter) 
instead of normalized 6th moment. This reduces file size by an order of magnitude. Unzipping of the file in 
the primary Makefile is no longer needed (file in run/ is no longer gzipped).

2. Improve flexibility by including Kogan (2013, JAS) option for cloud autoconversion/accretion. Addition of 
Kogan (2013) as an option for warm rain (set by iparam=4).

3. Minor bug fixes. 
   * Set rain-ice collection parameter to -99 when rain is not present, instead of 0 (since this is in log space) 
     --> non-answer changing
   * Fix to the double log bug for ice-rain collection
   * Minor LT2 fixes --> replace a few instances of cgp1 with cgp2, and array cgp1 and cgp2 in both i_Fr and 
      i_rhor instead of just rho_r.
   * Update to 2 of the 3 lookup table files and offline lookup table creation codes and corresponding changes 
      to the lookup table filenames in primary Makefile. 

LIST OF MODIFIED FILES: 
M phys/module_mp_p3.F
M ./Makefile 
M run/p3_lookupTable_1.dat-3momI_v5.1.7 --> lookup table file
M run/p3_lookupTable_2.dat-2momI_v5.2.3 --> lookup table file
M run/create_p3_lookupTable_1-2mom_v5.3.f90 --> offline lookup table generation code
M run/create_p3_lookupTable_1-3mom_v5.1.7.f90 --> offline lookup table generation code
M run/create_p3_lookupTable_2-v5.2.3.f90 --> offline lookup table generation code 
D run/p3_lookupTable_1.dat-5.3-3momI.gz 

TESTS CONDUCTED: 
1. The mods fix the problems identified. 
2. Extensive testing in idealized supercell case and real cases in Canadian GEM model.
3. Jenkins tests are all pass.

RELEASE NOTE: Updates to 4 options of the P3 scheme (MP options 50, 51, 52, 53). The 3-moment lookup table file size is reduced by an order of magnitude (file names are changed for the lookup tables).  Improved flexibility has been added by including Kogan (2013, JAS) option for cloud autoconversion/accretion. A few minor bug fixes to harden the code have been added. References are:

Morrison, H., and J. A. Milbrandt, 2015: Parameterization of cloud microphysics based on the prediction of bulk particle properties. Part 1: Scheme description and idealized tests. J. Atmos. Sci., 72, 287-311.

Milbrandt, J. A., and H. Morrison, 2016: Parameterization of cloud microphysics based on the prediction of bulk particle properties. Part 3: Introduction of multiple free categories. J. Atmos. Sci., 73, 975-995.

Milbrandt, J. A., H. Morrison, D. T. Dawson II, and M. Paukert, 2021: A triple-moment representation of ice in the Predicted Particle Properties (P3) microphysics scheme. J. Atmos. Sci., 439-458.

Makefile

@@ -116,10 +116,6 @@ wrf : framework_only
 	if [ $(ESMF_COUPLING) -eq 1 ] ; then \
 	  ( cd main ; $(MAKE) RLFLAGS="$(RLFLAGS)" MODULE_DIRS="$(ALL_MODULES)" SOLVER=em em_wrf_SST_ESMF ) ; \
 	fi
-	if [ ! -f run/p3_lookupTable_1.dat-5.3-3momI ] ; then \
-	  ( cd run ; cp p3_lookupTable_1.dat-5.3-3momI.gz hold.gz ; \
-	    gunzip hold.gz ; mv hold p3_lookupTable_1.dat-5.3-3momI ) ; \
-	fi
 	@echo "build started:   $(START_OF_COMPILE)"
 	@echo "build completed:" `date`
 
@@ -574,8 +570,8 @@ em_real : wrf
                ln -sf ../../run/eclipse_besselian_elements.dat . ;     \
                ln -sf ../../run/CCN_ACTIVATE.BIN . ;                   \
                ln -sf ../../run/p3_lookupTable_1.dat-5.3-2momI . ;     \
-               ln -sf ../../run/p3_lookupTable_1.dat-5.3-3momI . ;     \
-               ln -sf ../../run/p3_lookupTable_2.dat-2momI_v5.2.2 . ;  \
+               ln -sf ../../run/p3_lookupTable_1.dat-3momI_v5.1.7 . ;     \
+               ln -sf ../../run/p3_lookupTable_2.dat-2momI_v5.2.3 . ;  \
                ln -sf ../../run/HLC.TBL . ;                            \
                ln -sf ../../run/wind-turbine-1.tbl . ;                 \
                ln -sf ../../run/ishmael-gamma-tab.bin . ;              \
@@ -656,8 +652,8 @@ em_real : wrf
              ln -sf ../../run/bulkradii.asc_s_0_03_0_9 . ;          \
              ln -sf ../../run/CCN_ACTIVATE.BIN . ;                  \
              ln -sf ../../run/p3_lookupTable_1.dat-5.3-2momI . ;    \
-             ln -sf ../../run/p3_lookupTable_1.dat-5.3-3momI . ;    \
-             ln -sf ../../run/p3_lookupTable_2.dat-2momI_v5.2.2 . ; \
+             ln -sf ../../run/p3_lookupTable_1.dat-3momI_v5.1.7 . ;    \
+             ln -sf ../../run/p3_lookupTable_2.dat-2momI_v5.2.3 . ; \
              ln -sf ../../run/HLC.TBL . ;                           \
              ln -sf ../../run/wind-turbine-1.tbl . ;                \
              ln -sf ../../run/ishmael-gamma-tab.bin . ;             \

run/create_p3_lookupTable_1.f90-v5.3 → run/create_p3_lookupTable_1-2mom_v5.3.f90

@@ -13,8 +13,8 @@ PROGRAM create_p3_lookuptable_1
 ! All other parameter settings are linked uniquely to the version number.
 !
 !--------------------------------------------------------------------------------------
-! Version:       5.1.7_b04 + mods
-! Last modified: 2021-FEBRUARY
+! Version:       5.3
+! Last modified: 2022-JAN
 !______________________________________________________________________________________
 
 !______________________________________________________________________________________