Mesoscale Parallel Experiment Change Log

Experiment Name:  North American NEMS/NMMB test

When Who What
2014/03/08/00 W.-S. Wu
  • Updates to the GSI analysis:
    1. Turn on the surface obs in file nam_regional_convinfo.txt_cld
    2. Turn off type 250 which was on by mistake when updating satwnd subtypes.
    3. Put back gross error for surface winds to 5 m/s as in current operational NDAS; 6 m/s is used in the global.
  • 2014/03/06/00 W.-S. Wu, Melchior, Levine, Y. Ling
  • Begin using the new prepbufr report types 192, 193, 194, and 195 (and the corresponding 200-series wind component) in the GSI analysis:
    1. 192 - sfc reports from synop land missing station and msl pressure
    2. 193 - sfc reports from metar missing station and msl pressure
    3. 194 - sfc reports from ship/cman/buoy missing station and msl pressure
    4. 195 - sfc reports from mesonet missing station and msl pressure
  • 2014/03/03/18 Rogers, Ferrier
  • The gravity wave/mountain blocking parameter cleffamp was changed from 2.244 to 1.0
  • 2014/02/16/00 Ferrier
  • Fixes to the microphyiscs changes made on 1/30/2014:
    1. The original jan19 version can result in NSImax=0 if there is no ice present, and that in turn could incorrectly lead to no ice present with NInuclei=0.
    2. Fix made that set NSImax to be no less than NSI_max=250.e3 (250 per liter), so that NInuclei would be set to that value. This should lead to cloud ice now appearing at very cold temperatures. Note that NSImax can be >>250.e3 when there is a lot of ice present (like at the tops of convection); the value for NSI_max(=250.e3) just establishes a lower limit.
  • 2014/02/11/18 S. Liu, W.-S. Wu, Ferrier
  • Turn off separate microphysics species advection
  • GSI analysis / NDAS changes
    1. Turn on the diurnal reject lists and wbinuselist from RTMA
    2. Turn on cloud analysis and diabatic digital filter initialization with radar-derived temperature tendencies during the NDAS tm09, tm06, and tm03 analysis/forecast
  • 2014/02/10/18 Y. Wu, Ferrier
  • In the land-surface physics, the accumulated snow (array ACSNOW) calculation was changed to denote that only the frozen portion of the precipitation should contribute to the input surface snowfall.
  • 2014/01/30/18 Ferrier, Aligo
  • Microphysics changes ("jan19 version")
    1. Turn on separate microphysics species advection
    2. The maximum number concentration of ice (NLImax) is a function of the rime factor (RF) and temperature.
      1. For conditions when ice densities are high (RF>=10), NLImax=1.e3, which allows for mean ice diameters to exceed 1 mm in size. This so-called "hail mode", which allows the mean ice diameters to exceed 1 mm in size, only occurs in conditions associated with intense convection.
      2. For conditions when ice densities are not high (RF<10), NLImax is 10 L-1 at 0C and decreases to 5 L-1 at <=-40C. When the mean diameters of ice particles reach 1 mm in size, then the criterion is relaxed and the number concentrations of ice (NLICE) can exceed NLImax. This change helps with promoting the development of stratiform precipitation.
    3. Never adjust RF values when the number concentrations of ice (NLICE) exceed the maximum value (NLImax). The previous code would artificially increase the rime factor (ice densities) so that NLICE=NLImax, but this is no longer justified when the mass-weighted rime factor is being advected. This change also helps with promoting the development of stratiform precipitation.
    4. Ice deposition does not change rime factor values (RF) of ice when RF>10 at temperatures >-40C. This change preserves helps with preserving the presence of high-density ice particles in intense convective cores.
    5. Limit the increase in ice fall speeds due to reduced air resistance at low pressures (high altitudes) to not exceed a factor of 1.5 times the fall speeds at a reference pressure of 1000 hPa. This change prevents ice particles from falling too quickly in the upper troposphere.
    6. At cold temperatures (<-30C), the maximum number concentrations of small ice crystals (NSImax) is calculated based on assuming 10% of the total ice content is due to cloud ice. This change allows for higher amounts of cloud ice to exist at the tops of intense onvective cores, which can then gradually be converted into snow as parcels move away from areas of intense convection.
  • 2014/01/27/18 Ferrier, H.-M. Lin
  • The RRTM radiation was changed to revert back to the solar constant value of 1366 W/m2, instead of 1360.8 W/m2 which was being used inadvertently
  • 2014/01/23/18 Rogers, Janjic
  • Increased the amount of non-zero residual deformation applied in the top 3 model layers in the stratosphere by changing this equation in modules_DYNAMICS_ROUTINES:




    where l=the model level. This change eliminated the instability in the top model layers that caused the failures in the 00z 1/3 and 00z 1/4 NAMX runs in the Alaska nest.

  • 2014/01/22/00 Y. Zhu
  • Fixed a typo in the GSI scripts which caused the satbiaspc files to have all zero values and for the satbias files to evolve slowly.
  • 2014/01/09/06 Rogers, Katz
  • Terminate the use of GLERL water temperatures over the Great Lakes in the NDAS SST update due to a persistent warm bias in eastern Lake Superior. The surface water temperatures from the RTG_SST_HR will be used over the Great Lakes.
  • 2014/01/07/06 Rogers, Gayno
  • NAMB/NDASB will no longer use the AFWA snow depth analysis. If the IMS snow cover analysis indicates no snow cover, it will removes now. If IMS snow cover analysis indicates snow cover where the NDASB has none, it will set the snow depth to 5 cm. Otherwise, the snow depth will be cycled.
  • 2013/12/13/00 Ferrier, H-M. Lin
  • Replace the parameterization of effective radius of (cloud) ice from Wyser (1998, J. Clim) with that from McFarquhar & Heymsfield (1996, J. Atmos. Sci.). Effect will be to reduce upper tropospheric warm bias caused by the Wyser scheme with the model top at 2 mb.
  • 2013/12/06/18 Wu, Yanqiu Zhu
  • Turn on cloud detection channels for monitored instruments
  • Turn off known bad channels:
    1. NOAA-19 AMSU-A Ch 7
    2. NOAA-19 MHS Ch 3
  • Begin use of new RTMA (to be implemented in January 2014) reject lists as the primary source; use operational RTMA reject lists as a backup
  • 2013/12/05/18 Gayno, Rogers
  • Restarted the full cycling of NDAS land states from the operational NDAS due to the parallel land states being degraded starting ~ 6/1/2013 when the NMMB model was changed to output latitide/longitude in radians instead of degrees.
  • 2013/11/25/12 Ferrier, H-M. Lin, Janjic
  • Provide the full time step (DT) rather than integer of the time step (DT_INT) into the RRTM radiation driver, and provide the double-precision version of the full sigma array into the RRTM's subroutine rad_initialize (affected only diagnostic fields).
  • 2013/11/08/18 Ferrier, H.-M. Lin, Yihua Wu
  • Physics changes:
    1. Minimum effective radius for cloud droplets varies from 10 microns at 20C to 15 microns at 0C
    2. Number concentration of cloud droplets varies from 200 cm^-3 at 20C to 50 cm^-3 at 0C.
    3. Treat all ice as cloud ice in RRTM's progcld2. The effective radius for cloud ice is parameterized following Wyser (1998 J. Climate), in which the smaller cloud ice that is assumed should lead to slightly warmer (improved) temperatures aloft
    4. Build in tiny amounts of cloud within RRTM driver when no cloud is present and RH>95%.
    5. Used the "sep17" version of the modified fer_hires microphysics with RHgrd=98%.
    6. Changes to the JSFC surface layer, increasing ZTmax1 & ZTmax2 from 1.0 to 9.0.
    7. Added the snow cover array (SNOWC, the fraction of the grid covered by snow) that is calculated in the Noah LSM. The maximum fraction is capped at 0.98 (98% coverage) rather than 1.00 (100%) as in the repository code. The SNOWC array is provided as input to the RRTM radiation for the snow cover, where it is calculated following Marshall et al. (1994; =SNOW/[SNOW+70] for SNOW depth in mm).
    8. The fraction of frozen precipitation is partitioned into liquid and frozen contributions for accumulation in the land surface model. The liquid portion goes into the soil, the frozen portion contributes to the snow pack. Rime factor is passed into the Noah LSM and used to increase the density of the incoming frozen precipitation in order to reduce the depth of snow accumulations in marginal winter conditions with complex precipitation types. Also, the upper limit of snow density is increased from 0.4 g/cm^3 to 0.9 g/cm^3.
  • 2013/11/05/18 Wu, Yanqiu Zhu
  • GSI analysis changes:
    1. Turn on new satellite bias correction.
    2. Upgrade GSI to include GFS spectral IO modification which will improve wall time and reduce memory requirement; needed for global ensemble with higher resolution.
    3. Monitor GOES hourly visible winds
  • 2013/09/16/18 Ferrier, H.-M. Lin, Pyle
  • Radiation/microphysics changes:
    1. Increase long wave cloud emissivities by 1.5x
    2. Increase minimum cloud droplet radius from 5 to to 15 microns
  • Bug fix in NPS code : the index value LPT2 that specifies the interface between the terrain following and isobaric portions of the vertical coordinate was being recomputed incorrectly in the output routine, and that incorrect value was written to NPS output. This index value is used within the horizontal diffusion of the NMMB model to determine which levels need to check the coordinate slope to see if the model surface is flat enough to allow for lateral diffusion. The isobaric levels are flat by definition, and thus don't need to be checked. The end result of the bug is that NPS-initialized runs of NMMB have treated all levels as though they are in the isobaric portion of the hybrid coordinate, and thus wrongly have been applying lateral diffusion everywhere.
  • NDASB cycled land-states/satellite bias corrections restarted from NDASX
  • 2013/08/29/18 Wu
  • New version of the GSI analysis (r29871) with the following:
    1. Redo the use of radiosonde significant level data
    2. Changes to use CRTM2.1.3
    3. Assorted bug fixes
    4. Add endian-independent open mechanism
  • 2013/08/20/18 Ferrier
  • Begin use of the reengineered (but scientifically equivalent) version of the gravity wave drag/mounatin blocking code. With this version, cleffamp=2.244 is used which is roughly equivalent to the value of 10 used in the older code.
  • 2013/07/25/12 Rogers
  • Parallel offically moved to WCOSS.
  • 2013/07/11/00 Ferrier
  • Modified this expression (used to calculate the denominator of the scorer parameter below reference level) in the gravity wave drag/mountain blocking code:

    TEM = MAX(VELCO(I,K)*VELCO(I,K), .0001)


    TEM = MAX(VELCO(I,K)*VELCO(I,K), .1)


    SCR = BNV2(I,K) / TEM

    where BNV2 is the Brunt-Vaisala frequency. This change was made to make this expression consistent with the GFS version.

  • 2013/06/17/18 Ferrier
  • Parameter "hpmax" in module_GWD.F90 changed from 200 to 2400 to be consistent with the GFS version of the code. It is used to calculate the Froude Number of the flow (variable FR),

    FR = BNV * ULOI(I) * min(HPRIME(J),hpmax)

    where hprime is the standard deviation of the terrain height, BNV is Brunt-Vaisala frequency, and ULOI is 1/Umean, the mean low-level flow.

  • 2013/06/17/00 Ferrier, Vasic, Rogers
  • Bug fix in module_GWD.F90 : this expression "cleff = 0.5E-5 * SQRT(FLOAT(IMX)/192.0)" changed to "cleff = 0.5E-5 / SQRT(FLOAT(IMX)/192.0)" to be consistent with the GFS version
  • 2013/06/15/00 Ferrier, Wu
  • Turned off raob level enhancement in GSI analysis
  • Increased cleffamp (gravity wave parameterization fudge factor) from 1.0 to 10.0
  • 2013/04/24/00 Ferrier, Rogers
  • Change parameter cdmb in GWDMB scheme from 0.25 to 0.1
  • 2013/04/09/00 Ferrier, Rogers
  • Change parameter cdmb in GWDMB scheme from 0.5 to 0.25
  • 2013/03/26/00 Ferrier, Rogers
  • Change parameter cdmb in GWDMB scheme from 0 to 0.5
  • Change the elevated critical threshold for the onset of condensation (RHgrd) from 0.95 to 0.98
  • Restart cycled land-states and satelllite bias corrections from NDASX
  • 2013/03/14/00 Ferrier, Tassone
  • Restore sigfac value in gravity wave drag scheme/mountain blocking scheme (GWDMB) to 3
  • Set parameter cdmb=0 to turn off mountain blocking in the GWDMB code.
  • Revert back to the cloud fractions before the 3/9/2013 change in terms of the functional dependence of cloud fractions with cloud mixing ratios. But the March 9 change that allows for lower thresholds of nonzero cloud fractions is kept, which allow for cloud fractions as small as 0.1%.
  • Increase z0clim roughness length values for these IGBP vegetation types:
    1. Evergreen Needleleaf Forest
    2. Deciduous Broadleaf Forest
    3. Mixed Forests
    4. Croplands
    5. Cropland/natural vegetation mosaic
  • Restart cycled land-states and satelllite bias corrections from NDASX
  • 2013/03/09/00 Ferrier
  • Turned back on impact of aerosols in the RRTM scheme
  • Change the elevated critical threshold for the onset of condensation (RHgrd) to be 0.95 from 1.0
  • Increased cloud fractions for thin (low mixing ratio) clouds, reduced cloud/no cloud thresholds
  • 2013/03/06/00 Wu
  • Begun ingest of the GDAS ozone field in the GSI analysis for radiance assimilation enhancement
  • 2013/02/25/12 Ferrier
  • Turned off impact of aerosols in the RRTM scheme
  • 2013/02/21/00 Ferrier, Rogers
  • Changed sigfac parameter in GWD/MB scheme from 3.0 to 1.5
  • 2013/02/19/12 Ferrier, Rogers
  • Changed GWD parameter dpmin (height in cb of lower reference level) from 5 to 0 (equivalent to the first model level above ground or the PBL height, whichever is higher)
  • 2013/02/15/18 Janjic
  • Bug fix in gravity wave drag: the critical pressure level (RLOLEV), which has been set to 500 in the model configure file, should be in centibars, so it has been changed to 50.
  • 2013/02/01/00 Y. Wu, Ek
  • Bug fixed in the setup of soil parameters in the Bulk Urban Parameterization in Noah LSM to avoid too high LH fluxes and 2-m dewpoints in urban areas, done by setting SMCDRY=0.4
  • 2013/01/30/00 Ferrier
  • Change to the microphysics:
    1. Decreased NLImax (the maximum number of large ice particles) from 20.e3 (20 per liter) to 10.e3 (20 per liter).
  • 2013/01/29/18 Janjic, Wu
  • Changes to the BMJ convective scheme
    1. An error in the LCL layer in the triggering part of the code was fixed
    2. CAPE is checked instead of entropy in the triggering part of the code
    3. Redundant constant definitions through argument list and use are removed
    4. Weak momentum transport by shallow convection
    5. Combination of old and new shallow convection
  • Changes to the GSI aanlysis
    1. Turn on variational QC scheme in GSI code
    2. Assimilate Meteosat 10 wind subtypes with different data thinning
  • 2012/12/11/00 Wu
  • Turn on use of GOES-15 radiances in the GSI analysis
  • All parallel NDAS cycles restarted with satellite bias correction files spun up with GOES-15 information
  • 2012/12/05/00 Wu
  • GSI analysis upgraded to the same version (r20439) that is running in NAMZ on Zeus
  • 2012/11/21/18 Ferrier, Wu
  • Added subroutine CLTEND from the WRF NMM to the NMMB solver code so that temperature changes from moist cloud processes (convection, grid-scale microphysics) are incrementally applied at each dynamics time step.
  • Begin use of new VAD winds in GSI analysis
  • 2012/11/08/18 Wu
  • GSI analysis changes : use GPS bending angle observations
    1. Replace refractivity with bending angle from satellites GRACE A, COSMIC 1-6, METOP A
    2. Turn on bending angle obs from TSX, SAC-C, C/NOFS
    3. Set cut off height to 30 km (12mb)
  • 2012/10/02/00 Ferrier, H-M. Lin, Hou
  • Added an extra layer above the top of the model domain within the internal RRTM radiation driver for the NMMB (grrad_nmmb.f) to account for the radiative effects of ozone and other trace gases (CO2, N2O, CH4, O2, CO, and various CFCs) from the model top to the top of the atmosphere. A fixed profile with height is assumed for the trace gases.
  • 2012/09/27/18 Ferrier, H-M. Lin
  • Changes to the RRTM radiation
    1. Use global monthly CO2 concentrations that are used in the GFS
    2. Removed diurnal variation of albedo
    3. Code changed to add in small amounts of convective clouds into the input cloud calculations, but only if all of the following conditions are met:
      1. The maximum condensate mixing ratio is < QWmax.
      2. Only shallow convection is present, do not apply to deep convection.
      3. Only apply if the depth of shallow convection is between CU_DEEP_MIN (50 hPa) and CU_DEEP_MAX (200 hPa).
      4. Convective precipitation rate must be <0.01 mm/h.
  • Restarted NDASB land states and satellite bias corrections from NDASX
  • 2012/07/24/18 Ferrier, H-M. Lin
  • Begin use of July 3rd version of RRTM
  • Installed a fix for the calculation of the time-averaged mean of the cosine of the zenith angle in the RRTM radiation scheme.
  • Restarted NDASB land states and satellite bias corrections from NDASX
  • 2012/07/05/12 Janjic, Rogers
  • Reduce resolution factor for DSP's ("fres') from 0.85 to 0.75 in BMJ_DEV convection
  • Restarted NDASB cycle from NDASX soil states and NDASB satellite bias corrections
  • 2012/06/07/00 Wu, Rogers
  • Begin use of EnKF version of the GSI analysis, using ensemble forecasts from the global EnKF system to compute background error covariances.
  • Restarted NDASB soil states and satellite bias corrections from NDASX.
  • 2012/05/22/18 Wu
  • Begin use of mesonet wind observation reject list from RTMA.
  • 2012/04/24/00 Ferrier, H-M. Lin, Hou
  • Revert back to 28 March version of cloud fraction formulation in RRTM radiation (radiation_clouds.f routine).
  • 2012/04/13/18 Ferrier, H-M. Lin, Hou
  • Retuned version of cloud fraction formulation in RRTM radiation with a modified radiation_clouds.f routine put into NAMB (April 13th version).
  • 2012/04/12/18 Ferrier, H-M. Lin, Hou
  • Revert back to NAM cloud fraction formulation in RRTM radiation with a modified radiation_clouds.f routine (April 11th version).
  • 2012/04/06/18 Pyle
  • Changes to the NPS preprocessing code:
    1. When GFS spectral files are used (as is the case in starting the NDAS), the specific humidity rather than the RH field now is used as the moisture source.
    2. The RH processing is modified to prevent setting a floor value of 1% in the RH --> mixing ratio conversion routine. Now a floor value of 0.001% is utilized.
  • 2012/04/05/00 Y. Lin, Rogers
  • Discovered that the NDASB was not assimilating precipitation due to the precip assimilation flag being hardwired to false in the NDASB forecast script. This is the reason the long-term budget adjustment file for NDASB was corrupted, because the input NDASB land-sfc precip accumulation (LSPA) was undefined.
  • Because of the previous problem, the cycled land-states for NDASB were restarted from operational NDAS.
  • Starting with the ops NDAS long-term budget adjustment file (which has not been updated since 1/23/2012 due to the termination of the CPC 1/8th degree precip analysis), the NDASB long-term precip budget adjustment was turned back on, now using the CCPA product.
  • 2012/03/31/00 Ferrier, H-M. Lin, Hou
  • Begin use of a newer modified version of radiation_clouds.f routine (March 28th version).
  • 2012/03/28/12 Y. Lin
  • Turned off long-term budget adjustment for the precipitation analyses used to update soil states in the NDAS forecast due to issues with the now defunct 1/8 deg CPC precip analysis.
  • 2012/03/27/18 Ferrier
  • Pulled out modified version of radiation_clouds.f routine (from Mar 19) and revert back to repository code; removed smoothing of vertical profiles of heating/cooling rates from the 3/24/2012 changes.
  • 2012/03/24/00 Ferrier, H-M. Lin, Hou
  • The following changes were made to the RRTM radiation:
    1. Use the dynamic albedos from the land surface model
    2. Smoothing of vertical profiles of heating/cooling rates
    3. Uses cloud fractions implemented in October 2011 in the NAM (see here for a description of the change)
    4. Correctly separate cloud ice and snow mixing ratios from the input total ice mixing ratios
    5. Calculates cloud-layer (e.g., low/middle/high) and broadband radiation from the prognostic cloud code rather than from the old NAM/GFDL driver (19 Mar version of radiation_clouds.f)
  • 2012/02/02/18 Ferrier, H-M. Lin, Hou
  • Switched from GFDL to RRTM radiation physics.
  • 2012/01/27/18 Wu
  • Begin use of raob level enhancement in GSI analysis.
  • 2011/12/20/00 Ferrier, Rogers
  • In the gravity wave drag parameterization, changed the parameter SIGFAC from 0 to 3. SIGFAC scales the effects of subgrid-scale orography (SIGFAC*standard deviation of the terrain height).
  • 2011/11/22/12 Ferrier
  • Changes to the microphysics:
    1. Decreased NLImax (the maximum number of large ice particles) from 50.e3 (50 per liter) to 20.e3 (20 per liter).
    2. Fixed a bug in the algorithm that calculates the number of large ice particles and the mean ice-particle diameter. It will also alter the rime factor (ice particle density) under special, rare circumstances when the mean diameter and number concentration of the large ice are set to maximum values.
  • 2011/11/07/18 Ferrier, Rogers
  • Removed 4x increase in diffusion of specific humidity and cloud water mixing ratio
  • Restarted cycled NDAS from ops NAM land-states and satellite bias corrections
  • 2011/10/24/18 Rogers
  • Restarted cycled NDAS from ops NAM land-states and satellite bias corrections
  • 2011/10/04/18 Janjic, Pyle
  • Restored smag2 (Smagorinsky constant for 2nd order diffusion) to the original value of 0.4
  • Restored codamp (divergence damping constant) to the original value of 9.0
  • Increased resolution factor for DSP's ("fres') increased to 0.85
  • NDAS cycle (soil states and satellite bias corrections) restarted from NAMX control
  • 2011/09/20/12 Janjic, Pyle, Black
  • Prior to outputting the full non-hydrostatic vertical velocity, add back the barotropic (external mode) component, which was removed some time ago to alleviate a noise problem in the CONUS nest domain
  • 2011/09/15/00 Pyle
  • Reduced smag2 (Smagorinsky constant for 2nd order diffusion) from 0.4 to 0.1
  • Reduced codamp (divergence damping constant) from 9.0 to 3.0
  • 2011/09/02/00 Rogers, Yanqiu Zhu
  • Changed GSI analysis scripts to use operational NDAS/NAM prepbufr files, they were instead using test files from Dennis Keyser with the NRL Aircraft QC applied
  • NDAS cycle (soil states and satellite bias corrections) restarted from NAMX control
  • 2011/08/29/18 Janjic
  • Resolution factor for DSP's ("fres') increased from 0.5 to 0.75
  • NDAS cycle (soil states and satellite bias corrections) restarted from NAMX control
  • 2011/08/08/12 Ferrier, Janjic
  • Start of new NMMB parallel with the following changes to the convective parameterization
    1. Run with experimental BMJ scheme ("bmj_dev")
    2. Reduce resolution factor for DSP's ("fres') from 1.0 to 0.5