[3880] | 1 | MODULE YOM_YGFL |
---|
| 2 | |
---|
| 3 | USE PARKIND1 , ONLY : JPIM, JPRB |
---|
| 4 | |
---|
| 5 | IMPLICIT NONE |
---|
| 6 | SAVE |
---|
| 7 | |
---|
| 8 | !------------------------------------------------------------------------- |
---|
| 9 | ! Contains the descriptors of GFL arrays |
---|
| 10 | !------------------------------------------------------------------------- |
---|
| 11 | |
---|
| 12 | ! JPGFL : Max number of GFL fields |
---|
| 13 | ! JPNAMED_GFL : Number of currently pre-defined components of GFL |
---|
| 14 | ! JPGHG : Number of greenhouse gas fields |
---|
| 15 | ! JPGRG : Number of reactive gas fields |
---|
| 16 | ! JPCHEM : Number of chemical species |
---|
| 17 | ! JPAERO : Number of active aerosol fields |
---|
| 18 | ! JPAEROUT: Number of output aerosol fields |
---|
| 19 | ! JPUVP : Number of output from UV processor |
---|
| 20 | ! JPTRAC : Number of tracers for diagnostics |
---|
| 21 | ! JPERA40 : Number of ERA40 diagnostic fields |
---|
| 22 | ! JPCH4S : Number of added fields related to methane |
---|
| 23 | ! JPNOGW : Number of diagnostic fields for NORO GWD SCHEME |
---|
| 24 | ! JPSLDIA : Number of SL dynamics diagnostic fields |
---|
| 25 | ! JPCHEM_ASSIM : Maximum number of assimilated of chemical species |
---|
| 26 | !------------------------------------------------------------------------- |
---|
| 27 | |
---|
| 28 | INTEGER(KIND=JPIM), PARAMETER :: JPGFL=2163 |
---|
| 29 | INTEGER(KIND=JPIM), PARAMETER :: JPNAMED_GFL=27 |
---|
| 30 | INTEGER(KIND=JPIM), PARAMETER :: JPGHG=3 |
---|
| 31 | INTEGER(KIND=JPIM), PARAMETER :: JPTRAC=10 |
---|
| 32 | INTEGER(KIND=JPIM), PARAMETER :: JPCHEM=66 |
---|
| 33 | INTEGER(KIND=JPIM), PARAMETER :: JPGRG=5 |
---|
| 34 | INTEGER(KIND=JPIM), PARAMETER :: JPCHEM_ASSIM=5 |
---|
| 35 | INTEGER(KIND=JPIM), PARAMETER :: JPAERO=16 |
---|
| 36 | INTEGER(KIND=JPIM), PARAMETER :: JPFORC=1800 |
---|
| 37 | INTEGER(KIND=JPIM), PARAMETER :: JPERA40=14 |
---|
| 38 | INTEGER(KIND=JPIM), PARAMETER :: JPSLDIA=7 |
---|
| 39 | INTEGER(KIND=JPIM), PARAMETER :: JPEZDIAG=50 |
---|
| 40 | INTEGER(KIND=JPIM), PARAMETER :: JPCH4S=2 |
---|
| 41 | INTEGER(KIND=JPIM), PARAMETER :: JPNOGW=2 |
---|
| 42 | INTEGER(KIND=JPIM), PARAMETER :: JPAEROUT=17 |
---|
| 43 | INTEGER(KIND=JPIM), PARAMETER :: JPUVP=2 |
---|
| 44 | INTEGER(KIND=JPIM), PARAMETER :: JPPHYS=8 |
---|
| 45 | INTEGER(KIND=JPIM), PARAMETER :: GRIB_CODE_GFL_PHYS=81 ! AJGDB hopefully harmless |
---|
| 46 | |
---|
| 47 | TYPE TYPE_GFL_COMP ! Individual field descriptor |
---|
| 48 | |
---|
| 49 | SEQUENCE ! daand: necessary to avoid memory corruption with gfortran 4.3.3 |
---|
| 50 | |
---|
| 51 | CHARACTER(LEN=16) :: CNAME = '' ! ARPEGE field name |
---|
| 52 | INTEGER(KIND=JPIM) :: IGRBCODE = -999 ! GRIB code |
---|
| 53 | LOGICAL :: LADV = .FALSE. ! Field advected or not |
---|
| 54 | LOGICAL :: LADV5 = .FALSE. ! Field advected without wind increments |
---|
| 55 | LOGICAL :: LTDIABLIN = .FALSE. ! Diabatic tendency is interpolated by lin. int. |
---|
| 56 | LOGICAL :: LHORTURB = .FALSE. ! Horizontal part affected by 3D turbulence |
---|
| 57 | INTEGER(KIND=JPIM) :: NREQIN = 0 ! 1 if field requiered in input, 0 if not, -1 if initialised |
---|
| 58 | ! with a reference value REFVALI |
---|
| 59 | LOGICAL :: LREQOUT = .FALSE. ! T if field requiered in output |
---|
| 60 | LOGICAL :: LGPINGP = .TRUE. ! GP field input as GP |
---|
| 61 | LOGICAL :: LGP = .FALSE. ! Field exists and of grid-point type |
---|
| 62 | LOGICAL :: LSP = .FALSE. ! Field exists and of spectral type |
---|
| 63 | LOGICAL :: LCDERS = .FALSE. ! Derivatives required (spectral only) |
---|
| 64 | LOGICAL :: LACTIVE = .FALSE. ! Field in use |
---|
| 65 | LOGICAL :: LTHERMACT = .FALSE. ! Field thermodynamically active |
---|
| 66 | REAL(KIND=JPRB) :: R = 0.0_JPRB |
---|
| 67 | REAL(KIND=JPRB) :: RCP = 0.0_JPRB |
---|
| 68 | LOGICAL :: LT9 = .FALSE. ! Field in t-dt GFL |
---|
| 69 | LOGICAL :: LT1 = .FALSE. ! Field in t+dt GFL |
---|
| 70 | LOGICAL :: LT5 = .FALSE. ! Field in trajectory GFL |
---|
| 71 | LOGICAL :: LPHY = .FALSE. ! Field in physics GFL |
---|
| 72 | LOGICAL :: LPT = .FALSE. ! Field in PC phy. tend. GFL (GFLPT) |
---|
| 73 | LOGICAL :: LTRAJIO = .FALSE. ! Field written to and from trajectory structure |
---|
| 74 | LOGICAL :: LDIAG = .FALSE. ! Field is "diagnostic" at t; e.g. cloud fraction |
---|
| 75 | LOGICAL :: LPC = .FALSE. ! Field in predictor/corrector time stepping (GFLPC) |
---|
| 76 | REAL(KIND=JPRB) :: REFVALI = 0.0_JPRB ! Reference value for init, used in case NREQIN==-1 |
---|
| 77 | ! LAM specific attributes (Arome/Aladin) |
---|
| 78 | LOGICAL :: LADJUST0 = .FALSE. ! True if field is thermodynamically adjusted at t |
---|
| 79 | ! (immediatly after inverse spectral transforms) |
---|
| 80 | LOGICAL :: LADJUST1 = .FALSE. ! True if field is thermodynamically adjusted at t+dt |
---|
| 81 | ! (after SL interpolations and NL residuals) |
---|
| 82 | INTEGER(KIND=JPIM) :: NCOUPLING = 0 ! 1 if field is coupled by Davies relaxation, 0 if not, |
---|
| 83 | ! -1 if coupled with reference value for coupling REFVALC |
---|
| 84 | REAL(KIND=JPRB) :: REFVALC = 0.0_JPRB ! Reference value for coupling, used in case NCOUPLING==-1 |
---|
| 85 | LOGICAL :: LBIPER = .FALSE. ! True if field must be biperiodised inside the transforms |
---|
| 86 | ! End LAM specific attributes (Arome/Aladin) |
---|
| 87 | CHARACTER(LEN=12) :: CSLINT = '' ! S.L interpolaion "type" |
---|
| 88 | INTEGER(KIND=JPIM) :: MP = -99999999 ! Basic field "pointer" |
---|
| 89 | INTEGER(KIND=JPIM) :: MPL = -99999999 ! zonal derivative "pointer" |
---|
| 90 | INTEGER(KIND=JPIM) :: MPM = -99999999 ! Meridional derivative "pointer" |
---|
| 91 | INTEGER(KIND=JPIM) :: MP9 = -99999999 ! Basic field "pointer" t-dt |
---|
| 92 | INTEGER(KIND=JPIM) :: MP9_PH = -99999999 ! Basic field "pointer" for Physics |
---|
| 93 | INTEGER(KIND=JPIM) :: MP1 = -99999999 ! Basic field "pointer" t+dt |
---|
| 94 | INTEGER(KIND=JPIM) :: MP5 = -99999999 ! Basic field "pointer" trajectory |
---|
| 95 | INTEGER(KIND=JPIM) :: MP5L = -99999999 ! zonal derivative "pointer" trajectory |
---|
| 96 | INTEGER(KIND=JPIM) :: MP5M = -99999999 ! Meridional derivative "pointer" trajectory |
---|
| 97 | INTEGER(KIND=JPIM) :: MPSLP = -99999999 ! Basic field "pointer" physics |
---|
| 98 | INTEGER(KIND=JPIM) :: MPSP = -99999999 ! Basic field "pointer" spectral space |
---|
| 99 | INTEGER(KIND=JPIM) :: MP_SPL = -99999999 ! Basic field "pointer" spline interpolation |
---|
| 100 | INTEGER(KIND=JPIM) :: MP_SL1 = -99999999 ! Basic field "pointer" in SLBUF1 |
---|
| 101 | INTEGER(KIND=JPIM) :: MP_SLX = -99999999 ! Basic field "pointer" in SLBUF1 for CPG_PT |
---|
| 102 | INTEGER(KIND=JPIM) :: MPPT = -99999999 ! Physics tendency "pointer" |
---|
| 103 | INTEGER(KIND=JPIM) :: MPPC = -99999999 ! Predictor/corrector auxiliary array "pointer" |
---|
| 104 | |
---|
| 105 | ! daand: INTFLEX attributes |
---|
| 106 | LOGICAL :: LWATER ! TRUE for water species |
---|
| 107 | LOGICAL :: LPRECIP ! TRUE for precipitating water species |
---|
| 108 | REAL(KIND=JPRB) :: RLZER ! Latent heat change at 0K |
---|
| 109 | |
---|
| 110 | ! gems nl ext |
---|
| 111 | INTEGER(KIND=JPIM) :: NCOUPLO4 ! Coupled to CTM by OASIS4 intefrace |
---|
| 112 | LOGICAL :: LASSIM ! use as Control Variable (either monitored or assimilated) |
---|
| 113 | INTEGER(KIND=JPIM) :: IGRIBDV ! GRIB code of deposition velocity |
---|
| 114 | INTEGER(KIND=JPIM) :: IGRIBTC ! GRIB code of Total Column |
---|
| 115 | INTEGER(KIND=JPIM) :: IGRIBSFC ! GRIB code of Surface Flux |
---|
| 116 | LOGICAL :: LDIFF ! Diffusion on |
---|
| 117 | LOGICAL :: LCONV ! Convection on |
---|
| 118 | REAL(KIND=JPRB) :: RMOLMASS ! Molar Mass |
---|
| 119 | REAL(KIND=JPRB) :: REFOLD ! Efolding decay time |
---|
| 120 | REAL(KIND=JPRB) :: HENRYA ! Henry constant a |
---|
| 121 | REAL(KIND=JPRB) :: HENRYB ! Henry constant b |
---|
| 122 | LOGICAL :: LNEGFIX ! Cut off negative values in sugridug an |
---|
| 123 | LOGICAL :: LMASSFIX ! Correct mass error of sl advection in gpmodel (if LTRCMFIX) |
---|
| 124 | TYPE(TYPE_GFL_COMP),POINTER :: PREVIOUS ! Pointer to previously def. field |
---|
| 125 | |
---|
| 126 | END TYPE TYPE_GFL_COMP |
---|
| 127 | |
---|
| 128 | TYPE TYPE_GFL_NAML ! Individual field descriptor for namelist input |
---|
| 129 | |
---|
| 130 | SEQUENCE ! daand: necessary to avoid memory corruption with gfortran 4.3.3 |
---|
| 131 | |
---|
| 132 | CHARACTER(LEN=16) :: CNAME ! ARPEGE field name |
---|
| 133 | INTEGER(KIND=JPIM) :: IGRBCODE ! GRIB code |
---|
| 134 | INTEGER(KIND=JPIM) :: NREQIN ! 1 if field required in input, 0 if not, -1 if initialised |
---|
| 135 | ! with a reference value REFVALI |
---|
| 136 | REAL(KIND=JPRB) :: REFVALI ! Reference value for initialisation, used in case NREQIN==-1 |
---|
| 137 | LOGICAL :: LREQOUT ! T if field requiered in output |
---|
| 138 | LOGICAL :: LGPINGP ! GP field input as GP |
---|
| 139 | LOGICAL :: LGP ! Field exists and of grid-point type |
---|
| 140 | LOGICAL :: LSP ! Field exists and of spectral type |
---|
| 141 | LOGICAL :: LCDERS ! Derivatives required (spectral only) |
---|
| 142 | LOGICAL :: LT9 ! Field in t-dt GFL |
---|
| 143 | LOGICAL :: LT1 ! Field in t+dt GFL |
---|
| 144 | LOGICAL :: LT5 ! Field in trajectory GFL |
---|
| 145 | LOGICAL :: LPHY ! Field with physics tendencies GFL |
---|
| 146 | LOGICAL :: LPT ! Field in PC physics tendency GFLPT |
---|
| 147 | LOGICAL :: LTRAJIO ! Field written to and from trajectory structure |
---|
| 148 | LOGICAL :: LDIAG ! Field is "diagnostic" at t; e.g. cloud fraction |
---|
| 149 | LOGICAL :: LPC ! Field in predictor/corrector time stepping GFLPC |
---|
| 150 | LOGICAL :: LADV ! Field advected or not |
---|
| 151 | LOGICAL :: LADV5 ! Field advected without wind increments |
---|
| 152 | LOGICAL :: LINTLIN ! Linear interpolation for field |
---|
| 153 | LOGICAL :: LTDIABLIN ! Diabatic tendency is interpolated by linear int. |
---|
| 154 | LOGICAL :: LHORTURB ! Horizontal part affected by 3D turbulence |
---|
| 155 | LOGICAL :: LQM ! quasi-monotonous interpolation for field |
---|
| 156 | LOGICAL :: LQMH ! quasi-monotonous interpolation in horizontal for field |
---|
| 157 | LOGICAL :: LQM3D ! quasi-monotone interpolation applied directly in 3 dimensions |
---|
| 158 | LOGICAL :: LSLHD ! Semi-lagrangian horizontal diffusion used for field |
---|
| 159 | LOGICAL :: LCOMAD ! COMAD weights used for SL interpolation of field |
---|
| 160 | LOGICAL :: LHV ! Hermite vertical interpolation used for field (only ozone sofar) |
---|
| 161 | LOGICAL :: LVSPLIP ! vertical spline interpolation used for field (only ozone sofar) |
---|
| 162 | INTEGER(KIND=JPIM) :: NCOUPLING ! 1 if field is coupled by Davies relaxation, 0 if not, |
---|
| 163 | ! -1 if coupled with reference value for coupling REFVALC |
---|
| 164 | REAL(KIND=JPRB) :: REFVALC ! Reference value for coupling, used in case |
---|
| 165 | ! NCOUPLING==-1 |
---|
| 166 | ! gems nl ext |
---|
| 167 | INTEGER(KIND=JPIM) :: NCOUPLO4 ! Coupled to CTM by OASIS4 intefrace =1 input,=2 in&output,=-1 none |
---|
| 168 | LOGICAL :: LASSIM ! use as Control Variable (either monitored or assimilated) |
---|
| 169 | INTEGER(KIND=JPIM) :: IGRIBDV ! GRIB code of deposition velocity |
---|
| 170 | INTEGER(KIND=JPIM) :: IGRIBTC ! GRIB code of Total Column |
---|
| 171 | INTEGER(KIND=JPIM) :: IGRIBSFC ! GRIB code of Surface Flux |
---|
| 172 | LOGICAL :: LDIFF ! Diffusion on |
---|
| 173 | LOGICAL :: LCONV ! Convection on |
---|
| 174 | LOGICAL :: LNEGFIX ! Cut off negative values in sugridug and callpar |
---|
| 175 | LOGICAL :: LMASSFIX ! Correct mass error of sl advection in gpmodel (if LTRCMFIX) |
---|
| 176 | REAL(KIND=JPRB) :: RMOLMASS ! Molar Mass |
---|
| 177 | REAL(KIND=JPRB) :: REFOLD ! Efolding decay time |
---|
| 178 | REAL(KIND=JPRB) :: HENRYA ! Henry constant a |
---|
| 179 | REAL(KIND=JPRB) :: HENRYB ! Henry constant b |
---|
| 180 | |
---|
| 181 | END TYPE TYPE_GFL_NAML |
---|
| 182 | |
---|
| 183 | !------------------------------------------------------------------------- |
---|
| 184 | ! Derived types for describing the GFL structure. |
---|
| 185 | !------------------------------------------------------------------------- |
---|
| 186 | ! Modifications: |
---|
| 187 | ! 03/07/09 C. Fischer - add Arome/Aladin attributes |
---|
| 188 | ! 03/10/01 C. Moussy - add Arome/Aladin attributes coupling |
---|
| 189 | ! 03/10/31 M. Tudor - add physics tendencies for predictor-corrector |
---|
| 190 | ! 05/10/10 J. Haseler - switch for I/O to trajectory structure |
---|
| 191 | ! 2004-Nov F. Vana - update of CSLINT attribute |
---|
| 192 | ! 20-Feb-2005 Vivoda - 3TL Eul PC scheme (GFLPC) |
---|
| 193 | ! 07/06/27 E. Holm - TL/AD advection without wind increments LADV5 |
---|
| 194 | ! 12/04/08 J. Flemming - GFL attribute extention for GEMS |
---|
| 195 | ! 22-Feb-11 F. Vana - LTDIABLIN and LHORTURB |
---|
| 196 | ! spring 2011 ECMWF - LINTLIN |
---|
| 197 | ! Nov. 2013 - LCOMAD |
---|
| 198 | ! 2013-11, D. Degrauwe - INTFLEX attributes |
---|
| 199 | |
---|
| 200 | TYPE TYPE_GFLD |
---|
| 201 | |
---|
| 202 | SEQUENCE ! daand: necessary to avoid memory corruption with gfortran 4.3.3 |
---|
| 203 | |
---|
| 204 | ! Overall descriptor,dimensioning etc. |
---|
| 205 | INTEGER(KIND=JPIM) :: NUMFLDS = 0 ! Number of GFL fields |
---|
| 206 | INTEGER(KIND=JPIM) :: NDERS = 0 ! Number of horizontal derivatives fields |
---|
| 207 | INTEGER(KIND=JPIM) :: NUMSPFLDS = 0 ! Number of spectrally represented GFL fields |
---|
| 208 | INTEGER(KIND=JPIM) :: NUMGPFLDS = 0 ! Number of grid-point GFL fields |
---|
| 209 | INTEGER(KIND=JPIM) :: NUMFLDS9 = 0 ! Number of GFL fields in (t-dt) part |
---|
| 210 | INTEGER(KIND=JPIM) :: NUMFLDS1 = 0 ! Number of GFL fields in (t+dt) array |
---|
| 211 | INTEGER(KIND=JPIM) :: NUMSPFLDS1 = 0 ! Number of spectrally represented GFL fields (t+dt) |
---|
| 212 | INTEGER(KIND=JPIM) :: NUMFLDS5 = 0 ! Number of GFL fields (trajectory) |
---|
| 213 | INTEGER(KIND=JPIM) :: NUMFLDSPHY = 0 ! Number of GFL fields (phys.) |
---|
| 214 | INTEGER(KIND=JPIM) :: NUMFLDS_SPL = 0 ! Number of GFL fields (S.L. spline interpolation) |
---|
| 215 | INTEGER(KIND=JPIM) :: NUMFLDS_SL1 = 0 ! Number of GFL fields in S.L. buffer 1 |
---|
| 216 | INTEGER(KIND=JPIM) :: NUMFLDSPC = 0 ! Number of GFL fields (predictor/corrector) |
---|
| 217 | INTEGER(KIND=JPIM) :: NDIM = 0 ! Dimension of main array holding GFL fields(GFL) |
---|
| 218 | INTEGER(KIND=JPIM) :: NUMFLDSPT = 0 ! Number of GFL fields (phy. tend.) |
---|
| 219 | INTEGER(KIND=JPIM) :: NDIM0 = 0 ! Dimension of t0 part of GFL |
---|
| 220 | INTEGER(KIND=JPIM) :: NDIM9 = 0 ! Dimension of t-dt part of GFL |
---|
| 221 | INTEGER(KIND=JPIM) :: NDIM1 = 0 ! Dimension of t+dt array (GFLT1) |
---|
| 222 | INTEGER(KIND=JPIM) :: NDIM5 = 0 ! Dimension of traj. GFL array (GFL5) |
---|
| 223 | INTEGER(KIND=JPIM) :: NDIMSLP = 0 ! Diminsion of S.L. phys. GFL array (GFLSLP) |
---|
| 224 | INTEGER(KIND=JPIM) :: NDIM_SPL = 0 ! Dim. of arrays holding GFL fields (S.L.spline int.) |
---|
| 225 | INTEGER(KIND=JPIM) :: NDIMPT = 0 ! Dimension of phy. tend. GFL array (GFLPT) |
---|
| 226 | INTEGER(KIND=JPIM) :: NDIMPC = 0 ! Dimension of iterative scheme auxiliary array (GFLPC) |
---|
| 227 | |
---|
| 228 | INTEGER(KIND=JPIM) :: NGFL_EXT |
---|
| 229 | INTEGER(KIND=JPIM) :: NGFL_FORC |
---|
| 230 | INTEGER(KIND=JPIM) :: NGFL_EZDIAG |
---|
| 231 | INTEGER(KIND=JPIM) :: NGHG |
---|
| 232 | INTEGER(KIND=JPIM) :: NTRAC |
---|
| 233 | INTEGER(KIND=JPIM) :: NGRG |
---|
| 234 | INTEGER(KIND=JPIM) :: NGRG_CPLO4 |
---|
| 235 | INTEGER(KIND=JPIM) :: NGRG_ASSIM |
---|
| 236 | INTEGER(KIND=JPIM) :: NAERO |
---|
| 237 | INTEGER(KIND=JPIM) :: NACTAERO |
---|
| 238 | INTEGER(KIND=JPIM) :: NDDHAERO |
---|
| 239 | INTEGER(KIND=JPIM) :: NERA40 |
---|
| 240 | INTEGER(KIND=JPIM) :: NNOGW |
---|
| 241 | INTEGER(KIND=JPIM) :: NAEROUT |
---|
| 242 | INTEGER(KIND=JPIM) :: NUVP |
---|
| 243 | INTEGER(KIND=JPIM) :: NSLDIA |
---|
| 244 | INTEGER(KIND=JPIM) :: NSLDIAGP |
---|
| 245 | INTEGER(KIND=JPIM) :: NGFL_PHYS |
---|
| 246 | LOGICAL :: LCO2SFC |
---|
| 247 | LOGICAL :: LCH4SFC |
---|
| 248 | LOGICAL :: LAEROSFC |
---|
| 249 | LOGICAL :: LFIRE |
---|
| 250 | LOGICAL :: LAERODIU |
---|
| 251 | LOGICAL :: LTRCMFIX ! Activates tracer mass fixer |
---|
| 252 | LOGICAL :: LTRCMFIX_PS ! Adjust pressure to conserve dry mass in mass fixer calculations |
---|
| 253 | LOGICAL :: LAEROUT |
---|
| 254 | LOGICAL :: LUVPOUT |
---|
| 255 | LOGICAL :: LCHEM |
---|
| 256 | |
---|
| 257 | INTEGER(KIND=JPIM) :: NGEMS ! The total number of "GEMS" fields. |
---|
| 258 | INTEGER(KIND=JPIM) :: NCHEM |
---|
| 259 | INTEGER(KIND=JPIM) :: NCHEM_ASSIM |
---|
| 260 | INTEGER(KIND=JPIM) :: NCHEM_FLX |
---|
| 261 | INTEGER(KIND=JPIM) :: NCHEM_DV |
---|
| 262 | INTEGER(KIND=JPIM) :: NCHEM_TC |
---|
| 263 | INTEGER(KIND=JPIM) :: NCHEM_SCV |
---|
| 264 | |
---|
| 265 | ! ------------------------------------------------------------------ |
---|
| 266 | ! Mass fixers |
---|
| 267 | ! ------------------------------------------------------------------ |
---|
| 268 | INTEGER(KIND=JPIM) :: NNEGAFIX ! Num of fields to apply -ve fixer |
---|
| 269 | INTEGER(KIND=JPIM) :: NOPTNEGFIX ! 1: simple negative fixer (reset to 0) |
---|
| 270 | ! 2: reset to local minimum |
---|
| 271 | |
---|
| 272 | LOGICAL :: LQM3DCONS ! Bermejo & Staniforth quasi-monotone limiter with improved |
---|
| 273 | ! conservation option. When true, applied to all GFL s.t. LQM3D=true |
---|
| 274 | LOGICAL :: LADVNEGFIX ! Activates negative fixer for advection |
---|
| 275 | LOGICAL :: LTRCMFBC ! Activate Bermejo & Conde if true |
---|
| 276 | LOGICAL :: LTRCMFPR ! Activate Priestley algorithm if true |
---|
| 277 | LOGICAL :: LTRCMFMG ! Activate Mac Gregor's algorithm if true |
---|
| 278 | LOGICAL :: LEXTRADF ! Extra diagnostics |
---|
| 279 | |
---|
| 280 | |
---|
| 281 | INTEGER(KIND=JPIM) :: NFLDSFIX ! Number of fields to be fixed |
---|
| 282 | INTEGER(KIND=JPIM) :: NOPTMFIX ! Bermejo & Conde fixer option for calculating its weight |
---|
| 283 | INTEGER(KIND=JPIM) :: NOPTVFE ! Use Vertical FE in calculation of column mass total |
---|
| 284 | INTEGER(KIND=JPIM) :: NPMFIX ! Parameter used in weight calculation |
---|
| 285 | INTEGER(KIND=JPIM) :: NMFDIAGLEV ! Determines global diagnostic output level for fixer: |
---|
| 286 | ! 0 - nothing, 1 - norms printed, 2 - norms + monotonicity |
---|
| 287 | INTEGER(KIND=JPIM) :: NMFIXFLDS(JPNAMED_GFL+JPGHG+JPGRG+JPCHEM+JPAERO+JPTRAC) |
---|
| 288 | ! Index of fields to be corrected by mass fixers |
---|
| 289 | INTEGER(KIND=JPIM) :: NNEGFLDS(JPNAMED_GFL+JPGHG+JPGRG+JPCHEM+JPAERO+JPTRAC) |
---|
| 290 | ! Index of fields to be corrected by SL -ve fixer |
---|
| 291 | REAL(KIND=JPRB) :: ZMFIXEPS ! Threshold for mass fixing scheme |
---|
| 292 | |
---|
| 293 | TYPE(TYPE_GFL_COMP) :: YCOMP(JPGFL) ! General descriptor of all components |
---|
| 294 | |
---|
| 295 | TYPE(TYPE_GFL_COMP),POINTER :: YQ => NULL() ! Specific humidity |
---|
| 296 | TYPE(TYPE_GFL_COMP),POINTER :: YI => NULL() ! Ice water |
---|
| 297 | TYPE(TYPE_GFL_COMP),POINTER :: YL => NULL() ! Liquid water |
---|
| 298 | TYPE(TYPE_GFL_COMP),POINTER :: YLCONV => NULL() ! Liquid water (CONV. PART) |
---|
| 299 | TYPE(TYPE_GFL_COMP),POINTER :: YICONV => NULL() ! Ice water (CONV. PART) |
---|
| 300 | TYPE(TYPE_GFL_COMP),POINTER :: YRCONV => NULL() ! Rain (CONV. PART) |
---|
| 301 | TYPE(TYPE_GFL_COMP),POINTER :: YSCONV => NULL() ! Snow (CONV. PART) |
---|
| 302 | TYPE(TYPE_GFL_COMP),POINTER :: YIRAD => NULL() ! Radiative cloud Ice water |
---|
| 303 | TYPE(TYPE_GFL_COMP),POINTER :: YLRAD => NULL() ! Radiative cloud Liquid water |
---|
| 304 | TYPE(TYPE_GFL_COMP),POINTER :: YS => NULL() ! Snow |
---|
| 305 | TYPE(TYPE_GFL_COMP),POINTER :: YR => NULL() ! Rain |
---|
| 306 | TYPE(TYPE_GFL_COMP),POINTER :: YG => NULL() ! Graupel |
---|
| 307 | TYPE(TYPE_GFL_COMP),POINTER :: YH => NULL() ! Hail |
---|
| 308 | TYPE(TYPE_GFL_COMP),POINTER :: YTKE => NULL() ! Turbulent Kinetic Energy |
---|
| 309 | TYPE(TYPE_GFL_COMP),POINTER :: YTTE => NULL() ! Turbulent Total Energy |
---|
| 310 | TYPE(TYPE_GFL_COMP),POINTER :: YEFB1 => NULL() ! First variable EFB scheme |
---|
| 311 | TYPE(TYPE_GFL_COMP),POINTER :: YEFB2 => NULL() ! Second variable EFB scheme |
---|
| 312 | TYPE(TYPE_GFL_COMP),POINTER :: YEFB3 => NULL() ! Third variable EFB scheme |
---|
| 313 | TYPE(TYPE_GFL_COMP),POINTER :: YA => NULL() ! Cloud fraction |
---|
| 314 | TYPE(TYPE_GFL_COMP),POINTER :: YO3 => NULL() ! Ozone |
---|
| 315 | TYPE(TYPE_GFL_COMP),POINTER :: YSRC => NULL() ! Second-order flux for AROME s'rc'/2Sigma_s2 multiplied by Lambda_3 |
---|
| 316 | TYPE(TYPE_GFL_COMP),POINTER :: YMXL => NULL() ! Prognostic mixing length |
---|
| 317 | TYPE(TYPE_GFL_COMP),POINTER :: YSCC2 => NULL() ! Saturation deficit^2 for Tompkins |
---|
| 318 | TYPE(TYPE_GFL_COMP),POINTER :: YGCCA => NULL() ! Skewness for Tompkins |
---|
| 319 | TYPE(TYPE_GFL_COMP),POINTER :: YCPF => NULL() ! Convective precipitation flux |
---|
| 320 | TYPE(TYPE_GFL_COMP),POINTER :: YSPF => NULL() ! Stratiform precipitation flux |
---|
| 321 | TYPE(TYPE_GFL_COMP),POINTER :: YCVGQ => NULL() ! Moisture Convergence for french physics |
---|
| 322 | TYPE(TYPE_GFL_COMP),POINTER :: YQVA => NULL() ! total humidity variation |
---|
| 323 | TYPE(TYPE_GFL_COMP),POINTER :: YGHG(:) => NULL() ! Greenhouse Gases |
---|
| 324 | TYPE(TYPE_GFL_COMP),POINTER :: YGRG(:) => NULL() ! Reactive Gases |
---|
| 325 | TYPE(TYPE_GFL_COMP),POINTER :: YCHEM(:) => NULL() ! Chemistry |
---|
| 326 | TYPE(TYPE_GFL_COMP),POINTER :: YGRGTEND(:) => NULL() ! Reactive Gases Tendecies |
---|
| 327 | TYPE(TYPE_GFL_COMP),POINTER :: YAERO(:) => NULL() ! Aerosols |
---|
| 328 | TYPE(TYPE_GFL_COMP),POINTER :: YTRAC(:) => NULL() ! tracers for diagnostics |
---|
| 329 | TYPE(TYPE_GFL_COMP),POINTER :: YLRCH4 => NULL() ! CH4 loss rate (instantaneous field) |
---|
| 330 | TYPE(TYPE_GFL_COMP),POINTER :: YCH4S => NULL() ! CH4 atmospheric sink (accumulated field) |
---|
| 331 | TYPE(TYPE_GFL_COMP),POINTER :: YFORC(:) => NULL() ! large scale forcing |
---|
| 332 | TYPE(TYPE_GFL_COMP),POINTER :: YEZDIAG(:) => NULL() ! easy diagnostics |
---|
| 333 | TYPE(TYPE_GFL_COMP),POINTER :: YERA40(:) => NULL() ! ERA40 diagnostic fields |
---|
| 334 | TYPE(TYPE_GFL_COMP),POINTER :: YNOGW(:) => NULL() ! NORO GWD SCHEME |
---|
| 335 | TYPE(TYPE_GFL_COMP),POINTER :: YSLDIA(:) => NULL() ! SL dynamics diagnostics |
---|
| 336 | TYPE(TYPE_GFL_COMP),POINTER :: YAEROUT(:) => NULL() ! Aerosol outputs |
---|
| 337 | TYPE(TYPE_GFL_COMP),POINTER :: YUVP(:) => NULL() ! UV-processor output |
---|
| 338 | TYPE(TYPE_GFL_COMP),POINTER :: YPHYS(:) => NULL() ! PHYS output |
---|
| 339 | |
---|
| 340 | |
---|
| 341 | TYPE(TYPE_GFL_COMP),POINTER :: YSDSAT => NULL() ! Standard Deviation of the |
---|
| 342 | ! SATuration Depression (Sigma_s) |
---|
| 343 | TYPE(TYPE_GFL_COMP),POINTER :: YCVV => NULL() ! Convective Vertical Velocity |
---|
| 344 | TYPE(TYPE_GFL_COMP),POINTER :: YRKTH => NULL() ! Rasch-Kristjansson H tendency |
---|
| 345 | TYPE(TYPE_GFL_COMP),POINTER :: YRKTQV => NULL() ! Rasch-Kristjansson Qv tendency |
---|
| 346 | TYPE(TYPE_GFL_COMP),POINTER :: YRKTQC => NULL() ! Rasch-Kristjansson Qc tendency |
---|
| 347 | |
---|
| 348 | ! Prognostic convection variables: add 6 named components |
---|
| 349 | TYPE(TYPE_GFL_COMP),POINTER :: YUOM => NULL() ! Updraught vert velocity |
---|
| 350 | TYPE(TYPE_GFL_COMP),POINTER :: YUAL => NULL() ! Updraught mesh fraction |
---|
| 351 | TYPE(TYPE_GFL_COMP),POINTER :: YDOM => NULL() ! Downdraught vert velocity |
---|
| 352 | TYPE(TYPE_GFL_COMP),POINTER :: YDAL => NULL() ! Downdraught mesh fraction |
---|
| 353 | TYPE(TYPE_GFL_COMP),POINTER :: YUEN => NULL() ! Updraught entrainment |
---|
| 354 | TYPE(TYPE_GFL_COMP),POINTER :: YUNEBH => NULL() ! pseudo-historic convective |
---|
| 355 | |
---|
| 356 | ! Extra fields |
---|
| 357 | |
---|
| 358 | TYPE(TYPE_GFL_COMP),POINTER :: YEXT(:) => NULL() ! Extra fields |
---|
| 359 | |
---|
| 360 | TYPE(TYPE_GFL_NAML) :: YQ_NL ! Specific humidity |
---|
| 361 | TYPE(TYPE_GFL_NAML) :: YI_NL ! Ice water |
---|
| 362 | TYPE(TYPE_GFL_NAML) :: YL_NL ! Liquid water |
---|
| 363 | TYPE(TYPE_GFL_NAML) :: YLCONV_NL ! Liquid water (CONV. PART) |
---|
| 364 | TYPE(TYPE_GFL_NAML) :: YICONV_NL ! Ice water (CONV. PART) |
---|
| 365 | TYPE(TYPE_GFL_NAML) :: YRCONV_NL ! Rain (CONV. PART) |
---|
| 366 | TYPE(TYPE_GFL_NAML) :: YSCONV_NL ! Snow (CONV. PART) |
---|
| 367 | TYPE(TYPE_GFL_NAML) :: YIRAD_NL ! Radiative cloud Ice water |
---|
| 368 | TYPE(TYPE_GFL_NAML) :: YLRAD_NL ! Radiative cloud Liquid water |
---|
| 369 | TYPE(TYPE_GFL_NAML) :: YS_NL ! Snow |
---|
| 370 | TYPE(TYPE_GFL_NAML) :: YR_NL ! Rain |
---|
| 371 | TYPE(TYPE_GFL_NAML) :: YG_NL ! Graupels |
---|
| 372 | TYPE(TYPE_GFL_NAML) :: YH_NL ! Hail |
---|
| 373 | TYPE(TYPE_GFL_NAML) :: YTKE_NL ! Turbulent Kinetic Energy |
---|
| 374 | TYPE(TYPE_GFL_NAML) :: YTTE_NL ! Turbulent Total Energy |
---|
| 375 | TYPE(TYPE_GFL_NAML) :: YEFB1_NL ! First variable EFB scheme |
---|
| 376 | TYPE(TYPE_GFL_NAML) :: YEFB2_NL ! Second variable EFB scheme |
---|
| 377 | TYPE(TYPE_GFL_NAML) :: YEFB3_NL ! Third variable EFB scheme |
---|
| 378 | TYPE(TYPE_GFL_NAML) :: YA_NL ! Cloud fraction |
---|
| 379 | TYPE(TYPE_GFL_NAML) :: YO3_NL ! Ozone |
---|
| 380 | TYPE(TYPE_GFL_NAML) :: YSRC_NL ! Second-order flux for AROME |
---|
| 381 | ! s'rc'/2Sigma_s2 |
---|
| 382 | ! multiplied by Lambda_3 |
---|
| 383 | TYPE(TYPE_GFL_NAML) :: YMXL_NL ! Prognostic mixing length |
---|
| 384 | TYPE(TYPE_GFL_NAML) :: YSCC2_NL ! Saturation deficit^2 for Tompkins |
---|
| 385 | TYPE(TYPE_GFL_NAML) :: YGCCA_NL ! Skewness for Tompkins |
---|
| 386 | TYPE(TYPE_GFL_NAML) :: YCPF_NL ! Convective precipitation flux |
---|
| 387 | TYPE(TYPE_GFL_NAML) :: YSPF_NL ! Stratiform precipitation flux |
---|
| 388 | TYPE(TYPE_GFL_NAML) :: YCVGQ_NL ! Moisture Convergence for french physics |
---|
| 389 | TYPE(TYPE_GFL_NAML) :: YQVA_NL ! Total humidity variation |
---|
| 390 | |
---|
| 391 | TYPE(TYPE_GFL_NAML) :: YGHG_NL(JPGHG) ! Greenhouse Gases |
---|
| 392 | TYPE(TYPE_GFL_NAML) :: YGRG_NL(JPGRG) ! Reactive Gases |
---|
| 393 | TYPE(TYPE_GFL_NAML) :: YCHEM_NL(JPCHEM) ! Chemical species |
---|
| 394 | TYPE(TYPE_GFL_NAML) :: YGRGTEND_NL(JPGRG) ! Reactive Gases Tendecies |
---|
| 395 | TYPE(TYPE_GFL_NAML) :: YAERO_NL(JPAERO) ! Aerosol fields |
---|
| 396 | TYPE(TYPE_GFL_NAML) :: YTRAC_NL(JPTRAC) ! Tracers for diagnostics |
---|
| 397 | TYPE(TYPE_GFL_NAML) :: YERA40_NL(JPERA40) ! ERA40 diagnostic fields |
---|
| 398 | TYPE(TYPE_GFL_NAML) :: YNOGW_NL(JPNOGW) ! NORO GWD SCHEME |
---|
| 399 | TYPE(TYPE_GFL_NAML) :: YSLDIA_NL(JPSLDIA) ! SL dynamics diagnostics |
---|
| 400 | TYPE(TYPE_GFL_NAML) :: YLRCH4_NL ! CH4 loss rate |
---|
| 401 | TYPE(TYPE_GFL_NAML) :: YCH4S_NL ! CH4 atmospheric sink |
---|
| 402 | TYPE(TYPE_GFL_NAML) :: YAEROUT_NL(JPAEROUT) ! Aerosol outputs |
---|
| 403 | TYPE(TYPE_GFL_NAML) :: YUVP_NL(JPUVP) ! UV-processor outputs |
---|
| 404 | TYPE(TYPE_GFL_NAML) :: YRKTH_NL ! Rasch-Kristjansson H tendency |
---|
| 405 | TYPE(TYPE_GFL_NAML) :: YRKTQV_NL ! Rasch-Kristjansson Qv tendency |
---|
| 406 | TYPE(TYPE_GFL_NAML) :: YRKTQC_NL ! Rasch-Kristjansson Qc tendency |
---|
| 407 | TYPE(TYPE_GFL_NAML) :: YPHYS_NL(JPPHYS) ! PHYS outputs |
---|
| 408 | |
---|
| 409 | ! Extra fields |
---|
| 410 | TYPE(TYPE_GFL_NAML) :: YSDSAT_NL ! Standard Deviation of the |
---|
| 411 | ! SATuration Depression (Sigma_s) |
---|
| 412 | TYPE(TYPE_GFL_NAML) :: YCVV_NL ! Convective Vertical Velocity |
---|
| 413 | TYPE(TYPE_GFL_NAML) :: YFORC_NL(JPFORC) ! Forcing precursor |
---|
| 414 | TYPE(TYPE_GFL_NAML) :: YEZDIAG_NL(JPEZDIAG) ! Easy diagnostics |
---|
| 415 | TYPE(TYPE_GFL_NAML) :: YEXT_NL(JPGFL-JPNAMED_GFL-JPGHG-JPGRG-JPFORC-JPEZDIAG-JPAERO-JPTRAC-JPERA40-& |
---|
| 416 | & JPNOGW-JPSLDIA-JPCH4S-JPAEROUT-JPUVP-JPCHEM-JPPHYS) ! Extra fields |
---|
| 417 | |
---|
| 418 | ! Prognostic convection variables: 6 more namelist components |
---|
| 419 | TYPE(TYPE_GFL_NAML) :: YUOM_NL ! Updraught vert velocity |
---|
| 420 | TYPE(TYPE_GFL_NAML) :: YUAL_NL ! Updraught mesh fraction |
---|
| 421 | TYPE(TYPE_GFL_NAML) :: YDOM_NL ! Downdraught vert velocity |
---|
| 422 | TYPE(TYPE_GFL_NAML) :: YDAL_NL ! Downdraught mesh fraction |
---|
| 423 | TYPE(TYPE_GFL_NAML) :: YUEN_NL ! Updraught entrainment |
---|
| 424 | TYPE(TYPE_GFL_NAML) :: YUNEBH_NL ! Pseudi Hist Conv cloud fraction |
---|
| 425 | |
---|
| 426 | END TYPE TYPE_GFLD |
---|
| 427 | |
---|
| 428 | ! GFL general descriptor |
---|
| 429 | TYPE(TYPE_GFLD), POINTER :: YGFL => NULL() |
---|
| 430 | |
---|
| 431 | END MODULE YOM_YGFL |
---|