readmp.F90

Go to the documentation of this file.

 
 
  program readmp
 
  implicit none
 
  character(255)     :: file
  character(8)       :: subset
  character(1)       :: go
  integer, parameter :: lunit = 20
  integer            :: idate,ireadmg,ireadsb,i4dy
  logical            :: exist
 
! get the filename to open and read
 
  call get_command_argument(1,file); file=trim(adjustl(file))
  if (file == '') then
     print *, 'Usage: readmp <bufrfile> will print reports one at a time'
     call exit(2)
  endif
  inquire(file=file,exist=exist)
  if (.not.exist) then
     print *,trim(file)//' does not exist'
     call exit(3)
  endif
  call get_command_argument(2,go); go=trim(adjustl(go)) ! this for testing
  open(lunit,file=file,form='unformatted')
 
! open the file to bufr and dump the subsets to standard outout one at a time
 
  call openbf(lunit,'IN',lunit)
  do while(ireadmg(lunit,subset,idate)==0)
     do while(ireadsb(lunit)==0)
        print*,'message date=',i4dy(idate)
        call ufdump(lunit,6)
        if(go/='q') read(5,'(a)') go
        if(go=='q') stop
     enddo
  enddo
 
  end program readmp