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.

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

"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 code structure was put in the trunk in June 2010. w1 has been changed back to 0.9 as intended.

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

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

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

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.

1. ACARS humidity
2. "WDSATR" : superob'ed Scatterometer winds over oceans
3. HIRS and AMSU-A radiances from NOAA-19
4. IASI radiances from METOP
5. AMSU-A radiances from AQUA

• IM = 954
• JM = 835
• DX = .126
• DY = .108
• WBD = -60.039
• SBD = -45.036
• PHI0 = 54 (North latitude)
• LAM0 = -106 (West longitude)
• Southwest corner point = -7.491214752N 144.1336823W
• Northeast corner point = 44.53958130N 14.80149937E

1. Changed an aspect of the boundary wind flux correction to match what is done in the NPS. Only the winds that can impact the mass flux in/out of the domain are matched to the source data. On the W/E boundaries, only the U component is touched; on the N/S boundaries, only the V component is touched
2. Changed the limits of inclusion on the source wind data. All input data that falls over PT to PSFC is considered. Previously only data falling between the mid-layer pressures at the top and bottom would be used.

1. Enforce a minimum number concentration for large ice particles (1 per liter) at all temperatures, rather than not enforcing it at >0C in the operational code where a fixed mean diameter is assumed.
2. When initiating ice in water saturated conditions, growth rates are calculated by combining the number concentrations of ice initiation from eq. (2.4) of Meyers et al. (JAM, 1992) with estimated mass growth rates of ice crystals using the calculations of Miller and Young (JAS, 1979). The Miller-Young calculations were in cgs units (g per unit time), so the values in the model need to be multiplied by 1e.-3 to convert to mks units (kg per unit time) in subroutine MY_GROWTH_RATES.
3. A dummy variable representing the mean mass of large ice particles (snow) was instead the first-guess mean diameter of the ice particles. The bug occurs only when the the first-guess number concentrations of large ice particles are outside of a predefined range between NLImin an

1. Land-sfc precipitation accumulation (LSPA) is not yet output from the NDAS forecast, so the operational NDAS long-term precipitation budget is used to adjust the precipitation analysis. LSPA added to output and began recalculaution of long-term precipitation budget at 18z 8/17/2009/
2. The operational NAM 00z 12-36 h precip is used for LSPA outside of the CONUS (OCONUS). Switched to 00z NAMB 12-36 h precip at 12z 9/28/2009.
3. 19 June 2007 NAM Change: 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 (Added this change at 12z 7/30/2009)
4. The NDAS forecast runs with the same amount of divergence damping as the NAM forecast; in operations, WRF-NMM NDAS forecast has 5x higher divergence damping than NAM forecast. (NOTE : NDASB forecast runs w/5x damoing as of 12z 9/28/2009)
5. Two model changes in the 16 December 2008 NAM change package:
   1. The PBL/turbulance schemes were modified to mix each hydrometeor species in the vertical.
   2. To apply vertical diffusion for separate water species, the model was changed so that (a) it can apply vertical diffusion to an arbitrary number of species, (b) the counter gradient option can be applied to some or all of the species if desired, and (c) option to set to zero some or all of the surface fluxes is also made available.