Implementing changes to the initial conditions. This commit contains code changes related to the modified initial conditions for convergence test simulations. These code changes can be activated by setting l_modify_ic_with_cubic_int = .true. in the namelist on a case-by-case basis.
Along with this option, the sounding profiles are also modified for the BOMEX, RICO, DYCOMS2_RF02 and Wangara cases.
-- For BOMEX, RICO, Wangara cases, we add more height levels in the original sounding profiles so that the cubic spline interpolation produces consistent profiles with those with linear interpolation
-- For DYCOMS2_RF02, instead of using the formulations in the code to derive the initial condition profiles (which will results in the grid-spacing dependent initial condition when we refine grid), we construct a sounding profiles (still use the same formulas as in the src/sounding.F90) on a high-resolution grid (refine the standard grid by a factor of 2^7), then save the profile in dycoms2_rf02_sounding.in. In this way, the model initiliazation will always read the same sounding profile for initialization when user refine the vertical model grid.
This commit most definitely does not change any bits, but commit 67878ef was BIT_CHANGING for the DYCOMS-II RF02 family of cases, RICO (and RICO SILHS), BOMEX, and Wangara.