oi_acsolw.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
SUBROUTINE oi_acsolw (KST,KNBPT,&
!-----------------------------------------------------------------------
! - INPUT  1D
   parg,pd2,pws,piveg,psab,&
! - INPUT  LOGIQUE
   odhmt,&
! - OUTPUT 1D .
   pwfc,pwpmx,pwsat,pwsmx,pwwilt)  
  

!**** *ACSOLW  * - DETERMINATION DES CONTENUS EN EAU CARACTERISTIQUES DU SOL

!     Sujet.
!     ------

!     - ROUTINE DE CALCUL INTERMEDIAIRE.
!       DETERMINATION DES CONTENUS EN EAU CARACTERISTIQUES DU SOL.

!    SUR MER      (PITM=0. , PIVEG=NTVMER    )
!       wwilt=0.        wfc=1.        wsat=1.
!    SUR BANQUISE (PITM=1. , PIVEG=NTVGLA*1.1) (MODE CLIMAT SEULEMENT)
!       wwilt=0.        wfc=1.        wsat=1.
!    SUR GLACIER  (PITM=1. , PIVEG=NTVGLA    )
!       wwilt=f1(arg)   wfc=f2(arg)   wsat=wfc
!    SUR TERRE    (PITM=1. , AUTRES PIVEG    )
!       wwilt=f1(arg)   wfc=f2(arg)   wsat=f3(sab)

!**   Interface.
!     ----------
!        *CALL* *ACSOLW*

!-----------------------------------------------------------------------
! WARNING: THE ENGLISH VERSION OF VARIABLES' NAMES IS TO BE READ IN THE
!          "APLPAR" CODE.
!-----------------------------------------------------------------------

! -   ARGUMENTS D'ENTREE.
!     -------------------

! - NOM DES PARAMETRES DE DIMENSIONNEMENT DE LA PHYSIQUE.

! KST      : INDICE DE DEPART DES BOUCLES VECTORISEES SUR L'HORIZONT..
! KNBPT      : INDICE DE FIN DES BOUCLES VECTORISEES SUR L'HORIZONTALE.


! - NOM DES CLES LOGIQUES

! LDHMT      : CALCULS REDUITS (POUR L'ANALYSE OU FULL-POS)
!              SEULS LES TABLEAUX (*) SONT UTILISES/CALCULES EN E/S

! - NOM DES VARIABLES DE LA PHYSIQUE (PAR ORDRE ALPHABETIQUE DANS CHAQUE
!   CATEGORIE).

! - 1D (GEOGRAPHIQUE) .

! PARG (*)   : POURCENTAGE D'ARGILE DANS LA MAILLE.
! PD2        : EPAISSEUR DU RESERVOIR PROFOND.
! PWS (*)    : CONTENU EN EAU SUPERFICIEL (pour identification NATURE)
! PIVEG      : TYPE DE SURFACE (MER, BANQUISE OU GLACIER, TERRE).
! PSAB       : POURCENTAGE DE SABLE DANS LA MAILLE.

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

! -   ARGUMENTS DE SORTIE.
!     --------------------

! - NOM DES VARIABLES DE LA PHYSIQUE (PAR ORDRE ALPHABETIQUE DANS CHAQUE
!   CATEGORIE).

! - 1D (DIAGNOSTIQUE) .

! PWFC   (*) : TENEUR EN EAU A LA CAPACITE AUX CHAMPS.
! PWPMX      : TENEUR EN EAU MAXIMALE DU RESERVOIR PROFOND.
! PWSAT      : TENEUR EN EAU A LA SATURATION
! PWSMX      : TENEUR EN EAU MAXIMALE DU RESERVOIR DE SURFACE.
! PWWILT (*) : TENEUR EN EAU CORRESPONDANT AU POINT DE FLETRISSEMENT.

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

! -   ARGUMENTS IMPLICITES.
!     ---------------------

! COMMON/YOMPHY1/
! the common that contains the parameters for the land surface part of the 
! physics of the model
!-----------------------------------------------------------------------

!     Externes.
!     ---------

!     Methode.
!     --------

!     Auteur.
!     -------
!        99-02, D. Giard, d'apres ACSOL

!     Modifications.
!     --------------
!-----------------------------------------------------------------------
!
!  
USE modd_assim, ONLY : ntvgla, xg1wsat, xg2wsat, xewfc, xgwfc, xewwilt, xgwwilt, xrd1, xreps2
USE modd_csts,  ONLY : xrholw
USE modd_surf_par, ONLY : xundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!  
REAL    ,INTENT(IN)    :: parg(knbpt)
REAL    ,INTENT(IN)    :: pd2(knbpt)
REAL    ,INTENT(IN)    :: pws(knbpt)
REAL    ,INTENT(IN)    :: piveg(knbpt) 
REAL    ,INTENT(IN)    :: psab(knbpt)
!
LOGICAL, INTENT(IN)    :: odhmt
!
REAL    ,INTENT(OUT)   :: pwfc(knbpt) 
REAL    ,INTENT(OUT)   :: pwpmx(knbpt) 
REAL    ,INTENT(OUT)   :: pwsat(knbpt) 
REAL    ,INTENT(OUT)   :: pwsmx(knbpt) 
REAL    ,INTENT(OUT)   :: pwwilt(knbpt)
!
REAL    :: zarg, zsab
INTEGER    :: knbpt, kst
INTEGER    :: jrof
!
REAL(KIND=JPRB) :: zhook_handle
!
!  ** initialization of the parameters of the land surface part of 
!  ** the physics of the model
!
!     ------------------------------------------------------------------
!     I - CALCULS REDUITS DANS LE CAS LDHMT=.TRUE.
!         ----------------------------------------
!         REDUCED COMPUTATIONS IN CASE LDHMT=.T.
!         --------------------------------------
!
! ** WFC , WWILT **
!
IF (lhook) CALL dr_hook('OI_ACSOLW',0,zhook_handle)
!
DO jrof=kst,knbpt
  IF ( pws(jrof) == xundef ) THEN
    pwwilt(jrof) = 0.0
    pwfc(jrof) = 1.0
  ELSE
    zarg = max(xreps2,parg(jrof))
    pwwilt(jrof) = xgwwilt*(zarg**xewwilt)
    pwfc(jrof) = xgwfc  *(zarg**xewfc)
  ENDIF
ENDDO
!
IF (odhmt .AND. lhook) CALL dr_hook('OI_ACSOLW',1,zhook_handle)
IF (odhmt) RETURN
!
!     ------------------------------------------------------------------
!     II - CALCUL DES AUTRES CHAMPS
!          ------------------------
!          COMPUTING OTHER FIELDS
!          ----------------------
!
DO jrof=kst,knbpt
!
! ** CAS DE LA BANQUISE EN MODE CLIMAT - SEA-ICE - (CLIMAT) **
!
  IF ( nint(10.*piveg(jrof))==(10*ntvgla+1) ) THEN
    pwwilt(jrof) = 0.0
    pwfc(jrof) = 1.0
  ENDIF
!
! ** WSAT **
!
  IF ( pws(jrof)==xundef .OR. nint(piveg(jrof))==ntvgla ) THEN
    pwsat(jrof) = pwfc(jrof)
  ELSE
    zsab = max(xreps2,psab(jrof))
    pwsat(jrof) = xg1wsat*zsab + xg2wsat
  ENDIF
!
! **  WSMX , WPMX **
!
  pwpmx(jrof) = pwsat(jrof)*pd2(jrof)*xrholw
  pwsmx(jrof) = pwsat(jrof)*xrd1*xrholw
!
ENDDO
!
IF (lhook) CALL dr_hook('OI_ACSOLW',1,zhook_handle)
!
END SUBROUTINE oi_acsolw