program main
!
! General routine to read MEDS ASCII
! Structured to allow code to be added to do counts, spit out 
! a special record, or whatever is needed.
!
  implicit none
! include 'gtspp.h'
  character*80 arg,input
  real*8,  parameter :: fildoubl = 9.9692099683868690e+36
  real,    parameter :: filfloat = 9.9692099683868690e+36

  integer :: no_stations
  logical :: if_parseobs
  real*8 :: time, day, ref_date_time
  integer :: yyyy, mm
  integer :: i
  real  :: lat_min,lat_max,lon_min,lon_max,vertical_min,vertical_max
  real*8 :: time_min,time_max

  if_parseobs = .FALSE.
  time = fildoubl
  yyyy = 1900
  mm = 1
  day = 1.0
  call greg2jd(mm,day,yyyy,ref_date_time)
  write(*,*) 'Reference_Data_Time=>',ref_date_time
  if (iargc() .eq. 0) then
    write(*,*)'Invalid command line option(s).'
    write(*,*)'Use readASCII.x input'
    stop
  endif
  do i = 1, iargc()
    call getarg(i,arg)
    select case(i)
      case (1)
        input = trim(arg)
      case (2)
        if (trim(arg) .NE. '') if_parseobs = .TRUE.
      case default
    end select
  end do
  write(*,*) 'Input=',input
  if (if_parseobs) then
    write(*,*) 'Parseobs is on.'
  else
    write(*,*) 'Parseobs is off.'
  end if

  call readASCII(input,if_parseobs,no_stations,lat_min,lat_max,&
    lon_min,lon_max,vertical_min,vertical_max,time_min,time_max)
  write(*,*)'no_stations=',no_stations
  write(*,*)'lat_min=',lat_min
  write(*,*)'lat_max=',lat_max
  write(*,*)'lon_min=',lon_min
  write(*,*)'lon_max=',lon_max
  write(*,*)'time_min=',time_min
  write(*,*)'time_max=',time_max
  write(*,*)'vertical_min=',vertical_min
  write(*,*)'vertical_max=',vertical_max
end program