yomgem.F90

Go to the documentation of this file.

MODULE yomgem

USE parkind1  ,ONLY : jpim     ,jprb

IMPLICIT NONE

SAVE

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

!*    * Number of grid points

!     NGPTOT   : Total number of grid columns on a PE
!     NGPTOT_CAP  : Size of grid points arrays for ALADIN
!     NGPTOTMX : Maximum number of grid columns on any of the PEs
!     NGPTOTG  : Total number of grid columns on the Globe
!     NGPTOTL(NPRGPNS,NPRGPEW)  : Total number of grid columns on on eacch PE

INTEGER(KIND=JPIM) :: ngptot
INTEGER(KIND=JPIM) :: ngptot_cap
INTEGER(KIND=JPIM) :: ngptotmx
INTEGER(KIND=JPIM) :: ngptotg
INTEGER(KIND=JPIM),ALLOCATABLE, TARGET :: ngptotl(:,:)

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

!*    * Defining the transformed sphere

!     RMUCEN : MU OF THE POLE OF STRETCHING
!     RLOCEN : LONGITUDE OF THE POLE OF STRETCHING
!     RSTRET : STRETCHING FACTOR
!     NSTTYP : 1 = POLE OF STRETCHING, POLE OF THE COLLOCATION GRID
!                 AT THE NORTHERN POLE OF THE REAL EARTH.
!              2 = THE POLE OF STRETCHING IS ANYWHERE ON THE REAL EARTH
!             AND ON THE EQUATOR OF THE COLLOCATION GRID ON THE MERIDIAN PI.
!                  THE EQUATOR OF THE COLLOCATION GRID IS TANGENT
!             TO A PARALLEL OF THE EARTH.

!     NHTYP  : 0 = regular grid
!            : 1 = number of points proportional to sqrt(1-mu**2)
!            : 2 = number of points read on namelist namrgri

!     RNLGINC: increment to get non-linear grid

!     R4JP    inverse de delta(teta) approche a l'ordre 1
!     RC2P1   RSTRET*RSTRET+1.
!     RC2M1   RSTRET*RSTRET-1.
!     RCOR0   COMPONENT (0,0) OF CORIOLIS
!     RCOR1   COMPONENT (0,1) OF CORIOLIS
!     RCOR2   COMPONENT (1,1) OF CORIOLIS

!     RCOLON(NGPTOT) cosine of longitude on transformed sphere
!     RSILON(NGPTOT)   sine        "             "         "
!     RINDX (NGPTOT) Longitude index
!     RINDY (NGPTOT) Latitude index
!     RATATH(NGPTOT) RA*TAN(THETA) on real sphere
!     RATATX(NGPTOT) Curvature term for LAM (for u eq.)

!     NLOEN(NDGSAG:NDGENG)  : number of active points on a parallel
!     NLOENG(NDGSAG:NDGENG) : global version of NLOEN
!     NMEN(NDGSAG:NDGENG)   : associated cut-off wave number
!     NMENTC(NDGSAG:NDGENG) : same as NMEN but for truncation NTCMAX.
!     NMENG(NDGSAG:NDGENG)   : global version of NMEN
!     NDGLU(0:MAX(NSMAX,NMSMAX)) : number of active points in an hemisphere
!                          for a given wave number m
!     NSTAGP(NGPTOT)     : start position of latitude data for boundary fields
!     NESTAGP(NGPTOT)    : start position of latitude data for boundary fields
!                          in extension zone (ALADIN).
!     NTSTAGP(NGPTOT)    : start position of latitude data for boundary fields
!                          in C+I+E zone (ALADIN).


!  REFLRHC :  reference length for critical relative humidity (cloud scheme)
!  REFLKUO, REFLCAPE, REFLRHC are ONLY to be used in the SETUP 
!  of TEQK and TEQC and TEQH

!  TEQK : ratio between REFLKUO and the model equivalent mesh size
!  TEQC : ratio between REFLCAPE and the model equivalent mesh size
!  TEQH : ratio between REFLRHC and the model equivalent mesh size

!  TEQK, TEQC, TEQH are to be used in the PHYSICS.

REAL(KIND=JPRB) :: rmucen
REAL(KIND=JPRB) :: rlocen
REAL(KIND=JPRB) :: rstret
INTEGER(KIND=JPIM) :: nsttyp
INTEGER(KIND=JPIM) :: nhtyp
REAL(KIND=JPRB) :: rnlginc
REAL(KIND=JPRB) :: r4jp
REAL(KIND=JPRB) :: rc2p1
REAL(KIND=JPRB) :: rc2m1
REAL(KIND=JPRB) :: rcor0
REAL(KIND=JPRB) :: rcor1
REAL(KIND=JPRB) :: rcor2
REAL(KIND=JPRB),ALLOCATABLE:: rcolon(:)
REAL(KIND=JPRB),ALLOCATABLE:: rsilon(:)
REAL(KIND=JPRB),ALLOCATABLE:: rindx(:)
REAL(KIND=JPRB),ALLOCATABLE:: rindy(:)
REAL(KIND=JPRB),ALLOCATABLE:: ratath(:)
REAL(KIND=JPRB),ALLOCATABLE:: ratatx(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nloen(:)
INTEGER(KIND=JPIM),ALLOCATABLE,TARGET :: nloeng(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nmen(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nmentc(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nmeng(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: ndglu(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nstagp(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nestagp(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: ntstagp(:)

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

!*    * Defining the transformed sphere: physics input

REAL(KIND=JPRB) :: reflrhc
REAL(KIND=JPRB) :: teqh
REAL(KIND=JPRB) :: reflkuo
REAL(KIND=JPRB) :: reflcape
REAL(KIND=JPRB) :: teqk
REAL(KIND=JPRB) :: teqc

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

!*    * DEFINING THE VERTICAL COORDINATE

!     VP00  : REFERENCE PRESSURE FOR DEFINING VERTICAL COORDINATE
!     VALH  : (0:NFLEVG)
!     VBH   : (0:NFLEVG) : B of the vertical coordinate
!     VETAH : (0:NFLEVG) ; VERTICAL COORDINATE = VALH+VBH
!     VETAF : (0:NFLEVG+1) ; VERTICAL COORDINATE ON LAYERS.
!     VCUICO: is used to compute denominators of weights
!             for semi-Lagrangian vertical interpolations
!             applied to full-level variables.
!     VCUICOH:is used to compute denominators of weights 
!             for semi-Lagrangian vertical interpolations
!             applied to half-level variables.
!     VRLEVX: REAL(NRLEVX)
!     NVAUTF: NVAUTF(VRLEVX*eta) is the number of the layer (full level)
!             immediately above "eta", and is bounded by 1 and nflevg-1.
!     NVAUTH: NVAUTH(VRLEVX*eta) is the number of the interlayer (half level)
!             immediately above "eta", and is bounded by 0 and nflevg-1.
!     VAH   : (0:NFLEVG) ;  =VALH*VP00
!     VC    : (NFLEVG)   ;  =VAH(J)*VBH(J-1)-VAH(J-1)*VBH(J)
!     VDELB : (NFLEVG)   ;  =VBH(J)-VBH(J-1)
!     VDELA : (NFLEVG)   ;  =VAH(J)-VAH(J-1)
!     VAF   : like VAH but at full levels.
!     VBF   : like VBH but at full levels.
!     VRDETAH: 1/[Delta eta]
!     TOPPRES: REFERENCE "EVANESCENT" PRESSURE
!              TOPPRES allows to solve some calculations of singularities
!              when the top pressure of the model is zero (for ex. in
!              GPPREF, GPXYB, SUNHBMAT).

!     WE HAVE THEN FOR THE HALF LEVEL PRESSURE : VAH + VBH*(SURFACE PRESSURE)

!     NOTE THAT THE HALF LEVEL VALUE AT K+.5 IS VXXX(K)
!     (THE FULL LEVEL VALUES ARE FROM 1 TO NFLEVG)

REAL(KIND=JPRB) :: vp00
REAL(KIND=JPRB),ALLOCATABLE:: valh(:)
REAL(KIND=JPRB),ALLOCATABLE:: vbh(:)
REAL(KIND=JPRB),ALLOCATABLE:: vetah(:)
REAL(KIND=JPRB),ALLOCATABLE:: vetaf(:)
REAL(KIND=JPRB),ALLOCATABLE:: vcuico(:,:)
REAL(KIND=JPRB),ALLOCATABLE:: vcuicoh(:,:)
REAL(KIND=JPRB) :: vrlevx
INTEGER(KIND=JPIM),ALLOCATABLE:: nvautf(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nvauth(:)
REAL(KIND=JPRB),ALLOCATABLE:: vah(:)
REAL(KIND=JPRB),ALLOCATABLE:: vc(:)
REAL(KIND=JPRB),ALLOCATABLE:: vdelb(:)
REAL(KIND=JPRB),ALLOCATABLE:: vdela(:)
REAL(KIND=JPRB),ALLOCATABLE:: vaf(:)
REAL(KIND=JPRB),ALLOCATABLE:: vbf(:)
REAL(KIND=JPRB),ALLOCATABLE:: vrdetah(:)
REAL(KIND=JPRB) :: toppres

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

!*    * Miscellaneous

!     NBEEGP : ???
!     NBNEGP : ???

INTEGER(KIND=JPIM) :: nbeegp
INTEGER(KIND=JPIM) :: nbnegp

!     ------------------------------------------------------------------
!$OMP THREADPRIVATE(nbeegp,nbnegp,ngptot,ngptot_cap,ngptotg,ngptotmx,nhtyp,nsttyp,r4jp,rc2m1,rc2p1,rcor0)
!$OMP THREADPRIVATE(rcor1,rcor2,reflcape,reflkuo,reflrhc,rlocen,rmucen,rnlginc,rstret,teqc,teqh,teqk,toppres,vp00,vrlevx)
!$OMP THREADPRIVATE(ndglu,nestagp,ngptotl,nloen,nloeng,nmen,nmeng,nmentc,nstagp,ntstagp,nvautf,nvauth,ratath,ratatx)
!$OMP THREADPRIVATE(rcolon,rindx,rindy,rsilon,vaf,vah,valh,vbf,vbh,vc,vcuico,vcuicoh,vdela,vdelb,vetaf,vetah,vrdetah)
END MODULE yomgem