Mesoscale Parallel Experiment Change Log

Experiment Name:  North American NEMS/NMMB test

When Who What
2014/04/19/00 H.-M. Lin, Ferrier
  • Corrected a bug that was introduced when this change from 1/6/2014 was made: "If RH>95% and no condensate is present, then a small amount of subgrid condensate will be added in the form of liquid water no matter what the air temperature is, even if it is very cold when it should be in the form of ice". The bug fix assumes that the subgrid condensate is the form of liquid water at >=-30C and ice at colder temperatures.
  • 2014/03/28/18 H.-M. Lin, Ferrier
  • Modified the RRTM radiation driver to reduce the output cloud cover in association with subgrid clouds, which occur when tiny amounts of cloud are added to the input cloud arrays used by the radiation in marginal conditions where the relative humidity exceeds 95% and no grid-scale condensate is present. This change only affects the output cloud fractions diagnosed from the model and does not change any of the other forecast fields.
  • 2014/03/18/00 Aligo, Ferrier
  • Separate microphysics species advection was turned back on in all nests except Alaska.
  • 2014/03/14/00 Rogers, Ferrier
  • The gravity wave/mountain blocking parameter cleffamp was changed from 2.244 to 1.0
  • 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 Ferrier, Aligo, S. Liu, W.-S. Wu
  • Microphysics changes ("jan19 version")
    1. 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.
    2. 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.
    3. 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.
    4. 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.
    5. 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.
  • Fixes to the jan19 microphyiscs changes:
    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.
  • 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
  • NDASX land states and satellite bias corrections restarted from NDASB
  • 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 Rogers, Gayno
  • NAMX/NDASX will no longer use the AFWA snow depth analysis. If the IMS snow cover analysis indicates no snow cover, it will remove snow. If IMS snow cover analysis indicates snow cover where the NDASX has none, it will set the snow depth to 5 cm. Otherwise, the snow depth will be cycled.
  • NDASX land states and satellite bias corrections restarted from NDASB
  • 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:

    defc=0.005*(3-l+1.)/3.

    to

    defc=0.015*(3-l+1.)/3.

    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.

  • Changed the Smagorinsky constant for 2nd order diffusion for the NAM nests back to 0.3
  • 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/06/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 McFarquhar & Heymsfield (1996, J. Atmos. Sci.)
    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.
    9. 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).
  • NDASX cycled land-states and satellite bias corrections restarted from NDASB
  • 2014/01/05/00 Rogers
  • Increase Smagorinsky constant for 2nd order diffusion from 0.3 back to the original value to 0.4 for the parallel NAM nests.
  • 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/13/00 Ferrier, Aligo
  • Reduce Smagorinsky constant for 2nd order diffusion from 0.4 to 0.3 for the parallel NAM nests.
  • 2013/11/12/06 Rogers
  • Began running parallel DGEX off of the NAMX 78-h forecast
  • 2013/11/08/18 Ferrier, H.-M. Lin
  • Radiation/microphysics changes:
    1. Increase long wave cloud emissivities by 1.5x
    2. Increase minimum cloud droplet radius from 5 to to 15 microns
  • 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/26/12 Rogers
  • Turn off nest domains for the 12z NAMX run due to poor throughput on the WCOSS development machine.
  • 2013/09/16/18 Pyle
  • 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.
  • 2013/09/16/12 Ferrier
  • Begin use of the reengineered (but scientifically equivalent) version of the gravity wave drag/mountain blocking code. With this version, the following parameter settings are used for domains that run gravity wave drag/mountain blocking:
    1. 12 km North American parent domain : cleffamp=2.244, cdmb=0.1
    2. Alaska 6 km nest domain : cleffamp=1.0, cdmb=0.5
  • 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/06/00 Rogers
  • Increased the divergence damping parameter for the parallel CONUS, Alaska, Hawaii, and Puerto Rico nests from 9 to 12 so they are the same as the settings in production. They were changed in production on 3/13/2012 to address non-fatal, jet streak related stability issues in the CONUS nest in January 2012.
  • 2013/07/30/12 Ferrier
  • 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
  • Increased cleffamp (gravity wave parameterization fudge factor) from 1.0 to 10.0
  • Change parameter cdmb in GWDMB scheme from 0.25 to 0.1
  • 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.

  • 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)

    to

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

    where

    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/07/25/12 Rogers
  • Parallel offically moved to WCOSS; Hawaii, Puerto Rico, and Fire Weather nest added in 00z/12z NAMX runs.
  • 2013/06/15/00 Wu
  • Turned off raob level enhancement in GSI analysis
  • 2013/05/28/18 Ferrier, Janjic, 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)
  • 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/04/24/00 Ferrier, Tassone
  • Change the elevated critical threshold for the onset of condensation (RHgrd) to be 0.98
  • Change in radiation to allow 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
  • 2013/03/13/12 Rogers
  • Changed SIGFAC from 1.5 to 0 the parallel Alaska nest
  • 2013/03/06/00 Wu
  • Begun ingest of the GDAS ozone field in the GSI analysis for radiance assimilation enhancement
  • 2013/02/15/18 Wu
  • Changes to the GSI aanlysis
    1. Turn on variational QC scheme in GSI code
    2. Assimilate Meteosat 10 wind subtypes with different data thinning
  • 2013/02/13/00 Rogers, Janjic
  • Changed SIGFAC from 3 to 1.5 in the parallel Alaska nest
  • 2013/02/10/00 Rogers
  • NAMX forecast now running a CONUS and Alaska nest with these changes:
    1. No parameterized convection in the parallel CONUS nest
    2. SIGFAC=3 in the gravity wave drag formulation in the parallel Alaska nest
  • 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/00 Rogers
  • Use of GOES-15 radiances in the GSI analysis for NAMX was never turned on as indicated in the 12/11/2012 entry below. They will be used in NAMX as of this cycle, and the NDASX cycle will be restarted with bias correction files spun up for GOES-15 from NDASB.
  • 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 Janjic
  • 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
  • 2012/09/27/18 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/14/00 Ferrier
  • Restored aerosol impact in RRTM radiation
  • 2012/08/14/18 Ferrier
  • Removed diurnal variation of albedo in the RRTM radiation
  • 2012/07/31/18 Ferrier
  • Begin use of global monthly CO2 concentrations used in the GFS
  • Turned off aerosol impact in RRTM radiation
  • 2012/07/24/18 Ferrier, H-M. Lin
  • The RRTM code was 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.
  • The July 3rd version of RRTM and the fix for the calculation of the time-averaged mean of the cosine of the zenith angle in the RRTM radiation scheme was installed in the NAMB parallel.
  • 2012/07/09/18 Ferrier
  • Installed a fix for the calculation of the time-averaged mean of the cosine of the zenith angle in the RRTM radiation scheme.
  • 2012/07/05/12 Ferrier, H-M. Lin
  • Begin use of July 3rd version of RRTM
  • Reduction of resolution factor for DSP's ("fres') from 0.85 to 0.75 put into NAMB parallel
  • 2012/06/07/00 Janjic, Wu
  • Use of global EnKF members to compute background error covariances in the regional GSI was put into the NAMB control parallel
  • Reduce resolution factor for DSP's ("fres') from 0.85 to 0.75
  • 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/21/12 Wu
  • First run of NAM parallel with EnKF version of the GSI analysis, using ensemble forecasts from the global EnKF system to compute background error covariances; with NAMB as the control run.
  • NAMX running with all NAMB changes described here.