MODULE YOMPHY1

USE PARKIND1  ,ONLY : JPIM     ,JPRB

IMPLICIT NONE

SAVE

!*
!     ------------------------------------------------------------------
!     CONSTANTES PHYSIQUES REGLABLES INTERVENANT DANS LA DEFINITION DU
!     SOL, DE LA NEIGE OU DE LA VEGETATION :
!       ALB1    : COEFF. POUR LE SCHEMA (LVGSN).
!               : COEFF. FOR THE SNOW SCHEME (LVGSN).
!       ALB2    : COEFF. POUR LE SCHEMA (LVGSN).
!               : COEFF. FOR THE SNOW SCHEME (LVGSN).
!       ALBGLA  : ALBEDO BANQUISE (SANS NEIGE).
!               : SNOW FREE SEA-ICE ALBEDO.
!       ALBMAX  : ALBEDO MAXIMUM DE LA NEIGE.
!               : MAXIMUM SNOW ALBEDO.
!       ALBMER  : ALBEDO MER.
!               : OPEN SEA ALBEDO.
!       ALBMED  : ALBEDO MER DIFFUS.
!               : OPEN SEA DIFFUSE ALBEDO.
!       ALBMIN  : ALBEDO MINIMUM DE LA NEIGE.
!               : MINIMUM SNOW ALBEDO.
!       ALCRIN  : ALBEDO DE LA NEIGE.
!               : SNOW ALBEDO.
!       ALRCN1  : PREMIERE CONSTANTE POUR L'INFLUENCE NEIGEUSE SUR Z0.
!               : FIRST CONSTANT FOR THE SNOW'S INFLUENCE ON Z0.
!       ALRCN2  : DEUXIEME CONSTANTE POUR L'INFLUENCE NEIGEUSE SUR Z0.
!               : SECOND CONSTANT FOR THE SNOW'S INFLUENCE ON Z0.
!       EA      : EXPOSANT DANS LE SCHEMA SOL-VEGETATION
!               : EXPONENT IN THE SOIL-VEGETATION SCHEME
!       EC2REF  : IDEM
!       EMCRIN  : EMISSIVITE DE LA NEIGE.
!               : SNOW EMISSIVITY.
!       EMMGLA  : EMISSIVITE DE LA BANQUISE
!               : SEA-ICE EMISSIVITY.
!       EMMMER  : EMISSIVITE DE LA MER
!               : OPEN SEA EMISSIVITY.
!       EWFC    : EXPOSANT DANS LE SCHEMA SOL-VEGETATION
!               : EXPONENT IN THE SOIL-VEGETATION SCHEME
!       EWWILT  : IDEM
!       GA      : COEFFICIENT DE REGRESSION DANS LE SCHEMA SOL-VEGETATION
!               : REGRESSION COEFFICIENT IN THE SOIL-VEGETATION SCHEME
!       GC1     : IDEM.
!       GC1S1, GC1S2, GC1S3, GC1S4 : COEFFICIENT POUR LE C1 PHASE VAPEUR
!                                  : COEFFICIENT FOR C1 VAPOUR PHASE
!       GTSVAP  : =0 PTS VARIABLE POUR LE CALCUL DU MAX DE C1. SI < 0 -GTSVAP
!                 REPRESENTE LA TEMPERATURE CONSTANTE EN DEGRE CELSIUS
!                 =0 PTS VARIABLE FOR THE COMPUTATION OF THE MAX OF C1. IF
!                 < 0 -GTSVAP IS CONSTANT TEMPERATURE IN DEGRE
!       GC1Y1   : VALEUR DU C1 POUR WS=0 POUR LA PHASE VAPEUR DE C1
!               : VALUE OF C1 FOR WS=0 FOR VAPOUR PHASE OF C1
!       GVEGMX,GLAIMX,GNEIMX,GWPIMX,GCGEL: COEFF POUR EAU GELEE DU SOL PROFOND.
!       GVEGMX,GLAIMX,GNEIMX,GWPIMX,GCGEL: COEFF FOR ICE IN DEEP SOIL.
!       GVEGMXS,GLAIMXS,GNEIMXS,GCGELS: COEFF POUR EAU GELEE DE SURFACE.
!       GVEGMXS,GLAIMXS,GNEIMXS,GCGELS: COEFF FOR ICE IN THE SUPERFICIAL
!               RESERVOIR.
!       GC2     : REGRESSION COEFFICIENT IN THE SOIL-VEGETATION SCHEME
!       GC2REF  : IDEM.
!       GC3     : IDEM.
!       GCZ0H   : IDEM.
!       GC31    : COEFFICIENT POUR LE DRAINAGE PROFOND (SOL-VEGETATION)
!               : COEFFICIENT FOR THE DEEP SOIL DRAINAGE (SOIL-VEGETATION)
!       GC32    : IDEM
!       GCONV   : FACTEUR DE CHANGEMENT D'UNITE (MASSE VOL. DE L'EAU LIQUIDE)
!       GF1     : COEFFICIENT DE LA FONCTION F1 (SOL-VEGETATION).
!               : COEFFICIENT FOR THE FUNCTION F1 (SOIL-VEGETATION).
!       GF3(18) : COEFFICIENT DE LA FONCTION F3 (SOL-VEGETATION)
!               : COEFFICIENT FOR THE FUNCTION F3 (SOIL-VEGETATION)
!       GF4(18) : COEFFICIENT DE LA FONCTION F4 (SOL-VEGETATION)
!               : COEFFICIENT FOR THE FUNCTION F4 (SOIL-VEGETATION)
!       GWFC    : COEFFICIENT DE REGRESSION DANS LE SCHEMA SOL-VEGETATION
!               : REGRESSION COEFFICIENT IN THE SOIL-VEGETATION SCHEME
!       GWWILT  : IDEM.
!       G1B     : IDEM.
!       G1CGSAT : IDEM.
!       G1C1SAT : IDEM.
!       G1P     : COEFFICIENT DE REGRESSION DANS LE SCHEMA SOL-VEGETATION
!               : REGRESSION COEFFICIENT IN THE SOIL-VEGETATION SCHEME
!       G1WSAT  : IDEM.
!       G2B     : IDEM.
!       G2CGSAT : IDEM.
!       G2C1SAT : IDEM.
!       G2P     : IDEM.
!       G2WSAT  : IDEM.
!       G3CGSAT : IDEM.
!       GSNC1   : IDEM (LSNV).
!       GSNC2   : IDEM (LSNV).
!       GWLEX   : (WL/WLMX)**GWLEX POUR LE CALCUL DU COEF. D'HALSTEAD.
!               : (WL/WLMX)**GWLEX FOR THE HALSTEAD COEF.
!       GWLMX   : COEFFICIENT POUR LE RESERVOIR D'INTERCEPTION MAX.
!               : COEFFICIENT FOR INTERCEPTION WATER CONTENT.
!       HSOL    : INVERSE DE L'INERTIE DU SOL EN SURFACE.
!               : INVERSE OF THE SURFACE SOIL INERTIA.
!       HSOLIWR : INVERSE DE LA QUANTITE CRITIQUE EN EAU (CALCUL DE HSOL).
!               : INVERSE OF CRITICAL SOIL WATER CONTENT (HSOL COMPUTATION).
!       HSOLIT0 : INVERSE DE LA LARGEUR CRITIQUE EN TEMPERATURE (CALCUL DE HSOL).
!               : INVERSE OF TEMPERATURE CRITICAL DEVIATION (HSOL COMPUTATION).
!       NTVGLA  : INDICE DE VEGETATION SUR GLACE.
!               : VEGETATION INDEX OVER ICE.
!       NTVMER  : INDICE DE VEGETATION SUR MER.
!               : VEGETATION INDEX OVER SEA.
!       OMTPRO  : 2 PI SUR LA PER. DE RAPP. DE LA TEMPERATURE PROFONDE.
!               : 2 PI DIVIDED BY THE TIME CONSTANT OF DEEP TEMPERATURE.
!       OMWPRO  : 1.   SUR LA PER. DE RAPP. DU CONTENU EN EAU PROFOND.
!               : 1.   DIVIDED BY THE TIME CONSTANT OF DEEP SOIL WATER.
!       RC1MAX  : VALEUR MAXIMALE DU COEFFICIENT HYDRIQUE G1=C1/RD1
!               : MAXIMUM VALUE FOR THE HYDRIC COEFFICIENT G1
!       RCTGLA  : COEFFICIENT THERMIQUE DE LA GLACE
!               : ICE THERMAL COEFFICIENT.
!       RCGMAX  : COEFFICIENT THERMIQUE DU SOL NU MAXIMUM
!               : MAXIMUM THERMAL BARE GROUND COEFFICIENT.
!       RCTVEG  : COEFFICIENT THERMIQUE DE LA VEGETATION.
!          (18) : VEGETATION THERMAL COEFFICIENT.
!       RD1     : EPAISSEUR DU RESERVOIR SUPERFICIEL.
!               : UPPER RESERVOIR DEPTH.
!       RD2GLA  : EPAISSEUR DU RESERVOIR BANQUISE.
!               : SEA-ICE RESERVOIR DEPTH.
!       RD2MER  : EPAISSEUR DU RESERVOIR MER.
!               : OPEN SEA RESERVOIR DEPTH.
!       RGL(18) : VALEUR LIMITE DU RAYONNEMENT GLOBAL DEPENDANT DU TYPE
!                 DE VEGETATION. RGL INTERVIENT DANS LE CALCUL DE RSTO.
!               : VEGETATION RADIATIVE COEFFICIENT LIMITING GLOBAL
!                 RADIATION. RGL IS USED IN THE COMPUTATION OF RSTO.
!       RHOMAX  : DENSITE MAXIMUM DE LA NEIGE.
!               : MAXIMUM SNOW DENSITY.
!       RHOMIN  : DENSITE MINIMUM DE LA NEIGE.
!               : MINIMUM SNOW DENSITY.
!       RLAIMX  : LAI MAX POUR LE SCHEMA (LVGSN).
!               : LAI MAX FOR THE SNOW SCHEME (LVGSN).
!       RLAI    : LAI A PARTIR DUQUEL ON REDUIT LA FRACTION DE NEIGE (LVGSN).
!               : LAI FROM WHICH THE SNOW FRACTION IS REDUCED (LVGSN).
!       RSMAX   : RESISTANCE STOMATIQUE MAXIMALE.
!               : MAXIMUM STOMATAL RESISTANCE
!       RTINER  : RAPPORT DES INERTIES THERMIQUES PROFONDEUR/SURFACE.
!               : RATIO OF THE DEEP/SURFACE THERMAL INERTIAS.
!       RZ0GLA  : LONGUEUR RUGOSITE DYNAMIQUE DE LA BANQUISE
!               : DYNAMICAL ROUGHNESS LENGTH OF SEA-ICE
!       RZ0MER  : LONGUEUR RUGOSITE DYNAMIQUE DE LA MER
!               : DYNAMICAL ROUGHNESS LENGTH OF SEA
!       RZHZ0G  : RAPPORT LONGUEUR DE RUGOSITE THERMIQUE SUR LONGUEUR
!                 RUGOSITE DYNAMIQUE DE LA BANQUISE
!               : RATIO OF THERMAL ROUGHNESS LENGTH ON DYNAMICAL 
!                 ROUGHNESS LENGTH OF SEA-ICE
!       RZHZ0M  : RAPPORT LONGUEUR DE RUGOSITE THERMIQUE SUR LONGUEUR
!                 RUGOSITE DYNAMIQUE DE LA MER
!               : RATIO OF THERMAL ROUGHNESS LENGTH ON DYNAMICAL 
!                 ROUGHNESS LENGTH OF SEA
!       SODELX(0:9): DISCRETISATION VERTICALE DU SOL (MAXI 10 COUCHES)
!                : SOIL VERTICAL DISCRETIZATION (MAX 10 LAYERS)
!       TMERGL  : TEMPERATURE DE FONTE DE LA GLACE DE MER.
!               : MELTING TEMPERATURE OF FLOATING ICE.
!       TOEXP   : TAUX DE DECROISSANCE EXPONENTIELLE DE L'ALBEDO ET
!                 ET DE LA DENSITE DE LA NEIGE EN PERIODE DE FONTE.
!               : RATE OF EXPONENTIAL DECREASE IN SNOW ALBEDO AND DENSITY
!                 DURING SNOW MELT.
!       TOLIN   : TAUX DE DECROISSANCE LINEAIRE DE L'ALBEDO DE LA NEIGE
!                 EN PERIODE D'ACCUMULATION.
!               : RATE OF LINEAR DECREASE IN SNOW ALBEDO
!                 DURING SNOW ACCUMULATION.
!       TREF4(18): TEMPERATURE SEUIL DE LA FONCTION F4 (SOL-VEGETATION)
!                : THRESHOLD TEMPERATURE IN F4 FUNCTION (SOIL-VEGETATION)
!       WCRIN   : QUANTITE CRITIQUE POUR LA NEIGE MI-COUVRANTE.
!               : CRITICAL TRANSITION VALUE FOR SNOW DEPTH (HALF-COVER).
!       WCRINC  : QUANTITE CRITIQUE POUR LA NEIGE MI-COUVRANTE.
!                 (CALCULS THERMIQUES)
!               : CRITICAL TRANSITION VALUE FOR SNOW DEPTH (HALF-COVER).
!                 (THERMIC CALCULATIONS)
!       WCRING  : QUANTITE CRITIQUE POUR LA NEIGE MI-COUVRANTE.
!                 (CALCULS RADIATIFS ET EVAPORATION)
!               : CRITICAL TRANSITION VALUE FOR SNOW DEPTH (HALF-COVER).
!                 (RADIATIVE CALCULATIONS AND EVAPORATION)
!       WNEW    : PRECIPITIONS DE NEIGE MINIMUM (MM) POUR RENOUVELER
!                 L'ALBEDO DE LA NEIGE.
!               : MINIMUM SNOWFALL (MM) TO RENEW THE SNOW ALBEDO.
!       WPMX    : VALEUR MAXIMUN POUR LE CONTENU EN EAU EN PROFONDEUR.
!               : MAXIMUM VALUE FOR DEEP SOIL WATER CONTENT.
!       WSMX    : VALEUR MAXIMUN POUR LE CONTENU EN EAU DE SURFACE.
!               : MAXIMUM VALUE FOR SURFACE SOIL WATER CONTENT.
!       XCRINR  : COEFFICIENT DANS LE FACTEUR ZNR DETERMINANT
!                 L'INFLUENCE DU RELIEF SOUS-MAILLE
!                 SUR LES FRACTIONS DE NEIGE
!               : COEFFICIENT IN THE FACTOR ZNR DETERMINING
!                 THE INFLUENCE OF SUBGRID OROGRAPHY
!                 ON THE SNOW COVER FRACTIONS
!       XCRINV  : COEFFICIENT DE LA FRACTION PNEIJV DE NEIGE RECOUVRANT
!                 LA VEGETATION.
!               : COEFFICIENT FOR THE FRACTION PNEIJV OF SNOW COVERING
!                 THE CANOPY.
!       LIMC    : BORNE SUPERIEURE POUR CT
!               : UPPER BOUND FOR CT
!       LIMW    : CONTROLE DE W/Wwilt POUR CG,C1,C2
!               : CHECKING W/Wwilt
!       LC1VAP  : PHASE VAPEUR POUR LE C1
!               : VAPOUR PHASE FOR C1

REAL(KIND=JPRB) :: GF3(18)
REAL(KIND=JPRB) :: GF4(18)
REAL(KIND=JPRB) :: TREF4(18)
REAL(KIND=JPRB) :: RCTVEG(18)
REAL(KIND=JPRB) :: RGL(18)
REAL(KIND=JPRB) :: SODELX(0:9)
REAL(KIND=JPRB) :: GCZ0H(0:3,4)
REAL(KIND=JPRB) :: ALBGLA
REAL(KIND=JPRB) :: ALBMAX
REAL(KIND=JPRB) :: ALBMER
REAL(KIND=JPRB) :: ALBMED
REAL(KIND=JPRB) :: ALBMIN
REAL(KIND=JPRB) :: ALCRIN
REAL(KIND=JPRB) :: ALRCN1
REAL(KIND=JPRB) :: ALRCN2
REAL(KIND=JPRB) :: EA
REAL(KIND=JPRB) :: EC2REF
REAL(KIND=JPRB) :: EMCRIN
REAL(KIND=JPRB) :: EMMGLA
REAL(KIND=JPRB) :: EMMMER
REAL(KIND=JPRB) :: EWFC
REAL(KIND=JPRB) :: EWWILT
REAL(KIND=JPRB) :: GA
REAL(KIND=JPRB) :: GC1
REAL(KIND=JPRB) :: GC1S1
REAL(KIND=JPRB) :: GC1S2
REAL(KIND=JPRB) :: GC1S3
REAL(KIND=JPRB) :: GC1S4
REAL(KIND=JPRB) :: GC1Y1
REAL(KIND=JPRB) :: GTSVAP
REAL(KIND=JPRB) :: GVEGMX
REAL(KIND=JPRB) :: GLAIMX
REAL(KIND=JPRB) :: GNEIMX
REAL(KIND=JPRB) :: GWPIMX
REAL(KIND=JPRB) :: GCGEL
REAL(KIND=JPRB) :: GC2
REAL(KIND=JPRB) :: GC2REF
REAL(KIND=JPRB) :: GC3
REAL(KIND=JPRB) :: GC31
REAL(KIND=JPRB) :: GC32
REAL(KIND=JPRB) :: GCONV
REAL(KIND=JPRB) :: GF1
REAL(KIND=JPRB) :: GWFC
REAL(KIND=JPRB) :: GWLEX
REAL(KIND=JPRB) :: GWLMX
REAL(KIND=JPRB) :: GWWILT
REAL(KIND=JPRB) :: G1B
REAL(KIND=JPRB) :: G1CGSAT
REAL(KIND=JPRB) :: G1C1SAT
REAL(KIND=JPRB) :: G1P
REAL(KIND=JPRB) :: G1WSAT
REAL(KIND=JPRB) :: G2B
REAL(KIND=JPRB) :: G2CGSAT
REAL(KIND=JPRB) :: G2C1SAT
REAL(KIND=JPRB) :: G2P
REAL(KIND=JPRB) :: G2WSAT
REAL(KIND=JPRB) :: G3CGSAT
REAL(KIND=JPRB) :: GSNC1
REAL(KIND=JPRB) :: GSNC2
REAL(KIND=JPRB) :: HSOL
REAL(KIND=JPRB) :: HSOLIWR
REAL(KIND=JPRB) :: HSOLIT0
REAL(KIND=JPRB) :: OMTPRO
REAL(KIND=JPRB) :: OMWPRO
REAL(KIND=JPRB) :: RC1MAX
REAL(KIND=JPRB) :: RCTGLA
REAL(KIND=JPRB) :: RCGMAX
REAL(KIND=JPRB) :: RD1
REAL(KIND=JPRB) :: RD2GLA
REAL(KIND=JPRB) :: RD2MER
REAL(KIND=JPRB) :: RHOMAX
REAL(KIND=JPRB) :: RHOMIN
REAL(KIND=JPRB) :: RSMAX
REAL(KIND=JPRB) :: RTINER
REAL(KIND=JPRB) :: RZ0GLA
REAL(KIND=JPRB) :: RZ0MER
REAL(KIND=JPRB) :: RZHZ0G
REAL(KIND=JPRB) :: RZHZ0M
REAL(KIND=JPRB) :: RZHGLA
REAL(KIND=JPRB) :: RZHMER
REAL(KIND=JPRB) :: TMERGL
REAL(KIND=JPRB) :: TOEXP
REAL(KIND=JPRB) :: TOLIN
REAL(KIND=JPRB) :: WCRIN
REAL(KIND=JPRB) :: WCRINC
REAL(KIND=JPRB) :: WCRING
REAL(KIND=JPRB) :: WNEW
REAL(KIND=JPRB) :: WPMX
REAL(KIND=JPRB) :: WSMX
REAL(KIND=JPRB) :: XCRINR
REAL(KIND=JPRB) :: XCRINV
LOGICAL :: LIMC
LOGICAL :: LIMW
LOGICAL :: LC1VAP
INTEGER(KIND=JPIM) :: NTVGLA
INTEGER(KIND=JPIM) :: NTVMER
REAL(KIND=JPRB) :: GCGELS
REAL(KIND=JPRB) :: GVEGMXS
REAL(KIND=JPRB) :: GLAIMXS
REAL(KIND=JPRB) :: GNEIMXS
REAL(KIND=JPRB) :: ALB1
REAL(KIND=JPRB) :: ALB2
REAL(KIND=JPRB) :: RLAIMX
REAL(KIND=JPRB) :: RLAI
INTEGER(KIND=JPIM) :: NCHSP
!     ------------------------------------------------------------------
!$OMP THREADPRIVATE(alb1,alb2,albgla,albmax,albmed,albmer,albmin,alcrin,alrcn1,alrcn2,ea,ec2ref,emcrin,emmgla)
!$OMP THREADPRIVATE(emmmer,ewfc,ewwilt,g1b,g1c1sat,g1cgsat,g1p,g1wsat,g2b,g2c1sat,g2cgsat,g2p,g2wsat,g3cgsat)
!$OMP THREADPRIVATE(ga,gc1,gc1s1,gc1s2,gc1s3,gc1s4,gc1y1,gc2,gc2ref,gc3,gc31,gc32,gcgel,gcgels,gconv,gcz0h,gf1)
!$OMP THREADPRIVATE(gf3,gf4,glaimx,glaimxs,gneimx,gneimxs,gsnc1,gsnc2,gtsvap,gvegmx,gvegmxs,gwfc,gwlex,gwlmx)
!$OMP THREADPRIVATE(gwpimx,gwwilt,hsol,hsolit0,hsoliwr,lc1vap,limc,limw,nchsp,ntvgla,ntvmer,omtpro,omwpro,rc1max)
!$OMP THREADPRIVATE(rcgmax,rctgla,rctveg,rd1,rd2gla,rd2mer,rgl,rhomax,rhomin,rlai,rlaimx,rsmax,rtiner,rz0gla)
!$OMP THREADPRIVATE(rz0mer,rzhgla,rzhmer,rzhz0g,rzhz0m,sodelx,tmergl,toexp,tolin,tref4,wcrin,wcrinc,wcring)
!$OMP THREADPRIVATE(wnew,wpmx,wsmx,xcrinr,xcrinv)
END MODULE YOMPHY1