Add a new namelist variable, pmin, to &ungrib to specify the minimum pressure level...
commit5f9c79a6b4f659cccace1bb960259b83ed7eae9f
authorJim Bresch <bresch@ucar.edu>
Tue, 5 Jun 2018 22:31:13 +0000 (5 16:31 -0600)
committerMichael Duda <duda@ucar.edu>
Thu, 7 Jun 2018 21:35:07 +0000 (7 15:35 -0600)
tree7b541be6f85cb78ae8575f73938edea9871db30c
parentd237af17a2aef1be13c46d4d148011a7ee5d5d92
Add a new namelist variable, pmin, to &ungrib to specify the minimum pressure level to process.

Some new NCEP output adds fields as high as 40 Pa. This is too high for typical
wrf applications. Because users will disagree about what is 'too high', they can
now specify their desired minimum pressure level. Setting pmin works for both
grib1 (EC) and grib2 files. Default is set to 100 Pa (1 mb or about 48 km).
For reference 1000 Pa is about 31 km, 40 Pa is 70 km, and 2000 Pa is 27 km.
It is recommended that pmin be equal to or smaller than p_top_requested in the
wrf namelist. Setting pmin=0. is the same as the current behavior which is to process
every pressure level (level code=100) in the grib file.

Using pmin has the added benefit of reducing wall clock time as well as the size
of the intermediate files since stratospheric levels can now be skipped.

Tested on GFS/FV3 and ERA-I files.
namelist.wps.all_options
ungrib/src/rd_grib1.F
ungrib/src/rd_grib2.F
ungrib/src/read_namelist.F
ungrib/src/ungrib.F