This change log starts from 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)
WPS
Switch to 2.2.1 WPS code with the bug-fix for winds.
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.
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.
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)
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.
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)
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.
cold-start RR is switched to WPSV301 and WRFV301 (version 3.01)
As a result of this change, the RR1h cycle is now once again included
in the RAOB verification results available on the web.
-WRM
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
Cold-start RR - switch to Goddard shortwave radiation (ra_sw_physics=2)
Cold-start RR - TKE is added to WRF model output.
WRFPOST - switch to RUC-like visibility algorithm Implemented in both cold-start and cycled RR.
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.
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
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
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.
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)
Turn on 6th order diffusion:
diff_6th_opt = 2,
diff_6th_factor = 0.12
Turn on fractional sea ice option:
fractional_seaice = 1,
seaice_threshold = 271.4
Change in dyn_em/start_em.F to call physics initialization only once
- this saves ~500s of run time.
Changed diagnostics of T2 in phys/modelu_sf_sfcdiag.F
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.
More changes to module_sf_sfcdiags_ruclsm.F.
2-m diagnostics changes are implemented also in chemRR
Changes in WRFPOST: esveral constants are adjusted to match WRF;
Using TH array for potential temperature.
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.
HRRR-CHEM is switched to use WRF 3.1.1 code from /whome/smirnova/RR-svn/trunk/WRF311chem
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).
WRFPOST - fixed bug in the previous change (INITPOST.f)
WRFPOST - switch to revision 106
WRF - switch to WRF code from revision #3982
Cycled RR Change log after Nov 2008 (and some of the change log before Nov08)