NCEP reanalysis for TRMM Kwajex domain for the period (01 July - 30 Sept., 1999)
The data are provided here as ascii files so that they can be easily
read on any system. To read the data:
Edit the appropriate FORTRAN program to serve your purposes.
Compile, then run the appropriate program on your machine.
Two programs are provided to read the analysis fields:
read_kwj_up - read kwajex upper-air data (listed below)
read_kwj_sl - read kwajex 1 level grids (listed below)
Description of data:
(1) Upper-Air Data
domain
f(ilon,ilat) ilon=7, jlat=5
lon = 157.5E + (i-1)*2.5 (i=1,ilon) or 157.5E to 172.5E
lat = 2.5N - (j-1)*2.5 (j=1,jlat) or 2.5N to 12.5N
pressure levels in mbs (1000, 925, 850, 700, 600, 500, 400, 300,
250, 200, 150, 100, 70, 50)
z - geopotential height (m) at 14 pressure levels
t - temperature (C) at 14 pressure levels
u - zonal wind component (m/s) at 14 pressure levels
v - meridional wind component (m/s) at 14 pressure levels
q - specific humidity (g/kg) at 14 pressure levels
w - vertical motion (mb/hr) at 12 pressure levels
r - relative humidity (percent) at 8 pressure levels and tropopause
NOTE: q was computed from the relative humidity field which only is given
up to 300 mb, above this level q is computed by assuming that the
dew point (Td) = T - 30 C.
(2) Single-level Data
These fields are provided over the KWAJEX domain on a qaussian
grid. The program to read this data also contains the software to
interpolate the data onto a regular grid (e.g., 1D in latitude and
longitude).
f(ilon,ilat) ilon=9, jlat=8
lon = 157.5E + (i-1)*1.875 (i=1,ilon) or 157.5E to 172.5E
lat is in gaussian latitude going from 0.9524N to 14.2855N
lat(j),j=1,jlat) /0.9524, 2.8571, 4.7618, 6.6666,
8.5713, 10.4760, 12.3808, 14.2855/
p - surface pressure (m)
t - surface temperature (C)
u - surface zonal wind component (m/s)
v - surface meridional wind component (m/s)
q - surface specific humidity (g/kg)
lh - surface latent heat flux (W/m^2)
sh - surface sensible heat flux (W/m^2)
pr - surface rainrate (mm/day)
cf - cloud fraction (%)
lwds - longwave down at sfc (W/m^2)
lwus - longwave up at sfc (W/m^2)
lwut - longwave down at toa (W/m^2)
swds - shortwave down at sfc (W/m^2)
swdt - shortwave down at toa (W/m^2)
swus - shortwave up at sfc (W/m^2)
swut - shortwave up at toa (W/m^2)
rnet - Qr_net (W/m^2) = lwus - lwds + swus - swds + swdt - swut - lwut
Questions should be directed to Paul Ciesielski
(paulc@atmos.colostate.edu).