site map - email us - search  
 

Back to Example List


Run MM5 with North American Regional Reanalysis (NARR) data.

This instructions are valid for NARR data obtained from NCAR (http://dss.ucar.edu/pub/narr).

When downloading data from this site, one needs to download the 3D, sfc, flx and FIXED files. Note that the FIXED fields have a datestamp of 1979-11-08_00.

pregrid:

Run pregrid twice.

First with the following Vtable (use FIXED input data and date stamp 1979-11-08_00)

GRIB | Level| Level| Level| REGRID  | REGRID | REGRID                           |
Code | Code | 1    | 2    | Name    | Units  | Description                      |
-----+------+------+------+---------+--------+----------------------------------+
81   | 1    | 0    |      | LANDSEA |        | Land/Sea flag                    |
7    | 1    | 0    |      | SOILHGT | m      | Terrain field of source analysis |
224  | 1    | 0    |      | SOILCAT |        | Dominant soil type category      |
225  | 1    | 0    |      | VEGCAT  |        | Dominant land use category       |
-----+------+------+------+---------+--------+----------------------------------+

RENAME the intermediate file NARRSFC

Run pregrid again - this time for your case date, using the 3D, sfc and flx files as input and the following Vtable:

GRIB|Level|Level|Level| REGRID  | REGRID| REGRID                                  |
Code|Code |1    |2    | Name    | Units | Description                             |
----+-----+-----+-----+---------+-------+-----------------------------------------+
11  |100  | *   |     | T       | K     | Temperature                             |
33  |100  | *   |     | U       | m s-1 | U                                       |
34  |100  | *   |     | V       | m s-1 | V                                       |
51  |100  | *   |     | SPECHUMD|       | Specific Humidity                       |
52  |100  | *   |     | RH      |       | Relative Humidity                       |
7   |100  | *   |     | HGT     | m     | Height                                  |
11  |105  | 2   |     | T       | K     | Temperature                             |
52  |105  | 2   |     | RH      |       | Relative Humidity                       |
33  |105  | 10  |     | U       | m s-1 | U                                       |
34  |105  | 10  |     | V       | m s-1 | V                                       |
1   |1    | 0   |     | PSFC    | Pa    | Surface Pressure                        |
2   |102  | 0   |     | PMSL    | Pa    | Sea-level Pressure                      |
65  |1    | 0   |     | WEASD   | kg m-2| Water Equivalent of Accum Snow Depth    |
11  |1    | 0   |     | SKINTEMP| K     | Sea-Surface Temperature                 |
223 |1    | 0   |     | CANOPYM | kg m-2| Plant Canopy Surface Water              |
85  |112  | 0   | 10  | SOILT010| K     | T of 0-10 cm ground layer               |
85  |112  | 10  | 40  | SOILT040| K     | T of 10-40 cm ground layer              |
85  |112  | 40  |100  | SOILT100| K     | T of 40-100 cm ground layer             |
85  |112  | 100 |200  | SOILT200| K     | T of 100-200 cm ground layer            |
144 |112  | 0   | 10  | SOILM010|       | Soil Moisture of 0-10 cm ground layer   |
144 |112  | 10  | 40  | SOILM040|       | Soil Moisture of 10-40 cm ground layer  |
144 |112  | 40  |100  | SOILM100|       | Soil Moisture of 40-100 cm ground layer |
144 |112  | 100 |200  | SOILM200|       | Soil Moisture of 100-200 cm ground layer|
----+-----+-----+-----+---------+-------+-----------------------------------------+

regridder:

MAKE the following CHANGE to the regridder source code:

REGRID/regridder/src/regridder.F, COMMENT the CALL on line 358 OUT. RE-compile.

IF ( ( wind_llflag .EQ. 3 ) .OR. ( wind_llflag .EQ. 5 ) ) THEN
   CALL proc_map_to_met_winds ( wind_xlonc , wind_cone_fac , wind_xlatc )
END IF

By commenting the CALL to proc_map_to_met_winds, out, you will prevent the double rotation to earth coordinates, as NARR data is already in earth coordinates. Just REMEMBER to add this call back in for any other data you are using, else the winds for other datasets will be wrong.

Change the namelist.input file, so that NARRSFC is read in during the execution of this program (MAKE sure the "/" is from after terrain_file_name, to after constants_full_name).

root                          = '../pregrid/FILE'
terrain_file_name             = '../../TERRAIN/TERRAIN_DOMAIN1'
constants_full_name           = '../pregrid/NARRSFC' /

 

 

UCAR Privacy Policy - Terms of Use - ©2004UCAR