yomphy1.F90

Go to the documentation of this file.

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