external/io_grib1/test_grib1_routines.F

   1
   2 PROGRAM test_grib1_routines
   3   IMPLICIT NONE
   4
   5   CHARACTER, DIMENSION(:), POINTER :: fileindex
   6   REAL     , DIMENSION(:), POINTER :: data
   7   INTEGER :: fid
   8   INTEGER :: err
   9   INTEGER :: ret
  10   INTEGER :: size
  11   INTEGER :: index
  12   INTEGER :: istat
  13
  14   INTEGER :: parmid
  15   CHARACTER(40) :: datestr
  16   INTEGER :: leveltype
  17   INTEGER :: level1
  18   INTEGER :: level2
  19   INTEGER :: fcsttime1
  20   INTEGER :: fcsttime2
  21   character(200) :: value
  22   character(20) :: strval
  23   integer :: test
  24   character(50) :: form
  25   integer :: NumberTimes
  26   character (19) :: Time
  27   character (19), pointer  :: Times(:)
  28   integer :: ierr
  29   integer :: numcols, numrows
  30   integer :: center, subcenter, parmtbl
  31
  32
  33   parmid=33
  34   datestr = '2005041412'
  35   leveltype = 119
  36   level1 = 9965
  37   level2 = -HUGE(1)
  38   fcsttime1 = 0
  39   fcsttime2 = -HUGE(1)
  40   center = 250
  41   subcenter = 2
  42   parmtbl = 200
  43
  44   print *,'about to call GET_FILEINDEX_SIZE'
  45   CALL GET_FILEINDEX_SIZE(size)
  46   ALLOCATE(fileindex(1:size), STAT=istat)
  47
  48   print *,'about to call ALLOC_INDEX_FILE'
  49   CALL ALLOC_INDEX_FILE(fileindex)
  50   print *,'about to call OPEN_FILE'
  51   CALL OPEN_FILE('test.grb','r',fid,err)
  52   print *,'about to call INDEX_FILE'
  53   CALL INDEX_FILE(fid,fileindex)
  54
  55   print *,'about to call GET_GRIB_INDEX'
  56   CALL GET_GRIB_INDEX(fileindex, center, subcenter, parmtbl, &
  57        parmid,trim(datestr),leveltype,level1,level2, &
  58        fcsttime1,fcsttime2,index)
  59   print *,'got grib index: ',index
  60
  61   print *,'about to call GET_METADATA_VALUE'
  62   CALL GET_METADATA_VALUE(fileindex, 'GRIB_GRID_ID', "none", "none", &
  63        Value, istat)
  64
  65   print *,'about to call GET_NUM_TIMES'
  66   CALL GET_NUM_TIMES(fileindex, NumberTimes)
  67   print *,'found ',NumberTimes,' times'
  68
  69   ALLOCATE(Times(1:NumberTimes), STAT=ierr)
  70   print *,'about to call GET_TIME'
  71   CALL GET_TIME(fileindex,1,Time)
  72   print *,'Time: ',Time
  73
  74   print *,'about to call GET_SIZEOF_GRID'
  75   CALL GET_SIZEOF_GRID(fileindex,index,numcols,numrows)
  76   allocate(data(1:numcols*numrows))
  77
  78   print *,'about to call READ_GRIB'
  79   CALL READ_GRIB(fileindex,fid,index,data)
  80   print *,'data(20): ',data(20)
  81
  82   deallocate(data)
  83
  84   print *,'about to call FREE_INDEX_FILE'
  85   CALL FREE_INDEX_FILE(fileindex)
  86
  87   print *,'program completed'
  88
  89 END PROGRAM