yomfpc.F90

Go to the documentation of this file.

MODULE yomfpc

USE parkind1  ,ONLY : jpim     ,jprb

USE parfpos, ONLY : jpos2df  ,jpos3df  ,jpos3h   ,jpos3p   ,jpos3pv  , &
 & jpos3s   ,jpos3th  ,jposcfu  ,jposdir  ,jposlen  ,jposdom  , &
 & jpossgp  ,jposxfu  

IMPLICIT NONE

SAVE

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

!*    Scientific and technical variables

! === TECHNICAL VARIABLES ===

!     CFPDIR : prefix (path) for the output files
!     CFPIDEN: identificator of the output files
!     CFPFMT : format of the output files ('MODEL','GAUSS','LELAM' or 'LALON')
!     CFPDOM, C1FPDOM : names of the subdomains
!     CFP3DF, C1FP3DF : names 3D dynamics fields to compute
!     CFP3DFS,C1FP3DFS: names 3D dynamics fields to compute on model levels
!     CFP3DFP,C1FP3DFP: names 3D dynamics fields to compute on Pressure levels
!     CFP3DFH,C1FP3DFH: names 3D dynamics fields to compute on H levels
!     CFP3DFT,C1FP3DFT: names 3D dynamics fields to compute on THETA levels
!     CFP3DFV,C1FP3DFV: names 3D dynamics fields to compute on PV levels
!     CFP2DF, C1FP2DF : names 2D dynamics fields to compute
!     CFPPHY, C1FPPHY : names of physical fields to be post-processed
!     CFPCFU, C1FPCFU : names of cumulated fluxes fields to be post-processed
!     CFPXFU, C1FPXFU : names of instantaneous fluxes fields to be post-proc
!     Remark: C1FPXXX=CFPXXX(1)//CFPXXX(2)//...//CFPXXX(JPOSXXX)
!      for XXX=DOM,3DF,2DF,PHY,CFU,XFU.

!     Variables used for ECMWF 

!     MFP3DFS: Gribcodes of  3D dynamics fields to compute on model levels
!     MFP3DFH: Gribcodes of  3D dynamics fields to compute on H levels
!     MFP3DFT: Gribcodes of  3D dynamics fields to compute on THETA levels
!     MFP3DFV: Gribcodes of  3D dynamics fields to compute on PV levels
!     MFP3DFP: Gribcodes of  3D dynamics fields to compute on P levels
!     MFP2DF : Gribcodes of  2D dynamics fields to compute
!     MFPPHY : Gribcodes of  physical fields to processed 

!     RFP3P  : post-processing pressure levels
!     RFP3H  : post-processing height (above orography) levels
!     RFP3TH : post-processing potential temperature levels
!     RFP3PV : post-processing potential vorticity levels
!     NRFP3S : post-processing eta levels (CONF. 927 only)

!     LFPCNT : control varying output variables according to time step
!     LTRACEFP: trace for Full-POS
!     NFPGRIB : level of GRIB coding in output file ARPEGE/ALADIN.
!               0 : no packing at all
!               1 : standart GRIB encoding
!               2 : modified GRIB encoding 

!     NFPDOM : useful dimension of CFPDOM
!     NFP3DF : useful dimension of CFP3DF
!     NFP3DFS: useful dimension of CFP3DFS
!     NFP3DFH: useful dimension of CFP3DFH
!     NFP3DFT: useful dimension of CFP3DFT
!     NFP3DFV: useful dimension of CFP3DFV
!     NFP3DFP: useful dimension of CFP3DFP
!     NFP2DF : useful dimension of CFP2DF
!     NFPPHY : useful dimension of CFPPHY
!     NFPCFU : useful dimension of CFPCFU
!     NFPXFU : useful dimension of CFPXFU
!     NFP3P  : useful dimension of RFP3P
!     NFP3H  : useful dimension of RFP3H
!     NFP3TH : useful dimension of RFP3TH
!     NFP3PV : useful dimension of RFP3PV
!     NFP3S  : useful dimension of NRFP3S
!     NFPXLEV: max. number of pp. levels
!     MFP3DYN : max. number of 3D-dynamic fields neeeded for LFPART2
!     MFP2DYN : max. number of 2D-dynamic fields neeeded for LFPART2
!     NFPDPHY : maximum number of physical fields

!     The next keys force variables to be primitive for the post-processing :
!     LFPNHPD : .TRUE. if pp of NH pressure departure
!     LFPNHVD : .TRUE. if pp of NH vertical divergence
!     LFPNHVW : .TRUE. if pp of NH true vertical velocity

!     LFPLOSP : .TRUE. = Fill Ps array with Log(Ps)
!     LFPRH100 : .TRUE. to convert relative humidity in percent

! === SCIENTIFIC VARIABLES ===

!     LFPSPEC: =.T. if pp. dyn. fields are written out as spectral coefficients
!              =.F. if pp. dyn. fields are written out as grid point values.
!     LFITP  : =1 if pp. fields on P levels should be fited ; =0 otherwise
!     LFITT  : =1 if pp. fields on THETA levels should be fited ; =0 otherwise
!     LFITV  : =1 if pp. fields on PV levels should be fited ; =0 otherwise

!     NFPCLI : usage level for climatology
!              =0 no climatology
!              =1 orography and land-sea mask of output only
!              =2 all available climatological fields of the current month
!              =3 shifting mean from the climatological fields of the current
!                 month to the ones of the closest month

!     NFPINDYN : type of interpolations for dynamical fields
!     NFPINPHY : type of interpolations for  physical fields
!                4 for bilinear, 12 for quadratic

!     LFPQ   : =.TRUE. if specific humidity is interpolated
!            : =.FALSE. if relative humidity is interpolated
!     LASQ   : =.TRUE. if ASQ set to 80% saturation

!     WSXI   : max. surface moisture in input
!     WDXI   : max. deep soil moisture in input
!     WSXO   : max. surface moisture in output
!     WDXO   : max. deep soil moisture in output

!     FPBL   : Critical Thickness of PBL
!     RFPCORR: Critical orography difference for correcting surface
!              temperature through standart profile.
!     RFPCSAB: Critical sand percentage difference for computing relative 
!              soil moisture in ISBA
!     RFPCD2 : Critical soil depth difference for computing relative
!              soil moisture in ISBA
!     RFPVCAP: Minimum pressure of model level to provide an equatorial
!              cap in the computation of variables on constant PV surfaces
!     NFPLNPR : 0 => conventional formulation of (delta P) : ln(P(l)/P(l-1))
!               1 => formulation of (delta P) used in non hydrostatic model,
!                    i.e. (P(l)-P(l-1))/SQRT(P(l)*P(l-1))
!     LFPMOIS: Month allowed for climatology usage : 
!               .F. => month of the model (forecast)
!               .T. => month of the file
!     NFPINCR : Bogussing indicator ; =0 for no bogussing ; =1 for bogussing
!     NFPLAKE : To overwrite created lakes or islands by specific data :
!               0 => do not overwrite
!              -1 => overwrite with rhoughly interpolated data
!              +1 => overwrite with climatology
!     NFPCAPE : Kind of computation for CAPE & CIN : 
!               1 => from bottom model layer
!               2 => from the most unstable layer
!               3 => from mto standart height (2 meters) as recomputed values
!               4 => from mto standart height (2 meters) out of fluxes
!                    (used for analysis)
!     NFPSURFEX : Subcontract surface fields to SURFEX
!               0 => no subcontract
!               1 => transform native arp/ald surface fields to surfex fields
!                    and write out by surfex
!     NFPMASK   : number of masks for the interpolation of surface fields
!                 0 => no mask
!                 1 => land-sea mask
!                 2 => land mask, sea mask
!     LMOCONVAR  : if true mocon is computed from N-1 level parameters (for varpack use)

CHARACTER (LEN = JPOSDIR) ::  cfpdir
CHARACTER (LEN = 5) ::  cfpfmt
CHARACTER (LEN = 8) ::  cfpiden
CHARACTER (LEN = JPOSLEN) ::  cfpdom(jposdom)
CHARACTER (LEN = 12) ::  cfp3df(jpos3df)
CHARACTER (LEN = 12) ::  cfp3dfs(jpos3df)
CHARACTER (LEN = 12) ::  cfp3dfp(jpos3df)
CHARACTER (LEN = 12) ::  cfp3dfh(jpos3df)
CHARACTER (LEN = 12) ::  cfp3dft(jpos3df)
CHARACTER (LEN = 12) ::  cfp3dfv(jpos3df)
CHARACTER (LEN = 16) ::  cfp2df(jpos2df)
CHARACTER (LEN = 16) ::  cfpphy(jpossgp)
CHARACTER (LEN = 16) ::  cfpcfu(jposcfu)
CHARACTER (LEN = 16) ::  cfpxfu(jposxfu)
CHARACTER (LEN = 12*JPOS3DF) :: c1fp3df
CHARACTER (LEN = 12*JPOS3DF) :: c1fp3dfs
CHARACTER (LEN = 12*JPOS3DF) :: c1fp3dfp
CHARACTER (LEN = 12*JPOS3DF) :: c1fp3dfh
CHARACTER (LEN = 12*JPOS3DF) :: c1fp3dft
CHARACTER (LEN = 12*JPOS3DF) :: c1fp3dfv
CHARACTER (LEN = 16*JPOS2DF) :: c1fp2df
CHARACTER (LEN = 16*JPOSSGP) :: c1fpphy
CHARACTER (LEN = 16*JPOSCFU) :: c1fpcfu
CHARACTER (LEN = 16*JPOSXFU) :: c1fpxfu
CHARACTER (LEN = JPOSLEN*JPOSDOM) :: c1fpdom
INTEGER(KIND=JPIM) :: mfp3dfs(jpos3df)
INTEGER(KIND=JPIM) :: mfp3dfp(jpos3df)
INTEGER(KIND=JPIM) :: mfp3dfh(jpos3df)
INTEGER(KIND=JPIM) :: mfp3dft(jpos3df)
INTEGER(KIND=JPIM) :: mfp3dfv(jpos3df)
INTEGER(KIND=JPIM) :: mfp2df(jpos2df)
INTEGER(KIND=JPIM) :: mfpphy(jpossgp)

REAL(KIND=JPRB) :: rfp3p(jpos3p)
REAL(KIND=JPRB) :: rfp3h(jpos3h)
REAL(KIND=JPRB) :: rfp3th(jpos3th)
REAL(KIND=JPRB) :: rfp3pv(jpos3pv)
INTEGER(KIND=JPIM) :: nrfp3s(jpos3s)

INTEGER(KIND=JPIM) :: nfpdom
INTEGER(KIND=JPIM) :: nfp3df
INTEGER(KIND=JPIM) :: nfp3dfs
INTEGER(KIND=JPIM) :: nfp3dfp
INTEGER(KIND=JPIM) :: nfp3dfh
INTEGER(KIND=JPIM) :: nfp3dft
INTEGER(KIND=JPIM) :: nfp3dfv
INTEGER(KIND=JPIM) :: nfp2df
INTEGER(KIND=JPIM) :: nfpphy
INTEGER(KIND=JPIM) :: nfpcfu
INTEGER(KIND=JPIM) :: nfpxfu
INTEGER(KIND=JPIM) :: nfp3p
INTEGER(KIND=JPIM) :: nfp3h
INTEGER(KIND=JPIM) :: nfp3th
INTEGER(KIND=JPIM) :: nfp3pv
INTEGER(KIND=JPIM) :: nfp3s
INTEGER(KIND=JPIM) :: nfpgrib
INTEGER(KIND=JPIM) :: nfpcli
INTEGER(KIND=JPIM) :: mfp3dyn
INTEGER(KIND=JPIM) :: mfp2dyn
INTEGER(KIND=JPIM) :: nfplnpr
INTEGER(KIND=JPIM) :: nfpincr
INTEGER(KIND=JPIM) :: nfpindyn
INTEGER(KIND=JPIM) :: nfpinphy
INTEGER(KIND=JPIM) :: nfpcape
REAL(KIND=JPRB) :: wsxi
REAL(KIND=JPRB) :: wdxi
REAL(KIND=JPRB) :: wsxo
REAL(KIND=JPRB) :: wdxo
REAL(KIND=JPRB) :: fpbl
REAL(KIND=JPRB) :: rfpcorr
REAL(KIND=JPRB) :: rfpcsab
REAL(KIND=JPRB) :: rfpcd2
REAL(KIND=JPRB) :: rfpvcap
LOGICAL :: lfpspec
LOGICAL :: lfitp
LOGICAL :: lfitt
LOGICAL :: lfitv
LOGICAL :: lfpq
LOGICAL :: ltracefp
LOGICAL :: lfpcnt
LOGICAL :: lasq
LOGICAL :: lfpnhpd
LOGICAL :: lfpnhvd
LOGICAL :: lfpnhvw
LOGICAL :: lfpmois
LOGICAL :: lfplosp
LOGICAL :: lfprh100
LOGICAL :: lmoconvar
INTEGER(KIND=JPIM) :: nfplake
INTEGER(KIND=JPIM) :: nfpxlev
INTEGER(KIND=JPIM) :: nfpdphy
INTEGER(KIND=JPIM) :: nfpsurfex
INTEGER(KIND=JPIM) :: nfpmask
!     ------------------------------------------------------------------
!$OMP THREADPRIVATE(c1fp2df,c1fp3df,c1fp3dfh,c1fp3dfp,c1fp3dfs,c1fp3dft,c1fp3dfv,c1fpcfu,c1fpdom,c1fpphy,c1fpxfu)
!$OMP THREADPRIVATE(cfp2df,cfp3df,cfp3dfh,cfp3dfp,cfp3dfs,cfp3dft,cfp3dfv,cfpcfu,cfpdir,cfpdom,cfpfmt,cfpiden)
!$OMP THREADPRIVATE(cfpphy,cfpxfu,fpbl,lasq,lfitp,lfitt,lfitv,lfpcnt,lfplosp,lfpmois,lfpnhpd,lfpnhvd,lfpnhvw)
!$OMP THREADPRIVATE(lfpq,lfprh100,lfpspec,lmoconvar,ltracefp,mfp2df,mfp2dyn,mfp3dfh,mfp3dfp,mfp3dfs,mfp3dft)
!$OMP THREADPRIVATE(mfp3dfv,mfp3dyn,mfpphy,nfp2df,nfp3df,nfp3dfh,nfp3dfp,nfp3dfs,nfp3dft,nfp3dfv,nfp3h,nfp3p)
!$OMP THREADPRIVATE(nfp3pv,nfp3s,nfp3th,nfpcape,nfpcfu,nfpcli,nfpdom,nfpdphy,nfpgrib,nfpincr,nfpindyn,nfpinphy)
!$OMP THREADPRIVATE(nfplake,nfplnpr,nfpmask,nfpphy,nfpsurfex,nfpxfu,nfpxlev,nrfp3s,rfp3h,rfp3p,rfp3pv,rfp3th)
!$OMP THREADPRIVATE(rfpcd2,rfpcorr,rfpcsab,rfpvcap,wdxi,wdxo,wsxi,wsxo)
END MODULE yomfpc