- Removed all options and conditionals for 'dblMom_{x} = .F.'; i.e. all hydrometeor categories are 2-moment only
- Removed the passing in of both V{t*} and V{t-1}. Now, only V{t*} is passed in. Process rates are based on this time level; its value is then updated to V{t}, which is passed out. All tendency variables (V_TEND) have been removed. Note, this was relevant for GEM only. In WRF, the wrapper subroutine (mp_milbrandt2mom_driver' no longer passes copies of the prognostic variables (which was done to make the main MY2 subroutine modular between GEM and WRF).
- Diagnostic quantities (including h_CB, h_ML1, etc., VISx, and SLW) have been removed. These should be computed as model post-processing (outside of microphysics)
- 2D diagnostic arrays (slabs in physics; 3D in model) SS01, SS02, etc. have been consolidated into a single array SS(i,k,n).
- Maximum values VsMax, VgMax, DgMax have been modified (to 4 m/s, 6 m/s, and 5 mm).
- Indexing for vertical levels has been generalized, with an input switch from driving model to determine if level k=1 is at the model top (GEM) or surface (WRF).
- Modified procedure to check for inconsistency in moment passed in at beginning of scheme, with initialization of number concentrations if inconsistent and clipping to zero for small values in subsaturated conditions.
- Replaced formulation of saturation vapor mixing ratio from that of Teten (orig. GEM physics, FOQSA) to that of Flatau (accurate polynomial approximation).
- Added back sublimation/deposition of graupel and hail.
- Reformulation of CNis (conversion from ice to snow). New approach results in much less upper-level ice and more mid-level snow. This addresses a known bias of excessive anvil ice and weak stratiform regions in convective systems.
- Reformulation of CNsg (conversion of snow to graupel).
- Reformulation of CNgh (conversion of graupel to hail). The new approach attempts to convert the portion of the graupel size distrition with diameters greater than the Schumann-Ludlam size limit. Also, the value of Dh_min, the size threshold below which hail is converted to graupel, has been reduced from 5 mm to 1 mm. The result of these changes is to produce more hail.
- The collection efficiencies Ecg and Erg are now formulated based on Cober and List (1993, JAS).
- The rain drop breakup formulation has been modified. It is now more aggressive and results in smaller mean rain sizes and increased evaporation rates.
- The condensation/evaporation section has been recoded, primarily for clarity.
- Homogeneous freezing of cloud can now occur after condensation.
- Changes to sedimentation:
+ The box-Lagrangian sedimentation has been replaced by a simpler (and cheaper) Eulerian forward-in-time advection.
+ Optimization has been added to compute sedimentation only in columns with non-negligible mass contents.
Code has also been added which computes sedimentation only on specific levels. This is for future use, for model configurations with very high vertical resolution.
- Unused functions in the module 'my_fncs_mod' have been removed.
- In the WRF interface, the array flipping before and after the call to the microphysics has been removed. This is no longer necessary since the vertical levels have now been generalized.