ungrib/src/ngl/g2/gribmod.f

   1       module grib_mod
   2 !$$$  SUBPROGRAM DOCUMENTATION BLOCK
   3 !                .      .    .                                       .
   4 ! MODULE:    grib_mod
   5 !   PRGMMR: Gilbert         ORG: W/NP11    DATE: 2002-01-23
   6 !
   7 ! ABSTRACT: This Fortran Module contains the declaration
   8 !   of derived type gribfield.
   9 !   If variable gfld is declared of type gribfield
  10 !   ( i.e. TYPE(GRIBFIELD) :: GFLD ), it would have the following componenets:
  11 !
  12 !        gfld%version = GRIB edition number ( currently 2 )
  13 !        gfld%discipline = Message Discipline ( see Code Table 0.0 )
  14 !        gfld%idsect() = Contains the entries in the Identification
  15 !                        Section ( Section 1 )
  16 !                        This element is actually a pointer to an array
  17 !                        that holds the data.
  18 !            gfld%idsect(1)  = Identification of originating Centre
  19 !                                    ( see Common Code Table C-1 )
  20 !                             7 - US National Weather Service
  21 !            gfld%idsect(2)  = Identification of originating Sub-centre
  22 !            gfld%idsect(3)  = GRIB Master Tables Version Number
  23 !                                    ( see Code Table 1.0 )
  24 !                             0 - Experimental
  25 !                             1 - Initial operational version number
  26 !            gfld%idsect(4)  = GRIB Local Tables Version Number
  27 !                                    ( see Code Table 1.1 )
  28 !                             0     - Local tables not used
  29 !                             1-254 - Number of local tables version used
  30 !            gfld%idsect(5)  = Significance of Reference Time (Code Table 1.2)
  31 !                             0 - Analysis
  32 !                             1 - Start of forecast
  33 !                             2 - Verifying time of forecast
  34 !                             3 - Observation time
  35 !            gfld%idsect(6)  = Year ( 4 digits )
  36 !            gfld%idsect(7)  = Month
  37 !            gfld%idsect(8)  = Day
  38 !            gfld%idsect(9)  = Hour
  39 !            gfld%idsect(10)  = Minute
  40 !            gfld%idsect(11)  = Second
  41 !            gfld%idsect(12)  = Production status of processed data
  42 !                                    ( see Code Table 1.3 )
  43 !                              0 - Operational products
  44 !                              1 - Operational test products
  45 !                              2 - Research products
  46 !                              3 - Re-analysis products
  47 !            gfld%idsect(13)  = Type of processed data ( see Code Table 1.4 )
  48 !                              0  - Analysis products
  49 !                              1  - Forecast products
  50 !                              2  - Analysis and forecast products
  51 !                              3  - Control forecast products
  52 !                              4  - Perturbed forecast products
  53 !                              5  - Control and perturbed forecast products
  54 !                              6  - Processed satellite observations
  55 !                              7  - Processed radar observations
  56 !        gfld%idsectlen = Number of elements in gfld%idsect().
  57 !        gfld%local() = Pointer to character array containing contents
  58 !                       of Local Section 2, if included
  59 !        gfld%locallen = length of array gfld%local()
  60 !        gfld%ifldnum = field number within GRIB message
  61 !        gfld%griddef = Source of grid definition (see Code Table 3.0)
  62 !                      0 - Specified in Code table 3.1
  63 !                      1 - Predetermined grid Defined by originating centre
  64 !        gfld%ngrdpts = Number of grid points in the defined grid.
  65 !        gfld%numoct_opt = Number of octets needed for each
  66 !                          additional grid points definition.
  67 !                          Used to define number of
  68 !                          points in each row ( or column ) for
  69 !                          non-regular grids.
  70 !                          = 0, if using regular grid.
  71 !        gfld%interp_opt = Interpretation of list for optional points
  72 !                          definition.  (Code Table 3.11)
  73 !        gfld%igdtnum = Grid Definition Template Number (Code Table 3.1)
  74 !        gfld%igdtmpl() = Contains the data values for the specified Grid
  75 !                         Definition Template ( NN=gfld%igdtnum ).  Each
  76 !                         element of this integer array contains an entry (in
  77 !                         the order specified) of Grid Defintion Template 3.NN
  78 !                         This element is actually a pointer to an array
  79 !                         that holds the data.
  80 !        gfld%igdtlen = Number of elements in gfld%igdtmpl().  i.e. number of
  81 !                       entries in Grid Defintion Template 3.NN
  82 !                       ( NN=gfld%igdtnum ).
  83 !        gfld%list_opt() = (Used if gfld%numoct_opt .ne. 0)  This array
  84 !                          contains the number of grid points contained in
  85 !                          each row ( or column ).  (part of Section 3)
  86 !                          This element is actually a pointer to an array
  87 !                          that holds the data.  This pointer is nullified
  88 !                          if gfld%numoct_opt=0.
  89 !        gfld%num_opt = (Used if gfld%numoct_opt .ne. 0)  The number of entries
  90 !                       in array ideflist.  i.e. number of rows ( or columns )
  91 !                       for which optional grid points are defined.  This value
  92 !                       is set to zero, if gfld%numoct_opt=0.
  93 !        gfdl%ipdtnum = Product Definition Template Number (see Code Table 4.0)
  94 !        gfld%ipdtmpl() = Contains the data values for the specified Product
  95 !                         Definition Template ( N=gfdl%ipdtnum ).  Each element
  96 !                         of this integer array contains an entry (in the
  97 !                         order specified) of Product Defintion Template 4.N.
  98 !                         This element is actually a pointer to an array
  99 !                         that holds the data.
 100 !        gfld%ipdtlen = Number of elements in gfld%ipdtmpl().  i.e. number of
 101 !                       entries in Product Defintion Template 4.N
 102 !                       ( N=gfdl%ipdtnum ).
 103 !        gfld%coord_list() = Real array containing floating point values
 104 !                            intended to document the vertical discretisation
 105 !                            associated to model data on hybrid coordinate
 106 !                            vertical levels.  (part of Section 4)
 107 !                            This element is actually a pointer to an array
 108 !                            that holds the data.
 109 !        gfld%num_coord = number of values in array gfld%coord_list().
 110 !        gfld%ndpts = Number of data points unpacked and returned.
 111 !        gfld%idrtnum = Data Representation Template Number
 112 !                       ( see Code Table 5.0)
 113 !        gfld%idrtmpl() = Contains the data values for the specified Data
 114 !                         Representation Template ( N=gfld%idrtnum ).  Each
 115 !                         element of this integer array contains an entry
 116 !                         (in the order specified) of Product Defintion
 117 !                         Template 5.N.
 118 !                         This element is actually a pointer to an array
 119 !                         that holds the data.
 120 !        gfld%idrtlen = Number of elements in gfld%idrtmpl().  i.e. number
 121 !                       of entries in Data Representation Template 5.N
 122 !                       ( N=gfld%idrtnum ).
 123 !        gfld%unpacked = logical value indicating whether the bitmap and
 124 !                        data values were unpacked.  If false,
 125 !                        gfld%bmap and gfld%fld pointers are nullified.
 126 !        gfld%expanded = Logical value indicating whether the data field
 127 !                         was expanded to the grid in the case where a
 128 !                         bit-map is present.  If true, the data points in
 129 !                         gfld%fld match the grid points and zeros were
 130 !                         inserted at grid points where data was bit-mapped
 131 !                         out.  If false, the data values in gfld%fld were
 132 !                         not expanded to the grid and are just a consecutive
 133 !                         array of data points corresponding to each value of
 134 !                         "1" in gfld%bmap.
 135 !        gfld%ibmap = Bitmap indicator ( see Code Table 6.0 )
 136 !                     0 = bitmap applies and is included in Section 6.
 137 !                     1-253 = Predefined bitmap applies
 138 !                     254 = Previously defined bitmap applies to this field
 139 !                     255 = Bit map does not apply to this product.
 140 !        gfld%bmap() = Logical*1 array containing decoded bitmap,
 141 !                      if ibmap=0 or ibap=254.  Otherwise nullified.
 142 !                      This element is actually a pointer to an array
 143 !                      that holds the data.
 144 !        gfld%fld() = Array of gfld%ndpts unpacked data points.
 145 !                     This element is actually a pointer to an array
 146 !                     that holds the data.
 147 !
 148 !
 149 ! PROGRAM HISTORY LOG:
 150 ! 2002-01-23  Gilbert
 151 ! 2007-04-24  Vuong     -  Added GDT 3.204  Curvilinear Orthogonal Grids
 152 ! 2008-05-29  Vuong     -  Added GDT 3.32768  Rotate Lat/Lon E-grid
 153 ! 2009-02-17  Vuong     -  Allow negative scale factors and limits for
 154 !                          Templates 4.5 and 4.9
 155 ! 2009-12-14  Vuong     -  Fixed bug in routine getidx.f
 156 !                       -  Modified to increase length of seek(512)
 157 !                       -  Added Templates (Satellite Product) 4.31
 158 !                       -  Added Templates (ICAO WAFS) 4.15
 159 !
 160 ! USAGE:    use grib_mod
 161 !
 162 ! ATTRIBUTES:
 163 !   LANGUAGE: Fortran 90
 164 !   MACHINE:  IBM SP
 165 !
 166 !$$$
 167
 168       character(len=12) :: G2_VERSION="g2lib-1.2.2"
 169
 170       type gribfield
 171           integer :: version,discipline
 172           integer,pointer,dimension(:) :: idsect
 173           integer :: idsectlen
 174           character(len=1),pointer,dimension(:) :: local
 175           integer :: locallen
 176           integer :: ifldnum
 177           integer :: griddef,ngrdpts
 178           integer :: numoct_opt,interp_opt,num_opt
 179           integer,pointer,dimension(:) :: list_opt
 180           integer :: igdtnum,igdtlen
 181           integer,pointer,dimension(:) :: igdtmpl
 182           integer :: ipdtnum,ipdtlen
 183           integer,pointer,dimension(:) :: ipdtmpl
 184           integer :: num_coord
 185           real,pointer,dimension(:) :: coord_list
 186           integer :: ndpts,idrtnum,idrtlen
 187           integer,pointer,dimension(:) :: idrtmpl
 188           logical :: unpacked
 189           logical :: expanded
 190           integer :: ibmap
 191           logical*1,pointer,dimension(:) :: bmap
 192           real,pointer,dimension(:) :: fld
 193       end type gribfield
 194
 195       end module