program read_ssmi
c
c	Peter M. Woiceshyn	August 22, 2002
c       818-354-1450
c	Peter.M.Woiceshyn@jpl.nasa.gov
c
c   Sample program provides read statements to read in SSM/I data
c		of wind speed ( meters per second) as determined from
C               the F. Wentz SSMI Alogorithms.
c
c   The wind speed data are on a 1440 X 720 element grid with the origin at
c      -89.875 South latitude and 0.125 degrees East longitude (i.e.,
C      a one quarter by one quarter degree lat/lon grid resolution ranging from
C      -89.875 (S) to 89.875 (N) in latitude & from 0.125 to 359.875
C      in East longitude).  A num = 1 indicates good data and
C      a num = 0 indicates no data in the num array.
c
c   There are 12 months of data for the year 2000.
c
c   Note that the data for each month corresponds
c        to a calender month.
c
c  Declarations:

	REAL del_deg, ssmi_spd(1440,720)
	REAL lon(1440), lat(720)
C	
	INTEGER*2    lon2, lat2, num(1440,720)
c
	
	CHARACTER*3  c_mon(12)
	CHARACTER*2  yc, c_mo(12)
	CHARACTER*20 f_in
	CHARACTER*5  f_label
	CHARACTER*4  suff, f_lab

C----------------------------------------------------------------------C

	DATA   c_mon/'Jan','Feb','Mar','Apr','May','Jun',
     &	              'Jul','Aug','Sep','Oct','Nov','Dec'/

	DATA   c_mo/'01','02','03','04','05','06',
     &	            '07','08','09','10','11','12'/


C----------------------------------------------------------------------C

	     yc = '00'
	
	   suff = '.dat'

        f_label = c_mon(1)//yc  	! Jan00
          f_lab = yc//c_mo(1)	!  0001
  
           f_in = 'Ssmi'//f_lab//suff

C----------------------------------------------------------------------C

	     lon2 = 1440
	     lat2 =  720
  
	  del_deg = 0.25
  
	   PRINT *,' ' 
	   PRINT *,'    longitude array dimension =  ', lon2
	   PRINT *,'     latitude array dimension =  ', lat2
	   PRINT *,'            delta lat/lon deg = ', del_deg
	   PRINT *,' ' 

C  Construct lat & lon arrays

	DO jj = 0, lat2 - 1
	
                   lat(jj+1) = jj*del_deg - 90.0 + del_deg/2.0
	ENDdo

	DO ii = 0, lon2 - 1
	
                   lon(ii+1) = ii*del_deg + del_deg/2.0
	ENDdo
  
c
c
c	Lat range:	-89.875(S) to  89.875(N) degrees 
c
c	Lon range:	  0.125(E) to 359.875(E) degrees 
c
c
	   PRINT *,' ' 
	   PRINT *,'    Reading file: ',f_in
	   PRINT *,' ' 

	       OPEN(22,file=f_in,form='unformatted')

		      read(22) ssmi_spd
		      read(22) num

	       CLOSE(22)
	       
	DO j = 256, 257
	     DO i = 1, 5

	  WRITE (6,'(A5,I3,2x,A4,I3,2x,A6,F6.2,2X,A6,F5.2,2X,
     &               A6,F5.2,2x,A6,I1)')
     &       ' i = ', i,'j = ', j,
     &       'lat = ', lat(j),'lon = ', lon(i),
     &          'spd = ',ssmi_spd(i,j), 'num = ',num(i,j)
	       
	     ENDDO
	        WRITE (6,'(A2)') ' ' 
	ENDDO
	        WRITE (6,'(A2)') ' ' 
C----------------------------------------------------------------------C
	       
	STOP

	END

C----------------------------------------------------------------------C
C
C  Output from executing this program:
C
C host{pw}432> read_ssmi.x
C 
C     longitude array dimension = 1440
C      latitude array dimension =  720
C             delta lat/lon deg =    0.25
C  
C     Reading file: Ssmi0001.dat        
C  
C i =   1  j = 256  lat = -26.12  lon =  0.12  spd =  6.40  num =  1
C i =   2  j = 256  lat = -26.12  lon =  0.38  spd =  6.45  num =  1
C i =   3  j = 256  lat = -26.12  lon =  0.62  spd =  6.49  num =  1
C i =   4  j = 256  lat = -26.12  lon =  0.88  spd =  6.54  num =  1
C i =   5  j = 256  lat = -26.12  lon =  1.12  spd =  6.60  num =  1
  
C i =   1  j = 257  lat = -25.88  lon =  0.12  spd =  6.44  num =  1
C i =   2  j = 257  lat = -25.88  lon =  0.38  spd =  6.49  num =  1
C i =   3  j = 257  lat = -25.88  lon =  0.62  spd =  6.51  num =  1
C i =   4  j = 257  lat = -25.88  lon =  0.88  spd =  6.60  num =  1
C i =   5  j = 257  lat = -25.88  lon =  1.12  spd =  6.74  num =  1
C 
C host{pw}433> 
C
C----------------------------------------------------------------------C