c The Gridded S1 and S2 data files are in Unix Fortran standard binary 
c format (i.e., sequential access; note the ocean data are in an ASCII 
c file). Each file contains 5 records (or maps) of seasonal and annual mean 
c values (averaged over the 1976-1997 period).
c The grid is a 4 (lat) x 5 (lon) deg. grid with the 1st box centered at
c (177.5W, 88.0S) and lon increases eastward and lat increases northward.
c Note: the amplitude (in mb) is smoothed while the phase (local solar time 
c in hrs of the peak or the 1st peak for S2) is not. More details are given 
c in Dai and Wang (1999, J. Atmos. Sci. 56, 3874-3891). Sample codes to read 
c these two files:


       parameter(mi=72,mj=45,di=5.0,dj=4.0,aln1=-177.5,alt1=-88.)
       real amplitude(mi,mj),phase(mi,mj)          

       open(1,form='unformatted',file=S1FileNameHere) 
         

       do isea =1, 5   ! 1 --> DJF, 2--> MAM, 3--> JJA, 4--> SON, 5--> ANN
   
       read(1)amplitude,phase    ! ampl in mb unit, phase in hours LST

ccc    do something with the data here.

       end do
         
       end