external/io_grib1/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