yomgrb.F90

Go to the documentation of this file.

MODULE yomgrb

USE parkind1  ,ONLY : jpim     ,jprb
USE yom_ygfl  , ONLY : jpghg, jpgrg , jptrac

IMPLICIT NONE

SAVE

!     ------------------------------------------------------------------
!*    GRIB CODING DESCRIPTORS

! NLOCGRB - ECMWF LOCAL USAGE IDENTIFIER
! NTOTENS - TOTAL NUMBER OF FORCASTS IN ENSEMBLE
! NENSFNB - ENSAMBLE FORECAST NUMBER
! NCALVAL - cal/val number
! NSTREAM - EXPLICIT STREAM NUMBER (OR ZERO FOR DEFAULT STREAMS TO BE USED)
! NSYSTEM - FOR USE IN SEASONAL STREAM (DIFFERENT OPERATIONAL SYSTEMS)
! NMETHOD - FOR USE IN SEASONAL STREAM (DIFFERENT ENSEMBLES)
! NREFERENCE - FOR USE IN HINDCAST STREAM 
! NCONSENSUS - FOR MULTI_ANALYSIS STREAM (1 if consensus)
! NDWD    - FOR MULTI_ANALYSIS STREAM (1 if DWD analysis used)
! NMFR    - FOR MULTI_ANALYSIS STREAM
! NNCEP   - FOR MULTI_ANALYSIS STREAM
! NUKM   - FOR MULTI_ANALYSIS STREAM
! NSTEPLPP(II) - STEP FOR PREVIOUS P.P. OF VARIABLE WITH GRIB CODE II
! NSMAXNP - SUBTRUNCTATION FOR COMPLEX PACKING OF SPECTRAL COEFF
! NBITSSH - NUMBER OF BITS FOR PACKING SPECTRAL COEFF.
! NBITSGG - NUMBER OF BITS FOR PACKING GRID-POINT DATA
! NJDIAG  - SENSITIVITY DIAGNOSTIC NUMBER (1=>J1, 2=>J2, 3=>J3, 4=>J4)
! NJDOMAI - SENSITIVITY DIAGNOSTIC MASK REGION (0=>global, 1=>europe,
!                                               2=>N.H., 3=>S.H.)
! NLEG    - current VAREPS leg number (eg 1(2) for the T399(T255) part of a T399-T255 VAREPS)

! NGRB...  - GRIB codes according to ECMWF local Code Table 2

! NGRBTH   -   3 Potential Temperature
! NGRBALUVP-  15 MODIS albedo UV-vis parallel radiation
! NGRBALUVD-  16 MODIS albedo UV-vis diffuse radiation
! NGRBALNIP-  17 MODIS albedo Near-IR parallel radiation
! NGRBALNID-  18 MODIS albedo Near-IR diffuse radiation
! NGRBSPARC-  20 surface clear-sky PARadiation

! NGRBCVL  -  27 Low vegetation cover
! NGRBCVH  -  28 High vegetation cover
! NGRBTVL  -  29 Low vegetation type
! NGRBTVH  -  30 High vegetation type
! NGRBCI   -  31 Sea ice cover
! NGRBASN  -  32 Snow albedo
! NGRBRSN  -  33 Snow density
! NGRBSST  -  34 Sea surface temperature
! NGRBISTL1-  35 Ice surface temperature layer 1
! NGRBISTL2-  36 Ice surface temperature layer 2
! NGRBISTL3-  37 Ice surface temperature layer 3
! NGRBISTL4-  38 Ice surface temperature layer 4
! NGRBSWL1 -  39 Volumetric soil water content layer 1
! NGRBSWL2 -  40 Volumetric soil water content layer 2
! NGRBSWL3 -  41 Volumetric soil water content layer 3      
! NGRBSWL4 -  42 Volumetric soil water content layer 4 
! NGRBES   -  44 Evaporation of snow
! NGRBSMLT -  45 Snow melt
! NGRB10FG -  49 gust at 10 m level
! NGRBLSPF -  50 large scale precipitation fraction 

! NGRBMONT - 053 Montgomery Geopotential
! NGRBPTHPV- 054 Pressure on Theta and PV surfaces

! NGRBSPAR -  57 surface PARadiation
! NGRBSUVB -  58 surface UV-B radiation
! NGRBCAPE -  59 convect.avail.potential energy
! NGRBPV   -  60 Potential Vorticity

! NGRBSDFOR - 74 standard deviation of a filtered orography
! NGRBSF6  -  210185 SF6 - anthropogenic emissions

! NGRBTCLW -  78 Total column liquid water
! NGRBTCIW -  79 Total column ice water
! NGRBSPD  -  80 !! 80 and 81 extra grib code introduced to 
! NGRBSVD  -  81 !! introduce extra fields for NH

! NGRB082 to NGRB117 reserved for extra fields. Do not use for permanent post-processed fields

! NGRBTCTRAC(JPTRAC)  - 210184 Total column TRAC2: Sulfure hexaFluoride
!                     - 210183 Total column TRAC1: Radon

! Codes for 2D and 3D extra fields
! NGRBMINXTRA  to NGRBMAXXTRA

! NGRBEMIS - 124 Surface Longwave emissivity # has replaced NGRB212

! NGRBAT   - 127 Atmospheric tide
! NGRBBV   - 128 Budget values
! NGRBZ    - 129 Geopotential (at the surface orography)
! NGRBT    - 130 Temperature
! NGRBU    - 131 U-velocity
! NGRBV    - 132 V-velocity
! NGRBQ    - 133 Specific humidity
! NGRBSP   - 134 Surface pressure
! NGRBW    - 135 Vertical velocity
! NGRBTCW  - 136 Total column water
! NGRBTCWV - 137 Total column water vapour
! NGRBVO   - 138 Vorticity (relative)
! NGRBSTL1 - 139 Surface temperature level 1
! NGRBSD   - 141 Snow depth
! NGRBLSP  - 142 Large scale precipitation
! NGRBCP   - 143 Convective precipitation
! NGRBSF   - 144 Snow fall
! NGRBBLD  - 145 Boundary layer dissipation
! NGRBSSHF - 146 Surface sensible heat flux
! NGRBSLHF - 147 Surface latent heat flux
! NGRBCHAR - 148 Charnock parameter 
! NGRB149  - 149 Not used
! NGRB150  - 150 Not used
! NGRBMSL  - 151 Mean sea level pressure
! NGRBLNSP - 152 Log surface pressure
! NGRB153  - 153 Not used
! NGRB154  - 154 Not used
! NGRBD    - 155 Divergence
! NGRBGH   - 156 Height (geopotential)
! NGRBR    - 157 Relative humidity
! NGRBTSP  - 158 Tendency of surface pressure
! NGRBBLH  - 159 Boundary layer height
! NGRBSDOR - 160 Standard deviation of orography
! NGRBISOR - 161 Anisotropy of subgrid scale orography
! NGRBANOR - 162 Angle of subgrid scale orography
! NGRBSLOR - 163 Slope of subgrid scale orography
! NGRBTCC  - 164 Total cloud cover
! NGRB10U  - 165 10 metre u wind
! NGRB10V  - 166 10 metre v wind
! NGRB2T   - 167 2 metre temperature
! NGRB2D   - 168 2 metre dewpoint temperature
! NGRBSSRD - 169 Surface solar radiation downwards
! NGRBSTL2 - 170 Soil temperature level 2
! NGRBLSM  - 172 Land/sea mask
! NGRBSR   - 173 Surface roughness
! NGRBAL   - 174 Albedo
! NGRBSTRD - 175 Surface thermal radiation downwards
! NGRBSSR  - 176 Surface solar radiation
! NGRBSTR  - 177 Surface thermal radiation
! NGRBTSR  - 178 Top solar radiation
! NGRBTTR  - 179 Top thermal radiation
! NGRBEWSS - 180 U-stress
! NGRBNSSS - 181 V-stress
! NGRBE    - 182 Evaporation
! NGRBSTL3 - 183 Soil temperature level 3
! NGRBCCC  - 185 Convective cloud cocer
! NGRBLCC  - 186 Low cloud cover
! NGRBMCC  - 187 Medium cloud cover
! NGRBHCC  - 188 High cloud cover
! NGRBSUND - 189 Sunshine duration
! NGRBEWOV - 190 EW component of sub-grid scale orographic variance
! NGRBNSOV - 191 NS component of sub-grid scale orographic variance
! NGRBNWOV - 192 NWSE component of sub-grid scale orographic variance
! NGRBNEOV - 193 NESW component of sub-grid scale orographic variance
! NGRBTBT  - 194 Brightness temperature (K)
! NGRBLGWS - 195 Latitudinal component of gravity wave stress
! NGRBMGWS - 196 Meridional component of gravity wave stress
! NGRBGWD  - 197 Gravity wave dissipation
! NGRBSRC  - 198 Skin reservoir content
! NGRBVEG  - 199 Percentage of vegetation
! NGRBVSO  - 200 variance of sub-grid scale orogrophy
! NGRBMX2T - 201 Maximum temperature at 2m since last post-processing
! NGRBMN2T - 202 Minimum temperature at 2m since last post-processing
! NGRBO3   - 203 Ozone mixing ratio (EC prognostic ozone)
! NGRBPAW  - 204 Precipitation analysis weights
! NGRBRO   - 205 Runoff
! NGRBTCO3 - 206 Total column ozone 
! NGRB207  - 207 Not used
! NGRBTSRC - 208 Top solar radiation clear sky 
! NGRBTTRC - 209 Top thermal radiation clear sky
! NGRBSSRC - 210 Surface solar radiation clear sky
! NGRBSTRC - 211 Surface thermal radiation clear sky

! NGRBSTINC- 212 TOA incident solar radiation

!-- bunch of codes, confusing in their use ...
! NGRB214  - 214 PROFPROP.RMAX.EA
! NGRB215  - 215 SURFPROP.RMAX.EA
! NGRB216  - 216 RELAPROP.RMAX.EA
! NGRB217  - 217 PROFRESERV.EAU
! NGRB218  - 218 INTSURFTEMPERATU
! NGRB219  - 219 PROFTEMPERATURE

! NGRB221  - 221 Not used
! NGRB222  - 222 Not used
! NGRB223  - 223 Not used
! NGRB224  - 224 used in suafn1
! NGRB225  - 225 used in suafn1
! NGRB226  - 226 Not used
! NGRB227  - 227 used in suafn1
! NGRBTP   - 228 Total precipitation
! NGRBIEWS - 229 Intantaneous X-surface stress
! NGRBINSS - 230 Intantaneous Y-surface stress
! NGRBISHF - 231 Intantaneous surface heat flux
! NGRBIE   - 232 Intantaneous moisture flux (evaporation)
! NGRBLSRH - 234 Logarithm of surface roughness length for heat
! NGRBSKT  - 235 Skin temperature
! NGRBSTL4 - 236 Soil temperature level 4
! NGRBTSN  - 238 Temperature of snow layer
! NGRBCSF  - 239 Convective snow-fall
! NGRBLSF  - 240 Large scale snow-fall
! NGRB241  - 241 Not used
! NGRB242  - 242 Not used
! NGRBFAL  - 243 Forecast albedo
! NGRBFSR  - 244 Forecast surface roughness
! NGRBFLSR - 245 Forecast logarithm of surface roughness for heat
! NGRBCLWC - 246 Cloud liquid water content
! NGRBCIWC - 247 Cloud ice water content
! NGRBCC   - 248 Cloud cover



!-- aerosols in Table 210 -------------------------
! NGRBAERMR01 - 001 aerosol mixing ratio 1
! NGRBAERMR02 - 002 aerosol mixing ratio 2
! NGRBAERMR03 - 003 aerosol mixing ratio 3
! NGRBAERMR04 - 004 aerosol mixing ratio 4
! NGRBAERMR05 - 005 aerosol mixing ratio 5
! NGRBAERMR06 - 006 aerosol mixing ratio 6
! NGRBAERMR07 - 007 aerosol mixing ratio 7
! NGRBAERMR08 - 008 aerosol mixing ratio 8
! NGRBAERMR09 - 009 aerosol mixing ratio 9
! NGRBAERMR10 - 010 aerosol mixing ratio 10
! NGRBAERMR11 - 011 aerosol mixing ratio 11
! NGRBAERMR12 - 012 aerosol mixing ratio 12
! NGRBAERMR13 - 013 aerosol mixing ratio 13
! NGRBAERMR14 - 014 aerosol mixing ratio 14
! NGRBAERMR15 - 015 aerosol mixing ratio 15
! NGRBAERGN01 - 016 aerosol gain acc. 1                         2D
! NGRBAERGN02 - 017 aerosol gain acc. 2                         2D
! NGRBAERGN03 - 018 aerosol gain acc. 3                         2D
! NGRBAERGN04 - 019 aerosol gain acc. 4                         2D
! NGRBAERGN05 - 020 aerosol gain acc. 5                         2D
! NGRBAERGN06 - 021 aerosol gain acc. 6                         2D
! NGRBAERGN07 - 022 aerosol gain acc. 7                         2D
! NGRBAERGN08 - 023 aerosol gain acc. 8                         2D
! NGRBAERGN09 - 024 aerosol gain acc. 9                         2D
! NGRBAERGN10 - 025 aerosol gain acc. 10        2D
! NGRBAERGN11 - 026 aerosol gain acc. 11        2D
! NGRBAERGN12 - 027 aerosol gain acc. 12        2D
! NGRBAERGN13 - 028 aerosol gain acc. 13        2D
! NGRBAERGN14 - 029 aerosol gain acc. 14        2D
! NGRBAERGN15 - 030 aerosol gain acc. 15        2D
! NGRBAERLS01 - 031 aerosol loss acc. 1                         2D
! NGRBAERLS02 - 032 aerosol loss acc. 2                         2D
! NGRBAERLS03 - 033 aerosol loss acc. 3                         2D
! NGRBAERLS04 - 034 aerosol loss acc. 4                         2D
! NGRBAERLS05 - 035 aerosol loss acc. 5                         2D
! NGRBAERLS06 - 036 aerosol loss acc. 6                         2D
! NGRBAERLS07 - 037 aerosol loss acc. 7                         2D
! NGRBAERLS08 - 038 aerosol loss acc. 8                         2D
! NGRBAERLS09 - 039 aerosol loss acc. 9                         2D
! NGRBAERLS10 - 040 aerosol loss acc. 10        2D
! NGRBAERLS11 - 041 aerosol loss acc. 11        2D
! NGRBAERLS12 - 042 aerosol loss acc. 12        2D
! NGRBAERLS13 - 043 aerosol loss acc. 13        2D
! NGRBAERLS14 - 044 aerosol loss acc. 14        2D
! NGRBAERLS15 - 045 aerosol loss acc. 15        2D
! NGRBAERPR   - 046 aerosol precursor mixing ratio 
! NGRBAERSM   - 047 small aerosols mixing ratio 
! NGRBAERLG   - 048 large aerosols mixing ratio 
! NGRBAODPR   - 049 aerosol precursor opt.depth 2D
! NGRBAODSM   - 050 small aerosols opt. depth   2D
! NGRBAODLG   - 051 large aerosols opt. depth   2D
! NGRBAERDEP  - 052 dust emission potential clim2D 
! NGRBAERLTS  - 053 lifting threshold speed clim2D
! NGRBAERSCC  - 054 soli clay content       clim2D

!---------------------------------------------------

! NGRBGHG(JPGHG)   - 210061 GHG1: Carbon dioxide
!                  - 210062 GHG2: Methane
!                  - 210063 GHG3: Nitrous oxide
! NGRBTCGHG(JPGHG) - 210064 Total column GHG1: Carbon Dioxide
!                  - 210065 Total column GHG2: Methane
!                  - 210066 Total column GHG3: Nitrous Oxide

!---------------------------------------------------
! NGRBCO2O         - 210067 CO2 - ocean flux
! NGRBCO2B         - 210068 CO2 - biosphere flux
! NGRBCO2A         - 210069 CO2 - anthropogenic emissions
!
!---------------------------------------------------
! NGRBGRG(JPGRG)   - 210121 GRG1: Nitrogen dioxide
!                  - 210122 GRG2: Sulphur dioxide
!                  - 210123 GRG3: Carbon monoxide
!                  - 210124 GRG4: Formaldehyde

! NGRBTCGRG(JPGRG) - 210125 Total column GRG1: Nitrogen dioxide
!                  - 210126 Total column GRG2: Sulphur dioxide
!                  - 210127 Total column GRG3: Carbon monoxide
!                  - 210128 Total column GRG4: Formaldehyde

!                  - 210203 GRG5: Ozone
!                  - 210206 Total column GRG5: GEMS Ozone

! NGRBTRAC(JPTRAC)  - 210182 TRAC2: SF6
!                   - 210181 TRAC1: Radon


! NGRBSP2, NGRBSP3 - Grib codes for fields in SPA2 and SPA3
! NGRBGP2, NGRBGP3 - Grib codes for grid point fields

INTEGER(KIND=JPIM), PARAMETER :: jpn3sp=10
INTEGER(KIND=JPIM), PARAMETER :: jpn2sp=4

INTEGER(KIND=JPIM) :: nsec0(2)
INTEGER(KIND=JPIM) :: nsec1(2048)
INTEGER(KIND=JPIM) :: nsec2spp(22)
INTEGER(KIND=JPIM) :: nsec2spm(22)
INTEGER(KIND=JPIM) :: nsec3(2)
INTEGER(KIND=JPIM) :: nsec4(42)
INTEGER(KIND=JPIM) :: ngrbs3(0:jpn3sp)
INTEGER(KIND=JPIM) :: ngrbs2(0:jpn2sp)
INTEGER(KIND=JPIM) :: nsteplpp(255)
REAL(KIND=JPRB) :: rsec3(2)

INTEGER(KIND=JPIM),ALLOCATABLE:: nsec2gg(:)
REAL(KIND=JPRB),ALLOCATABLE:: rsec2(:)

INTEGER(KIND=JPIM) :: msec0(2)
INTEGER(KIND=JPIM) :: msec1(2048)
INTEGER(KIND=JPIM) :: msec2spp(22)
INTEGER(KIND=JPIM) :: msec2spm(22)
INTEGER(KIND=JPIM) :: msec3(2)
INTEGER(KIND=JPIM) :: msec4(42)
INTEGER(KIND=JPIM) :: mgrbs3(0:jpn3sp)
INTEGER(KIND=JPIM) :: mgrbs2(0:jpn2sp)
REAL(KIND=JPRB) :: ssec3(2)

INTEGER(KIND=JPIM),ALLOCATABLE:: msec2gg(:)
REAL(KIND=JPRB),ALLOCATABLE:: ssec2(:)

INTEGER(KIND=JPIM), ALLOCATABLE :: ngrbsp2(:), ngrbsp3(:)
INTEGER(KIND=JPIM), ALLOCATABLE :: ngrbgp2(:), ngrbgp3(:)

INTEGER(KIND=JPIM) :: nlocgrb
INTEGER(KIND=JPIM) :: ntotens
INTEGER(KIND=JPIM) :: nensfnb
INTEGER(KIND=JPIM) :: ncalval
INTEGER(KIND=JPIM) :: nstream
INTEGER(KIND=JPIM) :: nsystem
INTEGER(KIND=JPIM) :: nmethod
INTEGER(KIND=JPIM) :: nreference
INTEGER(KIND=JPIM) :: nconsensus
INTEGER(KIND=JPIM) :: ndwd
INTEGER(KIND=JPIM) :: nmfr
INTEGER(KIND=JPIM) :: nncep
INTEGER(KIND=JPIM) :: nukm
INTEGER(KIND=JPIM) :: nsmaxnp
INTEGER(KIND=JPIM) :: nbitssh
INTEGER(KIND=JPIM) :: nbitsgg
INTEGER(KIND=JPIM) :: njdiag
INTEGER(KIND=JPIM) :: njdomai
INTEGER(KIND=JPIM) :: njiter
INTEGER(KIND=JPIM) :: mlocgrb
INTEGER(KIND=JPIM) :: mtotens
INTEGER(KIND=JPIM) :: mensfnb
INTEGER(KIND=JPIM) :: msmaxnp
INTEGER(KIND=JPIM) :: mbitssh
INTEGER(KIND=JPIM) :: mbitsgg
INTEGER(KIND=JPIM) :: mjdiag
INTEGER(KIND=JPIM) :: mjdomai
INTEGER(KIND=JPIM) :: mjiter

INTEGER(KIND=JPIM) :: nleg

INTEGER(KIND=JPIM) :: ngrbth   =  3
INTEGER(KIND=JPIM) :: ngrbaluvp= 15
INTEGER(KIND=JPIM) :: ngrbaluvd= 16
INTEGER(KIND=JPIM) :: ngrbalnip= 17
INTEGER(KIND=JPIM) :: ngrbalnid= 18
INTEGER(KIND=JPIM) :: ngrbsparc= 20
INTEGER(KIND=JPIM) :: ngrb21   = 21
INTEGER(KIND=JPIM) :: ngrb22   = 22
INTEGER(KIND=JPIM) :: ngrb23   = 23
INTEGER(KIND=JPIM) :: ngrbcvl  = 27
INTEGER(KIND=JPIM) :: ngrbcvh  = 28
INTEGER(KIND=JPIM) :: ngrbtvl  = 29
INTEGER(KIND=JPIM) :: ngrbtvh  = 30
INTEGER(KIND=JPIM) :: ngrbci   = 31
INTEGER(KIND=JPIM) :: ngrbasn  = 32
INTEGER(KIND=JPIM) :: ngrbrsn  = 33
INTEGER(KIND=JPIM) :: ngrbsst  = 34
INTEGER(KIND=JPIM) :: ngrbistl1= 35
INTEGER(KIND=JPIM) :: ngrbistl2= 36
INTEGER(KIND=JPIM) :: ngrbistl3= 37
INTEGER(KIND=JPIM) :: ngrbistl4= 38
INTEGER(KIND=JPIM) :: ngrbswl1 = 39
INTEGER(KIND=JPIM) :: ngrbswl2 = 40
INTEGER(KIND=JPIM) :: ngrbswl3 = 41
INTEGER(KIND=JPIM) :: ngrbswl4 = 42
INTEGER(KIND=JPIM) :: ngrbes   = 44
INTEGER(KIND=JPIM) :: ngrbsmlt = 45
INTEGER(KIND=JPIM) :: ngrb10fg = 49
INTEGER(KIND=JPIM) :: ngrblspf = 50

INTEGER(KIND=JPIM) :: ngrbmont = 53
INTEGER(KIND=JPIM) :: ngrbpthpv= 54

INTEGER(KIND=JPIM) :: ngrbspar = 57
INTEGER(KIND=JPIM) :: ngrbsuvb = 58
INTEGER(KIND=JPIM) :: ngrbcape = 59
INTEGER(KIND=JPIM) :: ngrbpv   = 60

INTEGER(KIND=JPIM) :: ngrbsdfor= 74
INTEGER(KIND=JPIM) :: ngrbtclw = 78
INTEGER(KIND=JPIM) :: ngrbtciw = 79
! LARPEGE
INTEGER(KIND=JPIM) :: ngrbspd = 80
INTEGER(KIND=JPIM) :: ngrbsvd = 81
! LECMWF
INTEGER(KIND=JPIM) :: ngrb080 = 80
INTEGER(KIND=JPIM) :: ngrb081 = 81

INTEGER(KIND=JPIM) :: ngrbsf6  = 210185
INTEGER(KIND=JPIM) :: ngrbminxtra  = 082
INTEGER(KIND=JPIM) :: ngrbmaxxtra  = 117

INTEGER(KIND=JPIM) :: ngrb082  = 082
INTEGER(KIND=JPIM) :: ngrb083  = 083
INTEGER(KIND=JPIM) :: ngrb084  = 084
INTEGER(KIND=JPIM) :: ngrb085  = 085
INTEGER(KIND=JPIM) :: ngrb086  = 086
INTEGER(KIND=JPIM) :: ngrb087  = 087
INTEGER(KIND=JPIM) :: ngrb088  = 088
INTEGER(KIND=JPIM) :: ngrb089  = 089
INTEGER(KIND=JPIM) :: ngrb090  = 090
INTEGER(KIND=JPIM) :: ngrb091  = 091
INTEGER(KIND=JPIM) :: ngrb092  = 092
INTEGER(KIND=JPIM) :: ngrb093  = 093
INTEGER(KIND=JPIM) :: ngrb094  = 094
INTEGER(KIND=JPIM) :: ngrb095  = 095
INTEGER(KIND=JPIM) :: ngrb096  = 096
INTEGER(KIND=JPIM) :: ngrb097  = 097
INTEGER(KIND=JPIM) :: ngrb098  = 098
INTEGER(KIND=JPIM) :: ngrb099  = 099
INTEGER(KIND=JPIM) :: ngrb100  = 100
INTEGER(KIND=JPIM) :: ngrb101  = 101
INTEGER(KIND=JPIM) :: ngrb102  = 102
INTEGER(KIND=JPIM) :: ngrb103  = 103
INTEGER(KIND=JPIM) :: ngrb104  = 104
INTEGER(KIND=JPIM) :: ngrb105  = 105
INTEGER(KIND=JPIM) :: ngrb106  = 106
INTEGER(KIND=JPIM) :: ngrb107  = 107
INTEGER(KIND=JPIM) :: ngrb108  = 108
INTEGER(KIND=JPIM) :: ngrb109  = 109
INTEGER(KIND=JPIM) :: ngrb110  = 110
INTEGER(KIND=JPIM) :: ngrb111  = 111
INTEGER(KIND=JPIM) :: ngrb112  = 112
INTEGER(KIND=JPIM) :: ngrb113  = 113
INTEGER(KIND=JPIM) :: ngrb114  = 114
INTEGER(KIND=JPIM) :: ngrb115  = 115
INTEGER(KIND=JPIM) :: ngrb116  = 116
INTEGER(KIND=JPIM) :: ngrb117  = 117

INTEGER(KIND=JPIM) :: ngrb118  = 118
INTEGER(KIND=JPIM) :: ngrb119  = 119
INTEGER(KIND=JPIM) :: ngrb120  = 120

INTEGER(KIND=JPIM) :: ngrbemis = 124

INTEGER(KIND=JPIM) :: ngrbat   = 127
INTEGER(KIND=JPIM) :: ngrbbv   = 128
INTEGER(KIND=JPIM) :: ngrbz    = 129
INTEGER(KIND=JPIM) :: ngrbt    = 130
INTEGER(KIND=JPIM) :: ngrbu    = 131
INTEGER(KIND=JPIM) :: ngrbv    = 132
INTEGER(KIND=JPIM) :: ngrbq    = 133
INTEGER(KIND=JPIM) :: ngrbsp   = 134
INTEGER(KIND=JPIM) :: ngrbw    = 135
INTEGER(KIND=JPIM) :: ngrbtcw  = 136
INTEGER(KIND=JPIM) :: ngrbtcwv = 137
INTEGER(KIND=JPIM) :: ngrbvo   = 138
INTEGER(KIND=JPIM) :: ngrbstl1 = 139
INTEGER(KIND=JPIM) :: ngrbsd   = 141
INTEGER(KIND=JPIM) :: ngrblsp  = 142
INTEGER(KIND=JPIM) :: ngrbcp   = 143
INTEGER(KIND=JPIM) :: ngrbsf   = 144
INTEGER(KIND=JPIM) :: ngrbbld  = 145
INTEGER(KIND=JPIM) :: ngrbsshf = 146
INTEGER(KIND=JPIM) :: ngrbslhf = 147
INTEGER(KIND=JPIM) :: ngrbchar = 148
INTEGER(KIND=JPIM) :: ngrb149  = 149
INTEGER(KIND=JPIM) :: ngrb150  = 150
INTEGER(KIND=JPIM) :: ngrbmsl  = 151
INTEGER(KIND=JPIM) :: ngrblnsp = 152
INTEGER(KIND=JPIM) :: ngrb153  = 153
INTEGER(KIND=JPIM) :: ngrb154  = 154
INTEGER(KIND=JPIM) :: ngrbd    = 155
INTEGER(KIND=JPIM) :: ngrbgh   = 156
INTEGER(KIND=JPIM) :: ngrbr    = 157
INTEGER(KIND=JPIM) :: ngrbtsp  = 158
INTEGER(KIND=JPIM) :: ngrbblh  = 159
INTEGER(KIND=JPIM) :: ngrbsdor = 160
INTEGER(KIND=JPIM) :: ngrbisor = 161
INTEGER(KIND=JPIM) :: ngrbanor = 162
INTEGER(KIND=JPIM) :: ngrbslor = 163
INTEGER(KIND=JPIM) :: ngrbtcc  = 164
INTEGER(KIND=JPIM) :: ngrb10u  = 165
INTEGER(KIND=JPIM) :: ngrb10v  = 166
INTEGER(KIND=JPIM) :: ngrb2t   = 167
INTEGER(KIND=JPIM) :: ngrb2d   = 168
INTEGER(KIND=JPIM) :: ngrbssrd = 169
INTEGER(KIND=JPIM) :: ngrbstl2 = 170
INTEGER(KIND=JPIM) :: ngrblsm  = 172
INTEGER(KIND=JPIM) :: ngrbsr   = 173
INTEGER(KIND=JPIM) :: ngrbal   = 174
INTEGER(KIND=JPIM) :: ngrbstrd = 175
INTEGER(KIND=JPIM) :: ngrbssr  = 176
INTEGER(KIND=JPIM) :: ngrbstr  = 177
INTEGER(KIND=JPIM) :: ngrbtsr  = 178
INTEGER(KIND=JPIM) :: ngrbttr  = 179
INTEGER(KIND=JPIM) :: ngrbewss = 180
INTEGER(KIND=JPIM) :: ngrbnsss = 181
INTEGER(KIND=JPIM) :: ngrbe    = 182
INTEGER(KIND=JPIM) :: ngrbstl3 = 183
INTEGER(KIND=JPIM) :: ngrbccc  = 185
INTEGER(KIND=JPIM) :: ngrblcc  = 186
INTEGER(KIND=JPIM) :: ngrbmcc  = 187
INTEGER(KIND=JPIM) :: ngrbhcc  = 188
INTEGER(KIND=JPIM) :: ngrbsund = 189
INTEGER(KIND=JPIM) :: ngrbewov = 190
INTEGER(KIND=JPIM) :: ngrbnsov = 191
INTEGER(KIND=JPIM) :: ngrbnwov = 192
INTEGER(KIND=JPIM) :: ngrbneov = 193
INTEGER(KIND=JPIM) :: ngrbtbt  = 194
INTEGER(KIND=JPIM) :: ngrblgws = 195
INTEGER(KIND=JPIM) :: ngrbmgws = 196
INTEGER(KIND=JPIM) :: ngrbgwd  = 197
INTEGER(KIND=JPIM) :: ngrbsrc  = 198
INTEGER(KIND=JPIM) :: ngrbveg  = 199
INTEGER(KIND=JPIM) :: ngrbvso  = 200
INTEGER(KIND=JPIM) :: ngrbmx2t = 201
INTEGER(KIND=JPIM) :: ngrbmn2t = 202
INTEGER(KIND=JPIM) :: ngrbo3   = 203
INTEGER(KIND=JPIM) :: ngrbpaw  = 204
INTEGER(KIND=JPIM) :: ngrbro   = 205
INTEGER(KIND=JPIM) :: ngrbtco3 = 206
INTEGER(KIND=JPIM) :: ngrb207  = 207
INTEGER(KIND=JPIM) :: ngrbtsrc = 208
INTEGER(KIND=JPIM) :: ngrbttrc = 209
INTEGER(KIND=JPIM) :: ngrbssrc = 210
INTEGER(KIND=JPIM) :: ngrbstrc = 211

INTEGER(KIND=JPIM) :: ngrb214  = 214
INTEGER(KIND=JPIM) :: ngrb215  = 215
INTEGER(KIND=JPIM) :: ngrb216  = 216
INTEGER(KIND=JPIM) :: ngrb217  = 217
INTEGER(KIND=JPIM) :: ngrb218  = 218
INTEGER(KIND=JPIM) :: ngrb219  = 219

INTEGER(KIND=JPIM) :: ngrbstinc= 212
INTEGER(KIND=JPIM) :: ngrbvimd = 213

INTEGER(KIND=JPIM) :: ngrb222  = 222
INTEGER(KIND=JPIM) :: ngrb223  = 223
INTEGER(KIND=JPIM) :: ngrb224  = 224
INTEGER(KIND=JPIM) :: ngrb225  = 225
INTEGER(KIND=JPIM) :: ngrb226  = 226
INTEGER(KIND=JPIM) :: ngrb227  = 227
INTEGER(KIND=JPIM) :: ngrbtp   = 228
INTEGER(KIND=JPIM) :: ngrbiews = 229
INTEGER(KIND=JPIM) :: ngrbinss = 230
INTEGER(KIND=JPIM) :: ngrbishf = 231
INTEGER(KIND=JPIM) :: ngrbie   = 232
INTEGER(KIND=JPIM) :: ngrblsrh = 234
INTEGER(KIND=JPIM) :: ngrbskt  = 235
INTEGER(KIND=JPIM) :: ngrbstl4 = 236
INTEGER(KIND=JPIM) :: ngrbtsn  = 238
INTEGER(KIND=JPIM) :: ngrbcsf  = 239
INTEGER(KIND=JPIM) :: ngrblsf  = 240
INTEGER(KIND=JPIM) :: ngrb241  = 241
INTEGER(KIND=JPIM) :: ngrb242  = 242
INTEGER(KIND=JPIM) :: ngrbfal  = 243
INTEGER(KIND=JPIM) :: ngrbfsr  = 244
INTEGER(KIND=JPIM) :: ngrbflsr = 245
INTEGER(KIND=JPIM) :: ngrbclwc = 246
INTEGER(KIND=JPIM) :: ngrbciwc = 247
INTEGER(KIND=JPIM) :: ngrbcc   = 248
INTEGER(KIND=JPIM) :: ngrb249  = 249
INTEGER(KIND=JPIM) :: ngrb250  = 250
INTEGER(KIND=JPIM) :: ngrb251  = 251
INTEGER(KIND=JPIM) :: ngrb252  = 252
INTEGER(KIND=JPIM) :: ngrb253  = 253
INTEGER(KIND=JPIM) :: ngrb254  = 254
INTEGER(KIND=JPIM) :: ngrb255  = 255


!-- aerosols -- Table 210 --------------------------
INTEGER(KIND=JPIM) :: ngrbaermr01=210001
INTEGER(KIND=JPIM) :: ngrbaermr02=210002
INTEGER(KIND=JPIM) :: ngrbaermr03=210003
INTEGER(KIND=JPIM) :: ngrbaermr04=210004
INTEGER(KIND=JPIM) :: ngrbaermr05=210005
INTEGER(KIND=JPIM) :: ngrbaermr06=210006
INTEGER(KIND=JPIM) :: ngrbaermr07=210007
INTEGER(KIND=JPIM) :: ngrbaermr08=210008
INTEGER(KIND=JPIM) :: ngrbaermr09=210009
INTEGER(KIND=JPIM) :: ngrbaermr10=210010
INTEGER(KIND=JPIM) :: ngrbaermr11=210011
INTEGER(KIND=JPIM) :: ngrbaermr12=210012
INTEGER(KIND=JPIM) :: ngrbaermr13=210013
INTEGER(KIND=JPIM) :: ngrbaermr14=210014
INTEGER(KIND=JPIM) :: ngrbaermr15=210015

INTEGER(KIND=JPIM) :: ngrbaergn01=210016
INTEGER(KIND=JPIM) :: ngrbaergn02=210017
INTEGER(KIND=JPIM) :: ngrbaergn03=210018
INTEGER(KIND=JPIM) :: ngrbaergn04=210019
INTEGER(KIND=JPIM) :: ngrbaergn05=210020
INTEGER(KIND=JPIM) :: ngrbaergn06=210021
INTEGER(KIND=JPIM) :: ngrbaergn07=210022
INTEGER(KIND=JPIM) :: ngrbaergn08=210023
INTEGER(KIND=JPIM) :: ngrbaergn09=210024
INTEGER(KIND=JPIM) :: ngrbaergn10=210025
INTEGER(KIND=JPIM) :: ngrbaergn11=210026
INTEGER(KIND=JPIM) :: ngrbaergn12=210027
INTEGER(KIND=JPIM) :: ngrbaergn13=210028
INTEGER(KIND=JPIM) :: ngrbaergn14=210029
INTEGER(KIND=JPIM) :: ngrbaergn15=210030

INTEGER(KIND=JPIM) :: ngrbaerls01=210031
INTEGER(KIND=JPIM) :: ngrbaerls02=210032
INTEGER(KIND=JPIM) :: ngrbaerls03=210033
INTEGER(KIND=JPIM) :: ngrbaerls04=210034
INTEGER(KIND=JPIM) :: ngrbaerls05=210035
INTEGER(KIND=JPIM) :: ngrbaerls06=210036
INTEGER(KIND=JPIM) :: ngrbaerls07=210037
INTEGER(KIND=JPIM) :: ngrbaerls08=210038
INTEGER(KIND=JPIM) :: ngrbaerls09=210039
INTEGER(KIND=JPIM) :: ngrbaerls10=210040
INTEGER(KIND=JPIM) :: ngrbaerls11=210041
INTEGER(KIND=JPIM) :: ngrbaerls12=210042
INTEGER(KIND=JPIM) :: ngrbaerls13=210043
INTEGER(KIND=JPIM) :: ngrbaerls14=210044
INTEGER(KIND=JPIM) :: ngrbaerls15=210045

INTEGER(KIND=JPIM) :: ngrbaerpr  =210046
INTEGER(KIND=JPIM) :: ngrbaersm  =210047
INTEGER(KIND=JPIM) :: ngrbaerlg  =210048
INTEGER(KIND=JPIM) :: ngrbaodpr  =210049
INTEGER(KIND=JPIM) :: ngrbaodsm  =210050
INTEGER(KIND=JPIM) :: ngrbaodlg  =210051
INTEGER(KIND=JPIM) :: ngrbaerdep =210052
INTEGER(KIND=JPIM) :: ngrbaerlts =210053
INTEGER(KIND=JPIM) :: ngrbaerscc =210054
!---------------------------------------------------
INTEGER(KIND=JPIM) :: ngrbco2o = 210067
INTEGER(KIND=JPIM) :: ngrbco2b = 210068
INTEGER(KIND=JPIM) :: ngrbco2a = 210069

!---------------------------------------------------

INTEGER(KIND=JPIM), DIMENSION(JPGHG) :: ngrbghg = (/&
 & 210061, 210062, 210063/)
integer(kind=jpim), DIMENSION(JPTRAC) :: ngrbtrac = (/&
 & 210181, 210182/)

INTEGER(KIND=JPIM), DIMENSION(JPGHG) :: ngrbtcghg = (/&
 & 210064, 210065, 210066/)
integer(kind=jpim), DIMENSION(JPTRAC) :: ngrbtctrac = (/&
 & 210183, 210184/)

INTEGER(KIND=JPIM), DIMENSION(JPGRG) :: ngrbgrg = (/&
 & 210121, 210122, 210123, 210124, 210203/)

INTEGER(KIND=JPIM), DIMENSION(JPGRG) :: ngrbtcgrg = (/&
 & 210125, 210126, 210127, 210128, 210206/)

!     ------------------------------------------------------------------
!$OMP THREADPRIVATE(mbitsgg,mbitssh,mensfnb,mgrbs2,mgrbs3,mjdiag,mjdomai,mjiter,mlocgrb,msec0,msec1,msec2spm,msec2spp)
!$OMP THREADPRIVATE(msec3,msec4,msmaxnp,mtotens,nbitsgg,nbitssh,ncalval,nconsensus,ndwd,nensfnb,ngrb080,ngrb081)
!$OMP THREADPRIVATE(ngrb082,ngrb083,ngrb084,ngrb085,ngrb086,ngrb087,ngrb088,ngrb089,ngrb090,ngrb091,ngrb092,ngrb093)
!$OMP THREADPRIVATE(ngrb094,ngrb095,ngrb096,ngrb097,ngrb098,ngrb099,ngrb100,ngrb101,ngrb102,ngrb103,ngrb104,ngrb105)
!$OMP THREADPRIVATE(ngrb106,ngrb107,ngrb108,ngrb109,ngrb10fg,ngrb10u,ngrb10v,ngrb110,ngrb111,ngrb112,ngrb113,ngrb114)
!$OMP THREADPRIVATE(ngrb115,ngrb116,ngrb117,ngrb118,ngrb119,ngrb120,ngrb149,ngrb150,ngrb153,ngrb154,ngrb207,ngrb21)
!$OMP THREADPRIVATE(ngrb214,ngrb215,ngrb216,ngrb217,ngrb218,ngrb219,ngrb22,ngrb222,ngrb223,ngrb224,ngrb225,ngrb226)
!$OMP THREADPRIVATE(ngrb227,ngrb23,ngrb241,ngrb242,ngrb249,ngrb250,ngrb251,ngrb252,ngrb253,ngrb254,ngrb255,ngrb2d)
!$OMP THREADPRIVATE(ngrb2t,ngrbaerdep,ngrbaergn01,ngrbaergn02,ngrbaergn03,ngrbaergn04,ngrbaergn05,ngrbaergn06)
!$OMP THREADPRIVATE(ngrbaergn07,ngrbaergn08,ngrbaergn09,ngrbaergn10,ngrbaergn11,ngrbaergn12,ngrbaergn13,ngrbaergn14)
!$OMP THREADPRIVATE(ngrbaergn15,ngrbaerlg,ngrbaerls01,ngrbaerls02,ngrbaerls03,ngrbaerls04,ngrbaerls05,ngrbaerls06)
!$OMP THREADPRIVATE(ngrbaerls07,ngrbaerls08,ngrbaerls09,ngrbaerls10,ngrbaerls11,ngrbaerls12,ngrbaerls13,ngrbaerls14)
!$OMP THREADPRIVATE(ngrbaerls15,ngrbaerlts,ngrbaermr01,ngrbaermr02,ngrbaermr03,ngrbaermr04,ngrbaermr05,ngrbaermr06)
!$OMP THREADPRIVATE(ngrbaermr07,ngrbaermr08,ngrbaermr09,ngrbaermr10,ngrbaermr11,ngrbaermr12,ngrbaermr13,ngrbaermr14)
!$OMP THREADPRIVATE(ngrbaermr15,ngrbaerpr,ngrbaerscc,ngrbaersm,ngrbal,ngrbalnid,ngrbalnip,ngrbaluvd,ngrbaluvp,ngrbanor)
!$OMP THREADPRIVATE(ngrbaodlg,ngrbaodpr,ngrbaodsm,ngrbasn,ngrbat,ngrbbld,ngrbblh,ngrbbv,ngrbcape,ngrbcc,ngrbccc)
!$OMP THREADPRIVATE(ngrbchar,ngrbci,ngrbciwc,ngrbclwc,ngrbco2a,ngrbco2b,ngrbco2o,ngrbcp,ngrbcsf,ngrbcvh,ngrbcvl,ngrbd)
!$OMP THREADPRIVATE(ngrbe,ngrbemis,ngrbes,ngrbewov,ngrbewss,ngrbfal,ngrbflsr,ngrbfsr,ngrbgh,ngrbghg,ngrbgrg,ngrbgwd)
!$OMP THREADPRIVATE(ngrbhcc,ngrbie,ngrbiews,ngrbinss,ngrbishf,ngrbisor,ngrbistl1,ngrbistl2,ngrbistl3,ngrbistl4,ngrblcc)
!$OMP THREADPRIVATE(ngrblgws,ngrblnsp,ngrblsf,ngrblsm,ngrblsp,ngrblspf,ngrblsrh,ngrbmaxxtra,ngrbmcc,ngrbmgws)
!$OMP THREADPRIVATE(ngrbminxtra,ngrbmn2t,ngrbmont,ngrbmsl,ngrbmx2t,ngrbneov,ngrbnsov,ngrbnsss,ngrbnwov,ngrbo3,ngrbpaw)
!$OMP THREADPRIVATE(ngrbpthpv,ngrbpv,ngrbq,ngrbr,ngrbro,ngrbrsn,ngrbs2,ngrbs3,ngrbsd,ngrbsdfor,ngrbsdor,ngrbsf,ngrbsf6)
!$OMP THREADPRIVATE(ngrbskt,ngrbslhf,ngrbslor,ngrbsmlt,ngrbsp,ngrbspar,ngrbsparc,ngrbspd,ngrbsr,ngrbsrc,ngrbsshf,ngrbssr)
!$OMP THREADPRIVATE(ngrbssrc,ngrbssrd,ngrbsst,ngrbstinc,ngrbstl1,ngrbstl2,ngrbstl3,ngrbstl4,ngrbstr,ngrbstrc,ngrbstrd)
!$OMP THREADPRIVATE(ngrbsund,ngrbsuvb,ngrbsvd,ngrbswl1,ngrbswl2,ngrbswl3,ngrbswl4,ngrbt,ngrbtbt,ngrbtcc,ngrbtcghg)
!$OMP THREADPRIVATE(ngrbtcgrg,ngrbtciw,ngrbtclw,ngrbtco3,ngrbtctrac,ngrbtcw,ngrbtcwv,ngrbth,ngrbtp,ngrbtrac,ngrbtsn)
!$OMP THREADPRIVATE(ngrbtsp,ngrbtsr,ngrbtsrc,ngrbttr,ngrbttrc,ngrbtvh,ngrbtvl,ngrbu,ngrbv,ngrbveg,ngrbvimd,ngrbvo)
!$OMP THREADPRIVATE(ngrbvso,ngrbw,ngrbz,njdiag,njdomai,njiter,nleg,nlocgrb,nmethod,nmfr,nncep,nreference,nsec0,nsec1)
!$OMP THREADPRIVATE(nsec2spm,nsec2spp,nsec3,nsec4,nsmaxnp,nsteplpp,nstream,nsystem,ntotens,nukm,rsec3,ssec3)
!$OMP THREADPRIVATE(msec2gg,ngrbgp2,ngrbgp3,ngrbsp2,ngrbsp3,nsec2gg,rsec2,ssec2)
END MODULE yomgrb