program read_ssmi
c
c	Peter M. Woiceshyn	Feb 6, 1998
c
c   Sample program provides read statements to read in SSM/I F11 data
c		of wind speed ( meters per second) as determined from
C               F. Wentz Alogorithms.
c
c   The wind speed data are on a 1080 X 540 element grid with the origin at
c      -89.83333 South latitude and 0.16667 degrees East longitude (i.e.,
C      a one third by one third degree lat/lon grid resolution ranging from
C      -89.83333 (S) to 89.83333 (N) in latitude & from 0.16667 to 359.833333
C      in East longitude).
c
c   There are 12 months of data for 1995.
c
c   Note that the data for each month corresponds
c        to a calender month.
c
c  Declarations:

	REAL del_deg,ssmi_spd(1080,540)
	REAL lon(1080), lat(540)
C	
	INTEGER*2    lon2, lat2
c
	
	CHARACTER*3  c_mon(12)
	CHARACTER*2  yc
	CHARACTER*20 f_in
	CHARACTER*5  f_label
	CHARACTER*4  suff

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

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


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

	     yc = '95'
	
	   suff = '.dat'

        f_label = c_mon(1)//yc	! Jan95
  
           f_in = 'SsmiF11'//f_label//suff

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

	     lon2 = 1080
	     lat2 =  540
  
	  del_deg = 1.0/3.0
  
	   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.83333(S) to  89.83333(N) degrees 
c
c	Lon range:	  0.16667(E) to 359.83333(E) degrees 
c
c
	   PRINT *,' ' 
	   PRINT *,'    Reading file: ',f_in
	   PRINT *,' ' 

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

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

	  WRITE (6,'(A7,I5,2x,A7,I5,2x,A7,F8.2,2X,A7,F8.2,2X,A7,F8.2)') 
     &       '   i = ', i,'   j = ', j,
     &       ' lat = ', lat(j),' lon = ', lon(i),
     &          ' spd = ',ssmi_spd(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{wd21pw}432> read_ssmi.x
C 
C    longitude array dimension =    1080
C     latitude array dimension =    540
C            delta lat/lon deg =    0.333333
C 
C 
C    Reading file: SsmiF11Jan95.dat    
C 
C   i =     1     j =   256   lat =    -4.83   lon =     0.17   spd =     5.76
C   i =     2     j =   256   lat =    -4.83   lon =     0.50   spd =     6.06
C   i =     3     j =   256   lat =    -4.83   lon =     0.83   spd =     5.83
C   i =     4     j =   256   lat =    -4.83   lon =     1.17   spd =     5.93
C   i =     5     j =   256   lat =    -4.83   lon =     1.50   spd =     6.04
C  
C   i =     1     j =   257   lat =    -4.50   lon =     0.17   spd =     5.83
C   i =     2     j =   257   lat =    -4.50   lon =     0.50   spd =     5.88
C   i =     3     j =   257   lat =    -4.50   lon =     0.83   spd =     5.76
C   i =     4     j =   257   lat =    -4.50   lon =     1.17   spd =     6.03
C   i =     5     j =   257   lat =    -4.50   lon =     1.50   spd =     6.25
C  
C  
C host{wd21pw}433> 
C
C----------------------------------------------------------------------C