Change Log for 13-km development Rapid Refresh version at NOAA/ESRL/GSD

---

Main RR homepage

This change log starts from Dec 2007.

Starting September 2018 new changes are on the top.


1 may 2020 , 00 UTC

=========================================================================
RR-cold-2
Modifications to snow thermal conductivity and snow albedo for fresh snow.

9 April 2020 , 00 UTC

=========================================================================
RR-cold
WRFV3.9 - top of the master as of 20 Feb 2020.
Lake model is off.
=========================================================================
RR-cold-2
WRFV3.9 - top of rucfpe branch as of 7 April 2020.
1. Changes to MYNN
2. Changed to RUC LSM
- Cleaned up RUC LSM (land and ice are separate)
- New thresholds for snow model (land)
- New formulation for snow thermal conductivity as a function
of snow density over land (Sturm et al., 2017).
Lake model is off.
=========================================================================
RR-cold-3
WRFV3.9snow - top of master as of 3 Oct 2019.
New formulation for snow thermal conductivity as a function
of snow density (Sturm et al., 2017).
Lake model is off.
=========================================================================
RR-cold-4
WRFV3.9 - top of the master as of 3 Oct 2019.
Changed VEGPARM.TBL
Lake model is off.
=========================================================================

4 Oct 2019 , 12 UTC


Transition from theia to hera.


RR-cold - WRFV3.9 - top of the master as of 3 oct 19
Lake model is off.


RR-cold-2 - WRFV3.9snow, same as in RR-cold, but with the new formulation for snow thermal conductivity as a function
of snow density.
Lake model is off.


RR-cold-3 - WRFV3.9 as in RR-cold, but with old formulation of QVG initialization - WRFV3.9oldqvg.
Lake model is off.


RR-cold-4 - WRFV3.9
Lake model is off.
Changed VEGPARM.TBL

18 Jul 2019 , 00 UTC


RR-cold-4 - WRFV3.9
Changed VEGPARM.TBL also for wooded savanah: RS=125

18 Jul 2019 , 00 UTC


RR-cold-4 - WRFV3.9
Changed VEGPARM.TBL also for wooded savanah: RS=125

10 Jul 2019 , 00 UTC


RR-cold-3 - WRFV3.9 - control for RR-cold-4
9 soil model levels, lake is off

RR-cold-4 - WRFV3.9
Changed VEGPARM.TBL for grassland and broadleaf deciduous forest
1. IFOR is changed from 5 to 7 for grassland. This will turn on seasonal variation of Z0 in the grassland area based on the greenness.
(Similar to cropland and cropland/natural vegetation mosaic types).
2. RS for deciduous bradleaf forests is increased from 100. to 150. (similar to other forest types).
For now turned back to 100, because the smaller is the number the larger the transpiration.
Also, emissivity is increased from 0.93 to 0.94 as for other forest types. For now turned back to 100.
3. IFOR is changed from 3 to 2 for broadleaf deciduous forests. This will increase the rooting depth and increase evapotranspiration.
For now, set ifor=2 for all forests.

15 Jun 2019 , 00 UTC


RR-cold-3 - WRFV3.9 - control for RR-cold-4
9 soil model levels, lake is off

RR-cold-4 - WRFV3.9
6 soil model levels (4 layers, like in Noah LSM), lake is off.

12 Jun 2019 , 00 UTC


RR-cold - WRFV3.9 - control for RR-cold-2
Lake model is off.

RR-cold-2 - WRFV3.9snow with the new formulation for snow thermal conductivity as a function
of snow density.
Lake model is off.

1 Feb 2019 , 00 UTC


RR-cold-3 - WRFV3.9. Lake model is on. New geo* file with 15" landuse with Lakes.
Old soil data.
RR-cold-4 - WRFV3.9. Lake model is on. New geo* file with 15" landuse with Lakes.
New soil data.
Explore effect of new soil types.

8 Nov 2018 , 00 UTC


RR-cold - WRFV3.9 - change in lake model from Ayumi (GLERL), Zt is not equal to Zm.
Lake model is off.

RR-cold-2 - WRFV3.9_soilnew with changed run/SOILPARM.TBL
Wilting points for clay soils are reduced (recommendation from Johannes)
Lake model is off.
Explore effect of reduced wilting points.

6 Sep 2018 , 00 UTC


RR-cold - WRFV3.9 - change in lake model from Ayumi (GLERL), Zt is not equal to Zm.
Changes in RUC LSM except for the fc=max(qmin,ref*0.25).
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation:
use fc=max(qmin,ref*0.25) instead of fc=max(qmin,ref*0.5)
Explore only effect from fc change.

RR-cold-3 - WRFV3.9, change in lake model from Ayumi, changes in RUC LSM.
Lake model is on. Soil moisture adjustment is on.
RR-cold-4 - WRFV3.9, change in lake model from Ayumi, changes in RUC LSM.
Lake model is on. Soil moisture adjustment is off.
Effect from soil moisture adjustment.

=====================================================================================

5 Dec 2007

REAL, wrf.nl
Fix to sea-ice problem. Started using new sigma values in wrf.nl with lowest level now at 0.999 (previous 0.997). (Stan for Tanya)

10 Dec 2007

WPS
Switch to 2.2.1 WPS code with the bug-fix for winds.

11 Dec 2007

RR sigma levels
To avoid vertical instability the lowest two sigma levels in cold start RR are changed from 0.999, 0.995 to 0.998, 0.994. No change in HRRR3 runs.

6 Feb 2008

RR - new version of WRF code
The code is switched to beta-v3 for both WPS and WRF. WPS is build to be able to read GRIB2 data from GFS model. For SST - GRIB1 data is still used because the Vtable.SST needs to be updated for use with GRIB2 data. The DFI is turned on using new Hans's code for digital filter.

17 April 2008

RR - switch to WPSV3, WRFV3 and WPPV3.
The code is switched to v3 for both WPS and WRF. The only minor change (in WRFV3) was made to module_diagnostics.F to provide the noise level output in the units we are used to, and with different frequencies for DFI and forecast.
- DFI is still on.
- new Grell cumulus scheme - option 5
- ptop is changed to 10 mb (from 50 mb)

22 April 2008, 00 UTC

RR - change in WRFV3
- Updated phys/module_sf_ruclsm.F. Melting limits (Egglston et al. 1971) are introduced. Result - slower melting rate, especially in cases of warm advection over snow, when differences between skin temperature and first-level atmospheric temperature are significant. Also, the skin temperature is not limited to 273K under such conditions. Similar change to RUCLSM has been made in RUC, and it improved temperature forecast over snow.

29 April 2008

HRRR3 - switch to WPSV3, WRFV3 implemented in RR.
Several changes were made in WRFV3 to provide VIL output (Registry/Registy.EM, dyn_em/solve_em.F, namelist.input)

15 July 2008, 12 UTC

RR cold start - changes in WRFV3
New namelist variables are added to Registry.EM:
- mphdt - frequency of microphysics calls;
- time_step_dfi (DFI could have different time step);
- use_adaptive_time_step_dfi = .false. (no adaptive time step in DFI).
Changes to the code:
share/dfi.F - added reading in new namelist variables, and resetting time step to noraml after DFI step is done.
dyn_em/solve_em.F - added capability to call microphysics driver not every time step.
Starting 12z 15 July 08 these changes are included into cold-start RR.
Namelist variables are set to: mphdt=1, time_step_dfi=60, adaptive time step can vary between 60 and 100 s. (Only one run for 12z had diff_6th_opt=1, after than 6th order diffusion is turned off again) Also, quilt I/O is turned on.

7 August 2008, 12 UTC

cold-start RR is switched to WPSV301 and WRFV301 (version 3.01)

22 August 2008, 12 UTC

For the 1h cycle, sounding generation changed to correctly read the grib files, which evidently changed format around 12 UTC 8 August 2008. Previously, the grib files incorrectly identified each level as 'hybrid lev n'. Now, the grib levels are correctly identified as 'sigma=0.00nn' where nn is the level number.

As a result of this change, the RR1h cycle is now once again included in the RAOB verification results available on the web.
-WRM

25 August 2008, 12 UTC

RR 1h cycles uspdated to use WRFV3 and new GSI. Special features of GSI and ARW :

GSI: March 2008 NAM version with the following changes:
1) fixed lots of IO problems related to the use of ARW
2) moisture analysis (use qoption = 2 right now)
background error variance for qoption =1 (compute_derived.f90 from NCEP)
turned off relation between temperature and moisture (tv_to_tsen.f90)
qoption=2 works for ARW
3) parallel cloud analysis and cloud observation IO
4) turn off soil moisture and skin temperature updating.
5) normal wall time for GSI is 720 s with 64 processors
6) updating SST at 01Z

ARW: Version 3:
0) DFI is used (option 3, 40 min backward and 40 min forward)
1) use quilt to reduce IO time (8 processors)
2) combined wrf_input and wrfout
3) have boundary condition problems for dfi in last three cycles. right now use extra boundary fi les to bypass the problem.
4) normal wall time for ARW is 3200 s with 120 processors
5) use short DFI integration time step : time_step_dfi = 40,
6) microphysics step: mphdat = 1
7) Use new terrain that was smoothed by Joe in south boundary. The run started at 2008080723
8)Extend relax boundary zone for ARW forecast from 5 to 10 grid


21 October 2008, 12 UTC

Cold-start RR - switch to Goddard shortwave radiation (ra_sw_physics=2)

14 November 2008, 12 UTC

Cold-start RR - TKE is added to WRF model output.

15 December 2008, 12 UTC

WRFPOST - switch to RUC-like visibility algorithm Implemented in both cold-start and cycled RR.

25 December 2008, 12 UTC

WRFPOST - switch to RUC-like ceiling/cloud top algorithm To include convective clouds after a small change in Georg Grell's G3 scheme was made.

6 February 2009, 12 UTC

Added cycling of surface variables (snow water equivalent, snow depth, snow density, soil temperature and moisture, skin temperature)

9 February 2009 (17 Feb in cycled RR), 12 UTC

WRFPOST
- added option to output GRAUPEL on model and pressure surfaces;
- changed SLP computation using RUC-like procedure, it includes smoothing of temperature at 700 mb used in SLP computation and and also the final SLP (as done in RUC);
- changed computation of relative humidity on the pressure surfaces using RUC algorithm

4 March 2009 (17 Feb in cycled RR), 12 UTC

WRF model code - change in share/dfi.F subroutine not to use vertical damping
in the DFI backward integration.

Cold-start RR - switch to updated code, use damp_opt=3 with dampcoeff=0.05
to replace damp_opt=1 and damcoeff=0.02

19 March 2009 , 00 UTC

WRF model code - changes in RUCLSM to include parameterization of sea-ice/snow Changes to METGRID.TBL to allow initial snow be non-zero on the sea ice.

15 May 2009 , 12 UTC

Switch to 3.1 version of WRF code. The length of DFI integrations is reduced to 20 minutes from 40 min with the previous version of WRF code (3.1 version crashes on 26th minute of DFI backward integration)

13 June 2009 , 00 UTC

Turn on 6th order diffusion:
diff_6th_opt = 2,
diff_6th_factor = 0.12

26 June 2009 , 12 UTC

Turn on fractional sea ice option:
fractional_seaice = 1,
seaice_threshold = 271.4

22 July 2009 , 00 UTC

Change in dyn_em/start_em.F to call physics initialization only once
- this saves ~500s of run time.

23 July 2009 , 00 UTC

Changed diagnostics of T2 in phys/modelu_sf_sfcdiag.F


5 August 2009 , 00 UTC

Changes to phys/module_surface_driver.F, module_sf_ruclsm.F and addition of
new module - module_sf_sfcdiags_ruclsm.F. All these changes are made to ensure
consistency between different physics modules of WRF in interpretation of fluxes, surface
variables and diagnostics of 2-m values for temperature and moisture. Additional checks
are added to provide consistency bewtween 2-m and 1st level temperature
and moisture, also a saturation check on Q2 is added.

7 August 2009 , 00 UTC

More changes to module_sf_sfcdiags_ruclsm.F.
2-m diagnostics changes are implemented also in chemRR


14 August 2009 , 00 UTC

Changes in WRFPOST: esveral constants are adjusted to match WRF;
Using TH array for potential temperature.


8 October 2009 , 21 UTC

Fixed bug in definition of CONFLX in module_sf_ruclsm.F from thickness
between first full sigma level (array DZ8W) and the surface to half of this value.
This bug fix was made in cold-start RR, RR-chem, HRRR-chem and HRRR-RR.

13 October 2009 , 22 UTC

HRRR-CHEM is switched to use WRF 3.1.1 code from /whome/smirnova/RR-svn/trunk/WRF311chem


14 October 2009 , 22 UTC

WRFPOST - Use first level U, V, T and Q for 2-m T,Q and 10-m U, V. (INITPOST.f)
Also, switch to RUC-type dew point computation (SURFCE.f).

15 October 2009 , 18 UTC

WRFPOST - fixed bug in the previous change (INITPOST.f)


21 October 2009 , 18 UTC

WRFPOST - switch to revision 106


22 October 2009 , 12 UTC

WRF - switch to WRF code based on revision #3982 from NCAR WRF Repository


15 January 2010 , 00 UTC

Turned off 6th order diffusion and moist_adv_opt is changed from positive definite (1)
to monotonic (2)


29 January 2010 , 18 UTC

HRRR-CHEM.


Switch to 3.2 version of WRF code (from 25 Jan 2010).


12 March 2010 , 12 UTC

RR-cold-start
Switch to 3.2 version of WRF code (from 11 March 2010).
Constant BC in DFI.
Also, UNIPOST has replaced WRFPOST, and BUCKET and SFCEXTRACT tasks
in the XML are replaced with one task - COPYGB.

9 April 2010 , 12 UTC

RR-cold-start, RR-CHEM and HRRR-CHEM
Switch to official release 3.2 of version of WRF code.
Constant BC in DFI.
RR-CHEM - use shallow convection and radiation feedback


24 April 2010 , 12 UTC

RR-cold-start, RR-CHEM and HRRR-CHEM
Use tiling (numtiles = 4)in model runs which saves ~20% of run time
RR-CHEM - shallow convection is turned off because it causes NaNs in PM2.5 field


4 May 2010 , 18 UTC

RR-cold-start, RR-CHEM and HRRR-CHEM
Changed treatment of soil moisture variable when RUC LSM is used
Before the change wrfout* had soil moisture minus residual when RUC LSM is used.
With the change RUC LSM still has the same treatment of soil moisture variable internally
but wrfout* files will have total volumetric soil moisture at the model levels
This change will make it easier to compare soil moisture variable between different LSMs.
Graphics of RUC LSM soil moisture is more representative of actual soil moisture.
Changes include:
WPSV3.2 - ungrib/src/rprr.F (if initialized from RUC model add QMIN to soil moisture)
WRFV3.2 - dyn_em/module_initialize_real.F (don't subtract QMIN for RUC LSM initialized from Noah LSM
and don't add QMIN for Noah LSM initialized from RUC LSM)
phys/module_sf_ruclsm.F (subtract QMIN at the beginning and add at the end)
NO EFFECT on model performance.

26 May 2010 , 12 UTC

HRRR-CHEM
Switch to Unipost. Bucket and sfcextract steps are removed from XML.


8 June 2010 , 18 UTC

HRRR-CHEM
Use of mp_tend_lim=0.1 to prevent vertical CFL violations and crashes.


11 June 2010 , 00 UTC

RR-chem
Included computation of GD_CLOUD to WRF3.2 chem code.


11 June 2010 , 12 UTC

cold-start RR
Add shallow convection - ishallow=1.


2 August 2010 , 12 UTC

cold-start RR, RR-chem, HRRR-chem
WPS, WRF and UPP codes are recompiled in the new OS and run now on hcompnew nodes.


14 September 2010 , 00 UTC

RR-chem
From now on RR-CHEM will use GSI to assimilate aerosol


15 September 2010 , 18 UTC

RR-chem
From now on HRRR-CHEM will use GSI to assimilate aerosol


21 September 2010 , 00 UTC

HRRR-chem
The namelist parameter chemdt is changed from 3 to 30 min (for GOCART)


1 october 2010 , 18 UTC

RR-chem, HRRR-chem
Emissivity bug fix from Georg (for GOCART)


1 october 2010 , 18 UTC

RR-cold, RR-chem - WRF rev. 232
Bug fix in Grell's cumulusi scheme.
Changes in DFI code, compute QV after DFI from initial RH and filtered T
(Before QV was filtered and had only saturation check based on initial clouds.)

15 october 2010 , 00 UTC

RR-chem
Cold-started RR-chem, no ABBA fire data due to problems with the new format


15 october 2010 , 18 UTC

HRRR-chem
Cold-started HRRR-chem, no ABBA fire data due to problems with the new format


20 october 2010 , 12 UTC

cold-start RR
Backed off a change to recompute QV from initial RH in DFI, added QVG to
filtered variablesin DFI.

28 october 2010 , 12 UTC

started parallel cold-start RR on rotated lat/lon


29 october 2010 , 12 UTC

cold-start RR-RLL, cold-start RR, RR-chem, hrrr-chem
Bug fix for rotated lat/lon grid, and other changes in dynamic part from 3.2.1 WRF

2 November 2010 , 00 UTC

cold-start RR-RLL, cold-start RR
Updated /phys/module_bl_myjpbl.F and phys/module_sf_myjsfc.F to 3.2.1 version.

13 November 2010 , 00 UTC

cold-start RR-RLL, cold-start RR, RR-chem, HRRR-chem
Switch to use interp_type=2 (in LOG pressure) instead of interp_type=1 (in pressure)

3 December 2010 , 00 UTC

cold-start RR
Switch back to use interp_type=1.

8 December 2010 , 00 UTC

RR-chem
Switch back to use interp_type=1. (still no rebalancing)

11 December 2010 , 00 UTC

cold-start RR
Switch to use wrf_cntrl.parm from Geoff.

28 December 2010 , 00 UTC

cold-start RR
Switch to use Milbrandt microphysics.

8 June 2011 , 12 UTC

cold-start RR
Use revision 3.3 of WRF and WPS code. Modifications to RUC LSM to use mosaic approach to landuse properties. Switch form USGS to MODIS landuse classification.

8 June 2011 , 12 UTC

HRRR-chem
Use revision 3.3 of WRFchem and WPS code.

27 June 2011 , 12 UTC

cold-start RR
Modifications to RUC LSM to use mosaic approach also to soil properties.

28 June 2011 , 00 UTC

HRRR-chem
Fixed problem with the map projection which was introduced when switching to new version of prep_chem_sources.inp. Now the fires have correct locations.

20 October 2011 , 00 UTC

cold-start RR
Switch to WRF and WPS 3.3.1 merged with the RR modifications.

21 December 2011 , 00 UTC

cold-start RR
Turn on ishallow scheme (after bug fix from Georg).

20 June 2012 , 00 UTC

HRRR-chem
HRRR-chem is transferred from jet to zeus. It is now initialized from Rapid Refresh instead of RUC. Everything else remains the same as in the old HRRR-chem.

29 August 2012 , 12 UTC

RR-cold
Switch WPS and WRF to revision 3.4.1 with RR modifications merged in.

19 October 2012 , 12 UTC

RR-cold and HRRR-chem
Switch to 9 levels in soil.

19 November 2013 , 18 UTC

RR-cold
Switch WPS and WRF to version 3.5.1.

26 November 2013 , 18 UTC

RR-cold
Discontinued RR-cold-start using Lambert Conformal projection.
Started two parallel cold-start runs using rotated lat/lon projection on RAP dpmain.
The first run (RLL) is primary, and the second one (RLL2) is set to test different WRF options.
RLL - MYNN pbl
RLL2 - MYJ

27 November 2013 , 18 UTC

RR-cold
RLL - MYNN pbl, smooth_cg_topo = .true.
RLL2 - MYNN pbl, smooth_cg_topo = .false.

21 February 2014 , 00 UTC

RR-cold
RLL - snow density change (weighted average from snow, graupel and ice) - hypsometric_opt = 2 - swint_opt = 1 - smooth_cg_topo = .true.
RLL2 - same namelist option as in RLL, but thickness of top snow layer is changed from 3 cm to 5 cm.

21 February 2014 , 00 UTC

RR-cold RLL2 - change in 2-m T diagnostics (threshold coefficient is increased from 1.e-5 to 3.e-3) - John's formulations for snow, ice and graupel density

28 February 2014 , 00 UTC

RR-cold
RLL - old 2-m T diagnostics
RLL2 - uses new 2-m T diagnostics
Controlled test of the change to 2-m T diagnostics.

1 March 2014 , 00 UTC

RR-cold
RLL - old vegfrac and lai
RLL2 - MODIS greenness and LAI
Controlled test of MODIS LAI and grenness

25 March 2014 , 18 UTC

RR-cold
RLL - MODIS vegfrac and lai with smoothed topography
RLL2 - MODIS greenness and LAI without topo smoothing
Effects of topo smoothing

8 April 2014 , 12 UTC

RR-cold
RLL - MODIS vegfrac and lai with smoothed topography
RLL2 - MODIS greenness and LAI without topo smoothing plus
Joe's changes for shcu and rad. feedback.
Effects of topo smoothing

20 May 2014 , 18 UTC

HRRR-chem
Fixed reading of fire data. Set dust_opt=0, and turned off GSIchem.

30 May 2014 , 12 UTC

RR-cold
RLL - WRFV3.5.1 with MODIS vegfrac and lai.
RLL2 - WRFV3.6 with aerosol Thompson microphysics, aer_opt = 2
and Lake model.
Effects of WRFV3.6 physics

25 July 2014 , 00 UTC

RR-cold
RLL - WRFV3.6 with aerosol Thompson microphysics, aer_opt = 2
and Lake model.
RLL2 - WRFV3.6 with aerosol Thompson microphysics, aer_opt = 2
Lake model and MYNN changes.
Effects of MYNN changes in WRFV3.6

25 July 2014 , 00 UTC

RR-cold
RLL2 - RUC LSM has a change for cropland.
Also, reduced wilting point in SOILPARM.TBL.
Effects of MYNN and RUC LSM changes in WRFV3.6

23 October 2014 , 00 UTC

RR-cold
RLL - same as RLL2, but using new RUC LSM
Effects of RUC LSM changes in WRFV3.6

11 November 2014 , 00 UTC

RR-cold
RLL2 - use 0.25 degree GFS data

6 April 2015 , 12 UTC

RR-cold - use wrf3.7 version of RUC LSM with modified snhei_crit and snth
- revision 1116

28 January 2016 , 12 UTC

RR-cold - use WRF3.7.1 with old RUC LSM, 15" MODIS land-use
RR-cold-2 - use WRFV3.7.1 with new RUC LSM, 15" MODIS land-use
Effects of RUC LSM changes in WRFV3.7.1

13 April 2016 , 00 UTC


RR-cold-2 - Lawrence formulation for water intercepted by the canopy.
Effect: canopy evaporation should be reduced, cool and moist biases
in the afternoon should be also reduced.

5 May 2016 , 00 UTC

RR-cold - use WRF3.7.1 with old RUC LSM, 30" MODIS land-use
RR-cold-2 - use WRFV3.7.1 with new RUC LSM, 15" MODIS land-use
Both models will run on the large North America domain.
Effects of RUC LSM changes in WRFV3.7.1 and of 15" MODIS land-use.

9 May 2016 , 12 UTC

RR-cold - use 30" MODIS land-use
RR-cold-2 - use 15" MODIS land-use
Both models will run on the large North America domain.
Effects of 15" MODIS land-use.

7 June 2016 , 12 UTC

RR-cold - use MODIS veg.frac climatology
RR-cold-2 - real-time VIIRS GVF
Effects of real-time greenness

6 July 2016 , 12 UTC

RR-cold - replace HGT with the one with fixed heights.
RR-cold-2 - replace HGT with the one with fixed heights. Still uses real-time VIIRS GVF
Effects of real-time greenness

14 July 2016 , 12 UTC

RR-cold - no change
RR-cold-2 - aer_opt=3 (aerosol interacts with radiation)
Effects of real-time greenness and aer_opt=3 in RR-cold-2.

28 July 2016 , 12 UTC

RR-cold - real-time greenness, aer_opt=3
RR-cold-2 - real-time greenness, aer_opt=3 and LAI dependence in vegfrac
Effects of LAI dependence in vegfrac

5 August 2016 , 00 UTC

RR-cold-3 - hybrid vertical coordinate, etac=0.3
RR-cold-4 - original vertical coordinate, the same WRF code from Dave Gill (merged with GSD 3.7.1+).
Effect of hybrid vertical coordinate.

5 August 2016 , 00 UTC

RR-cold-3 - hybrid vertical coordinate, etac=0.1
RR-cold-4 - original vertical coordinate, the same WRF code from Dave Gill (merged with GSD 3.7.1+).
Effect of hybrid vertical coordinate.

19 August 2016 , 00 UTC

RR-cold - real-time greenness, aer_opt=3
RR-cold-2 - same, but with soil resisitance to evaporation.

23 August 2016 , 00 UTC

RR-cold-3 - hybrid vertical coordinate, use etac=0.2
RR-cold-4 - original vertical coordinate, the same WRF code from Dave Gill (merged with GSD 3.7.1+).
Effect of hybrid vertical coordinate.

27 August 2016 , 00 UTC

RR-cold - real-time greenness, aer_opt=3
RR-cold-2 - same, but with soil resisitance to evaporation and solar function in transpiration.

7 September 2016 , 00 UTC

RR-cold-3 - hybrid vertical coordinate - bug fix from Dave Gill.
RR-cold-4 - original vertical coordinate, the same WRF code from Dave Gill (merged with GSD 3.7.1+).
Effect of hybrid vertical coordinate.

15 September 2016 , 00 UTC

RR-cold - real-time greenness, aer_opt=3
RR-cold-2 - same, but with soil resisitance to evaporation. Solar function in transpiration is switched from Nohlan&Planton to Mahrer&Avissar.

15 September 2016 , 00 UTC

RR-cold - real-time greenness, aer_opt=3
RR-cold-2 - same, but with soil resisitance to evaporation and September modificatgions to MYNN pbl.
Compbined effect of recent physcis changes.

6 Feb 2017 , 00 UTC

RR-cold - real-time greenness, aer_opt=2
RR-cold-2 - real-time greenness, aer_opt=3.
Effect from the change of aer_opt from 2 to 3.

15 February 2017 , 00 UTC

RR-cold-3 - new UPP with aerosol-dependent visibility.
RR-cold-4 - old UPP.
Testing of aerosol dependency in the UPP visibility computation.

24 February 2017 , 00 UTC

RR-cold-3 - new UPP with aerosol-dependent visibility.
RR-cold-4 - old UPP with changed base for RH visibility from 60 to 90 km.
Testing of aerosol dependency in the UPP visibility computation.

1 March 2017 , 00 UTC

RR-cold-3 - new UPP with aerosol-dependent visibility.

WRFV3.8_GSD_hybridcoord code from 28feb17.
Last commit from Jaymes Kenyon - 80909f5
Turned on the EDMF scheme in MYNN:
bl_mynn_tkebudget = 0,
bl_mynn_tkeadvect = .false.,
bl_mynn_cloudpdf = 2,
bl_mynn_mixlength = 2,
bl_mynn_mixqt = 0,
bl_mynn_edmf = 1,
bl_mynn_edmf_mom = 0,
bl_mynn_edmf_tke = 0,
bl_mynn_edmf_part = 0,
icloud_bl = 1,

RR-cold-4 - old UPP with changed base for RH visibility from 60 to 90 km.
Testing of aerosol dependency in the UPP visibility computation.

WRFV3.8_GSD_hybridcoord code from 28feb17.
Last commit from Jaymes Kenyon - 80909f5
Turned off EDMF scheme in MYNN:
bl_mynn_tkebudget = 0,
bl_mynn_tkeadvect = .false.,
bl_mynn_cloudpdf = 0,
bl_mynn_mixlength = 1,
bl_mynn_mixqt = 0,
bl_mynn_edmf = 0,
bl_mynn_edmf_mom = 1,
bl_mynn_edmf_tke = 0,
bl_mynn_edmf_part = 0,
icloud_bl = 1,



1 March 2017 , 00 UTC

RR-cold-3 - new UPP with aerosol-dependent visibility.
WRFV3.8_GSD_hybridcoord with bug fix for hybrid coordinate
hybrid_opt = 2
etac = 0.2
RR-cold-4 - old UPP.
no hybrid.



17 March 2017 , 00 UTC

RR-cold-3 - new UPP with aerosol-dependent visibility.
WRFV3.8_GSD_hybridcoord with bug fix for hybrid coordinate
hybrid_opt = 2
etac = 0.2
Turned on the EDMF scheme in MYNN:
bl_mynn_tkebudget = 0,
bl_mynn_tkeadvect = .false.,
bl_mynn_cloudpdf = 2,
bl_mynn_mixlength = 2,
bl_mynn_mixqt = 0,
bl_mynn_edmf = 1,
bl_mynn_edmf_mom = 0,
bl_mynn_edmf_tke = 0,
bl_mynn_edmf_part = 0,
icloud_bl = 1,

RR-cold-4 - old UPP.
WRFV3.8_GSD_hybridcoord with bug fix for hybrid coordinate
hybrid_opt = 2
etac = 0.2
Turned off EDMF scheme in MYNN:
bl_mynn_tkebudget = 0,
bl_mynn_tkeadvect = .false.,
bl_mynn_cloudpdf = 0,
bl_mynn_mixlength = 1,
bl_mynn_mixqt = 0,
bl_mynn_edmf = 0,
bl_mynn_edmf_mom = 1,
bl_mynn_edmf_tke = 0,
bl_mynn_edmf_part = 0,
icloud_bl = 1,



28 March 2017 , 00 UTC

RR-cold-3 - the same.
RR-cold-4 - one change in WRF namelist:
ishallow = 1


30 March 2017 , 12 UTC

RR-cold-3 - updated WRFV3.8 code to the top of the trunk
RR-cold-4 - updated WRFV3.8 code to the top of the trunk
RR-cold - updated WRFV3.8 code to the top of the trunk


4 April 2017 , 00 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on.
RR-cold-4 - WRFV3.8 code with hybrid option turned off.


11 April 2017 , 00 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on, diff_opt=1;
RR-cold-4 - WRFV3.8 code with hybrid option on, diff_opt=2.


25 April 2017 , 00 UTC

RR-cold - WRFV3.8 code with hybrid option on;
RR-cold-2 - WRFV3.9 code with hybrid option on.


2 May 2017 , 00 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on;
RR-cold-4 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport off.

Result: worse performance in RR-cold-4.


4 May 2017 , 12 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on,

doubled smoothing in UPP;
RR-cold-4 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on.

Result: no effect, will repeat the test when verification on isobaric surfaces is available.


6 May 2017 , 00 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on.
RR-cold-4 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on, 1 pass smoothing in topography,zdamp=8000.

Result: worse performance in RR-cold-4.


9 May 2017 , 00 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on.
RR-cold-4 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on, doubled smoothing in UPP.

Result: only slight improvements in upper-air verification from additional smoothing in UPP.


9 May 2017 , 00 UTC

RR-cold - WRFV3.8 code with hybrid option on, gwd_opt=0 (default),real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, gwd_opt=3 - gravity wave drag on, real-time GVF.

Result: slight improvement in low-level winds in Western US.


12 May 2017 , 12 UTC

RR-cold-3 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on.
RR-cold-4 - WRFV3.8 code with hybrid option on, diff_opt=2, momentum transport on,

UPP reads actual u10 and v10, rather that 1-st level winds.

Result: actual diagnosed u10 and v10 increase RMS error and positive bias.


17 May 2017 , 00 UTC

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=3 - gravity wave drag on, real-time GVF.
RR-cold-3 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF.
RR-cold-4 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF, UPP post-process actual u10 and v10, rather that 1-st level winds.

23 May 2017 , 12 UTC

WRFV3.8 code has an update from Joe to compute U10/V10 diagnostics using neutral-log.

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=3 - gravity wave drag on, real-time GVF.
RR-cold-3 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF.
RR-cold-4 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF, UPP post-process actual u10 and v10, rather that 1-st level winds.

9 June 2017 , 00 UTC

WRFV3.8 code has an update from Joe to compute U10/V10 diagnostics using neutral-log.

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0, real-time GVF, turned off solar factor in RUC LSM.

Result: very close, RR-cold slightly better, especially for the Western US.

24 June 2017 , 00 UTC

WRFV3.8 code has an update from Joe to compute U10/V10 diagnostics using neutral-log.

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and topography data.

Result:

28 June 2017 , 12 UTC

WRFV3.8 code has an update from Joe to compute U10/V10 diagnostics using neutral-log.

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and old topo data.

Result: improved 2-m T and Td in the Western US, slight improvement in 2-m Td in the Eastern US, mixed result for 2-m T in Eastern US: better biases,
but slightly worse RMSE.

6 July 2017 , 00 UTC

WRFV3.8 code has an update from Joe to compute U10/V10 diagnostics using neutral-log.

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the old soil and new topo data.
Result: improvements in 2-m, T, Td and 10-m wind over Western US with the new topography, and slightly worse 2-m T RMSE over Eastern US.
Overall, improvement in RR-cold-2.


RR-cold-3 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF.
RR-cold-4 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF.

11 July 2017 , 00 UTC

WRFV3.8 code has an update from Joe to compute U10/V10 diagnostics using neutral-log.

RR-cold - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, update in GWD from Joe, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data.
Result:

13 July 2017 , 00 UTC


RR-cold-3 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data.
RR-cold-4 - WRFV3.8 code with hybrid option on, update in GWD, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
adjusted soil moisture initialized from the Noah model.
Result: reduced cold and moist biases.

17 July 2017 , 00 UTC


RR-cold-3 - WRFV3.8 code with hybrid option on, gwd_opt=0, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data.
RR-cold-4 - WRFV3.8 code with hybrid option on, gwd_opt=0, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model.
Result: reduced cold/moist bias in the lower troposphere from soil moisture adjustment.

30 August 2017 , 00 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data, fractional sea ice.

30 November 2017 , 00 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data, fractional sea
ice. Added read in GFS canopy water.

8 December 2017 , 00 UTC


RR-cold-3 - WRFV3.8 code with hybrid option on, gwd_opt=0, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model.
RR-cold-4 - WRFV3.8 code with hybrid option on, gwd_opt=0, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model, MODIS albedo.
Result from MODIS albedo test: Slightly increased cold bias and RMSE at 0z with the MODIS albedo.

12 January 2018 , 12 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data, fractional sea
ice. Use of GFS canopy water. Added MYNN sfc change to compute separate exchange coeff. for snow and snow-free reagions.

14 January 2018 , 12 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data, fractional sea
ice. Use of GFS canopy water. Removed MYNN sfc change to compute separate exchange coeff. for snow and snow-free reagions - buggy verson.

25 January 2018 , 00 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0 (default), real-time GVF;
RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data, fractional sea
ice. Use of GFS canopy water.
Added MYNN sfc change to compute separate exchange coeff. for snow and snow-free regions.

16 March 2018 , 00 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0 (default), real-time GVF;
Added soil moisture adjustment when initialized from GFS with the Noah LSM.
Added change in RUC LSM to compute snow density and not to use snow mosaic below 271 K.


RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=0, real-time GVF, new geo* file with the new soil and new topo data, MODIS albedo.
Use of GFS canopy water. Added MYNN sfc change to compute separate exchange coeff. for snow and snow-free regions.
Added soil moisture adjustment when initialized from GFS with the Noah LSM.
Added change in RUC LSM to compute snow density and not to use snow mosaic below 271 K.
Removed fractional sea ice.
Effect of new soil classes, GMTED2010 topo, MODIS albbedo, and MYNN sfc change in RR-cold-2 compared to control RR-cold.


RR-cold-3 - WRFV3.8 code with hybrid option on, gwd_opt=0, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model.
Switched to MODIS albedo.


RR-cold-4 - WRFV3.8 code with hybrid option on, gwd_opt=0, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model, MODIS albedo.
Fractional sea ice.
Effect from fractional sea ice.

19 April 2018 , 00 UTC


RR-cold-3 - WRFV3.8 code with hybrid option on, gwd_opt=3, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model.
Switched to MODIS albedo. Added fractional sea ice.


RR-cold-4 - WRFV3.9 code with hybrid option on, gwd_opt=3, diff_opt=2, climate GVF, new geo* file with the new soil and new topo data,
modified adjustment of soil moisture initialized from the Noah model, MODIS albedo with albedo dependence on solar angle. Fractional sea ice.
WRFv3.9 versus WRFv3.8.1

26 April 2018 , 00 UTC


RR-cold-2 - same as RR-cold-4, but with alternative 2-m diagnostics.

9 May 2018 , 00 UTC


RR-cold-4 - flux=.false. in module_sf_sfcdiags_ruclsm.F
RR-cold-2 - same as RR-cold-4, but with alternative 2-m diagnostics, flux=.false.

14 May 2018 , 12 UTC


RR-cold-4 - flux=.false. in module_sf_sfcdiags_ruclsm.F
RR-cold-2 - same as RR-cold-4, flux=.false., changes in surface MYNN.

16 May 2018 , 00 UTC


RR-cold-4 - added Jimy's fix for Q2: Q2<1.05Qlev1.
RR-cold-2 - same as RR-cold-4, flux=.false., changes in surface MYNN,
added Jimy's fix for Q2: Q2<1.05Qlev1.


1 June 2018 , 00 UTC


RR-cold - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=3, real-time GVF;
soil moisture adjustment. Old geo* file.
This run is a control for RR-cold-2.
RR-cold-2 - WRFV3.8 code with hybrid option on, diff_opt=2, gwd_opt=3, real-time GVF;
soil moisture adjustment. New geo* file (MODIS albedo and snow albedo, new topo, new soils)
Fractional sea ice.
Check effect of new geo* file and fractional sea ice.

1 June 2018 , 00 UTC


RR-cold-3 - WRFV3.9_mynn_stan - change in MYNN from Stan, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN,
new geo* file, albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
RR-cold-4 - WRFV3.9, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN, new geo* file.
albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
This run is a control for RR-cold-3 with Stan's change in MYNN.


8 June 2018 , 00 UTC


RR-cold-3 - WRFV3.9_mynn_stan - gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN,
new geo* file, albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
Reverted change from Stan in Eq. 7. New change from Jaymes in Eq. 9 (0.225 -> 0.250).


22 June 2018 , 00 UTC


RR-cold - WRFV3.9 code
RR-cold-2 - WRFV3.9_soilnew code, plus new soil moisture adjustment and change in bare soil evaporation.


26 June 2018 , 00 UTC


RR-cold - WRFV3.9 code, switch to climatological GVF, because NASA hasn't updated the data since February 2018.
RR-cold-2 - WRFV3.9_soilnew code, plus new soil moisture adjustment and change in bare soil evaporation.
Switched to climatological GVF.


4 July 2018 , 00 UTC


RR-cold - WRFV3.9 code, switch to climatological GVF, because NASA hasn't updated the data since February 2018.
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation.
use fc=max(qmin,ref) instead of fc=max(qmin,ref*0.5)
Soil moisture adjustment is pulled out.
Switched to climatological GVF.
Result: warmer and drier during the day, this change increases warm/dry biases in RAP.


21 July 2018 , 00 UTC


RR-cold-3 - WRFV3.9_lake - lake model is on, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN,
new geo* file, albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
RR-cold-4 - WRFV3.9, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN, new geo* file.
albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
This run is a control for RR-cold-3 with the lake model.


25 July 2018 , 00 UTC


RR-cold-4 - WRFV3.9_lake - lake model is on, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN,
new geo* file, albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
RR-cold-3 - WRFV3.9, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN, new geo* file.
albedo dependence on solar zenith angle. Old 2-m diagnostics.
Fractiona sea ice.
Switched places, now RR_cold-3 is control for RR-cold-4 with lakes.


4 Aug 2018 , 00 UTC


RR-cold - WRFV3.9 code, switch to climatological GVF, because NASA hasn't updated the data since February 2018.
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation:
use fc=max(qmin,ref*0.25) instead of fc=max(qmin,ref*0.5)
Result: reduced warm/dry bias.


5 Aug 2018 , 00 UTC


RR-cold - WRFV3.9 code, switch to climatological GVF, because NASA hasn't updated the data since February 2018.
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation:
use fc=max(qmin,ref*0.25) instead of fc=max(qmin,ref*0.5)
Also, remove iteration in energy budget when Qatm > Qsurface, which increases warm night-time bias.
Result:

14 Aug 2018 , 00 UTC


RR-cold - WRFV3.9 code, switch to climatological GVF, because NASA hasn't updated the data since February 2018.
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation:
use fc=max(qmin,ref*0.25) instead of fc=max(qmin,ref*0.5)
Also, remove iteration in energy budget when Qatm > Qsurface, which increases warm night-time bias.
Added bug fix for LAItoday.
Result:

27 Aug 2018 , 12 UTC


RR-cold - WRFV3.9 top of the trunk with MYNN changes for subgrid-scale clouds.
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation:
use fc=max(qmin,ref*0.25) instead of fc=max(qmin,ref*0.5)
RR-cold-3 - WRFV3.9, top of the trunk with MYNN changes for subgrid-scale clouds.
gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN, new geo* file.
albedo dependence on solar zenith angle. Old 2-m diagnostics.
This run is a control for RR-cold-3 with the lake model.
RR-cold-4 - WRFV3.9 top of the trunk with MYNN changes for subgrid-scale clouds.
Lake model is on, gwd_opt=3, diff_opt=2, climate GVF, change in sfc MYNN,
new geo* file, albedo dependence on solar zenith angle. Old 2-m diagnostics.


30 Aug 2018, 12 UTC


All four cold-start RAPs are updated to the top of the trunk WRFV3.9 code compiled for PNETCDF.
Unipost (UPP) is also updated, and scripts are changed to read new wrfout* files in PnetCDF
format. UPP is not writing out *grib1 files now, but creates grib2 files directly without convert,

as in real time RAP-dev1.

6 Sep 2018 , 00 UTC


RR-cold - WRFV3.9 - change in lake model from Ayumi (GLERL), Zt is not equal to Zm.
Changes in RUC LSM except for the fc=max(qmin,ref*0.25).
RR-cold-2 - WRFV3.9_soilnew code with change in bare soil evaporation:
use fc=max(qmin,ref*0.25) instead of fc=max(qmin,ref*0.5)
Explore only effect from fc change.
RR-cold-3 - WRFV3.9, change in lake model from Ayumi, changes in RUC LSM.
Lake model is on. Soil moisture adjustment is on.
RR-cold-4 - WRFV3.9, change in lake model from Ayumi, changes in RUC LSM.
Lake model is on. Soil moisture adjustment is off.
Effect from soil moisture adjustment.

Cycled RR Change log after Nov 2008 (and some of the change log before Nov08)