binv.F90

Go to the documentation of this file.

 
 
PROGRAM binv
 
  parameter(maxsub=100)
 
  CHARACTER*255 file
  CHARACTER*8   subset
  CHARACTER*8   sub(maxsub)
  integer*8     ninv(3,maxsub)
  real*8        xsub, xmsg
  logical       exist
 
  DATA lunbf  /20/
 
  !-----------------------------------------------------------------------
  nmbyt(lunit)= iupvs01(lunit,'LENM')
  !-----------------------------------------------------------------------
 
  !  get filename
 
  narg=command_argument_count()
  if(narg/=1) then
     print *,'Usage: binv <bufrfile> will print bufrfile inventory by message type'
     call exit(2)
  endif
 
  call get_command_argument(1,file)
  file = trim(file)//char(0)
  inquire(file=file,exist=exist)
  if (.not.exist) then
     print *,trim(file)//' does not exist'
     call exit(3)
  endif
  open(lunbf,file=file,form='unformatted')
 
  ninv = 0
  nsub = 0
 
 
  !  COMPUTE AN MESSAGE INVENTORY BY SUBSETS
  !  ---------------------------------------
 
  CALL openbf(lunbf,'IN',lunbf)
  DO WHILE(ireadmg(lunbf,subset,idate)==0)
     isub = 0
     DO i=1,nsub
        IF(subset==sub(i)) isub = i
     ENDDO
     IF(isub==0) THEN
        IF(nsub+1>maxsub) CALL bort('NSUB TOO BIG')
        sub(nsub+1) = subset
        nsub = nsub+1
        isub = nsub
     ENDIF
     ninv(1,isub) = ninv(1,isub)+1
     ninv(2,isub) = ninv(2,isub)+nmsub(lunbf)
     ninv(3,isub) = ninv(3,isub)+nmbyt(lunbf)
  ENDDO
 
  !  PRINT THE INVENTORY
  !  -------------------
 
  print*
  print'(a4,6x,2(a10,4x),a11)','type','messages','subsets','bytes'
  print*
  DO j=1,nsub
     xmsg = ninv(1,j)
     xsub = ninv(2,j)
     print'(a8,2x,2(i10,4x),i11,4x,f8.2)',sub(j),(ninv(i,j),i=1,3),xsub/xmsg
     IF(j>1) THEN
        ninv(1,1) = ninv(1,1)+ninv(1,j)
        ninv(2,1) = ninv(2,1)+ninv(2,j)
        ninv(3,1) = ninv(3,1)+ninv(3,j)
     ENDIF
  ENDDO
 
  print'(a8,2x,2(i10,4x),i11,4x)','TOTAL   ',(ninv(i,1),i=1,3)
  print*
 
  stop
END PROGRAM binv