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