init_teb_veg_optionsn.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 init_teb_veg_options_n (&
                                    cht, dgmto, tgdo, tvg, &
                                   hprogram)
!#############################################################
!
!!****  *INIT_TEB_TEB_VEG_n* - routine to initialize ISBA
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      A. Lemonsu  *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    09/2009
!!      B. Decharme 07/2011 : read pgd+prep
!!      B. Decharme 04/2013 : delete CTOPREG option (never used)
!!                            water table / surface coupling
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
USE modd_ch_teb_n, ONLY : ch_teb_t
USE modd_diag_misc_teb_option_n, ONLY : diag_misc_teb_options_t
USE modd_teb_garden_option_n, ONLY : teb_garden_options_t
USE modd_teb_veg_n, ONLY : teb_veg_options_t
!
USE modd_read_namelist,   ONLY : lnam_read
!
USE modd_type_date_surf
USE modd_type_snow
!
!


USE modd_data_cover_par,  ONLY: nvegtype
USE modd_surf_par,        ONLY: xundef, nundef
!
USE modd_isba_par,        ONLY : xoptimgrid
!
USE modn_teb_n,           ONLY : xtstep
!
USE modi_default_isba
USE modi_default_ch_dep
USE modi_default_ch_bio_flux
USE modi_read_default_teb_veg_n
USE modi_read_teb_veg_conf_n
USE modi_get_luout
USE modi_read_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
!
!
TYPE(ch_teb_t), INTENT(INOUT) :: cht
TYPE(diag_misc_teb_options_t), INTENT(INOUT) :: dgmto
TYPE(teb_garden_options_t), INTENT(INOUT) :: tgdo
TYPE(teb_veg_options_t), INTENT(INOUT) :: tvg
!
 CHARACTER(LEN=6),                   INTENT(IN)  :: hprogram  ! program calling surf. schemes
!
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: iversion, ibugfix  ! surface version
INTEGER           :: iluout   ! unit of output listing file
INTEGER           :: iresp    ! Error code after redding
 CHARACTER(LEN=12) :: yrecfm   ! Name of the article to be read
 CHARACTER(LEN=4 ) :: ylvl
!
INTEGER :: jlayer ! loop counter on layers
!
REAL                              :: zout_tstep
 CHARACTER(LEN=3)                  :: yrain 
LOGICAL                           :: gcanopy_drag
LOGICAL                           :: gglacier
LOGICAL                           :: gflood
LOGICAL                           :: gwtd
LOGICAL                           :: gvegupd
LOGICAL                           :: gspinupcarbs
LOGICAL                           :: gspinupcarbw
REAL                              :: zspinmaxs
REAL                              :: zspinmaxw
REAL                              :: zco2_start
REAL                              :: zco2_end
INTEGER                           :: inbyearspins
INTEGER                           :: inbyearspinw
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!               Initialisation for IO
!
IF (lhook) CALL dr_hook('INIT_TEB_VEG_OPTIONS_N',0,zhook_handle)
 CALL get_luout(hprogram,iluout)
!
!               Other little things
!
dgmto%LSURF_DIAG_ALBEDO = .false.
!
IF (lnam_read) THEN
 !
 !*       1.     Defaults
 !               --------
 !
 !        1.1. Hard defaults
 !      
 !       Definition of default options for ISBA (in MODD_TEB_VEG_n)
 !       REM - TSTEP, OUT_TSTEP, CANOPY_DRAG are defined as local variables
 !             because they are already in init_teb.f90 (these options are 
 !             forced to the same values for TEB and urban green areas)
 !
 CALL default_isba(xtstep, zout_tstep,                         &
                   tvg%CROUGH, tvg%CRUNOFF, tvg%CALBEDO, tvg%CSCOND,           &
                   tvg%CC1DRY, tvg%CSOILFRZ, tvg%CDIFSFCOND, tvg%CSNOWRES,     &
                   tvg%CCPSURF, tvg%XCGMAX, tvg%XCDRAG, tvg%CKSAT, tvg%LSOC,       &
                   yrain, tvg%CHORT, gglacier, gcanopy_drag,       &
                   gvegupd, gspinupcarbs, gspinupcarbw,        & 
                   zspinmaxs, zspinmaxw, zco2_start, zco2_end, &
                   inbyearspins, inbyearspinw, tvg%LNITRO_DILU     )
 !
 CALL default_ch_bio_flux(cht%LCH_BIO_FLUX)
 !
ENDIF
!        1.2. Defaults from file header
!    
 CALL read_default_teb_veg_n(cht, tvg, &
                             hprogram)
!
 CALL read_teb_veg_conf_n(cht, tvg, &
                          hprogram)
!
!-------------------------------------------------------------------------------
tvg%CRESPSL = 'DEF'
!-------------------------------------------------------------------------------
!
!         Initialisation for IO
!
!*       2.     Definition of version
!               ---------------------
!
yrecfm='VERSION'
 CALL read_surf(&
                hprogram,yrecfm,iversion,iresp)
!
yrecfm='BUG'
 CALL read_surf(&
                hprogram,yrecfm,ibugfix,iresp)
!
!*       2.     Initialisation of ISBA options
!               ------------------------------
!
!
yrecfm='TWN_ISBA'
IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=3) yrecfm='GD_ISBA'
 CALL read_surf(&
                hprogram,yrecfm,tvg%CISBA,iresp)
!
IF (iversion>=7) THEN
  !
  yrecfm='TWN_PEDOTF'
  IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=3) yrecfm='GD_PEDOTF'
  CALL read_surf(&
                hprogram,yrecfm,tvg%CPEDOTF,iresp)
  !
ELSE
  tvg%CPEDOTF = 'CH78'
ENDIF
!
yrecfm='TWN_PHOTO'
IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=3) yrecfm='GD_PHOTO'
 CALL read_surf(&
                hprogram,yrecfm,tvg%CPHOTO,iresp)
!
yrecfm='TWN_LAYER'
IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=3) yrecfm='GD_LAYER'
 CALL read_surf(&
                hprogram,yrecfm,tgdo%NGROUND_LAYER,iresp)
!
!* new radiative transfert
!
IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=2) THEN
  !
  yrecfm='TWN_TR_ML'
  IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=3) yrecfm='GD_TR_ML'
  CALL read_surf(&
                hprogram,yrecfm,tvg%LTR_ML,iresp)
  !
ELSE 
  tvg%LTR_ML = .false.
ENDIF
!
!* Reference grid for DIF
!
IF(tvg%CISBA=='DIF') THEN
  ALLOCATE(tgdo%XSOILGRID(tgdo%NGROUND_LAYER))
  tgdo%XSOILGRID=xundef
  IF (iversion>=8) THEN
     DO jlayer=1,tgdo%NGROUND_LAYER
        WRITE(ylvl,'(I4)') jlayer
        yrecfm='GD_SGRID'//adjustl(ylvl(:len_trim(ylvl)))
        CALL read_surf(&
                hprogram,yrecfm,tgdo%XSOILGRID(jlayer),iresp)
     ENDDO
  ELSEIF (iversion==7 .AND. ibugfix>=2) THEN
    yrecfm='TWN_SOILGRID'
    IF (iversion>7 .OR. iversion==7 .AND. ibugfix>=3) yrecfm='GD_SOILGRID'
    CALL read_surf(&
                hprogram,yrecfm,tgdo%XSOILGRID,iresp,hdir='-')
  ELSE
    tgdo%XSOILGRID(1:tgdo%NGROUND_LAYER)=xoptimgrid(1:tgdo%NGROUND_LAYER)
  ENDIF
ELSE
  ALLOCATE(tgdo%XSOILGRID(0))
ENDIF
!
!* number of biomass pools
!
tvg%NNBIOMASS=1
IF (tvg%CPHOTO=='NIT') tvg%NNBIOMASS=3
!
!-------------------------------------------------------------------------------
!
!
IF (lhook) CALL dr_hook('INIT_TEB_VEG_OPTIONS_N',1,zhook_handle)
END SUBROUTINE init_teb_veg_options_n