read_nam_pgd_isba.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 read_nam_pgd_isba(HPROGRAM, KPATCH, KGROUND_LAYER,                         &
                                   HISBA, HPEDOTF, HPHOTO, OTR_ML, HALBEDO, PRM_PATCH,      &
                                   HCLAY, HCLAYFILETYPE, PUNIF_CLAY, OIMP_CLAY,             &
                                   HSAND, HSANDFILETYPE, PUNIF_SAND, OIMP_SAND,             &
                                   HSOC_TOP, HSOC_SUB, HSOCFILETYPE, PUNIF_SOC_TOP,         &
                                   PUNIF_SOC_SUB, OIMP_SOC, HCTI, HCTIFILETYPE, OIMP_CTI,   &
                                   HPERM, HPERMFILETYPE, PUNIF_PERM, OIMP_PERM, OMEB,       &          
                                   HRUNOFFB, HRUNOFFBFILETYPE, PUNIF_RUNOFFB,               &
                                   HWDRAIN,  HWDRAINFILETYPE , PUNIF_WDRAIN, PSOILGRID,     &
                                   HPH, HPHFILETYPE, PUNIF_PH, HFERT, HFERTFILETYPE,        &
                                   PUNIF_FERT      )  
!     ##############################################################
!
!!**** *READ_NAM_PGD_ISBA* reads namelist for ISBA
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!OTR_ML,                        !
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    01/2005
!!       2008 B. Decharme : uniform value of subgrid drainage coefficient
!!    12/2008 E. Martin   : files of data for subgrid drainage 
!!                          and subgridrunoff
!!    06/2009 B. Decharme : files of data for topographic index
!!    07/2012 B. Decharme : files of data for permafrost area and for SOC top and sub soil
!!    10/2014 P. Samuelsson: MEB
!!    10/2016 B. Decharme : bug surface/groundwater coupling 
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_surf_par, ONLY : xundef, nundef
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
!
USE mode_pos_surf
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
 CHARACTER(LEN=6),    INTENT(IN)    :: HPROGRAM      ! Type of program
INTEGER,             INTENT(OUT)   :: KPATCH        ! number of patches
INTEGER,             INTENT(OUT)   :: KGROUND_LAYER ! number of soil layers
 CHARACTER(LEN=3),    INTENT(OUT)   :: HISBA         ! ISBA option
 CHARACTER(LEN=4),    INTENT(OUT)   :: HPEDOTF       ! Pedo-transfert function for DIF
 CHARACTER(LEN=3),    INTENT(OUT)   :: HPHOTO        ! photosynthesis option
LOGICAL,             INTENT(OUT)   :: OTR_ML        ! new radiative transfert
 CHARACTER(LEN=4),   INTENT(OUT)    :: HALBEDO
REAL,                INTENT(OUT)   :: PRM_PATCH     ! threshold to remove little fractions of patches
 CHARACTER(LEN=28),   INTENT(OUT)   :: HSAND         ! file name for sand fraction
 CHARACTER(LEN=28),   INTENT(OUT)   :: HCLAY         ! file name for clay fraction
 CHARACTER(LEN=28),   INTENT(OUT)   :: HCTI          ! file name for topographic index
 CHARACTER(LEN=28),   INTENT(OUT)   :: HPERM         ! file name for permafrost distribution
 CHARACTER(LEN=28),   INTENT(OUT)   :: HRUNOFFB      ! file name for runoffb parameter
 CHARACTER(LEN=28),   INTENT(OUT)   :: HWDRAIN       ! file name for wdrain parameter
 CHARACTER(LEN=6),    INTENT(OUT)   :: HSANDFILETYPE ! sand data file type
 CHARACTER(LEN=6),    INTENT(OUT)   :: HCLAYFILETYPE ! clay data file type
 CHARACTER(LEN=6),    INTENT(OUT)   :: HCTIFILETYPE  ! topographic index data file type
 CHARACTER(LEN=6),    INTENT(OUT)   :: HPERMFILETYPE    ! permafrost distribution data file type
 CHARACTER(LEN=6),    INTENT(OUT)   :: HRUNOFFBFILETYPE ! subgrid runoff data file type
 CHARACTER(LEN=6),    INTENT(OUT)   :: HWDRAINFILETYPE  ! subgrid drainage data file type
REAL,                INTENT(OUT)   :: PUNIF_SAND    ! uniform value of sand fraction
REAL,                INTENT(OUT)   :: PUNIF_CLAY    ! uniform value of clay fraction
REAL,                INTENT(OUT)   :: PUNIF_RUNOFFB ! uniform value of subgrid runoff coefficient
REAL,                INTENT(OUT)   :: PUNIF_WDRAIN  ! uniform value of subgrid drainage coefficient
REAL,                INTENT(OUT)   :: PUNIF_PERM    ! uniform value of permafrost distribution
LOGICAL,             INTENT(OUT)   :: OIMP_SAND     ! Imposed values for Sand
LOGICAL,             INTENT(OUT)   :: OIMP_CLAY     ! Imposed values for Clay
LOGICAL,             INTENT(OUT)   :: OIMP_CTI      ! Imposed values for topographic index statistics
LOGICAL,             INTENT(OUT)   :: OMEB          ! MEB
LOGICAL,             INTENT(OUT)   :: OIMP_PERM     ! Imposed maps of permafrost distribution
 CHARACTER(LEN=28),   INTENT(OUT)   :: HSOC_TOP      ! file name for organic carbon
 CHARACTER(LEN=28),   INTENT(OUT)   :: HSOC_SUB      ! file name for organic carbon
 CHARACTER(LEN=6),    INTENT(OUT)   :: HSOCFILETYPE  ! organic carbon data file type
REAL,                INTENT(OUT)   :: PUNIF_SOC_TOP ! uniform value of organic carbon top soil (kg/m2)
REAL,                INTENT(OUT)   :: PUNIF_SOC_SUB ! uniform value of organic carbon sub soil (kg/m2)
LOGICAL,             INTENT(OUT)   :: OIMP_SOC      ! Imposed maps of organic carbon
REAL, DIMENSION(:),  INTENT(OUT)   :: PSOILGRID     ! Soil layer thickness for DIF
 CHARACTER(LEN=28),   INTENT(OUT)   :: HPH           ! file name for pH
 CHARACTER(LEN=28),   INTENT(OUT)   :: HFERT         ! file name for fertilisation rate
 CHARACTER(LEN=6),    INTENT(OUT)   :: HPHFILETYPE   ! pH data file type
 CHARACTER(LEN=6),    INTENT(OUT)   :: HFERTFILETYPE ! fertilisation data file type
REAL,                INTENT(OUT)   :: PUNIF_PH      ! uniform value of pH
REAL,                INTENT(OUT)   :: PUNIF_FERT    ! uniform value of fertilisation rate
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER                           :: ILUOUT    ! output listing logical unit
INTEGER                           :: ILUNAM    ! namelist file logical unit
LOGICAL                           :: GFOUND    ! flag when namelist is present
!
!*    0.3    Declaration of namelists
!            ------------------------
!
INTEGER                  :: NPATCH           ! number of patches
INTEGER                  :: NGROUND_LAYER    ! number of soil layers
 CHARACTER(LEN=3)         :: CISBA            ! ISBA option
 CHARACTER(LEN=4)         :: CPEDO_FUNCTION   ! Pedo-transfert function for DIF
 CHARACTER(LEN=3)         :: CPHOTO           ! photosynthesis option
LOGICAL                  :: LTR_ML           ! new radiative transfert
 CHARACTER(LEN=4)         :: CALBEDO
REAL                     :: XRM_PATCH        ! threshold to remove little fractions of patches
 CHARACTER(LEN=28)        :: YSAND            ! file name for sand fraction
 CHARACTER(LEN=28)        :: YCLAY            ! file name for clay fraction
 CHARACTER(LEN=28)        :: YCTI             ! file name for topographic index
 CHARACTER(LEN=28)        :: YPERM            ! file name for permafrost distribution
 CHARACTER(LEN=28)        :: YRUNOFFB         ! file name for runoffb parameter
 CHARACTER(LEN=28)        :: YWDRAIN          ! file name for wdrain parameter
 CHARACTER(LEN=28)        :: YPH              ! file name for pH
 CHARACTER(LEN=28)        :: YFERT            ! file name for fertilisation rate
 CHARACTER(LEN=6)         :: YSANDFILETYPE    ! sand data file type
 CHARACTER(LEN=6)         :: YCLAYFILETYPE    ! clay data file type
 CHARACTER(LEN=6)         :: YCTIFILETYPE     ! topographic index data file type
 CHARACTER(LEN=6)         :: YPERMFILETYPE    ! permafrost distribution data file type
 CHARACTER(LEN=6)         :: YRUNOFFBFILETYPE ! subgrid runoff data file type
 CHARACTER(LEN=6)         :: YWDRAINFILETYPE  ! subgrid drainage data file type
 CHARACTER(LEN=6)         :: YPHFILETYPE      ! pH data file type
 CHARACTER(LEN=6)         :: YFERTFILETYPE    ! fertilisation data file type
LOGICAL                  :: LIMP_SAND        ! Imposed maps of Sand from another PGD file
LOGICAL                  :: LIMP_CLAY        ! Imposed maps of Clay from another PGD file
LOGICAL                  :: LIMP_CTI         ! Imposed values for topographic index statistics from another PGD file
LOGICAL                  :: LMEB             ! MEB
LOGICAL                  :: LIMP_PERM        ! Imposed maps of permafrost distribution
REAL                     :: XUNIF_SAND    ! uniform value of sand fraction
REAL                     :: XUNIF_CLAY    ! uniform value of clay fraction
REAL                     :: XUNIF_RUNOFFB ! uniform value of subgrid runoff coefficient
REAL                     :: XUNIF_WDRAIN  ! uniform value of subgrid drainage coefficient
REAL                     :: XUNIF_PERM    ! uniform value of permafrost distribution
REAL                     :: XUNIF_PH      ! uniform value of pH
REAL                     :: XUNIF_FERT    ! uniform value of fertilisation rate
!
REAL, DIMENSION(150)     :: XSOILGRID     ! Soil layer thickness for DIF
!
 CHARACTER(LEN=28)        :: YSOC_TOP      ! file name for organic carbon expressed in kg/m2
 CHARACTER(LEN=28)        :: YSOC_SUB      ! file name for organic carbon expressed in kg/m2
 CHARACTER(LEN=6)         :: YSOCFILETYPE  ! organic carbon data file type
REAL                     :: XUNIF_SOC_TOP ! uniform value of organic carbon (kg/m2)
REAL                     :: XUNIF_SOC_SUB ! uniform value of organic carbon (kg/m2)
LOGICAL                  :: LIMP_SOC      ! Imposed maps of organic carbon
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
NAMELIST/nam_isba/ npatch, nground_layer, cisba, cpedo_function, cphoto,   &
                   ltr_ml, calbedo, xrm_patch, yclay, yclayfiletype, xunif_clay,  &
                   limp_clay, ysand, ysandfiletype, xunif_sand, limp_sand, &
                   ysoc_top, ysoc_sub, ysocfiletype, xunif_soc_top,        &
                   xunif_soc_sub, limp_soc, ycti, yctifiletype, limp_cti,  &
                   yperm, ypermfiletype, xunif_perm, limp_perm, lmeb,      &                   
                   yrunoffb, yrunoffbfiletype, xunif_runoffb,              &
                   ywdrain,  ywdrainfiletype,  xunif_wdrain, xsoilgrid,    &
                   yph, yphfiletype, xunif_ph, yfert, yfertfiletype,       &
                   xunif_fert   
!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations of defaults
!             ---------------------------
!
!#####################
IF (lhook) CALL dr_hook('READ_NAM_PGD_ISBA',0,zhook_handle)
npatch         = 1
nground_layer  = nundef
CISBA          = '3-L'
CPEDO_FUNCTION = 'CH78'
CPHOTO         = 'NON'
ltr_ml         = .false.
CALBEDO        = 'DRY '
xsoilgrid(:)   = xundef
xrm_patch      = 0.0
!#####################
!
xunif_clay       = 0.33
xunif_sand       = 0.33
xunif_soc_top    = xundef
xunif_soc_sub    = xundef
xunif_runoffb    = 0.5
xunif_wdrain     = 0.
xunif_perm       = xundef
xunif_ph         = xundef
xunif_fert       = xundef
!
yclay            = '                          '
ysand            = '                          '
ysoc_top         = '                          '
ysoc_sub         = '                          '
ycti             = '                          '
yperm            = '                          '
yrunoffb         = '                          '
ywdrain          = '                          '
yph              = '                          '
yfert            = '                          '
!
yclayfiletype    = '      '
ysandfiletype    = '      '
ysocfiletype     = '      '
yctifiletype     = '      '
ypermfiletype    = '      '
yrunoffbfiletype = '      '
ywdrainfiletype  = '      ' 
yphfiletype      = '      '
yphfiletype      = '      '
!
limp_clay        = .false.
limp_sand        = .false.
limp_soc         = .false.
limp_cti         = .false.
lmeb             = .false.
limp_perm        = .false.
!
 CALL get_luout(hprogram,iluout)
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of namelist
!             -------------------
!
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_ISBA',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_isba)
!
 CALL close_namelist(hprogram,ilunam)
!
!-------------------------------------------------------------------------------
!
kpatch           = npatch           ! number of patches
kground_layer    = nground_layer    ! number of soil layers
psoilgrid        = xsoilgrid        ! soil layer tickness for DIF
hisba            = cisba            ! ISBA option
hpedotf          = cpedo_function   ! Pedo-transfert function for DIF
hphoto           = cphoto           ! photosynthesis option
otr_ml           = ltr_ml           ! new radiative transfert
halbedo          = calbedo 
prm_patch        = xrm_patch        ! threshol to remove little fractions of patches
hsand            = ysand            ! file name for sand fraction
hclay            = yclay            ! file name for clay fraction
hsoc_top         = ysoc_top         ! file name for organic carbon
hsoc_sub         = ysoc_sub         ! file name for organic carbon
hcti             = ycti             ! file name for topographic index
hperm            = yperm            ! file name for permafrost distribution
hrunoffb         = yrunoffb         ! file name for subgrid runoff
hwdrain          = ywdrain          ! file name for subgrid drainage
hsandfiletype    = ysandfiletype    ! sand data file type
hclayfiletype    = yclayfiletype    ! clay data file type
hsocfiletype     = ysocfiletype     ! organic carbon data file type
hctifiletype     = yctifiletype     ! topographic index data file type
hpermfiletype    = ypermfiletype    ! permafrost distribution data file type
hrunoffbfiletype = yrunoffbfiletype ! subgrid runoff data file type
hwdrainfiletype  = ywdrainfiletype  ! subgrid drainage data file type
punif_sand       = xunif_sand       ! uniform value of sand fraction
punif_clay       = xunif_clay       ! uniform value of clay fraction
punif_soc_top    = xunif_soc_top    ! uniform value of organic carbon top soil
punif_soc_sub    = xunif_soc_sub    ! uniform value of organic carbon sub soil
punif_runoffb    = xunif_runoffb    ! uniform value of subgrid runoff coefficient
punif_wdrain     = xunif_wdrain     ! uniform value of subgrid drainage coefficient
punif_perm       = xunif_perm       ! uniform value of permafrost distribution
oimp_sand        = limp_sand        ! Imposed values for SAND
oimp_clay        = limp_clay        ! Imposed values for CLAY
oimp_soc         = limp_soc         ! Imposed values for organic carbon
oimp_cti         = limp_cti         ! Imposed values for topographic index statistics
oimp_perm        = limp_perm        ! Imposed values for permafrost distribution
omeb             = lmeb             ! MEB
!
hph           = yph           ! file name for pH value
hfert         = yfert         ! file name for fertilisation data
hphfiletype   = yphfiletype   ! pH data file type
hfertfiletype = yfertfiletype ! Fertilisation data file type
punif_ph      = xunif_ph      ! uniform value of pH
punif_fert    = xunif_fert    ! uniform value of fertilisation rate
!
IF (lhook) CALL dr_hook('READ_NAM_PGD_ISBA',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_nam_pgd_isba