----------------------
File name: 
----------------------
terra.s30.byyyymm01.eyyyymm01x.psuvt

where byyyymm = beginning year and month (e.g., b200101)
      eyyyymm = ending year and month (e.g., b200102)

The psuvt contains the following instantaneous fields:
PS = Surface pressure (mb), 2-D
U  = U-winds (m/s), 3-D
V  = V-winds (m/s), 3-D
T  = Temperature (K), 3-D

Each file contains every timestep within the beginning month plus the first
timestep of the next month. For example, total timesteps in
terra.s30.b20010101.e20010201x.psuvt is 31(days) x 4(timestep/day) + 1(first
timestep in 20010201) = 125.
      
----------------------
Data format:
----------------------
These are direct access binary files, 32 bits, big endian.

parameter (imx=144, jmx=91, lmx=30)   ! nl is version dependent,30 in GEOS-3
dimension PS(imx,jmx),U(imx,jmx,lmx),V(imx,jmx,lmx),T(imx,jmx,lmx)
do k=1,nstep
 read(unit) nt1,nt2,PS
 read(unit) nt1,nt2,U
 read(unit) nt1,nt2,V
 read(unit) nt1,nt2,T
enddo

Where imx = total number of longitudinal grid (144),
      jmx = total number of latitudinal grid (91),
      lmx = total number of vertical layers (version dependent),
      nt1 = yyyymmdd (year-month-day, e.g., 20010201),
      nt2 = hhmmss (hour-minute-second, e.g., 120000)
      nstep = total timesteps
Note: nt1 and nt2 are long integers.

To read using IDL:

imx=144 & jmx=91 & lmx=30
PS=fltarr(imx,jmx)
U =fltarr(imx,jmx,lmx)
V =fltarr(imx,jmx,lmx)
T =fltarr(imx,jmx,lmx)
head=lonarr(2)
for k=0,nstep-1 do begin
 readu,1,head,PS
 readu,1,head,U
 readu,1,head,V
 readu,1,head,T
endfor

--------------------------
!!!! IMPORTANT !!!!
Since there are some difference in default format of data saving, Fortran
compiler on your workstation may not be the same as that on the Cray where
these files are saved.  Please print out the nt1, nt2 for every timestep to 
be sure that the files are read correctly.

If you use fortran on a SGI machine, please use the following statement when
read the binary file:

open(unit,file='xyz',status='old',access='sequential',form='system')

If you use linux or PC, you have to convert the data to little endian.