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, 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,       &          
                                   hgw, hgwfiletype, punif_gw, oimp_gw,                     &          
                                   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
!----------------------------------------------------------------------------
!
!*    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
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)   :: hgw           ! file name for groundwater 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)   :: hgwfiletype      ! groundwater 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
REAL,                INTENT(OUT)   :: punif_gw      ! uniform value of groundwater 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
LOGICAL,             INTENT(OUT)   :: oimp_gw       ! 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
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)        :: ygw              ! file name for groundwater map
 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)         :: ygwfiletype      ! groundwater 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
LOGICAL                  :: limp_gw          ! Imposed maps of groundwater 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_gw      ! uniform groundwater 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, 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,      &                   
                   ygw, ygwfiletype, xunif_gw, limp_gw,                    &                   
                   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.
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_gw         = xundef
xunif_ph         = xundef
xunif_fert       = xundef
!
yclay            = '                          '
ysand            = '                          '
ysoc_top         = '                          '
ysoc_sub         = '                          '
ycti             = '                          '
yperm            = '                          '
ygw              = '                          '
yrunoffb         = '                          '
ywdrain          = '                          '
yph              = '                          '
yfert            = '                          '
!
yclayfiletype    = '      '
ysandfiletype    = '      '
ysocfiletype     = '      '
yctifiletype     = '      '
ypermfiletype    = '      '
ygwfiletype      = '      '
yrunoffbfiletype = '      '
ywdrainfiletype  = '      ' 
yphfiletype      = '      '
yphfiletype      = '      '
!
limp_clay        = .false.
limp_sand        = .false.
limp_soc         = .false.
limp_cti         = .false.
lmeb             = .false.
limp_perm        = .false.
limp_gw          = .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
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
hgw            = ygw                ! file name for groundwater 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
hgwfiletype    = ygwfiletype        ! groundwater 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
punif_gw       = xunif_gw           ! uniform value of groundwater 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
oimp_gw        = limp_gw            ! Imposed values for groundwater distribution
!
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