dtg = '2009070100'    ; YYYYMMDDHH  date time group on the model initialization
fcst_hr = '00'        ;  00, 01, 02, 03, 04, 05, 06

imax=240   ; Number of Longitudes
jmax=120   ; Number of Latitudes
kmax=74    ; Number of Pressures

; Set up Lat/Lon grid
dlon = 360.0 / imax
glon = findgen(imax) * dlon
GAUSSGRID, imax, glat

pl = fltarr(imax,jmax,kmax)
phi = fltarr(imax,jmax,kmax)
tt = fltarr(imax,jmax,kmax)
uu = fltarr(imax,jmax,kmax)
vv = fltarr(imax,jmax,kmax)

file = 'fld3d' + dtg + '-0000' + fcst_hr
print, file

openr, 1, file, /swap_if_little_endian
  readu, 1, pl  ; Pressure
  readu, 1, tt  ; Temperature
  readu, 1, phi ; Geop. Height
  readu, 1, uu  ; Zonal Wind
  readu, 1, vv  ; Meridional Wind
;stop
close, 1

phi = phi / 1000.  ; Heights in km

help, pl, tt, phi, uu, vv, glon, glat

END