| 1 | MODULE YOMFPC |
|---|
| 2 | |
|---|
| 3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
|---|
| 4 | |
|---|
| 5 | USE PARFPOS, ONLY : JPOS2DF ,JPOS3DF ,JPOS3H ,JPOS3P ,JPOS3PV , & |
|---|
| 6 | & JPOS3S ,JPOS3TH ,JPOSCFU ,JPOSDIR ,JPOSLEN ,JPOSDOM , & |
|---|
| 7 | & JPOSSGP ,JPOSXFU |
|---|
| 8 | |
|---|
| 9 | IMPLICIT NONE |
|---|
| 10 | |
|---|
| 11 | SAVE |
|---|
| 12 | |
|---|
| 13 | ! ------------------------------------------------------------------ |
|---|
| 14 | |
|---|
| 15 | !* Scientific and technical variables |
|---|
| 16 | |
|---|
| 17 | ! === TECHNICAL VARIABLES === |
|---|
| 18 | |
|---|
| 19 | ! CFPDIR : prefix (path) for the output files |
|---|
| 20 | ! CFPIDEN: identificator of the output files |
|---|
| 21 | ! CFPFMT : format of the output files ('MODEL','GAUSS','LELAM' or 'LALON') |
|---|
| 22 | ! CFPDOM, C1FPDOM : names of the subdomains |
|---|
| 23 | ! CFP3DF, C1FP3DF : names 3D dynamics fields to compute |
|---|
| 24 | ! CFP3DFS,C1FP3DFS: names 3D dynamics fields to compute on model levels |
|---|
| 25 | ! CFP3DFP,C1FP3DFP: names 3D dynamics fields to compute on Pressure levels |
|---|
| 26 | ! CFP3DFH,C1FP3DFH: names 3D dynamics fields to compute on H levels |
|---|
| 27 | ! CFP3DFT,C1FP3DFT: names 3D dynamics fields to compute on THETA levels |
|---|
| 28 | ! CFP3DFV,C1FP3DFV: names 3D dynamics fields to compute on PV levels |
|---|
| 29 | ! CFP2DF, C1FP2DF : names 2D dynamics fields to compute |
|---|
| 30 | ! CFPPHY, C1FPPHY : names of physical fields to be post-processed |
|---|
| 31 | ! CFPCFU, C1FPCFU : names of cumulated fluxes fields to be post-processed |
|---|
| 32 | ! CFPXFU, C1FPXFU : names of instantaneous fluxes fields to be post-proc |
|---|
| 33 | ! Remark: C1FPXXX=CFPXXX(1)//CFPXXX(2)//...//CFPXXX(JPOSXXX) |
|---|
| 34 | ! for XXX=DOM,3DF,2DF,PHY,CFU,XFU. |
|---|
| 35 | |
|---|
| 36 | ! Variables used for ECMWF |
|---|
| 37 | |
|---|
| 38 | ! MFP3DFS: Gribcodes of 3D dynamics fields to compute on model levels |
|---|
| 39 | ! MFP3DFH: Gribcodes of 3D dynamics fields to compute on H levels |
|---|
| 40 | ! MFP3DFT: Gribcodes of 3D dynamics fields to compute on THETA levels |
|---|
| 41 | ! MFP3DFV: Gribcodes of 3D dynamics fields to compute on PV levels |
|---|
| 42 | ! MFP3DFP: Gribcodes of 3D dynamics fields to compute on P levels |
|---|
| 43 | ! MFP2DF : Gribcodes of 2D dynamics fields to compute |
|---|
| 44 | ! MFPPHY : Gribcodes of physical fields to processed |
|---|
| 45 | |
|---|
| 46 | ! RFP3P : post-processing pressure levels |
|---|
| 47 | ! RFP3H : post-processing height (above orography) levels |
|---|
| 48 | ! RFP3TH : post-processing potential temperature levels |
|---|
| 49 | ! RFP3PV : post-processing potential vorticity levels |
|---|
| 50 | ! NRFP3S : post-processing eta levels (CONF. 927 only) |
|---|
| 51 | |
|---|
| 52 | ! LFPCNT : control varying output variables according to time step |
|---|
| 53 | ! LTRACEFP: trace for Full-POS |
|---|
| 54 | ! NFPGRIB : level of GRIB coding in output file ARPEGE/ALADIN. |
|---|
| 55 | ! 0 : no packing at all |
|---|
| 56 | ! 1 : standart GRIB encoding |
|---|
| 57 | ! 2 : modified GRIB encoding |
|---|
| 58 | |
|---|
| 59 | ! NFPDOM : useful dimension of CFPDOM |
|---|
| 60 | ! NFP3DF : useful dimension of CFP3DF |
|---|
| 61 | ! NFP3DFS: useful dimension of CFP3DFS |
|---|
| 62 | ! NFP3DFH: useful dimension of CFP3DFH |
|---|
| 63 | ! NFP3DFT: useful dimension of CFP3DFT |
|---|
| 64 | ! NFP3DFV: useful dimension of CFP3DFV |
|---|
| 65 | ! NFP3DFP: useful dimension of CFP3DFP |
|---|
| 66 | ! NFP2DF : useful dimension of CFP2DF |
|---|
| 67 | ! NFPPHY : useful dimension of CFPPHY |
|---|
| 68 | ! NFPCFU : useful dimension of CFPCFU |
|---|
| 69 | ! NFPXFU : useful dimension of CFPXFU |
|---|
| 70 | ! NFP3P : useful dimension of RFP3P |
|---|
| 71 | ! NFP3H : useful dimension of RFP3H |
|---|
| 72 | ! NFP3TH : useful dimension of RFP3TH |
|---|
| 73 | ! NFP3PV : useful dimension of RFP3PV |
|---|
| 74 | ! NFP3S : useful dimension of NRFP3S |
|---|
| 75 | ! NFPXLEV: max. number of pp. levels |
|---|
| 76 | ! MFP3DYN : max. number of 3D-dynamic fields neeeded for LFPART2 |
|---|
| 77 | ! MFP2DYN : max. number of 2D-dynamic fields neeeded for LFPART2 |
|---|
| 78 | ! NFPDPHY : maximum number of physical fields |
|---|
| 79 | |
|---|
| 80 | ! The next keys force variables to be primitive for the post-processing : |
|---|
| 81 | ! LFPNHPD : .TRUE. if pp of NH pressure departure |
|---|
| 82 | ! LFPNHVD : .TRUE. if pp of NH vertical divergence |
|---|
| 83 | ! LFPNHVW : .TRUE. if pp of NH true vertical velocity |
|---|
| 84 | |
|---|
| 85 | ! LFPLOSP : .TRUE. = Fill Ps array with Log(Ps) |
|---|
| 86 | ! LFPRH100 : .TRUE. to convert relative humidity in percent |
|---|
| 87 | |
|---|
| 88 | ! === SCIENTIFIC VARIABLES === |
|---|
| 89 | |
|---|
| 90 | ! LFPSPEC: =.T. if pp. dyn. fields are written out as spectral coefficients |
|---|
| 91 | ! =.F. if pp. dyn. fields are written out as grid point values. |
|---|
| 92 | ! LFITP : =1 if pp. fields on P levels should be fited ; =0 otherwise |
|---|
| 93 | ! LFITT : =1 if pp. fields on THETA levels should be fited ; =0 otherwise |
|---|
| 94 | ! LFITV : =1 if pp. fields on PV levels should be fited ; =0 otherwise |
|---|
| 95 | |
|---|
| 96 | ! NFPCLI : usage level for climatology |
|---|
| 97 | ! =0 no climatology |
|---|
| 98 | ! =1 orography and land-sea mask of output only |
|---|
| 99 | ! =2 all available climatological fields of the current month |
|---|
| 100 | ! =3 shifting mean from the climatological fields of the current |
|---|
| 101 | ! month to the ones of the closest month |
|---|
| 102 | |
|---|
| 103 | ! NFPINDYN : type of interpolations for dynamical fields |
|---|
| 104 | ! NFPINPHY : type of interpolations for physical fields |
|---|
| 105 | ! 4 for bilinear, 12 for quadratic |
|---|
| 106 | |
|---|
| 107 | ! LFPQ : =.TRUE. if specific humidity is interpolated |
|---|
| 108 | ! : =.FALSE. if relative humidity is interpolated |
|---|
| 109 | ! LASQ : =.TRUE. if ASQ set to 80% saturation |
|---|
| 110 | |
|---|
| 111 | ! WSXI : max. surface moisture in input |
|---|
| 112 | ! WDXI : max. deep soil moisture in input |
|---|
| 113 | ! WSXO : max. surface moisture in output |
|---|
| 114 | ! WDXO : max. deep soil moisture in output |
|---|
| 115 | |
|---|
| 116 | ! FPBL : Critical Thickness of PBL |
|---|
| 117 | ! RFPCORR: Critical orography difference for correcting surface |
|---|
| 118 | ! temperature through standart profile. |
|---|
| 119 | ! RFPCSAB: Critical sand percentage difference for computing relative |
|---|
| 120 | ! soil moisture in ISBA |
|---|
| 121 | ! RFPCD2 : Critical soil depth difference for computing relative |
|---|
| 122 | ! soil moisture in ISBA |
|---|
| 123 | ! RFPVCAP: Minimum pressure of model level to provide an equatorial |
|---|
| 124 | ! cap in the computation of variables on constant PV surfaces |
|---|
| 125 | ! NFPLNPR : 0 => conventional formulation of (delta P) : ln(P(l)/P(l-1)) |
|---|
| 126 | ! 1 => formulation of (delta P) used in non hydrostatic model, |
|---|
| 127 | ! i.e. (P(l)-P(l-1))/SQRT(P(l)*P(l-1)) |
|---|
| 128 | ! LFPMOIS: Month allowed for climatology usage : |
|---|
| 129 | ! .F. => month of the model (forecast) |
|---|
| 130 | ! .T. => month of the file |
|---|
| 131 | ! NFPINCR : Bogussing indicator ; =0 for no bogussing ; =1 for bogussing |
|---|
| 132 | ! NFPLAKE : To overwrite created lakes or islands by specific data : |
|---|
| 133 | ! 0 => do not overwrite |
|---|
| 134 | ! -1 => overwrite with rhoughly interpolated data |
|---|
| 135 | ! +1 => overwrite with climatology |
|---|
| 136 | ! NFPCAPE : Kind of computation for CAPE & CIN : |
|---|
| 137 | ! 1 => from bottom model layer |
|---|
| 138 | ! 2 => from the most unstable layer |
|---|
| 139 | ! 3 => from mto standart height (2 meters) as recomputed values |
|---|
| 140 | ! 4 => from mto standart height (2 meters) out of fluxes |
|---|
| 141 | ! (used for analysis) |
|---|
| 142 | ! NFPSURFEX : Subcontract surface fields to SURFEX |
|---|
| 143 | ! 0 => no subcontract |
|---|
| 144 | ! 1 => transform native arp/ald surface fields to surfex fields |
|---|
| 145 | ! and write out by surfex |
|---|
| 146 | ! NFPMASK : number of masks for the interpolation of surface fields |
|---|
| 147 | ! 0 => no mask |
|---|
| 148 | ! 1 => land-sea mask |
|---|
| 149 | ! 2 => land mask, sea mask |
|---|
| 150 | ! LMOCONVAR : if true mocon is computed from N-1 level parameters (for varpack use) |
|---|
| 151 | |
|---|
| 152 | CHARACTER (LEN = JPOSDIR) :: CFPDIR |
|---|
| 153 | CHARACTER (LEN = 5) :: CFPFMT |
|---|
| 154 | CHARACTER (LEN = 8) :: CFPIDEN |
|---|
| 155 | CHARACTER (LEN = JPOSLEN) :: CFPDOM(JPOSDOM) |
|---|
| 156 | CHARACTER (LEN = 12) :: CFP3DF(JPOS3DF) |
|---|
| 157 | CHARACTER (LEN = 12) :: CFP3DFS(JPOS3DF) |
|---|
| 158 | CHARACTER (LEN = 12) :: CFP3DFP(JPOS3DF) |
|---|
| 159 | CHARACTER (LEN = 12) :: CFP3DFH(JPOS3DF) |
|---|
| 160 | CHARACTER (LEN = 12) :: CFP3DFT(JPOS3DF) |
|---|
| 161 | CHARACTER (LEN = 12) :: CFP3DFV(JPOS3DF) |
|---|
| 162 | CHARACTER (LEN = 16) :: CFP2DF(JPOS2DF) |
|---|
| 163 | CHARACTER (LEN = 16) :: CFPPHY(JPOSSGP) |
|---|
| 164 | CHARACTER (LEN = 16) :: CFPCFU(JPOSCFU) |
|---|
| 165 | CHARACTER (LEN = 16) :: CFPXFU(JPOSXFU) |
|---|
| 166 | CHARACTER (LEN = 12*JPOS3DF) :: C1FP3DF |
|---|
| 167 | CHARACTER (LEN = 12*JPOS3DF) :: C1FP3DFS |
|---|
| 168 | CHARACTER (LEN = 12*JPOS3DF) :: C1FP3DFP |
|---|
| 169 | CHARACTER (LEN = 12*JPOS3DF) :: C1FP3DFH |
|---|
| 170 | CHARACTER (LEN = 12*JPOS3DF) :: C1FP3DFT |
|---|
| 171 | CHARACTER (LEN = 12*JPOS3DF) :: C1FP3DFV |
|---|
| 172 | CHARACTER (LEN = 16*JPOS2DF) :: C1FP2DF |
|---|
| 173 | CHARACTER (LEN = 16*JPOSSGP) :: C1FPPHY |
|---|
| 174 | CHARACTER (LEN = 16*JPOSCFU) :: C1FPCFU |
|---|
| 175 | CHARACTER (LEN = 16*JPOSXFU) :: C1FPXFU |
|---|
| 176 | CHARACTER (LEN = JPOSLEN*JPOSDOM) :: C1FPDOM |
|---|
| 177 | INTEGER(KIND=JPIM) :: MFP3DFS(JPOS3DF) |
|---|
| 178 | INTEGER(KIND=JPIM) :: MFP3DFP(JPOS3DF) |
|---|
| 179 | INTEGER(KIND=JPIM) :: MFP3DFH(JPOS3DF) |
|---|
| 180 | INTEGER(KIND=JPIM) :: MFP3DFT(JPOS3DF) |
|---|
| 181 | INTEGER(KIND=JPIM) :: MFP3DFV(JPOS3DF) |
|---|
| 182 | INTEGER(KIND=JPIM) :: MFP2DF(JPOS2DF) |
|---|
| 183 | INTEGER(KIND=JPIM) :: MFPPHY(JPOSSGP) |
|---|
| 184 | |
|---|
| 185 | REAL(KIND=JPRB) :: RFP3P(JPOS3P) |
|---|
| 186 | REAL(KIND=JPRB) :: RFP3H(JPOS3H) |
|---|
| 187 | REAL(KIND=JPRB) :: RFP3TH(JPOS3TH) |
|---|
| 188 | REAL(KIND=JPRB) :: RFP3PV(JPOS3PV) |
|---|
| 189 | INTEGER(KIND=JPIM) :: NRFP3S(JPOS3S) |
|---|
| 190 | |
|---|
| 191 | INTEGER(KIND=JPIM) :: NFPDOM |
|---|
| 192 | INTEGER(KIND=JPIM) :: NFP3DF |
|---|
| 193 | INTEGER(KIND=JPIM) :: NFP3DFS |
|---|
| 194 | INTEGER(KIND=JPIM) :: NFP3DFP |
|---|
| 195 | INTEGER(KIND=JPIM) :: NFP3DFH |
|---|
| 196 | INTEGER(KIND=JPIM) :: NFP3DFT |
|---|
| 197 | INTEGER(KIND=JPIM) :: NFP3DFV |
|---|
| 198 | INTEGER(KIND=JPIM) :: NFP2DF |
|---|
| 199 | INTEGER(KIND=JPIM) :: NFPPHY |
|---|
| 200 | INTEGER(KIND=JPIM) :: NFPCFU |
|---|
| 201 | INTEGER(KIND=JPIM) :: NFPXFU |
|---|
| 202 | INTEGER(KIND=JPIM) :: NFP3P |
|---|
| 203 | INTEGER(KIND=JPIM) :: NFP3H |
|---|
| 204 | INTEGER(KIND=JPIM) :: NFP3TH |
|---|
| 205 | INTEGER(KIND=JPIM) :: NFP3PV |
|---|
| 206 | INTEGER(KIND=JPIM) :: NFP3S |
|---|
| 207 | INTEGER(KIND=JPIM) :: NFPGRIB |
|---|
| 208 | INTEGER(KIND=JPIM) :: NFPCLI |
|---|
| 209 | INTEGER(KIND=JPIM) :: MFP3DYN |
|---|
| 210 | INTEGER(KIND=JPIM) :: MFP2DYN |
|---|
| 211 | INTEGER(KIND=JPIM) :: NFPLNPR |
|---|
| 212 | INTEGER(KIND=JPIM) :: NFPINCR |
|---|
| 213 | INTEGER(KIND=JPIM) :: NFPINDYN |
|---|
| 214 | INTEGER(KIND=JPIM) :: NFPINPHY |
|---|
| 215 | INTEGER(KIND=JPIM) :: NFPCAPE |
|---|
| 216 | REAL(KIND=JPRB) :: WSXI |
|---|
| 217 | REAL(KIND=JPRB) :: WDXI |
|---|
| 218 | REAL(KIND=JPRB) :: WSXO |
|---|
| 219 | REAL(KIND=JPRB) :: WDXO |
|---|
| 220 | REAL(KIND=JPRB) :: FPBL |
|---|
| 221 | REAL(KIND=JPRB) :: RFPCORR |
|---|
| 222 | REAL(KIND=JPRB) :: RFPCSAB |
|---|
| 223 | REAL(KIND=JPRB) :: RFPCD2 |
|---|
| 224 | REAL(KIND=JPRB) :: RFPVCAP |
|---|
| 225 | LOGICAL :: LFPSPEC |
|---|
| 226 | LOGICAL :: LFITP |
|---|
| 227 | LOGICAL :: LFITT |
|---|
| 228 | LOGICAL :: LFITV |
|---|
| 229 | LOGICAL :: LFPQ |
|---|
| 230 | LOGICAL :: LTRACEFP |
|---|
| 231 | LOGICAL :: LFPCNT |
|---|
| 232 | LOGICAL :: LASQ |
|---|
| 233 | LOGICAL :: LFPNHPD |
|---|
| 234 | LOGICAL :: LFPNHVD |
|---|
| 235 | LOGICAL :: LFPNHVW |
|---|
| 236 | LOGICAL :: LFPMOIS |
|---|
| 237 | LOGICAL :: LFPLOSP |
|---|
| 238 | LOGICAL :: LFPRH100 |
|---|
| 239 | LOGICAL :: LMOCONVAR |
|---|
| 240 | INTEGER(KIND=JPIM) :: NFPLAKE |
|---|
| 241 | INTEGER(KIND=JPIM) :: NFPXLEV |
|---|
| 242 | INTEGER(KIND=JPIM) :: NFPDPHY |
|---|
| 243 | INTEGER(KIND=JPIM) :: NFPSURFEX |
|---|
| 244 | INTEGER(KIND=JPIM) :: NFPMASK |
|---|
| 245 | ! ------------------------------------------------------------------ |
|---|
| 246 | !$OMP THREADPRIVATE(c1fp2df,c1fp3df,c1fp3dfh,c1fp3dfp,c1fp3dfs,c1fp3dft,c1fp3dfv,c1fpcfu,c1fpdom,c1fpphy,c1fpxfu) |
|---|
| 247 | !$OMP THREADPRIVATE(cfp2df,cfp3df,cfp3dfh,cfp3dfp,cfp3dfs,cfp3dft,cfp3dfv,cfpcfu,cfpdir,cfpdom,cfpfmt,cfpiden) |
|---|
| 248 | !$OMP THREADPRIVATE(cfpphy,cfpxfu,fpbl,lasq,lfitp,lfitt,lfitv,lfpcnt,lfplosp,lfpmois,lfpnhpd,lfpnhvd,lfpnhvw) |
|---|
| 249 | !$OMP THREADPRIVATE(lfpq,lfprh100,lfpspec,lmoconvar,ltracefp,mfp2df,mfp2dyn,mfp3dfh,mfp3dfp,mfp3dfs,mfp3dft) |
|---|
| 250 | !$OMP THREADPRIVATE(mfp3dfv,mfp3dyn,mfpphy,nfp2df,nfp3df,nfp3dfh,nfp3dfp,nfp3dfs,nfp3dft,nfp3dfv,nfp3h,nfp3p) |
|---|
| 251 | !$OMP THREADPRIVATE(nfp3pv,nfp3s,nfp3th,nfpcape,nfpcfu,nfpcli,nfpdom,nfpdphy,nfpgrib,nfpincr,nfpindyn,nfpinphy) |
|---|
| 252 | !$OMP THREADPRIVATE(nfplake,nfplnpr,nfpmask,nfpphy,nfpsurfex,nfpxfu,nfpxlev,nrfp3s,rfp3h,rfp3p,rfp3pv,rfp3th) |
|---|
| 253 | !$OMP THREADPRIVATE(rfpcd2,rfpcorr,rfpcsab,rfpvcap,wdxi,wdxo,wsxi,wsxo) |
|---|
| 254 | END MODULE YOMFPC |
|---|
