p2sig: conversion from pressure level to sigma level

p2sig converts pressure leveled grib data to sigma leveled spectral coefficients data. Usage is;
$ runscr/p2sig input(pressure) output(sigma) YYYY:MM:DD:HH:FF

Updates

  • Added ISFCP=2 option for using mean sea level pressure. (2008/11/12)
  • Added IINTP for interpolation option. (2008/11/12)

How to install

p2sig will be installed with the latest version of GSM and RSM. If you want to install it separately,
  1. $ cd gsm (or rsm) ; cvs update src/p2sig ; cvs update src/mtn ; cvs update configure-model
  2. $ configure-model (make Makefile)
  3. $ cd gsm/src/p2sig ; make
  4. $ cd gsm/src/mtn ; make
  5. $ cd gruns/runscr ; cvs co p2sig.in p2sig_era.in
  6. $ cd gruns ; configure-scr gsm (or rsm, rsim, etc)

Options (see runscr/p2sig)

  • STDPRS: choices of pressure levels.
  • INDMOIS: 0 Using specific humidity, 1 Using relative humidity.
  • ISFCP: 0 Using geopotential height data (without topography), 1 Using surface pressure and topography data, 2 Using mean sea level pressure only.
  • INST: 0 for NCEP-Reanalysis, 1 for ECMWF Reanalysis (ERA40).
  • IINTP: 0 for horizontal interpolation, 1 for using closest grid (when undefined values are filled in original data).

Notes

  • Only for global data.
  • Input data should be GRIB format and in a single file. (If not, some pre-processes should be needed.)
  • Both gaussian and lon-lat coordinates can be used. But all variables should be the same coordinate and the coordinate information should be correctly written in the input GRIB file.
  • GRIB ID sensitive: Default GRIB ID is set for NCEP convention. In case of using ECMWF data, you should change the codes.
  • Time stamp sensitive: If you want to change the time stamp of GRIB file, use libs/etc/force_grib_date_mon
    • Forecast time should be 00. If not, use libs/etc/force_grib_date_mon to change the time stamp of grib file.

Example

  • Using NCEP reanalysis (pgb) data
    • Download any pgb file (pgb.anl.YYYYMMDDHH.grib) from ECPC Data server.
    • runscr/p2sig pgb.anl.YYYYMMDDHH.grib test.sig YYYY:MM:DD:HH:00
  • Using ERA40 data
    • runscr/p2sig_era ../libs/con/era40_19900309.grib test_era.sig 1990:03:09:00:00