Land-Surface Models in MM5 Version 3
How to run OSU LSM
The introduction of the Oregon State University / NCEP Eta Land-Surface
Model (LSM) in MM5 (Chen and Dudhia 1999) has required that the
modeling system be enhanced to
allow for additional inputs and outputs.
LSM Requirements in Model Pre-Processing System
To operate with the LSM, MM5 requires several additional input
fields. The Version 3 Terrain program provides an annual-mean surface
temperature adjusted to model terrain
elevation, a monthly climatological vegetation fraction, dominant
soil type, and dominant vegetation type in each grid cell.
The Regrid program provides soil moisture and temperature at various
depths, water-equivalent snow depth, sea ice, and optionally canopy
moisture. These additional fields are currently available from the
NCEP/NCAR Reanalysis, and the Eta AWIP analyses (for the US only).
These input data are illustrated in Fig. D-1.
How to Set Program Switches to Run LSM
Terrain -
In terrain.deck, set namelist variable
LSMDATA = .T.,
and select values 1 or 2 for VEGTYPE. e.g.,
VEGTYPE = 1,
In V3.3 or earlier terrain.deck, select
set NewLandUseOnly = FALSE
set LandSurface = TRUE
and either
set VegType = USGS
or
set VegType = SiB
The rest of the terrain deck is the same.
Note that SiB data is only available for North America, but these
categories corresepond to those used in the Eta model's operational
LSM. However, either set can be used in MM5's version of the LSM.
Also, note that the SiB classification lacks an `urban' category.

These setups will make use of global datasets to create the following
additional fields on the model grid:
1. VEGFRCnn (nn=1,12): vegetation fraction monthly climatology
2. TEMPGRD: annual mean ground temp adjusted to model terrain elevation
3. SOILINDX: soil type (currently only 30" over US, 1 deg elsewhere)
REGRID -
The datasets that have required additional fields to run LSM in
MM5 are either Eta (AWIP or Eta212 grid) or NNRP data.
To get the NNRP data from NCAR archive, use either get_nnrp.deck
for batch Cray job, or get_nnrp.csh for running interactively from
pregrid/nnrp directory.
There isn't a deck or script to get NCAR archived AWIP data. But
you can do the following to get the data (Note the Eta dataset only
covers US. It starts May 1995, and may have missing periods):
1. Use Web browser to go to
ftp://ncardata.ucar.edu/datasets/ds609.2/inventories/eta.inv
and find out which Gxxxxx file contains the time period of your
interest. Should get both 3Danal and SFanal files.
2. On NCAR's Cray, type the following to get the dataset in non-cos-blocked
format:
msread -fBI Gxxxxx /DSS/Gxxxxx
(You may then ftp this file back to your local workstation to do
the rest. Note though these files are very BIG in size: the 3Danal
file is about 1Gb each, and SFanal is about 250 mb. If ftping big
file is a problem, do steps 3 and 4 on Cray, and ftp the files after
step 4.)
3. Type the following to obtain all file listed in the Gxxxxx file:
tar tvf Gxxxxx > tar.list
Or find out the files from
ftp://ncardata.ucar.edu/datasets/ds609.2/inventories/TARLIST
and click on the appropriate tarlist file.
4. Extract the tmXX files to use by typing, using G40001 (containing
upperair data) and G40006 (containing surface data) files as an
example:
tar -xvf G40001 9706_3Danal/97062400.AWIP3D00.tm00
and
tar -xvf G40006 9706_SFanal/97062400.AWIPSF00.tm00
Repeat the last two commands several times to obtain all time periods.
The tmXX files are to be used by pregrid program. The extracted
file for each time period is considerably smaller in file size:
about 5 Mb each for upperair data and 1.2 Mb for surface data. You
can ftp each file back to your workstation, or tar them up and then
ftp the file back. Note, if you would like to run REGRID on Cray
using AWIP dataset, you can either run the job interactively if
your domain size is not too big (Cray allows 32 Mw of memory for
interactive job only), or you need to modify the deck to get your
extracted files.
In pregrid.csh make sure you have either Eta (AWIP or Eta212 grid)
or NNRP data for the relevant dates and set SRCSOIL and either
set VTSOIL = ../grib.misc/Vtable.AWIPSOIL
or
set VTSOIL = ../grib.misc/Vtable.NNRPSOIL
For snow cover you will also want SRCSNOW and Vtable.AWIPSNOW, or
Vtable.NNRPSNOW for VTSNOW.
Eta data is only available over the US for recent years and has
40 km resolution. NNRP is available globally for the past 40 years
or more, but has about 2.5 deg resolution.
If you have other LSM data, check the Vtables above to see which
fields are required. The fields added to the standard meteorological
fields by setting SRCSOIL (and SRCSNOW) are:
1. SOILT010, SOILT200, SOILT400: Soil temp at various depths (10,
200, and 400 cm; unit: K)
2. SOILM010, SOILM200: Soil moisture at various depths (10 and 200
cm; unit: fraction)
3. SOILHGT: Analysis surface elevation (unit: m)
4. SEAICE: Sea-ice mask (optional; 0 or 1)
5. WEASD: Water-equivalent snow depth (optional; unit: kg m{-2})
If your LSM data contain fields at levels other than those listed,
you can still use them. For example, you may get soil temperature
at 10 and 200 cm only, or perhaps at other levels, you can modify
Vtable to extract these fields. Later when you run MM5, you can
use the new namelist options ISTLYR and ISMLYR (available in 3-2)
to define where your data are. The current default prediction levels
for soil temperature and moisture in MM5/LSM are bounded between
surface and 300 cm below the surface, and there are 4 levels for
each variable.
MM5 -
To use the LSM, set ISOIL=2 in configure.user prior to compilation.
IBLTYP=5 (the MRF PBL) must be used for now. It is the only one
coupled to the LSM.
To use input soil temperature and moisture at levels other than
the default, one can add namelist variables ISTLYR and ISMLYR in
LPARAM section to specify these levels. For example,
ISTLYR = 10,200,0,0,
ISMLYR = 10,200,0,0,
This shows that the input soil temperature only comes in at 10 and
200 cm levels. Note that one can only input up to 4 levels of soil
temperature and moisture. The prediction levels in current MM5/LSM
are 5, 25, 70, and 150 cm. The climatological deep soil temperature
generated in program Terrain is used as the lower boundary condition
for LSM, while open-boundary condition is used for soil moisture.
Additional LSM outputs from MM5 are:
1. SOIL Tn (n=1,4): The soil temp at all 4 soil levels.
2. SOIL Mn (n=1,4): The soil moisture at all 4 soil levels.
3. CANOPYM: Canopy moisture.
4. SNODPTH: Snow depth.
5. SFCRNOFF: Surface runoff accumulation.
6. UGDRNOFF: Underground runoff accumulation.
Note 1: If using NNRP soil moisture, you may want to check
init.F and initnest.F, where there is a correction for some known
biases. Comment GOTO 1001, and GOTO 2001 in
those routines respectively, if you want to use this correction.
Note 2: If you want diagnostic LSM prints at a grid point,
change the line in surfce.F (SURFCE.241) to NOOUT=0 and edit the
IF statements determining where and how often to output data to
unit 29.
|