When |
Who |
What |
2006/01/11/00 |
Ferrier, Jovic |
Change to BMJ code : At very cold temperatures (THETA=210 K) and pressures (2 mb),
saturation mixing ratios were set to zero rather than having infinite values.
|
2006/01/10/00 |
Gayno, Rogers, Wu |
Fixed an error in the script which reads in the SST, snow, and sea-ice analysis, and updates
them in the cycled WRF-NMM restart file at the start of the 00z NDAS. This error
prevented these fields from being updated with the latest analysis, so that the SST and
sea-ice used in each cycle was that from last date that the cycled NDAS was coldstarted
(00z 12/28 for both NAMX and NAMY), while snow would cycle on itself continuously.
Begin use of the latest (Dec 2005) release of the GSI analysis
|
2006/01/07/00 |
Ferrier |
Modified WRF-NMM code to output the correct skin temperature in the
array NMM_TSK, which is used by the GSI. Previously, this field reflected
the skin temperature at 00-h if RESTRT=.FALSE., which is only invoked
when the NMM cycle needs to be cold started.
|
2006/01/06/12 |
Ferrier |
With respect to the 2nd item in the 2005/11/17/00 change, the upward LW
at the surface is calculated from the lowest model layer temperature only
rather than the average of the skin and lowest model layer temperatures.
|
2005/12/28/00 |
Rogers |
The cycled NDAS for both NAMX and NAMY have been restarted from ops NAM
land states and ops GFS atmospheric states, to facilitate a direct comparison
between the full 12-km domain NAMX and the 1/2 sized 120km domain NAMY
|
2005/12/21/00 |
Mitchell |
Changed LAI (Leaf Area Index) in module_sf_lsm_nmm.F from 2 (ops NAM value) to 4.
|
2005/12/21/00 |
Black |
Two bug fixes:
- The array SH2O was inproperly dimensioned in its declaration in subroutine TURBL.
SH2O is a 3-D state array but in TURBL it is dimensioned as only 2-D. It is passed
from solve_nmm through TURBL and the surface driver into the land surface package.
SH2O is properly dimensioned in both the surface driver and the LSM package.
- Variable DZ was incorrectly calculated in subroutine CUCNVC. Two 1-D local
working arrays, namely T(K) and Q(K), were being used to generate DZ(I,K,J) but
T and Q are never being updated for each different I,J. DZ is subsequently passed
into the BMJ code where it is used to calculate the layer pressure thicknesses.
|
2005/12/20/00 |
Pyle, Janjic, Rogers, Wu |
Changes to z0 : enable a weak topography-based component and introduce
a component based on the standard deviation of the topography
Changes to SFCDIF:
- Reduction of ZT over elevated terrain is only active in stable surface layers.
- A baseline lower CZET-like value of 3.0 was imposed, regardless of surface elevation
- PBL height computation changed to make it consistent with MYJPBL code and post
Revert back to use of normalized RH (qoption=2) as the moisture analysis
variable in the GSI
Restarted cycle from ops GFS atmospheric states, ops NAM land states
|
2005/12/14/12 |
Pyle, Janjic, Rogers |
Switch back to new hybrid profile with concurrent modifications to the
the MKBND code (which makes WRF-NMM lateral boundary conditions from
GFS sigma coefficient files)
Restarted cycle from ops GFS atmospheric states, ops NAM land states
|
2005/12/13/00 |
Rogers |
Revert back to old hybrid coordinate; code changes to SI for new vertical profile
were inconsistent with the MKBND code (which makes WRF-NMM lateral boundary conditions from
GFS sigma coefficient files)
Restarted cycle from ops GFS atmospheric states, ops NAM land states
|
2005/12/9/00 |
Janjic, Pyle |
Switched to new hybrid vertical coordinate profile
Restarted cycle from ops GFS atmospheric states (necessary for hybrid profile change),
ops NAM land states
|
2005/12/5/12 |
Janjic |
Modified MYJPBL to eliminate PBL oscillations noticed by Frank Freedman
Added code in PROFQ2 so that when there is no q2 (i.e. when q2 is equal to epsq2)
el also drops to its floor value because el is a measure of spatial extension of q2.
|
2005/12/2/00 |
Janjic |
This change, put into the NAM NMM parallels on 10/20, was changed back since the original formulation was correct:
- Fixed an error in subroutine SFCDIF in which the
first guess z0T is used to compute z+z0T instead of the converged value of z0T obtained iteratively.
|
2005/11/24/00 |
Wu, Rogers |
In GSI set qoption=1, use pseudo relative humidity instead of normalized RH
as analysis variable.
|
2005/11/23/12 |
Ferrier |
Longwave temperature tendencies from the lowest two model layers are
averaged together in order to address concerns that the layers may be too thin
for the longwave scheme to provide accurate heating and cooling rates.
|
2005/11/18/00 |
Wu, Treadon |
Changes to the GSI analysis:
- Turn off multi-variate correlation between temperature and humidity above sigma level =0.15
- Install small bug fixes from official GSI directory
- Restart cycle from ops GFS atmospheric states and ops NAM land states
|
2005/11/17/00 |
Ferrier |
Because the GFDL longwave (LW) radiation calculates heating rates
and net radiative fluxes (rather than separately calculating the upward and
downward components), the downward (incoming) LW is calculated as the sum of
the net and the upward LW at the surface. The upward LW at the surface was
based on the ground (skin) temperature, but now this has been changed to be
the average of the skin and lowest model layer temperatures. The objective
is to slow the rate that downward LW decreases when stable PBL conditions
form in clear skies overnight.
|
2005/11/16/00 |
Ferrier, Rogers |
Cloud optical depths were reduced by the square of the cloud
fraction in order to increase solar transmission through clouds, but now
they are reduced by a factor proportional to the cloud fraction
(as is the case for the NAM).
Modified ozone profile to concentrate all ozone above 2 mb
(the model top pressure) into the top model layer instead of being
distributed through the atmosphere when the model top was at 50 mb; this is
identical to ozone profile in NAMXF2 parallel.
Cycled NDAS restarted from the NAMXF2 parallel
|
2005/11/15/00 |
Wu, Treadon |
Changes to the GSI analysis:
- Turned off AMSUA channels 1-5 and 15 for NOAA 15, 16 and 18.
These channels, which influence the temperature near the surface, were causing
analysis increments to be in the opposite direction from the raob increments in
Northern Canada (sfc emissivity issues over snow).
- Increased the thinning of satellite data (currently about one data point in each grid box).
The large amount of satellite data led to the raobs being
effectively thrown out in certain regions.
|
2005/11/10/00 |
Wu, Treadon |
In the GSI, code was added to turn off the moisture constraint (limq) when using qoption=2
(normalized RH) as the moisture analysis variable. Using this constraint caused a negative
feedback on the temperature analysis. While limq only affects the moisture field, the forcing
it generated negatively fed back through the moisture-temperature-pressure coupling when using
normalized RH as the moisture variable.
Cycled NDAS restarted from ops NAM.
|
2005/11/8/00 |
Pyle |
In the 00z 11/2 change, the topography-based modification to the ZILFC term was applied
regardless of the value of DTHV instead of only when DTHV > 0.5. SFCDIF was changed to apply
this modification in stable conditions when DTHV > 0.5, as was originally intended.
|
2005/11/4/12 |
Ferrier, Rogers |
Revert back to using ops NAM ozone distribution and restore the 3% reduction in the solar
constant
Started cycled NDAS from the NAMXF2 cycle
|
2005/11/4/12 |
Chuang, Ferrier, Ek, Gayno, Janjic, Lin, Rogers, Wong |
Details on previous changes from the repository version of the WRF-NMM
are listed in the
old NAMXF parallel change log and in the
current NAMXF2 change log
|