yomct0.F90

Go to the documentation of this file.

MODULE yomct0

USE parkind1  ,ONLY : jpim     ,jprb

IMPLICIT NONE

SAVE

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

!*    Control variables for the job - constant within job

!========== TECHNICAL SWITCHES ================================================

! ----- configuration:
! NCONF      : configuration of the job

!                0- 99 : 3-D integration job
!              100-199 : variational job
!              200-299 : 2-D integration job
!              300-349 : KALMAN filter                    (currently obsolete)
!              350-399 : predictability model             (currently obsolete)
!              400-499 : test of the adjoint
!              500-599 : test of the tangent linear model
!              600-699 : eigenvalue/vector solvers
!              700-799 : optimal interpolation
!              800-899 : sensitivity
!              900-999 : miscellaneous other configurations.

!                    1 : 3-D primitive equation model

!                  131 : incremental 4-D VAR/3-D VAR

!                  201 : shallow-water model
!                  202 : vorticity equation model

!                  401 : test of adjoint with 3-D P.E. model
!                  421 : test of adjoint with shallow-water model
!                  422 : test of adjoint with vorticity equation model

!                  501 : test of tangent linear with 3-D P.E. model
!                  521 : test of tangent linear with shallow-water model
!                  522 : test of tangent linear with vorticity equation model

!                  601 : eigenvalue/vector solver for 3-D P.E. model

!                  701 : optimal interpolation with CANARI

!                  801 : sensitivity with 3-D P.E. model

!                  901 : set up initial conditions (CPREP1)
!                  903 : set up initial conditions (CPREP5)
!                  911 : computes dilatation matrices for spectral space
!                  912 : computes rotation matrices for spectral space
!                  923 : initialisation of climatologic files
!                  931 : creation of an ARPEGE file containing the NESDIS SST (INCLITC).
!                  932 : interpolates the sea-ice concentration field from
!                        satellite data to the ARPEGE grid (CSEAICE).
!                  940 : corrects the dry surface pressure in an ARPEGE file
!                        to keep total mass constant compared to a reference
!                        ARPEGE file (CORMASS)
!                  951 : difference between 2 model states (CPREP2)

! ----- ???:
! NTASKS     : ??? (no comment provided, currently used in CANARI only).

! ----- variables linked with diagnostics and outputs:
! LALLOPR    : .T. = print information about all allocations/deallocations
! NPRINTLEV  : 0 = "basic prints'; 1 = "more prints"; 2 = "debug prints"

! ----- type of file used:
! LFDBOP     : .T. = fields data base utilized
! LGRBOP     : .T. = output in GRIB (not ARPEGE)
! LFBDAP     : .T. = diagnostics written on trajectory files
! LARPEGEF   : .T. = use ARPEGE files
! LARPEGEF_TRAJHR   : .T. = use ARPEGE files for high resolution trajectory
! LARPEGEF_TRAJBG   : .T. = use ARPEGE files for background
! LARPEGEF_RDGP_INIT     : .T. = use grid-point ARPEGE files
! LARPEGEF_RDGP_TRAJHR   : .T. = use grid-point ARPEGE files for HR trajectory
! LARPEGEF_RDGP_TRAJBG   : .T. = use grid-point ARPEGE files for background
! CNDISPP    : directory of display files

! ----- variables linked to post-processing:
! LFPOS      : .T. = use of Full-POS rather than POS
! CFPNCF     : name of Full-POS control file (pseudo-configuration 927)
! LFPART2    : .T. = second part of interpolations for changing geometry
!              (pseudo-configuration 927)
! CFDIRLST   : path of postprocessing listing file
! CNPPATH    : directory of postprocessing namelist files

! ----- use of transmission coefficients stored in Fourier or spectral space:
! LRETCFOU   : .T. = reading Fourier transmission coefficients on file
!              for simplified physics.
! LWRTCFOU   : .T. = writing Fourier transmission coefficients on file.
! LRFOUTCNORM: activates diagnostics on Fourier transmission coefficients.
!              - LRFOUTCNORM(1)=.T.: global statistics.
!              - LRFOUTCNORM(2)=.T.: statistics per latitude.
!              - LRFOUTCNORM(3)=.T.: statistics per layer.
!              - LRFOUTCNORM(4)=.T.: statistics per wavenumber.
! LRGPTCNORM : activates diagnostics on grid-point transmission coefficients.
!              - LRGPTCNORM(1)=.T.: global statistics.
!              - LRGPTCNORM(2)=.T.: statistics per latitude.
!              - LRGPTCNORM(3)=.T.: statistics per layer.

! ----- other variables:
! LNF        : .T. = start, .F. = restart
! LSMSSIG    : .T. = send signals to SMS (ECMWF supervisor)
! LOPDIS     : .T. = calls OPDIS
! NCYCLE     : cycle identifier
! CNMEXP     : name of the experiment
!              An experiment is identified by its name (16 characters)
!              and its cycle (typically same experiment but without a bug)
! CFCLASS    : class for use by FDB
! CTYPE      : type for use by FDB
! LBACKG     : ???
! LMINIM     : ???

!========== MODEL SWITCHES ====================================================
! Remark: this section containes some dynamical variables which cannot be
!         put in NAMDYN because they are used in routines called before SUDYN
!         and sometimes to allocate some arrays.

! ----- advection scheme, semi-Lagrangian scheme:
! LSLAG     : .TRUE. = semi-lagrangian on
! LTWOTL    : .TRUE. if two-time-level semi-Lagrangian scheme will be used
! LRFRIC    : .TRUE. = Rayleigh friction in zonal wind (for p < 9.9hPa)
! LVERCOR   : .T./.F.: thin layer hypothesis relaxed/applied.

! ----- vertical discretisation, vertical boundaries:
! LAPRXPK   : way of computing full-levels pressures in primitive equation
!             hydrostatic model.
!             .T.: full levels are computed by PK=(PK+1/2 + PK-1/2)*0.5
!             .F.: full levels are computed by a more complicated formula
!                  consistent with "alpha" in geopotential formula.
! LREGETA   : .T.: for the interlayer L, ETA(L)=L/NFLEVG
!             .F.: for the interlayer L, ETA(L)=A(L)/P0+B(L)
! LRUBC     : .T.: if radiational upper boundary condition

! ----- numbering of timesteps:
! NSTART    : first timestep of model
! NSTOP     : last timestep of model

! ----- quadrature:
! NQUAD     : 1 ====> GAUSS

! ----- type of equations:
! LNHDYN    : .T. if 3-D non-hydrostatic dynamics is active

! ----- way of computing initial state:
! N2DINI    : 1 initialization for 2D with Haurwitz wave
!           : 2 initialization for 2D with real fields
! N3DINI    : 1 initialization for 3D with standard atmosphere
!             (not available since cycle 12)
!           : 0 initialization for 3D with real fields

! ----- semi-implicit scheme:
! LSITRIC   : .T.: if "tridiagonal" solver is used for the vertically
!             coupled semi-implicit equations

! ----- control of variables which are transformed in spectral space:
! LSPRT     : .T.: if R*T/Rd "virtual temperature" as spectral variable

! ----- diagnostics and frequencies:
! NSPPR     : 0: no spectrum printed in spnorm; only global norms averaged
!             on the vertical are printed
!           : 1: no spectrum printed in spnorm; only global norms averaged
!             on the vertical and global norms on each layer are printed
!           : 2: both total wavenumber spectrum and zonal wavenumber spectrum
!             are printed
! NFRPOS    : frequency of post-processing events (time-steps)
! NFRISP    : frequency of isp (animation !) events (time-steps)
! NFRCO     : frequency of coupled fields (time-steps)
! NFRCORM   : mass correction frequency (>0 time-steps, <0 hours, 0 no
!             correction)
! NFRHIS    : frequency of history write_ups (time-steps)
! NFRMASSCON: frequency of mass conservation fixup (time-steps)
! NFRGDI    : frequency of grid-point space diagnostics
! NFRSDI    : frequency of spectral space diagnostics
! NFRDHFG   : write-up frequency of global DDH
! NFRDHFZ   : write-up frequency of zonal DDH
! NFRDHFD   : write-up frequency of "limited domain" DDH
! NFRDHP    : write-up frequency of DDH files
! NPOSTS    : array containing postprocessing steps
! NPISPS    : array containing isp (animation !) steps
! NHISTS    : array containing history write-up steps
! NMASSCONS : array containing mass conservation fixup steps
! NGDITS    : array containing grid point diagnostics steps
! NSDITS    : array containing spectral diagnostics steps
! NDHFGTS   : array containing write out steps for global DDH
! NDHFZTS   : array containing write out steps for zonal means DDH
! NDHFDTS   : array containing write out steps for limited areas DDH
! NDHPTS    : array containing write out steps for DDH
! Explanation for N[XXX]TS:
!             1) if N[XXX]TS(0)=0 action if MOD(JSTEP,NFR[XXX])=0
!             2) if N[XXX]TS(0)>0 N[XXX]TS(0) significant numbers in 
!                N[XXX]TS are then considered and:
!                action for JSTEP=N[XXX]TS(.)*NFR[XXX]
!             3) IF N[XXX]TS(0)<0
!                action for JSTEP=(N[XXX]TS(.)/DELTAT)*NFR[XXX]

!========== ASSIMILATION SWITCHES =============================================

! LNOBGON : .F. if eq. .T. no stop if unsuccessful reading of CMA files
! LCANARI : .T. = term to control French OI
! LCASIG  : .T. = interpolation of the model errors (French OI)
! LGUESS  : .T. = term of first guess included
! LOBS    : .T. = term of observations included
! LSIMOB  : .T. = if simulated observations
! LOBSC1  : .T. = term of observations included in configuration 1
! LSCREEN : .T. = observation screening for variational assimilation
! LSCREEN_OPENMP : .T. = 4DVAR screening runs in OpenMP-parallel mode over timeslots
! L_SPLIT_SCREEN .T. = to split screenng
! L_SCREEN_CALL: .T. = call to screening routine SCREEN
! LOBSREF : .T. = comparison to observation for the trajectory (NCONF=131)
! LIFSMIN : .T. = if running minimisation
! LIFSTRAJ: .T. = if running high resolution trajectory integration
! NCNTVAR : Definition of the control variable of a variational job.
!           = 1 ===> control variables are model variables
!           = 2 ===> control variables are normalized departures of
!                    model variables from the background field
!           = 3 ===> ..........
! LOLDPP  : .T. use "old" p.p. of T,Q,U and V
! NSTEPINI: Initial step in hours for the initial conditions
!           at the beginning of 4D-Var trajectory (usually 3 hours).
!           It is used to update the step while saving the FCs along
!           the first trajectory.
! NINTERPTRAJ : Interpolation method applied to increments
! NINTERPINCR : Interpolation method applied to increments
!           = 1 ===> Bi-linear interpolation (default)
!           = 2 ===> Bi-cubic Rinterpolation
!           = 3 ===> Conserving interpolation
! For Conserving interpolation the style is defined by:
! NINTERPTRAJORDER : order of interpolation for trajectories
! NINTERPINCRORDER : order of interpolation for increments
! NINTERPTRAJLIMIT : 0 if limiter is NOT used (default), 1 if used for trajectory interpolation
! NINTERPINCRLIMIT : 0 if limiter is NOT used (default), 1 if used for increment interpolation

!========== ALADIN SWITCHES ===================================================

! LELAM   : .T. = limited area model with coupling or fully biperiodic
!                 model (torus)
!           .F. = global model (sphere)
! LRPLANE : .T. = plane cartesian geometry
!           .F. = spherical geometry
! LTENC   : TRUE if tendency coupling of surface pressure
!           FALSE if no tend. coupling of surf. pressure
! LALLTC  : used together with LTENC when LTENC=.T.
!          - no meaning for quadratic tendency coupling, where just t1 coupling
!          is applied at every NEFRCL time step
!          - for lin. tendency coupling:
!          TRUE if tendency coupling of surf. pres. at every step
!          FALSE if tend. coupl., except at every NEFRCL time steps
!          when just t1 coupling
! RTENC : multiplier of EALFA in the tendency coupling scheme
!  for stability reasons (RTENC<=1. close to 1)

!========== AROME SWITCH ======================================================

! LAROME  : .T. = AROME limited area model

!========== ECMWF Single Column Model =========================================
! LSCMEC  : .T. = ECMWF Single Column Model
! LSFCFLX : .T. = forcing with surface fluxes (latent and sensible).
! REXTSHF : externally supplied sensible heat flux [W/m^2]
! REXTLHF : externally supplied latent   heat flux [W/m^2]
! LROUGH  : .T. = surface roughness length is externally specified 
! REXTZ0M : externally supplied roughness length for momentum [m]
! REXTZ0H : externally supplied roughness length for heat [m]

!========== DISTRIBUTED MEMORY SWITCHES =======================================

! LMPOFF  : .T. = switch off the message passing library initialisation
!                 requested for special cases LMESSP=.T. and NPROC=1
!         : .F. = (default) full message passing features
! NPROC   : Total number of processors requested for this run
! N_REGIONS_NS : Number of regions (LEQ_REGIONS=T) or NPRGPNS (LEQ_REGIONS=F)
! N_REGIONS_EW : Maximum number of partitions for all regions
! N_REGIONS : Number of partitions in each region
! NPRGPNS : Number of processors used during grid-point phase in North-South
!           direction (previously known as NPROCA)
! NPRGPEW : Number of processors used during grid-point phase in East-West
!           direction (previously known as NPROCB)
! NPRTRNS : Number of processors used during transform phase in North-South
!           direction (previously the same as NPROCA and now implemented
!           such that NPRTRNS=NPTRW)
! NPRTRM  : Number of processors used during spectral computations in
!           meridional wave direction (previously the same as NPROCA
!           and now implemented such that NPRTRM=NPTRW)
! NPRTRN  : Number of processors used during spectral computations in total
!           wave direction (previously the same as NPROCA and now implemented
!           such that NPRTRN=NPTRV)
! NPRTRW  : Number of processors used during transform phase in wave space
!           (previously the same as NPROCA)
! NPRTRV  : Number of processors used during transform phase in vertical
!           direction (previously known as NPROCB)
! LOUTPUT : .T. = diagnostic output requested on this PE
! NOUTPUT : 0 = No diagnostic output
!           1 = Only diagnostic output from PE1 ( default )
!           2 = Diagnostic output from all PEs into separate files
! LREFOUT : .T. compare to reference run
! LREFGEN : .T. to generate reference file
! LMPDIAG : .T. = extensive message passing diagnostic output requested

!========== PC SCHEMES ========================================================
!----------------------------------------------
! PC SCHEMES QUANTITIES CONSTANT DURING INTEGRATION
!----------------------------------------------

! LPC_FULL - full PC scheme switch (with reiterations of trajectories)

! LPC_OLD - iterative scheme without trajectories recalculations
!         - with iteration of D3 terms (Bubnova et all., 1995)

! LPC_NESC - non-extrapolating two-time level SL SI scheme
!          - X(t+dt/2) = X(t) during predictor

!========== FORCING SWITCH ====================================================
! LSFORC - switch to activate the large scale forcings in setup and cpg
!==============================================================================

! * Parameters:
INTEGER(KIND=JPIM), PARAMETER :: jpnpst=240

! * Technical switches:
INTEGER(KIND=JPIM) :: nconf
INTEGER(KIND=JPIM) :: ntasks ! ??????
LOGICAL :: lallopr
INTEGER(KIND=JPIM) :: nprintlev
LOGICAL :: lfdbop
LOGICAL :: lgrbop
LOGICAL :: lfbdap
LOGICAL :: larpegef
LOGICAL :: larpegef_trajhr
LOGICAL :: larpegef_trajbg
LOGICAL :: larpegef_rdgp_init
LOGICAL :: larpegef_rdgp_trajhr
LOGICAL :: larpegef_rdgp_trajbg
CHARACTER (LEN = 120) ::  cndispp
LOGICAL :: lfpos
CHARACTER (LEN = 6) ::  cfpncf
LOGICAL :: lfpart2
CHARACTER (LEN = 120) ::  cfdirlst
CHARACTER (LEN = 120) ::  cnppath
LOGICAL :: lretcfou
LOGICAL :: lwrtcfou
LOGICAL :: lrfoutcnorm(4)
LOGICAL :: lrgptcnorm(3)
LOGICAL :: lnf
LOGICAL :: lsmssig
LOGICAL :: lopdis
INTEGER(KIND=JPIM) :: ncycle=128
CHARACTER (LEN = 16) ::  cnmexp
CHARACTER (LEN = 2) ::  cfclass
CHARACTER (LEN = 2) ::  ctype
LOGICAL :: lbackg
LOGICAL :: lminim

! * Model switches:
LOGICAL :: lslag
LOGICAL :: ltwotl
LOGICAL :: lrfric
LOGICAL :: lvercor
LOGICAL :: laprxpk
LOGICAL :: lregeta
LOGICAL :: lrubc
INTEGER(KIND=JPIM) :: nstart
INTEGER(KIND=JPIM) :: nstop
INTEGER(KIND=JPIM) :: nquad
LOGICAL :: lnhdyn
INTEGER(KIND=JPIM) :: n2dini
INTEGER(KIND=JPIM) :: n3dini
LOGICAL :: lsitric
LOGICAL :: lsprt
INTEGER(KIND=JPIM) :: nsppr
INTEGER(KIND=JPIM) :: nfrpos
INTEGER(KIND=JPIM) :: nfrisp
INTEGER(KIND=JPIM) :: nfrcorm
INTEGER(KIND=JPIM) :: nfrco
INTEGER(KIND=JPIM) :: nfrhis
INTEGER(KIND=JPIM) :: nfrmasscon
INTEGER(KIND=JPIM) :: nfrgdi
INTEGER(KIND=JPIM) :: nfrsdi
INTEGER(KIND=JPIM) :: nfrdhfg
INTEGER(KIND=JPIM) :: nfrdhfz
INTEGER(KIND=JPIM) :: nfrdhfd
INTEGER(KIND=JPIM) :: nfrdhp
INTEGER(KIND=JPIM) :: nposts(0:jpnpst)
INTEGER(KIND=JPIM) :: npisps(0:jpnpst)
INTEGER(KIND=JPIM) :: nhists(0:jpnpst)
INTEGER(KIND=JPIM) :: nmasscons(0:jpnpst)
INTEGER(KIND=JPIM) :: ngdits(0:jpnpst)
INTEGER(KIND=JPIM) :: nsdits(0:jpnpst)
INTEGER(KIND=JPIM) :: ndhfgts(0:jpnpst)
INTEGER(KIND=JPIM) :: ndhfzts(0:jpnpst)
INTEGER(KIND=JPIM) :: ndhfdts(0:jpnpst)
INTEGER(KIND=JPIM) :: ndhpts(0:jpnpst)

! * Assimilation:
LOGICAL :: lnobgon
LOGICAL :: lcanari
LOGICAL :: lcasig
LOGICAL :: lguess
LOGICAL :: lobs
LOGICAL :: lsimob
LOGICAL :: lobsc1
LOGICAL :: lscreen
LOGICAL :: lscreen_openmp
LOGICAL :: l_split_screen
LOGICAL :: l_screen_call
LOGICAL :: lobsref
LOGICAL :: lifsmin
LOGICAL :: lifstraj
INTEGER(KIND=JPIM) :: ncntvar
LOGICAL :: loldpp
INTEGER(KIND=JPIM) :: nstepini
INTEGER(KIND=JPIM) :: ninterptraj
INTEGER(KIND=JPIM) :: ninterpincr
INTEGER(KIND=JPIM) :: ninterptrajlimit
INTEGER(KIND=JPIM) :: ninterpincrlimit
INTEGER(KIND=JPIM) :: ninterptrajorder
INTEGER(KIND=JPIM) :: ninterpincrorder

! * ALADIN:
LOGICAL :: lelam
LOGICAL :: lrplane
LOGICAL :: ltenc
LOGICAL :: lalltc
REAL(KIND=JPRB) :: rtenc

! * AROME:
LOGICAL :: larome

! * ECMWF Single Column Model:
LOGICAL :: lscmec
LOGICAL :: lsfcflx
REAL(KIND=JPRB) :: rextshf
REAL(KIND=JPRB) :: rextlhf
LOGICAL :: lrough
REAL(KIND=JPRB) :: rextz0m
REAL(KIND=JPRB) :: rextz0h

! * Distributed memory:
LOGICAL :: lmpoff
INTEGER(KIND=JPIM) :: nproc
INTEGER(KIND=JPIM) :: n_regions_ns
INTEGER(KIND=JPIM) :: n_regions_ew
INTEGER(KIND=JPIM),ALLOCATABLE :: n_regions(:)
INTEGER(KIND=JPIM) :: nprgpns
INTEGER(KIND=JPIM) :: nprgpew
INTEGER(KIND=JPIM) :: nprtrns
INTEGER(KIND=JPIM) :: nprtrm
INTEGER(KIND=JPIM) :: nprtrn
INTEGER(KIND=JPIM) :: nprtrw
INTEGER(KIND=JPIM) :: nprtrv
LOGICAL :: loutput
INTEGER(KIND=JPIM) :: noutput
LOGICAL :: lrefout
LOGICAL :: lrefgen
LOGICAL :: lmpdiag

! * PC schemes.
LOGICAL ::  lpc_full
LOGICAL ::  lpc_old
LOGICAL ::  lpc_nesc
!
! * FORCING
LOGICAL ::  lsforc
!     ------------------------------------------------------------------
!$OMP THREADPRIVATE(cfclass,cfdirlst,cfpncf,cndispp,cnmexp,cnppath,ctype,l_screen_call,l_split_screen,lallopr,lalltc)
!$OMP THREADPRIVATE(laprxpk,larome,larpegef,larpegef_rdgp_init,larpegef_rdgp_trajbg,larpegef_rdgp_trajhr)
!$OMP THREADPRIVATE(larpegef_trajbg,larpegef_trajhr,lbackg,lcanari,lcasig,lelam,lfbdap,lfdbop,lfpart2,lfpos)
!$OMP THREADPRIVATE(lgrbop,lguess,lifsmin,lifstraj,lminim,lmpdiag,lmpoff,lnf,lnhdyn,lnobgon,lobs,lobsc1,lobsref)
!$OMP THREADPRIVATE(loldpp,lopdis,loutput,lpc_full,lpc_nesc,lpc_old,lrefgen,lrefout,lregeta,lretcfou,lrfoutcnorm)
!$OMP THREADPRIVATE(lrfric,lrgptcnorm,lrough,lrplane,lrubc,lscmec,lscreen,lscreen_openmp,lsfcflx,lsforc,lsimob)
!$OMP THREADPRIVATE(lsitric,lslag,lsmssig,lsprt,ltenc,ltwotl,lvercor,lwrtcfou,n2dini,n3dini,n_regions_ew,n_regions_ns)
!$OMP THREADPRIVATE(ncntvar,nconf,ncycle,ndhfdts,ndhfgts,ndhfzts,ndhpts,nfrco,nfrcorm,nfrdhfd,nfrdhfg,nfrdhfz,nfrdhp)
!$OMP THREADPRIVATE(nfrgdi,nfrhis,nfrisp,nfrmasscon,nfrpos,nfrsdi,ngdits,nhists,ninterpincr,ninterpincrlimit)
!$OMP THREADPRIVATE(ninterpincrorder,ninterptraj,ninterptrajlimit,ninterptrajorder,nmasscons,noutput,npisps,nposts)
!$OMP THREADPRIVATE(nprgpew,nprgpns,nprintlev,nproc,nprtrm,nprtrn,nprtrns,nprtrv,nprtrw,nquad,nsdits,nsppr,nstart)
!$OMP THREADPRIVATE(nstepini,nstop,ntasks,rextlhf,rextshf,rextz0h,rextz0m,rtenc)
!$OMP THREADPRIVATE(n_regions)
END MODULE yomct0