PROGRAM VERIFY_GMI_v8

c	this program calls the fortran subroutines that read the 
c	daily and time-averaged GMI data files.  It is currently set
c	to write out the data in the verification files

c	remove or comment out sections you do not have files for
	
     
      CHARACTER*100  filename
      REAL*4, DIMENSION(1440,720,7,2):: gmi_data
      REAL*4, DIMENSION(1440,720,6)  :: avgd_data
        

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

c	uncomment and fix the path below - change to match your system
c	filename='your drive:\your directory\f35_20140519vN.N'  !GMI

	CALL READ_GMI_DAY(FILENAME,GMI_DATA,IEXIST)
	if(iexist.ne.0) stop

	write(*,*) 'gmt time, valid data range 0 to 1440 (in minutes)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,1,1)
	write(*,*) ' '
	write(*,*) 'sea surface temperature, valid data range -3 to 34.5 (degree C)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,2,1)
	write(*,*) ' '
	write(*,*) 'wind speed low frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,3,1)
	write(*,*) ' '
	write(*,*) 'wind speed medium frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,4,1)
	write(*,*) ' '
	write(*,*) 'water vapor, valid data range 0 to 75 (millimeters)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,5,1)
	write(*,*) ' '
	write(*,*) 'cloud, valid data range -0.05 to 2.45 (millimeters)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,6,1)
	write(*,*) ' '
	write(*,*) 'rain rate, valid data range 0 to 25 (millimeters/hour)'
	write(*,'(6f11.2)') gmi_data(170:175,274:278,7,1)

c
c
c
	write(*,*) '***********3-day******************'
c	uncomment and fix the path below  change to match your system
c	filename='your drive:\your directory\f35_20140519vN.N_d3d'	!gmi

	CALL READ_GMI_averaged(FILENAME,AVGD_DATA,IEXIST)
	if(iexist.ne.0) stop

	write(*,*) 'sea surface temperature, valid data range -3 to 34.5 (degree C)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,1)
	write(*,*) ' '
	write(*,*) 'wind speed low frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,2)
	write(*,*) ' '
	write(*,*) 'wind speed medium frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,3)
	write(*,*) ' '
	write(*,*) 'water vapor, valid data range 0 to 75 (millimeters)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,4)
	write(*,*) ' '
	write(*,*) 'cloud, valid data range -0.05 to 2.45 (millimeters)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,5)
	write(*,*) ' '
	write(*,*) 'rain rate, valid data range 0 to 25 (millimeters/hour)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,6)


	write(*,*) '***********weekly******************'
c	uncomment and fix the path  below
c	filename='your drive:\your directory\weeks\f35_20140524vN.N'  !GMI
	
	CALL READ_GMI_averaged(FILENAME,AVGD_DATA,IEXIST)
	if(iexist.ne.0) stop

	write(*,*) 'sea surface temperature, valid data range -3 to 34.5 (degree C)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,1)
	write(*,*) ' '
	write(*,*) 'wind speed low frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,2)
	write(*,*) ' '
	write(*,*) 'wind speed medium frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,3)
	write(*,*) ' '
	write(*,*) 'water vapor, valid data range 0 to 75 (millimeters)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,4)
	write(*,*) ' '
	write(*,*) 'cloud, valid data range -0.05 to 2.45 (millimeters)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,5)
	write(*,*) ' '
	write(*,*) 'rain rate, valid data range 0 to 25 (millimeters/hour)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,6)


	write(*,*) '***********month******************'
c	uncomment and fix the path below
c	filename='your drive:\your directory\f35_201405vN.N'  !GMI

	CALL READ_GMI_averaged(FILENAME,AVGD_DATA,IEXIST)
	if(iexist.ne.0) stop

	write(*,*) 'sea surface temperature, valid data range -3 to 34.5 (degree C)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,1)
	write(*,*) ' '
	write(*,*) 'wind speed low frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,2)
	write(*,*) ' '
	write(*,*) 'wind speed medium frequency, valid data range 0 to 50.0 (meters/second)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,3)
	write(*,*) ' '
	write(*,*) 'water vapor, valid data range 0 to 75 (millimeters)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,4)
	write(*,*) ' '
	write(*,*) 'cloud, valid data range -0.05 to 2.45 (millimeters)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,5)
	write(*,*) ' '
	write(*,*) 'rain rate, valid data range 0 to 25 (millimeters/hour)'
	write(*,'(6f11.2)') avgd_data(170:175,274:278,6)

	stop
	end

	include 'read_gmi_subroutines_v8.f'