When |
Who |
What |
2011/09/22/18 |
Janjic, Pyle, Black |
Prior to outputting the full non-hydrostatic vertical velocity, add back the
barotropic (external mode) component, which was earlier removed to alleviate a noise
problem in the CONUS nest domain run.
|
2011/08/09/18 |
Gayno, Rogers, X. Wang |
Due to compile glitch with the makefile for the parallel snow analysis
update code (nam_snow2mdl) this code was aborting in the NCO parallel. This means
that in the NCO parallel NDAS snow was continuously cycling on itself without an analysis update
since the last time that NCO restarted the production NDAS parallel from the EMC NDAS
parallel (which was ~mid June)
|
2011/07/14/18 |
Wu |
GSI analysis changed to not use mesonet surface pressure when 8th character of the station ID='X"
(meaning it is set to the standard atmosphere pressure or computed based on standard atmosphere).
|
2011/06/13/12 |
X. Wang, J. Zhu, Rogers |
NAMX development parallel run on devccs turned off; all NAMX verification and graphics jobs now pointing
to the NCO production NAM parallel, which officially started today. In the NCO production parallel all the NAM
nests will run 4x/day
|
2011/03/22/18 |
Rogers, Ferrier |
Revert back to the absorption coefficient values of 800 and 500, respectively,
for cloud water and ice in the radiation.
Restarted NDASX with cycled land-states and satellite bias corrections
from the NDASB parallel
|
2011/02/23/18 |
Rogers |
Turned off digital filter initialization
Restarted with cycled land-states from operational NDAS
|
2011/02/15/00 |
Ferrier |
Increase absorption coefficients for cloud water and ice in the
radiation to 1600 and 1000, respectively.
|
2011/02/09/00 |
Janjic, Rogers |
Changed parameter CODAMP (divergence damping constant) from 9 back to 12 for the
fire weather nest domains.
|
2011/02/05/18 |
Gayno, J. Wang |
Changed the NEMS code to read the base or snowfree albedo from the
NEMSIO input file created by NPS at the start of each NDAS run. The code was
prevously only reading the dynamic albedo, using it as the base albedo. This led
to erroneously high albedo over shallow/patchy snow cover.
|
2011/02/01/00 |
Janjic |
Changed parameter CODAMP (divergence damping constant) from 12 to 9 for the
all nest domains.
|
2011/02/01/00 |
S. Liu, Wu |
Changes to the GSI analysis:
- Assign a sub-type to level 2, level 2.5, and level3 radar data
- Limit use level 2.5 data in CONUS domain
- Increase observation error for level 3 data over CONUS domain.
|
2011/01/21/00 |
Ferrier |
A slight change is made to the new cloud fraction calculation
introduced into the NAMX in November 2010. Cloud fractions are set to 0
for QCLD (cloud water and total ice mixing ratios) <1.e-4 g/kg. A minimum
cloud fraction of 1% is assumed for QCLD=1.e-4 g/kg. Cloud fractions
increase by the square root of QCLD until a maximum cloud fraction of 100%
is reached when QCLD>=0.01 g/kg. In other words, the QCLD thresholds are
reduced by an order of magnitude compared to the November 2010 code.
|
2011/01/19/18 |
Janjic, Pyle |
To remove computational noise in the 1.333/1.5 km Fire Weather nest,
the NMMB model (in ./atmos/nmm/module_DYNAMICS_ROUTINES.F90) was changed to gradually reduce dw/dt
near the top of the model, by assuming attenuation of dw/dt of the form cos**2 in the
top 15 mb of the atmosphere, i.e. in the top 1.5% of the mass of the atmosphere.
The problem occurred in the fire weather run since it has high horizontal
resolution, combined with the high model top (2 mb), so the horizontal resolution near the
top becomes higher than the vertical resolution,
which is opposite to the situation lower in the atmosphere.
|
2011/01/06/00 |
Pyle, Rogers |
Control NDASB forecast now running with a digital filter with a +/- 40 minute window.
all nests)
|
2010/12/17/00 |
Gayno, Rogers |
Modified the land-sea masks for the CONUS/Alaska nests to resolve all lakes that are resolved by the
MMAB 1/12th degree RTG_SST; begin running surface update code to interpolate RTG_SST_HR analysis and IMS 4 km sea-ice
coverage to the CONUS 4 km and Alaska 6 km domains directly. The nest's snow cover is interpolated from the parent
12 km run, which is updated once/day and evolved during the cycled NDAS for 24-h in between analysis updates.
|
2010/12/09/12 |
Pyle |
NDAS forecast now running with a digital filter with a +/- 40 minute window.
Cloud fraction changes (see 11/22/2010 entry) implemented in the control run.
|
2010/11/27/00 |
Rogers |
When the MODIS-IGBP land-use definitions were implemented in the NAMB control run on 11/4,
the AVHRR weekly greenness fraction was inadvertently put into all NAM parallels at the same time.
The weekly greenness fraction data set was not intended to be implemented in the NAM-NEMS-NMMB change
package, so both NAMB, NAMX, and all nests will revert to using the operational monthly greenness
fraction as of this cycle.
|
2010/11/22/18 |
Ferrier, H-M. Lin |
Forecast cloud fractions associated with high, thin, cold cirrus clouds
were too high, and so a change was made in the grid-scale cloud fractions (Fc)
calculated from the GFDL radiation. Values of Fc were estimated from the total
relative humidity (RHtot) with respect to water (>0C) or ice (<0C), and this was
replaced by estimating Fc as a simple function of the cloud water and ice mixing
ratios. Cloud amounts were set to 1.0 in the internal GFDL calculations wherever
Fc>0 so as to preserve the long- and short-wave radiative fluxes and temperature
tendencies. This 9-slide PDF file
shows an example of the old and new Fc fields,
limitations of the old Fc calculation, and a description of the simple fix, which
is considered to be temporary until a more comprehensive approach is formulated.
Radiation code changed to calculate the mean value for the cosine of the zenith angle
every time step instead of every physics timestep
Enhanced diffusion for specific humidity and cloud water by 4x (equivalent to
setting smag2=0.8) was put into the NAMB control run.
|
2010/11/19/12 |
Black |
Nests turned back on, with one change: the number of boundary blending rows for the
Alaska nest was reduced from 5 to 3. With 5 blending rows, a 2:1 grid space ratio for the Alaska nest
to the parent, and only 2 halo rows for the parent task, the blending region was too wide. When Alaska
has its northern boundary extremely close to the northern boundary of a parent task then Alaska's
blending region can completely cross from one parent tasks' integration region through that parent
task's halo and into the next parent task's integration region, causing the code to hang. Reducing
the number of blending rows to 3 eliminates this problem.
|
2010/11/17/12 |
Rogers |
Nests temporaily turned off due to persistent code hangs.
|
2010/11/12/00 |
Janjic, Pyle, Rogers |
Changed parameter CODAMP (Divergence damping constant) from 9 to 12 for the
nest domains.
|
2010/11/10/18 |
Wong |
These IGBP roughness length values were changed:
- Mixed forests : changed from 0.563 to 0.8
- Grasslands : changed from 0.035 to 0.08
- Croplands : changed from 0.035 to 0.07
- Cropland/Natural Vegetation Mosaics : changed from 0.035 to 0.07
|
2010/11/09/18 |
Pyle, Janjic, Rogers |
This change put into NAMB on 4/29/10:
"Slight off-centering in time in the Crank-Nicholson scheme (parameter w1 in
/atmos/dyn/nmm/module_DYNAMICS_ROUTINES.F90 changed from 1.0 to 0.9). This change was made to eliminate
a failure in a test of a 4 km nested CONUS run, which was caused by large temperature oscillations in the
vertical over the Sierra Nevada mountains as a strong jet streak was impinging on the
mountains."
was inadvertently removed from the code (i.e., parameter w1 in subroutine ADV1
in ./atmos/dyn/nmm/module_DYNAMICS_ROUTINES.F90 reverted
back to 1.0), probably when the new NEMS code directory structure was put in the trunk in June 2010. w1 has been changed
back to 0.9 as intended.
|
2010/11/04/12 |
Ferrier, Janjic, Pyle, Wong, Ek |
MODIS_IGBP land-use definitions implemented in NAMB control run
Increase Smagorinsky constant for 2nd order diffusion (smag2) from 0.2 to 0.4 for
all nested runs
Enhance diffusion for specific humidity and cloud water by 4x (equivalent to
setting smag2=0.8) for the 12 km parent and all nests.
|
2010/11/03/18 |
Ferrier |
Bug fix in module_BL_MYJPBL: convert QSFC array to specific humidity as intended over land when CHKLOWQ=0
(95% RH).
|
2010/10/29/18 |
Wong |
The NOAH-LSM module with MODIS-IGBP land use has been changed by
linearly decreasing potential evaporation ETP to 0 when Bulk Richardson number Ri
is increasing from 0.9 to 1 (laminar flow), without imposing snow coverage
nor net radiative flux conditions. In subroutine NOPAC, for ETA being the
actual evapotranspiration, change the condition ETP>0.0 to ETP>0.000001
to avoid (ETA / ETP) being extremely large.
|
2010/10/20/18 |
Wu |
Changes to the GSI analysis:
- Update to new version (r9527), reduces CPU consumption and
now uses 2 parallel threads (~25 second speedup)
- New observation errors and retuned background errors
- Use of new data types : RASS virtual temperature, MAP winds, ASCAT winds
|
2010/10/09/00 |
Rogers, DiMego, Gayno |
Replace the 1/2 deg RTG_SST analysis with the high-resolution (1/12th degree) RTG_SST analysis.
|
2010/09/24/00 |
Vasic, Ferrier, Rogers |
Turned on gravity wave drag/mounatin blocking in the 6 km Alaska nest.
|
2010/09/21/18 |
Janjic |
Changed the routine that computes the pressure gradient force (pgforce):
zeroed contributions to the "divergence correction" term along lateral boundaries.
|
2010/09/21/18 |
Wong, Ek |
The following changes were made to the MODIS_IGBP values for minimum canopy
resistance (all in sec/m) for these land use types:
- Closed and open shrublands : old=170, new=225
- Grasslands : old=20, new=70
- Croplands : old=60, new=70
- Cropland/natural vegetable mosaic : old=40, new=70
|
2010/09/16/18 |
Rogers |
Both NDASB and NDASX restarted from satellite bias corrections from 12z 9/14 NDASX run
to ensure cleaner comparison of the MODIS_IGBP land use (NAMX) vs USGS land use (in NAMB)
|
2010/09/16/00 |
Rogers |
Turned off use of ob type 111 (synthetic tropical cyclone center surface
pressure and specific humidity profile) in the GSI analysis, due to it having a large impact that it
might obscure the comparison of the NAMB vs NAMX test of MODIS_IGBP vs USGS land use designations.
|
2010/09/14/18 |
Wu, Rogers |
Turn on use of ob type 111 (synthetic tropical cyclone center surface
pressure and specific humidity profile) in the GSI analysis.
|
2010/09/14/12 |
Ferrier, Ek, Gayno, Pyle, Rogers, Wong, Wu |
8/4/2010 microphysics changes put into NAMB control run.
Begin use of MODIS_IGBP land use definitions in NAMX and all nested runs.
|
2010/09/02/00 |
Parrish, Wu, Rogers |
Initial conditions for each nest are now created by a GSI analysis
on the nest domain using the NDAS first guess, instead of interpolating
the 12 km NAM analysis to the nest domain.
|
2010/08/30/00 |
Ferrier, Janjic, Vasic |
Begin using the BMJ_DEV convection option in the nested runs with these settings:
- fres=0.25 (resolution factor for dsp's)
- fr=1.00 (land factor for dsp's)
- fsl=0.75 (reduction factor for "slow" dsp's over land)
- fss=0.75 (reduction factor for "slow" dsp's over water)
|
2010/08/16/12 |
Pyle, Rogers |
Begin direct use of GDAS spectral coefficients in NPS (rather than a GDAS 1 deg lat/lon
grid) to create the first guess at the start of the NDAS.
|
2010/08/06/12 |
Black, Gayno, Pyle, Rogers |
3 km Puerto Rico
nested domain enlarged to cover all of Hispanola.
|
2010/08/04/12 |
Ferrier |
NAMX (and all 4 nests) turned back on with the following microphysics changes:
- Uses expanded rain look up tables, increasing the maximum size of raindrops from 0.45 to 1.0 mm.
- Parameterizes sedimentation of 50-micron ice crystals.
- New cloud water to rain autoconversion (PRAUT) scheme following work of Liu & Daum et al. This new
scheme no longer requires an autoconversion threshold, and the functional dependence of cloud water
mixing ratio is quite different from the original Kessler-based formulation.
- NCW increased from 100 to 200 cm^-3.
- T_ice_init=0 results in faster initiation of ice. In 1D column runs it produced higher maximum
surface rainfall rates.
- FLARGE2=0.07 based on tests by JMA visiting scientist Nakagawa in his global runs. The parameter
determines the number fraction of large precipitation ice particles.
- VEL_INC~VRIMEF**2 – even faster rimed ice fall speeds than in the HAIL microphysics.
- These last two changes make the code faster and/or easier to read.
- Uses SQRT intrinsic function more, slight code optimization.
- New GET_INDEXR function to get an estimate for the mean drop size in microns (INDEXR).
|
2010/07/26/12 |
Rogers |
Use of modified vertical level distribution with more layers in the stratosphere
implemented in NAMB control run; NAMX (and all 4 nests) turned off temporarily pending development
of next change to test.
|
2010/07/14/00 |
Wu, Pondeca |
Begin use of dynamic reject list (from RTMA) for mesonet data in the GSI analysis.
|
2010/07/12/12 |
Rogers |
Begin running all 4
nested domains in the NAMX parallel for both 00z and
12z cycles:
- 6 km Alaska domain
- 4 km CONUS domain
- 3 km Hawaii domain
- 3 km Puerto Rico domain
- Nests run to 60-h, with explicit convection
|
2010/05/20/00 |
Pyle, Rogers |
The NPS job that was making initial conditions for the nested runs was producing
spurious cloud water/ice fields because these variables were only output up to the 30 mb level in
the NAM analysis GRIB file. NPS was substituting other fields for the cloud water/ice values
above 30 mb. The change to add more vertical layers in the stratosphere on 5/18 exacerbated the
problem. This was fixed by outputting cloud water/ice/rain/snow up to 2 mb in the NAM analysis
GRIB file, as is done for all other variables.
|
2010/05/18/18 |
Pyle, Rogers, Wu |
Turned off use of global model ozone in the GSI analysis.
Use modified vertical level distribution with more layers in the stratosphere
|
2010/05/17/12 |
Wu, Rogers |
Turn on use of global model ozone in the GSI analysis.
Commence routine nested runs in the NAMX parallel:
- 00Z cycle : 6 km Alaska domain
- 12z cycle : 4 km CONUS domain
- Nests run to 60-h, and with explicit convection
The boundary condition treatment changed to specify 5 boundary rows (LNSH, LNSV=5) instead
of one row was put into the NAMB control run.
|
2010/04/29/12 |
Pyle |
Slight off-centering in time in the Crank-Nicholson scheme (parameter w1 in
/atmos/dyn/nmm/module_DYNAMICS_ROUTINES.F90 changed from 1.0 to 0.9). This change was made to eliminate
a failure in a test of a 4 km nested CONUS run, which was caused by large temperature oscillations in the
vertical over the Sierra Nevada mountains as a strong jet streak was impinging on the
mountains.
|
2010/04/05/18 |
Wu, Parrish |
Changes to the GSI analysis
- Update to new version (r7244)
- Turn on assimilation of GPS radio occultation (COSMIC) data
- GSI analysis will now update the first guess 2-m temperature/specific humidity and 10-m wind
|
2010/03/30/12 |
Pyle, Ferrier |
Boundary condition treatment changed to specify 5 boundary rows (LNSH, LNSV=5) instead
of one row.
Reduction of the absorbtion coefficients for water and ice by one-half to 800
and 500, respectively, was placed in the NAMB control run.
|
2010/03/18/18 |
Ferrier |
Bug fix in the radiation code to remove spikes in the vertical profile of shortwave heating
|
2010/03/09/00 |
Ferrier |
Longwave bug fixes (see 18z 2/26 entry) installed in NAMB control run
In the radiation parameterization, the absorbtion coefficients for water and ice have been halved to 800
and 500, respectively.
|
2010/02/26/18 |
Ferrier |
Bug fixes to the radiation code:
- Fixed error with EQUIVALENCE statements for SAVEd arrays which were zeroing out lookup tables
- Removed two modifications to longwave radiation made for the June 2006 WRF-in-NAM implementation,
that were removed from the ops NAM in March 2008:
- Remove averaging of longwave temperature tendencies from the lowest two model layer
- Compute the upward LW at the surface based on the ground (skin) temperature instead
of the average of the skin and lowest model layer temperatures.
Restored the absorbtion coefficients for water and ice to 1600 and 1000 as they are
in the ops NAM and the NAMB control run
|
2010/02/17/12 |
Rogers |
Due to an oversight, the long-term precipitation budget from the operational NDAS
has been used in every parallel NDAS run since 12z 2/2/10, instead of having it cycle independently.
This has been fixed as of this cycle.
|
2010/02/02/12 |
Rogers, Lin, Gayno |
The b-grid copygb utility is used to 1) interpolate the daily CPC precip analysis
to the b-grid (for the long-term precip budget used to bias-correct the precip analysis
assimilated in the NDAS) and 2) interpolate the merged Stage II/IV precipitation analysis used in the NDAS forecast
as forcing to the land-surface model. On 25 January, the b-grid copygb binary was changed on cirrus only, which
corrupted these interpolations to the b-grid. This lasted until 27 January, when the parallel runs
switched to stratus, where the copygb binary was unchanged. Because of this finding, it was decided
to restart the NDASX land-states and long-term precipitation budget from the operational NDAS.
|
2010/01/27/12 |
Ferrier |
In the radiation parameterization, the absorbtion coefficients for water and ice have been halved to 800
and 500, respectively.
NDASX land states and the satellite bias corrections for the GSI analysis were restarted from the NDASB control.
|
2010/01/09/00 |
Treadon |
Turned off assimilation of NOAA-19 AMSUA channel 8 radiances due to excessive noise.
|
2009/12/18/06 |
Janjic |
Removed the factor sqrt(2) from the denominator in all the definitions
of hdac* in ./src/nmm_libutil/module_CONTROL.F90. This will lead to increased horizontal
diffusion.
|
2009/12/11/00 |
Janjic, Pyle |
Changed equation for parameter "steep" in the tracer advection so that it equals the preferred
value of 0.96 if parameter idtadt (number of adjustment timesteps between tracer advection calls) is set to 2
|
2009/12/08/00 |
Wu, Keyser |
"WDSATR" : superob'ed Scatterometer winds over oceans, were not used in the assimilation
until this cycle because the prepbufr parm card did not set the observation error for this
ob type, so they were given a quality mark=9, and therefore not used in the GSI analysis.
|
2009/12/01/12 |
NMMB/NEMS Development Team |
First run of new NMMB parallel with tracer advection, parallel using with all
changes running in the NEMS/NMMB control parallel run (NAMB).
|