Mesoscale Parallel Experiment Change Log

Experiment Name:  North American WRF-NMM with various physics and/or analysis changes

When Who What
2007/12/10/18 Treadon, Wu
  • Turned off assimilation of GOES-11 and GOES-12 1x1 sounder radiances. On 4 December NESDIS stopped processing GOES-12 sounding/radiance data (GOES-E). On 5 December they replaced it with GOES-10 in the GOES-E slot. The last GOES-12 (satellite id 256) data was in 18Z 20071204 GDAS dump. When GOES-10 replaced GOES-12 in the GOES-E files (~ 18Z 20071205) it was assigned the GOES-11 satellite id 255 instead of the correct value of 254. As a result all GOES-E data are being processed as though it were GOES-11 data.

    Radiance monitoring plots show noisy increments are being added to the analysis field over GOES-10 domain. Furthermore, since the GOES-10 is being processed by the gsi as GOES-11 data, the GOES-10 data is corrupting the GOES-11 bias correction. Thus, the assimilation of the GOES-11 data is also being adversely impacted. For this reason assimilation of all GOES sounder radiances will be turned off until it is verified that the data are being correctly encoded.

  • 2007/11/28/18 Rogers
  • Due to system maintainence resulting in a 30-h outage on the development computer, NDASY was restarted from GDAS atmospheric states and NDASY land states
  • 2007/11/27/00 Liu, Treadon, Wu
  • In the GSI analysis, fixed an error in how skin temperature variances are applied at water points. The application of the variances over water points in subroutine bkgvar was not correctly updated.
  • 2007/11/16/12 Treadon, Wu
  • Turned off assimilation of NOAA-18 AMSU-A channel 9 due to increased radiometric noise. Change implemented in operations at 18z 11/16
  • 2007/11/13/18 Pyle, Gayno, Wu, Rogers
  • Minor modification to the model orography to ensure consistent smoothing at/near the model's lateral boundary. NDASY was restarted from GDAS atmospheric states and NDASY land states so that the new surface terrain field will be in the cycled model history file.
  • Fixed an error in the GSI analysis script (present since 10/23) which caused the global_angupdate code to fail, which prevented the satbias_angle file from being updated.
  • 2007/11/01/18 Treadon, Wu
  • Turned off AQUA AMSU-A channel 4 due to increased radiometric noise.
  • 2007/10/29/18 Parrish
  • A GSI error involving rotation of wind obs from earth coordinates to NAM grid coordinates was fixed. This only impacts a small number of observations (aircraft and sat winds) that fall close to the date line (180E).
  • 2007/10/24/18 Ferrier, Ek
  • Switched back to using the original version of the unified LSM, but corrected to remove a bug which induced spurious warming at the surface which was triggered by the presence of tiny amounts if ice (likely hail/graupel) reaching the surface.
  • Values in the F_RimeF array are set to 1.0 if no ice is present at a grid point. Before this change, the value of the array could be >1 if rimed ice was present in the previous time step, but then ice was no longer present in the current time step. The values of this array were not reset to 1, so that the "memory" of what happened in earlier times in the forecast affected this array.
  • 2007/10/23/18 Rogers
  • No 12z 10/23 run was made due to dew system maintainence, so the NDASEXP cycle was restarted from GDAS atmospheric states and NDASY land states
  • 2007/10/18/18 Ek, Ferrier, Rogers
  • In the WRF-NMM model, switched to using a version of the unified LSM which "mimics" the ops NMM-LSM (sf_surface_physics=99) land-surface physics module.
  • 2007/10/16/18 Ferrier, Rogers
  • Reset SIGFAC parameter back to zero in the gravity wave drag formulation.
  • 2007/09/26/12 Treadon, Wu
  • Turned off three AQUA AIRS channels (300, 843, 1783) due to increased noise.
  • 2007/09/24/18 Ek
  • Began running with Unified Land-Surface physics package (sf_surface_physics=2 in WRF-NMM namelist)
  • 2007/09/05/18 Ferrier
  • In the gravity wave drag, change SIGFAC, which scales the effects of subgrid-scale orography (SIGFAC*standard deviation of the terrain height) from 0 to 3.
  • 2007/09/04/18 Wu
  • Switch to using latest (August 2007) version of the GSI analysis
  • 2007/08/27/12 Janjic, Gayno, Pyle, Ferrier, Johnson
  • Due to open boundary conditions, the requirement for exact conservation of specific humidity, TKE, and cloud water is relaxed in the advection step. The exact conservation is still required in the antifiltering step following the advection step.
  • Begin use of "3x3" smoothed-desmoothed terrain
  • NDASY restarted from GDAS atmospheric states and NDASY land states
  • 2007/08/06/12 Ferrier, Black, Ek, Gayno, Rogers
  • First run with WRF-NMM V2.2 code with IJK array indexing
  • Revert back to using ops LSM physics (sf_surface_physics=99) instead of unified LSM (sf_surface_physics=2)
  • Begin use of gravity wave drag/mountain blocking. Run with SIGFAC=0: gravity wave drag responds only to ELEVMAX, the maximum elevation within the grid box only. No further inflation is done based on the standard deviation of the 30" heights in the grid box.
  • Use new terrain with ops (1x) smoothing, a smaller (more realistic) Great Salt Lake, a better depiction of the Channel Islands off the California coast, and spurious waterfalls removed.
  • NDASY (IJK version) started from GDAS atmospheric states and NDASY land states on 12z 1 August
  • 2007/07/31/12 Wu
  • Begin use of retuned background errors in the GSI analysis for improved obs fit to first guess.
  • In the GSI analysis' satinfo file, change the use flag for GOES ozone channel 9 from 1 to -1. With a value of 1, makes this channel used in the quality control for all of the other channels and thus, many other channels are rejected. This change will allow for more GOES 1x1 radiances to be used in the analysis
  • 2007/07/23/18 Wu, Rogers
  • Turned off QuikScat data in GSI due to negative impact on the first guess fits.
  • 2007/07/19/12 Wu, Rogers
  • Turn off strong balance constraint in GSI analysis
  • NDASY atmospheric states restarted from the ops NDAS
  • 2007/06/25/12 Wu, Parrish
  • Changes to the GSI analysis:
    1. New GSI version with ijk and ikj I/O and strong balance constraint
    2. use AIRS radiance data
    3. use GOES 1X1 radiance data
    4. QuikScat ocean surface winds
    5. SATWND (eumetsat and modis)
    6. MESONET uv (winds with uselist)
  • 2007/06/19/12 Ek, Wong, Mitchell, Janjic, Ferrier
  • These two changes were implemented into the ops NDAS/NAM on this date:
    1. Change in module_sf_myjsfc: under stable conditions, modify roughness length for heat so that it is only a function of surface-layer bulk Richardson number, removing the dependence on surface elevation. This change will tend to decrease surface heat flux and the 2-m warm bias.
    2. The minimum canopy resistance was increased for evergreen needleleaf forest (doubled from 125.0 to 250.0, units of s/m) and for mixed forest (from 125.0 to 150.0). This change reduced large latent heat fluxes and high 2-m dewpoint temperatures over the northwest in summer 2006 retrospective runs (esp. FVS subregions NWC and NMT). The change will also reduce moisture fluxes and 2-m dewpoint temperatures over the southeast.
  • 2007/06/08/18 Ek, Rogers
  • When the new unified Noah-LSM land surface module was put into this run on 3 May, the change to the roughness length for heat, put in on 24 April, was inadvertently left out and replaced with the ops version of phys/module_sf_myjsfc.F As of this cycle this oversight has been corrected.
  • 2007/05/03/12 Ek
  • Begin use of the new unified Noah-LSM land surface module (sf_surface_physics=2) in place of the operational NMM-LSM module (sf_surface_physics=99)
  • 2007/04/24/18 Ek, Wong, Mitchell, Janjic
  • Change in module_sf_myjsfc: under stable conditions, modify roughness length for heat so that it is only a function of surface-layer bulk Richardson number, removing the dependence on surface elevation. This change will tend to decrease surface heat flux and the 2-m warm bias.
  • 2007/04/09/12 Black
  • Fixed an error in the computation of latitude used for stratopsheric ozone concentrations; the error led to using climatological ozone values valid at the equator at all latitudes.
  • 2007/03/19/18 Ferrier, Ek
  • The minimum canopy resistance was increased for evergreen needleleaf forest (doubled from 125.0 to 250.0, units of s/m) and for mixed forest (from 125.0 to 150.0). This change reduced large latent heat fluxes and high 2-m dewpoint temperatures over the northwest in summer 2006 retrospective runs (esp. FVS subregions NWC and NMT). The change will also reduce moisture fluxes and 2-m dewpoint temperatures over the southeast.
  • 2007/03/16/18 Ferrier, Rogers
  • Remove the upper limit of 5 g/kg of total condensate when calculating parcel buoyancy (see 1/22/07 change).
  • Revert back to using the lateral boundary conditions created for the ops NDAS/NAM forecasts
  • NDASY atmospheric states restarted from the ops NDAS
  • 2007/03/13/00 Rogers
  • Starting with this cycle, the start time of the NAMY was pushed back early enough to guarantee that the 6-h old parallel GFS was used for lateral boundary conditions of the 84-h free forecast. Previously, the start time of the NAMY free forecast was such that the on-time GFS parallel was used.
  • 2007/03/08/12 Pyle, Rogers
  • Create lateral boundary conditions from the parallel GFS (with hybrid vertical coordinate) using the new version of the NAM_MKBND code.
  • 2007/02/02/12 Lin
  • Began use of the 12-36 h forecast precipitation from the 00z ops NAM run to fill in for the CONUS-based Stage II/IV analysis as a driver for NDAS soil moisture outside of the CONUS (OCONUS)
  • 2007/01/31/18 Rogers
  • Started running the tm12 and tm09 portions of the offtime (06z/18z) parallel NDASY to mimic the operational NDAS cycling configuration.
  • NDAS surface update script was changed so that NDASY snow cover/depth, sea-ice, and SST are updated with the latest analyses at the start of the 06z NDASY run instead of the 12z run. This mimics what is done in operations.
  • 2007/01/22/12 Ferrier
  • The December 19 changes to the operational NAM included the effects of full water loading (i.e., no precipitation sedimentation) in calculating parcel buoyancy. This change imposes an assumed upper limit of 5 g/kg of total condensate when calculating parcel buoyancy, and it will have the impact of increasing the intensity and depth of deep (precipitating) parameterized convection.
  • 2006/12/28/00 Ferrier
  • Revert back to ops NAM convective parameterization
  • Remove two modifications to longwave radiation made for the June 2006 WRF-in-NAM implementation
    1. Remove averaging of longwave temperature tendencies from the lowest two model layer
    2. Compute the upward LW at the surface was based on the ground (skin) temperature instead of the average of the skin and lowest model layer temperatures.
  • NDASY cycle restarted from ops NDAS
  • 2006/12/21/00 Ferrier
  • Changes to convective parameterization:
    1. Time scale of convection increased from 40 to 60 min
    2. Increased stability of deep convective profiles (STABDF,STABDS increased from 0.90 to 1.0). This change increases 700-mb temperatures, where the NAM has had a consistent cool bias.
    3. Slight modification to triggering algorithm to be more consistent when treating parcels that are initially saturated.
    4. Further code optimization to speed up execution (added some F90/95 stuff).
    5. Depth of shallow vs. deep convection was reduced from 0.2*Plow to 0.15*Plow (Plow is the midlayer pressure at the lowest level) based on sensitivity experiments.
    6. At swap points (shallow convection when test for deep convection failed) the shallow cloud top is calculated as a function of virtual temperature exceess and when the ambient relative humidity exceeds a threshold value based on a deficit saturation pressure of -40 hPa at cloud base. The deficit saturation pressure that is used is the one used from the failed deep convection calculation (DSPBK).
  • NDASY cycle restarted from ops NDAS
  • Page Last Modified: 12 April 2007