modn_isban.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.
!##################
MODULE modn_isba_n
!##################
!
!!****  *MODN_ISBA_n* - declaration of namelist NAM_ISBA_n
!!
!!    PURPOSE
!!    -------
!       The purpose of this module is to specify  the namelist NAM_ISBA_n
!     which concern the surface configuration of one nested  model.
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!       
!!    AUTHOR
!!    ------
!!      V. Masson    *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004                    
!!      Modified    08/2009 by B. Decharme : LSURF_BUDGETC for all tiles
!!      Modified by A.L. Gibelin, 04/2009: add carbon spinup
!!
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
 CHARACTER(LEN=4)  :: CC1DRY
 CHARACTER(LEN=4)  :: CSCOND
 CHARACTER(LEN=3)  :: CSOILFRZ
 CHARACTER(LEN=4)  :: CDIFSFCOND
 CHARACTER(LEN=3)  :: CSNOWRES
 CHARACTER(LEN=4)  :: CALBEDO
 CHARACTER(LEN=4)  :: CROUGH
 CHARACTER(LEN=3)  :: CCPSURF
 CHARACTER(LEN=4)  :: CRUNOFF
 CHARACTER(LEN=3)  :: CKSAT
LOGICAL           :: LSOC
 CHARACTER(LEN=3)  :: CRAIN
 CHARACTER(LEN=3)  :: CHORT         
REAL  :: XTSTEP
REAL  :: XOUT_TSTEP
REAL  :: XCGMAX
REAL  :: XCDRAG
REAL  :: XSPINMAXS
REAL  :: XSPINMAXW
REAL  :: XCO2_START
REAL  :: XCO2_END
INTEGER  :: NNBYEARSPINS
INTEGER  :: NNBYEARSPINW
LOGICAL  :: LPGD
LOGICAL  :: LSURF_EVAP_BUDGET
LOGICAL  :: LSURF_BUDGET
LOGICAL  :: LRAD_BUDGET
LOGICAL  :: LSURF_BUDGETC
LOGICAL  :: LRESET_BUDGETC
LOGICAL  :: LPATCH_BUDGET
LOGICAL  :: LSURF_MISC_BUDGET
LOGICAL  :: LSURF_DIAG_ALBEDO
LOGICAL  :: LSURF_MISC_DIF
LOGICAL  :: LWATER_BUDGET
REAL     :: XDIAG_TSTEP
REAL     :: XALBICE1, XALBICE2,XALBICE3
REAL     :: XRHOTHRESHOLD_ICE
REAL     :: XZ0ICEZ0SNOW
REAL     :: XVAGING_NOGLACIER,XVAGING_GLACIER
REAL     :: XPERCENTAGEPORE
INTEGER  :: N2M
LOGICAL  :: L2M_MIN_ZS
LOGICAL  :: LCOEF
LOGICAL  :: LSURF_VARS
LOGICAL  :: LCH_BIO_FLUX
LOGICAL  :: LCH_NO_FLUX
LOGICAL  :: LGLACIER
LOGICAL  :: LVEGUPD
LOGICAL  :: LNITRO_DILU
LOGICAL  :: LSPINUPCARBS
LOGICAL  :: LSPINUPCARBW
LOGICAL  :: LCANOPY_DRAG
LOGICAL  :: LPERTSURF
LOGICAL  :: LSNOWDRIFT
LOGICAL  :: LSNOWDRIFT_SUBLIM
LOGICAL  :: LSNOW_ABS_ZENITH
 CHARACTER(3) :: CSNOWMETAMO
 CHARACTER(3) :: CSNOWRAD
 CHARACTER(LEN=6)  :: CCH_DRY_DEP
 CHARACTER(LEN=28) :: CCHEM_SURF_FILE
!
namelist/nam_isban/cc1dry,cscond,csoilfrz,cdifsfcond,csnowres,calbedo,crough,ccpsurf, &
                     xtstep,xcgmax,xcdrag,lglacier,lcanopy_drag,lvegupd,lpertsurf
namelist/nam_isba_agsn/lnitro_dilu
namelist/nam_sgh_isban/crunoff,cksat,lsoc,crain,chort
namelist/nam_diag_isban/lpgd,lsurf_evap_budget,lsurf_misc_budget,lsurf_diag_albedo,lpatch_budget, &
                        lsurf_misc_dif, lwater_budget
namelist/nam_diag_surfn/n2m,l2m_min_zs,lsurf_budget,lrad_budget, &
                        lsurf_budgetc,lreset_budgetc,lcoef,lsurf_vars
namelist/nam_ch_controln/cchem_surf_file
namelist/nam_ch_isban/lch_bio_flux,cch_dry_dep,lch_no_flux
namelist/nam_spinup_carbn/lspinupcarbs,lspinupcarbw,xspinmaxs,xspinmaxw,nnbyearspins,&
        nnbyearspinw, xco2_start, xco2_end
namelist/nam_isba_snown/lsnowdrift,lsnowdrift_sublim, lsnow_abs_zenith, csnowmetamo, csnowrad
!
 CONTAINS
!
SUBROUTINE init_nam_isban (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBAN',0,zhook_handle)
  cc1dry = i%CC1DRY
  cscond = i%CSCOND
  csoilfrz = i%CSOILFRZ
  cdifsfcond = i%CDIFSFCOND
  csnowres = i%CSNOWRES
  calbedo = i%CALBEDO
  crough = i%CROUGH
  ccpsurf = i%CCPSURF
  xtstep = i%XTSTEP
  xout_tstep = i%XOUT_TSTEP
  xcgmax = i%XCGMAX
  xcdrag = i%XCDRAG
  lglacier = i%LGLACIER
  lvegupd = i%LVEGUPD  
  lcanopy_drag = i%LCANOPY_DRAG
  lpertsurf = i%LPERTSURF
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBAN',1,zhook_handle)
END SUBROUTINE init_nam_isban

SUBROUTINE init_nam_isba_agsn (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBA_AGSN',0,zhook_handle)
  lnitro_dilu = i%LNITRO_DILU
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBA_AGSN',1,zhook_handle)
END SUBROUTINE init_nam_isba_agsn

SUBROUTINE init_nam_sgh_isban (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_SGH_ISBAN',0,zhook_handle)
  crunoff = i%CRUNOFF
  cksat=i%CKSAT
  lsoc=i%LSOC
  crain=i%CRAIN
  chort=i%CHORT
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_SGH_ISBAN',1,zhook_handle)
END SUBROUTINE init_nam_sgh_isban

SUBROUTINE update_nam_isban (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_ISBAN',0,zhook_handle)
  i%CC1DRY = cc1dry
  i%CSCOND = cscond
  i%CSOILFRZ = csoilfrz
  i%CDIFSFCOND = cdifsfcond
  i%CSNOWRES = csnowres
  i%CALBEDO = calbedo
  i%CROUGH = crough
  i%CCPSURF = ccpsurf
  i%XTSTEP = xtstep
  i%XOUT_TSTEP = xout_tstep
  i%XCGMAX = xcgmax
  i%XCDRAG = xcdrag
  i%LGLACIER = lglacier
  i%LVEGUPD = lvegupd  
  i%LCANOPY_DRAG = lcanopy_drag
  i%LPERTSURF = lpertsurf
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_ISBAN',1,zhook_handle)
END SUBROUTINE update_nam_isban

SUBROUTINE update_nam_isba_agsn (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_ISBA_AGSN',0,zhook_handle)
  i%LNITRO_DILU = lnitro_dilu
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_ISBA_AGSN',1,zhook_handle)
END SUBROUTINE update_nam_isba_agsn

SUBROUTINE update_nam_sgh_isban (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_SGH_ISBAN',0,zhook_handle)
  i%CRUNOFF=crunoff
  i%CKSAT=cksat
  i%LSOC=lsoc
  i%CRAIN=crain
  i%CHORT=chort
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_SGH_ISBAN',1,zhook_handle)
END SUBROUTINE update_nam_sgh_isban

SUBROUTINE init_nam_diag_isban (DGEI, DGI, DGMI)
!
  USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
  USE modd_diag_isba_n, ONLY : diag_isba_t
  USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_t
!
  IMPLICIT NONE

!
  TYPE(diag_evap_isba_t), INTENT(INOUT) :: dgei
  TYPE(diag_isba_t), INTENT(INOUT) :: dgi
  TYPE(diag_misc_isba_t), INTENT(INOUT) :: dgmi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_DIAG_ISBAN',0,zhook_handle)
  lpgd = dgi%LPGD
  lsurf_evap_budget = dgei%LSURF_EVAP_BUDGET
  lwater_budget = dgei%LWATER_BUDGET
  lsurf_misc_budget = dgmi%LSURF_MISC_BUDGET
  lsurf_diag_albedo = dgmi%LSURF_DIAG_ALBEDO
  lsurf_misc_dif = dgmi%LSURF_MISC_DIF
  lpatch_budget = dgi%LPATCH_BUDGET
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_DIAG_ISBAN',1,zhook_handle)
END SUBROUTINE init_nam_diag_isban

SUBROUTINE update_nam_diag_isban (DGEI, DGI, DGMI)
!
  USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
  USE modd_diag_isba_n, ONLY : diag_isba_t
  USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_t
!
  IMPLICIT NONE

!
  TYPE(diag_evap_isba_t), INTENT(INOUT) :: dgei
  TYPE(diag_isba_t), INTENT(INOUT) :: dgi
  TYPE(diag_misc_isba_t), INTENT(INOUT) :: dgmi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_DIAG_ISBAN',0,zhook_handle)
  dgi%LPGD = lpgd
  dgei%LSURF_EVAP_BUDGET = lsurf_evap_budget
  dgei%LWATER_BUDGET = lwater_budget
  dgmi%LSURF_MISC_BUDGET = lsurf_misc_budget
  dgmi%LSURF_DIAG_ALBEDO = lsurf_diag_albedo
  dgmi%LSURF_MISC_DIF = lsurf_misc_dif
  dgi%LPATCH_BUDGET = lpatch_budget
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_DIAG_ISBAN',1,zhook_handle)
END SUBROUTINE update_nam_diag_isban

SUBROUTINE init_nam_diag_surfn (DGEI, DGI)
!
  USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
  USE modd_diag_isba_n, ONLY : diag_isba_t
!
  IMPLICIT NONE

!
  TYPE(diag_evap_isba_t), INTENT(INOUT) :: dgei
  TYPE(diag_isba_t), INTENT(INOUT) :: dgi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_DIAG_SURFN',0,zhook_handle)
  xdiag_tstep = dgi%XDIAG_TSTEP  
  n2m = dgi%N2M
  l2m_min_zs = dgi%L2M_MIN_ZS
  lsurf_budget = dgi%LSURF_BUDGET
  lrad_budget = dgi%LRAD_BUDGET
  lsurf_budgetc = dgei%LSURF_BUDGETC
  lreset_budgetc = dgei%LRESET_BUDGETC  
  lcoef = dgi%LCOEF
  lsurf_vars = dgi%LSURF_VARS
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_DIAG_SURFN',1,zhook_handle)
END SUBROUTINE init_nam_diag_surfn

SUBROUTINE update_nam_diag_surfn (DGEI, DGI)
!
  USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
  USE modd_diag_isba_n, ONLY : diag_isba_t
!
  IMPLICIT NONE

!
  TYPE(diag_evap_isba_t), INTENT(INOUT) :: dgei
  TYPE(diag_isba_t), INTENT(INOUT) :: dgi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_DIAG_SURFN',0,zhook_handle)
  dgi%XDIAG_TSTEP = xdiag_tstep  
  dgi%N2M = n2m
  dgi%L2M_MIN_ZS = l2m_min_zs
  dgi%LSURF_BUDGET = lsurf_budget
  dgi%LRAD_BUDGET = lrad_budget 
  dgei%LSURF_BUDGETC = lsurf_budgetc
  dgei%LRESET_BUDGETC = lreset_budgetc  
  dgi%LCOEF = lcoef
  dgi%LSURF_VARS = lsurf_vars
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_DIAG_SURFN',1,zhook_handle)
END SUBROUTINE update_nam_diag_surfn

SUBROUTINE init_nam_ch_controln (CHI)
!
  USE modd_ch_isba_n, ONLY : ch_isba_t
!
  IMPLICIT NONE

!
  TYPE(ch_isba_t), INTENT(INOUT) :: chi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_CH_CONTROLN',0,zhook_handle)
  cchem_surf_file = chi%CCHEM_SURF_FILE
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_CH_CONTROLN',1,zhook_handle)
END SUBROUTINE init_nam_ch_controln

SUBROUTINE update_nam_ch_controln (CHI)
!
  USE modd_ch_isba_n, ONLY : ch_isba_t
!
  IMPLICIT NONE

!
  TYPE(ch_isba_t), INTENT(INOUT) :: chi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_CH_CONTROLN',0,zhook_handle)
  chi%CCHEM_SURF_FILE = cchem_surf_file
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_CH_CONTROLN',1,zhook_handle)
END SUBROUTINE update_nam_ch_controln

SUBROUTINE init_nam_ch_isban (CHI)
!
  USE modd_ch_isba_n, ONLY : ch_isba_t
!
  IMPLICIT NONE

!
  TYPE(ch_isba_t), INTENT(INOUT) :: chi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_CH_ISBAN',0,zhook_handle)
  lch_bio_flux = chi%LCH_BIO_FLUX
  lch_no_flux = chi%LCH_NO_FLUX
  cch_dry_dep = chi%CCH_DRY_DEP
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_CH_ISBAN',1,zhook_handle)
END SUBROUTINE init_nam_ch_isban

SUBROUTINE update_nam_ch_isban (CHI)
!
  USE modd_ch_isba_n, ONLY : ch_isba_t
!
  IMPLICIT NONE

!
  TYPE(ch_isba_t), INTENT(INOUT) :: chi
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_CH_ISBAN',0,zhook_handle)
  chi%LCH_BIO_FLUX = lch_bio_flux
  chi%LCH_NO_FLUX = lch_no_flux
  chi%CCH_DRY_DEP = cch_dry_dep
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_CH_ISBAN',1,zhook_handle)
END SUBROUTINE update_nam_ch_isban

SUBROUTINE init_nam_spinup_carb_isban (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_SPINUP_CARB_ISBAN',0,zhook_handle)
  lspinupcarbs = i%LSPINUPCARBS
  lspinupcarbw = i%LSPINUPCARBW  
  xspinmaxs   = i%XSPINMAXS
  xspinmaxw   = i%XSPINMAXW  
  nnbyearspins = i%NNBYEARSPINS
  nnbyearspinw = i%NNBYEARSPINW
  xco2_start = i%XCO2_START
  xco2_end = i%XCO2_END
IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_SPINUP_CARB_ISBAN',1,zhook_handle)
END SUBROUTINE init_nam_spinup_carb_isban

SUBROUTINE update_nam_spinup_carb_isban (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_SPINUP_CARB_ISBAN',0,zhook_handle)
  i%LSPINUPCARBS = lspinupcarbs
  i%LSPINUPCARBW = lspinupcarbw  
  i%XSPINMAXS   = xspinmaxs
  i%XSPINMAXW   = xspinmaxw  
  i%NNBYEARSPINS = nnbyearspins
  i%NNBYEARSPINW = nnbyearspinw
  i%XCO2_START = xco2_start
  i%XCO2_END = xco2_end
IF (lhook) CALL dr_hook('MODN_ISBA_N:UPDATE_NAM_SPINUP_CARB_ISBAN',1,zhook_handle)
END SUBROUTINE update_nam_spinup_carb_isban

SUBROUTINE init_nam_isba_snown (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBA_SNOWn',0,zhook_handle)
  lsnowdrift = i%LSNOWDRIFT
  lsnowdrift_sublim =  i%LSNOWDRIFT_SUBLIM
  lsnow_abs_zenith = i%LSNOW_ABS_ZENITH
  csnowmetamo = i%CSNOWMETAMO
  csnowrad = i%CSNOWRAD  
  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBA_SNOWn',1,zhook_handle)

END SUBROUTINE init_nam_isba_snown

SUBROUTINE update_nam_isba_snown (I)
!
  USE modd_isba_n, ONLY : isba_t
!
  IMPLICIT NONE

!
  TYPE(isba_t), INTENT(INOUT) :: i
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBA_SNOWn',0,zhook_handle)
  i%LSNOWDRIFT = lsnowdrift
  i%LSNOWDRIFT_SUBLIM = lsnowdrift_sublim
  i%LSNOW_ABS_ZENITH = lsnow_abs_zenith
  i%CSNOWMETAMO = csnowmetamo
  i%CSNOWRAD=csnowrad  
  IF (lhook) CALL dr_hook('MODN_ISBA_N:INIT_NAM_ISBA_SNOWn',1,zhook_handle)

END SUBROUTINE update_nam_isba_snown

END MODULE modn_isba_n