test_grid_info.F90

1 
2 
3 PROGRAM test
4   IMPLICIT NONE
5 
6   CHARACTER, DIMENSION(:), POINTER :: grib_table_info
7   CHARACTER, DIMENSION(:), POINTER :: grid_info
8   INTEGER :: ret
9   INTEGER :: size
10   INTEGER :: index
11   INTEGER :: istat
12   
13 
14   CALL GET_GRIB1_TABLE_INFO_SIZE(size)
15   ALLOCATE(grib_table_info(1:size), STAT=istat)
16   CALL LOAD_GRIB1_TABLE_INFO("gribmap.txt",grib_table_info,ret)
17   print *,'ret: ',ret
18   print *,'again'
19   CALL GET_GRIB_PARAM (grib_table_info, "TSK", index);
20   print *,'got index: ',index
21 
22   CALL GET_GRID_INFO_SIZE(size)
23   ALLOCATE(grid_info(1:size), STAT=istat)
24   CALL LOAD_GRID_INFO("test", "200509081200", 1, &
25        1.0, 1.0, 1440.0, &
26        0, 240, 1, &
27        450, 250, 45.0, &
28        -100.0, 12.0, 12.0, -101.0, &
29        1, 30.0, &
30        60.0, grib_table_info, &
31        grid_info)
32 
33   print *,'here!'
34   CALL PRINT_GRID_INFO(grid_info)
35   CALL FREE_GRID_INFO(grid_info)
36 
37 print *,'here1'
38 
39 END PROGRAM
40