Mesoscale Parallel Experiment Change Log

Experiment Name:  North American NEMS/NMMB test

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:
    1. Assign a sub-type to level 2, level 2.5, and level3 radar data
    2. Limit use level 2.5 data in CONUS domain
    3. 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:
    1. Mixed forests : changed from 0.563 to 0.8
    2. Grasslands : changed from 0.035 to 0.08
    3. Croplands : changed from 0.035 to 0.07
    4. 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:
    1. Update to new version (r9527), reduces CPU consumption and now uses 2 parallel threads (~25 second speedup)
    2. New observation errors and retuned background errors
    3. 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:
    1. Closed and open shrublands : old=170, new=225
    2. Grasslands : old=20, new=70
    3. Croplands : old=60, new=70
    4. 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:
    1. fres=0.25 (resolution factor for dsp's)
    2. fr=1.00 (land factor for dsp's)
    3. fsl=0.75 (reduction factor for "slow" dsp's over land)
    4. 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:
    1. Uses expanded rain look up tables, increasing the maximum size of raindrops from 0.45 to 1.0 mm.
    2. Parameterizes sedimentation of 50-micron ice crystals.
    3. 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.
    4. NCW increased from 100 to 200 cm^-3.
    5. T_ice_init=0 results in faster initiation of ice. In 1D column runs it produced higher maximum surface rainfall rates.
    6. 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.
    7. VEL_INC~VRIMEF**2 – even faster rimed ice fall speeds than in the HAIL microphysics.
    8. These last two changes make the code faster and/or easier to read.
      1. Uses SQRT intrinsic function more, slight code optimization.
      2. 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:
    1. 6 km Alaska domain
    2. 4 km CONUS domain
    3. 3 km Hawaii domain
    4. 3 km Puerto Rico domain
    5. 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:
    1. 00Z cycle : 6 km Alaska domain
    2. 12z cycle : 4 km CONUS domain
    3. 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
    1. Update to new version (r7244)
    2. Turn on assimilation of GPS radio occultation (COSMIC) data
    3. 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:
    1. Fixed error with EQUIVALENCE statements for SAVEd arrays which were zeroing out lookup tables
    2. 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:
      1. Remove averaging of longwave temperature tendencies from the lowest two model layer
      2. 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).