pgd_isba_par.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 pgd_isba_par (DTCO, UG, U, USS, DTV, IO, S, KDIM, HPROGRAM)
!     ##############################################################
!
!!**** *PGD_ISBA_PAR* monitor for averaging and interpolations of cover fractions
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    10/12/97
!!
!!       Modified 08/12/05, P. Le Moigne: user defined fields
!!                 05/2012  R. Alkama   : 19 vegtypes rather than 12    
!!       Modified 02/2012,  P. Samuelsson: MEB
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_sso_n, ONLY : sso_t
USE modd_data_isba_n, ONLY : data_isba_t
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY : isba_s_t
!
USE modd_data_cover_par, ONLY : nvegtype
USE modd_surf_par,       ONLY : xundef
!
USE modd_agri,        ONLY : lagrip
!
USE modi_av_pgd
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
USE modi_ini_var_from_data
USE modi_extrapol_fields
!
USE mode_pos_surf
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
TYPE(surf_atm_t), INTENT(INOUT) :: U
TYPE(sso_t), INTENT(INOUT) :: USS
TYPE(data_isba_t), INTENT(INOUT) :: DTV
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(isba_s_t), INTENT(INOUT) :: S
!
INTEGER, INTENT(IN) :: KDIM
!
 CHARACTER(LEN=6),    INTENT(IN)    :: HPROGRAM     ! Type of program
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL, DIMENSION(KDIM,NVEGTYPE)     :: TEST,TEST2,TEST3
INTEGER               :: ILUOUT    ! output listing logical unit
INTEGER               :: ILUNAM    ! namelist file  logical unit
INTEGER               :: IHGROUND_LAYER ! Half number of NGROUND_LAYER
INTEGER               :: IIH       ! Ground layer counter
LOGICAL               :: GFOUND    ! true if namelist is found
!
INTEGER               :: JVEG, JL  ! loop counter on patch
LOGICAL               :: GPAR_STRESS   ! type of stress
!
INTEGER               :: ISIZE_LMEB_PATCH  ! Number of patches with MEB=true
!
!*    0.3    Declaration of namelists
!            ------------------------
!
INTEGER :: NTIME
INTEGER, PARAMETER :: NTIME_MAX    = 36
INTEGER, PARAMETER :: NGROUND_MAX  = 150
INTEGER, PARAMETER :: NVEGTYPE_MAX = 20
!
REAL, DIMENSION(NVEGTYPE_MAX)   :: XSTRESS   ! 1. if defensive /0. if offensive
!
REAL, DIMENSION(NVEGTYPE_MAX)             :: XUNIF_VEGTYPE    ! fractions of each vegtypes
!
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_VEG        ! vegetation fraction
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_LAI        ! LAI
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_Z0         ! roughness length of vegetation
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_EMIS       ! emissivity
!
REAL, DIMENSION(NVEGTYPE_MAX,NGROUND_MAX)   :: XUNIF_DG         ! soil depths
REAL, DIMENSION(NVEGTYPE_MAX,NGROUND_MAX)   :: XUNIF_ROOTFRAC   ! root fraction profiles
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_GROUND_DEPTH! ground depth
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_ROOT_DEPTH ! root depth
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_ROOT_EXTINCTION! root extinction parameter
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_ROOT_LIN   ! root linear parameter
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_DICE       ! soil ice depth for runoff
!
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_RSMIN      ! minimal stomatal resistance
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_GAMMA      ! gamma parameter
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_WRMAX_CF   ! coefficient for interception
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_RGL        !Rgl
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_CV         ! Cv
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_Z0_O_Z0H   ! ratio of roughness lengths
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_ALBNIR_VEG ! albedo of vegetation (near-infra-red)
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_ALBVIS_VEG ! albedo of vegetation (visible)
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_ALBUV_VEG  ! albedo of vegetation (UV)
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_ALBNIR_SOIL! albedo of soil (near-infra-red)
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_ALBVIS_SOIL! albedo of soil (visible)
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_ALBUV_SOIL ! albedo of soil (UV)
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_GMES       ! Gmes
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_BSLAI      ! Biomass over LAI
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_LAIMIN     ! minimum LAI
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_SEFOLD     ! Sefold
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_GC         ! Gc
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_DMAX       ! Dmax
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_F2I        ! F2I
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_H_TREE     ! height of trees
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_RE25       ! soil respiration parameter
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_CE_NITRO   ! CE for nitrogen
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_CF_NITRO   ! CF for nitrogen
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_CNA_NITRO  ! CNA for nitrogen
!
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_IRRIG
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_WATSUP
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_SEED_M
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_SEED_D
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_REAP_M
REAL, DIMENSION(NVEGTYPE_MAX)               :: XUNIF_REAP_D
!
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_Z0LITTER
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_GNDLITTER
REAL, DIMENSION(NVEGTYPE_MAX,NTIME_MAX)     :: XUNIF_H_VEG
!
LOGICAL, DIMENSION(NVEGTYPE_MAX)            :: LUNIF_STRESS     ! stress type
!
REAL, DIMENSION(NGROUND_MAX)   :: XUNIF_CONDSAT
REAL, DIMENSION(NGROUND_MAX)   :: XUNIF_MPOTSAT
REAL, DIMENSION(NGROUND_MAX)   :: XUNIF_BCOEF
REAL, DIMENSION(NGROUND_MAX)   :: XUNIF_WWILT
REAL, DIMENSION(NGROUND_MAX)   :: XUNIF_WFC
REAL, DIMENSION(NGROUND_MAX)   :: XUNIF_WSAT
!
! name of files containing data
!
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)           :: CFNAM_VEGTYPE    ! fractions of each vegtypes
!
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_VEG        ! vegetation fraction
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_LAI        ! LAI
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_Z0         ! roughness length
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_EMIS       ! emissivity
!
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NGROUND_MAX) :: CFNAM_DG         ! soil depth
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NGROUND_MAX) :: CFNAM_ROOTFRAC   ! root fraction profiles
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_GROUND_DEPTH! ground depth
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_ROOT_DEPTH ! root depth
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_ROOT_EXTINCTION! root extinction parameter
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_ROOT_LIN   ! root linear parameter
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_DICE       ! soil ice depth for runoff (m)
!
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_RSMIN      ! minimal stomatal resistance
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_GAMMA      ! gamma parameter
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_WRMAX_CF   ! coefficient for interception
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_RGL        ! Rgl
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_CV         ! Cv
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_Z0_O_Z0H   ! ratio of roughness lengths
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_ALBNIR_VEG ! albedo of vegetation (near-infra-red)
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_ALBVIS_VEG ! albedo of vegetation (visible)
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_ALBUV_VEG  ! albedo of vegetation (UV)
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_ALBNIR_SOIL! albedo of soil (near-infra-red)
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_ALBVIS_SOIL! albedo of soil (visible)
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_ALBUV_SOIL ! albedo of soil (UV)
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_GMES       ! Gmes
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_BSLAI      ! Biomass over LAI
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_LAIMIN     ! minimum LAI
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_SEFOLD     ! e-folding time for senesence
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_GC         ! cuticular conductance
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_DMAX       ! Dmax
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_F2I        ! F2I
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_H_TREE     ! height of trees
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_RE25       ! soil respiration parameter
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_CE_NITRO   ! CE for nitrogen
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_CF_NITRO   ! CF for nitrogen
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_CNA_NITRO  ! CNA for nitrogen
!
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_IRRIG  ! 
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_WATSUP  ! 
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_SEED_M  ! 
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_SEED_D  ! 
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_REAP_M  ! 
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX)             :: CFNAM_REAP_D  !  
!
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_Z0LITTER
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_GNDLITTER
 CHARACTER(LEN=28), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFNAM_H_VEG
!
 CHARACTER(LEN=28), DIMENSION(NGROUND_MAX)   :: CFNAM_CONDSAT
 CHARACTER(LEN=28), DIMENSION(NGROUND_MAX)   :: CFNAM_MPOTSAT
 CHARACTER(LEN=28), DIMENSION(NGROUND_MAX)   :: CFNAM_BCOEF
 CHARACTER(LEN=28), DIMENSION(NGROUND_MAX)   :: CFNAM_WWILT
 CHARACTER(LEN=28), DIMENSION(NGROUND_MAX)   :: CFNAM_WFC
 CHARACTER(LEN=28), DIMENSION(NGROUND_MAX)   :: CFNAM_WSAT
!
! types of file containing data
!
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)           :: CFTYP_VEGTYPE    ! fractions of each vegtypes
!
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_VEG        ! vegetation fraction
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_LAI        ! LAI
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_Z0         ! roughness length
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_EMIS       ! emissivity
!
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NGROUND_MAX) :: CFTYP_DG         ! soil depth
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NGROUND_MAX) :: CFTYP_ROOTFRAC   ! root fraction profiles
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_GROUND_DEPTH! ground depth
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_ROOT_DEPTH ! root depth
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_ROOT_EXTINCTION! root extinction parameter
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_ROOT_LIN   ! root linear parameter
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_DICE       ! soil ice depth for runoff
!
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_RSMIN      ! minimal stomatal resistance
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_GAMMA      ! gamma parameter
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_WRMAX_CF   ! coefficient for interception
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_RGL        ! Rgl
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_CV         ! Cv
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_Z0_O_Z0H   ! ratio of roughness lengths
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_ALBNIR_VEG ! albedo of vegetation (near-infra-red)
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_ALBVIS_VEG ! albedo of vegetation (visible)
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_ALBUV_VEG  ! albedo of vegetation (UV)
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_ALBNIR_SOIL! albedo of soil (near-infra-red)
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_ALBVIS_SOIL! albedo of soil (visible)
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_ALBUV_SOIL ! albedo of soil (UV)
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_GMES       ! Gmes
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_BSLAI      ! Biomass over LAI
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_LAIMIN     ! minimum LAI
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_SEFOLD     ! e-folding time for senesence
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_GC         ! cuticular conductance
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_DMAX       ! Dmax
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_F2I        ! F2I
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_H_TREE     ! height of trees
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_RE25       ! soil respiration parameter
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_CE_NITRO   ! CE for nitrogen
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_CF_NITRO   ! CF for nitrogen
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_CNA_NITRO  ! CNA for nitrogen
!
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_IRRIG 
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_WATSUP
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_SEED_M
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_SEED_D
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_REAP_M
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX)             :: CFTYP_REAP_D
!
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_Z0LITTER
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_GNDLITTER
 CHARACTER(LEN=6), DIMENSION(NVEGTYPE_MAX,NTIME_MAX)   :: CFTYP_H_VEG
!
 CHARACTER(LEN=6), DIMENSION(NGROUND_MAX)   :: CFTYP_CONDSAT
 CHARACTER(LEN=6), DIMENSION(NGROUND_MAX)   :: CFTYP_MPOTSAT
 CHARACTER(LEN=6), DIMENSION(NGROUND_MAX)   :: CFTYP_BCOEF
 CHARACTER(LEN=6), DIMENSION(NGROUND_MAX)   :: CFTYP_WWILT
 CHARACTER(LEN=6), DIMENSION(NGROUND_MAX)   :: CFTYP_WFC
 CHARACTER(LEN=6), DIMENSION(NGROUND_MAX)   :: CFTYP_WSAT
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
NAMELIST/nam_data_isba/ntime, xunif_vegtype, xunif_dg, xunif_rootfrac, xunif_dice,                  &
                         xunif_ground_depth, xunif_root_depth, xunif_root_extinction,               &
                         xunif_root_lin, xunif_lai, xunif_veg, xunif_z0, xunif_emis,                &
                         xunif_rsmin, xunif_gamma, xunif_wrmax_cf, xunif_rgl,                       &
                         xunif_cv, xunif_z0_o_z0h,                                                  &
                         xunif_albnir_veg, xunif_albvis_veg, xunif_albuv_veg,                       &
                         xunif_albnir_soil, xunif_albvis_soil, xunif_albuv_soil,                    &
                         xunif_gmes, xunif_bslai, xunif_laimin, xunif_sefold,                       &
                         xunif_gc, xunif_dmax, xunif_f2i, lunif_stress, xunif_h_tree, xunif_re25,   &
                         xunif_ce_nitro, xunif_cf_nitro, xunif_cna_nitro,                           &
                         xunif_irrig, xunif_watsup, xunif_seed_m, xunif_seed_d, xunif_reap_m, xunif_reap_d, &
                         cfnam_veg,cfnam_lai,cfnam_rsmin,cfnam_gamma,cfnam_wrmax_cf,                &
                         cfnam_rgl,cfnam_cv,cfnam_dg,cfnam_dice,cfnam_z0,cfnam_z0_o_z0h,            &
                         cfnam_albnir_veg,cfnam_albvis_veg,cfnam_albuv_veg,                         &
                         cfnam_albnir_soil,cfnam_albvis_soil,cfnam_albuv_soil,                      &
                         cfnam_emis,                                                                &
                         cfnam_vegtype,cfnam_rootfrac,                                              &
                         cfnam_ground_depth,cfnam_root_depth,cfnam_root_extinction,cfnam_root_lin,  &
                         cfnam_gmes,cfnam_bslai,cfnam_laimin,cfnam_sefold,cfnam_gc,                 &
                         cfnam_dmax,cfnam_f2i, cfnam_h_tree,cfnam_re25,                             &
                         cfnam_ce_nitro,cfnam_cf_nitro,cfnam_cna_nitro,                             &
                         cfnam_irrig, cfnam_watsup, cfnam_seed_m, cfnam_seed_d, cfnam_reap_m, cfnam_reap_d, &                         
                         cftyp_veg,cftyp_lai,cftyp_rsmin,cftyp_gamma,cftyp_wrmax_cf,                &
                         cftyp_rgl,cftyp_cv,cftyp_dg,cftyp_dice,cftyp_z0,cftyp_z0_o_z0h,            &
                         cftyp_albnir_veg,cftyp_albvis_veg,cftyp_albuv_veg,                         &
                         cftyp_albnir_soil,cftyp_albvis_soil,cftyp_albuv_soil,                      &
                         cftyp_emis,                                                                &
                         cftyp_vegtype,cftyp_rootfrac,                                              &
                         cftyp_ground_depth,cftyp_root_depth,cftyp_root_extinction,cftyp_root_lin,  &
                         cftyp_gmes,cftyp_bslai,cftyp_laimin,cftyp_sefold,cftyp_gc,                 &
                         cftyp_dmax,cftyp_f2i, cftyp_h_tree,cftyp_re25,                             &
                         cftyp_ce_nitro,cftyp_cf_nitro,cftyp_cna_nitro,                             &
                         cftyp_irrig, cftyp_watsup, cftyp_seed_m, cftyp_seed_d, cftyp_reap_m, cftyp_reap_d, &                     
                         xunif_z0litter, xunif_gndlitter, xunif_h_veg, cfnam_z0litter,              &
                         cfnam_gndlitter, cfnam_h_veg, cftyp_z0litter, cftyp_gndlitter, cftyp_h_veg, &
                         xunif_condsat, cfnam_condsat, cftyp_condsat, xunif_mpotsat, cfnam_mpotsat, &
                         cftyp_mpotsat, xunif_bcoef, cfnam_bcoef, cftyp_bcoef, xunif_wwilt, cfnam_wwilt, &
                         cftyp_wwilt, xunif_wfc, cfnam_wfc, cftyp_wfc, xunif_wsat, cfnam_wsat, cftyp_wsat

!-------------------------------------------------------------------------------
!
!*    1.      Initializations
!             ---------------
!
IF (lhook) CALL dr_hook('PGD_ISBA_PAR',0,zhook_handle)
!
IF (u%LECOSG) THEN
             !4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
  xstress = (/1.,1.,1.,0.,0.,0.,0.,0.,1.,1.,1.,0.,0.,0.,0.,1.,1.,0.,0.,0./)
ELSE
             !1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20
  xstress = (/1.,1.,1.,0.,1.,0.,1.,0.,1.,0.,0.,0.,0.,0.,1.,0.,1.,0.,0.,0./)
ENDIF
!
dtv%NVEGTYPE = nvegtype
!
ntime                 = 36
xunif_veg             = xundef ! vegetation fraction
xunif_lai             = xundef ! LAI
xunif_rsmin           = xundef ! minimal stomatal resistance
xunif_gamma           = xundef ! gamma parameter
xunif_wrmax_cf        = xundef ! coefficient for interception
xunif_rgl             = xundef ! Rgl
xunif_cv              = xundef ! Cv
xunif_dg              = xundef ! soil depths
xunif_dice            = xundef ! soil ice depths for runoff
xunif_z0              = xundef ! roughness length of vegetation
xunif_z0_o_z0h        = xundef ! ratio of roughness lengths
xunif_albnir_veg      = xundef ! albedo of vegetation (near-infra-red)
xunif_albvis_veg      = xundef ! albedo of vegetation (visible)
xunif_albuv_veg       = xundef ! albedo of vegetation (UV)
xunif_albnir_soil     = xundef ! albedo of soil (near-infra-red)
xunif_albvis_soil     = xundef ! albedo of soil (visible)
xunif_albuv_soil      = xundef ! albedo of soil (UV)
xunif_emis            = xundef ! emissivity of vegetation
xunif_vegtype         = xundef ! fractions of each vegtypes
xunif_rootfrac        = xundef ! root fraction profiles
xunif_ground_depth    = xundef ! ground depth
xunif_root_depth      = xundef ! root depth
xunif_root_extinction = xundef ! root extinction parameter
xunif_root_lin        = xundef ! root linear parameter
xunif_gmes            = xundef ! Gmes
xunif_bslai           = xundef ! Biomass over LAI
xunif_laimin          = xundef ! minimum LAI
xunif_sefold          = xundef ! Sefold
xunif_gc              = xundef ! Gc
xunif_dmax            = xundef ! Dmax
xunif_f2i             = xundef ! F2I
lunif_stress          = .true.! stress type
xunif_h_tree          = xundef ! height of trees
xunif_re25            = xundef ! soil respiration parameter
xunif_ce_nitro        = xundef ! CE for nitrogen
xunif_cf_nitro        = xundef ! CF for nitrogen
xunif_cna_nitro       = xundef ! CNA for nitrogen
!
xunif_irrig           = xundef
xunif_watsup          = xundef
xunif_seed_m          = xundef
xunif_seed_d          = xundef
xunif_reap_m          = xundef
xunif_reap_d          = xundef
!
xunif_z0litter        = xundef
xunif_gndlitter       = xundef
xunif_h_veg           = xundef
!
xunif_condsat         = xundef
xunif_mpotsat         = xundef
xunif_bcoef           = xundef
xunif_wwilt           = xundef
xunif_wfc             = xundef
xunif_wsat            = xundef
!
CFNAM_VEGTYPE (:)     = '                            '
!
CFNAM_VEG  (:,:)      = '                            '
CFNAM_LAI  (:,:)      = '                            '
CFNAM_Z0   (:,:)      = '                            '
CFNAM_EMIS (:,:)      = '                            '
!
CFNAM_DG       (:,:)  = '                            '
CFNAM_ROOTFRAC (:,:)  = '                            '
CFNAM_DICE     (:)    = '                            '
!
CFNAM_GROUND_DEPTH    (:) = '                            '
CFNAM_ROOT_DEPTH      (:) = '                            '
CFNAM_ROOT_EXTINCTION (:) = '                            '
CFNAM_ROOT_LIN        (:) = '                            '
!
CFNAM_RSMIN       (:) = '                            '
CFNAM_GAMMA       (:) = '                            '
CFNAM_WRMAX_CF    (:) = '                            '
CFNAM_RGL         (:) = '                            '
CFNAM_CV          (:) = '                            '
CFNAM_Z0_O_Z0H    (:) = '                            '
CFNAM_ALBNIR_VEG  (:,:) = '                            '
CFNAM_ALBVIS_VEG  (:,:) = '                            '
CFNAM_ALBUV_VEG   (:,:) = '                            '
CFNAM_ALBNIR_SOIL (:,:) = '                            '
CFNAM_ALBVIS_SOIL (:,:) = '                            '
CFNAM_ALBUV_SOIL  (:,:) = '                            '
CFNAM_GMES        (:) = '                            '
CFNAM_BSLAI       (:) = '                            '
CFNAM_LAIMIN      (:) = '                            '
CFNAM_SEFOLD      (:) = '                            '
CFNAM_GC          (:) = '                            '
CFNAM_DMAX        (:) = '                            '
CFNAM_F2I         (:) = '                            '
CFNAM_H_TREE      (:) = '                            '
CFNAM_RE25        (:) = '                            '
CFNAM_CE_NITRO    (:) = '                            '
CFNAM_CF_NITRO    (:) = '                            '
CFNAM_CNA_NITRO   (:) = '                            '
!
CFNAM_IRRIG       (:,:) = '                            '
CFNAM_WATSUP      (:,:) = '                            '
CFNAM_SEED_M      (:) = '                            '
CFNAM_SEED_D      (:) = '                            '
CFNAM_REAP_M      (:) = '                            '
CFNAM_REAP_D      (:) = '                            '
!
CFNAM_Z0LITTER    (:,:) = '                            '
CFNAM_GNDLITTER   (:,:) = '                            '
CFNAM_H_VEG       (:,:) = '                            '
!
CFNAM_CONDSAT     (:) = '                            '
CFNAM_MPOTSAT     (:) = '                            '
CFNAM_BCOEF       (:) = '                            '
CFNAM_WWILT       (:) = '                            '
CFNAM_WFC         (:) = '                            '
CFNAM_WSAT        (:) = '                            '
!
CFTYP_VEGTYPE (:)     = '      '
!
CFTYP_VEG  (:,:)      = '      '
CFTYP_LAI  (:,:)      = '      '
CFTYP_Z0   (:,:)      = '      '
CFTYP_EMIS (:,:)      = '      '
!
CFTYP_DG       (:,:)  = '      '
CFTYP_ROOTFRAC (:,:)  = '      '
CFTYP_DICE     (:)    = '      '
!
CFTYP_GROUND_DEPTH    (:) = '      '
CFTYP_ROOT_DEPTH      (:) = '      '
CFTYP_ROOT_EXTINCTION (:) = '      '
CFTYP_ROOT_LIN        (:) = '      '
!
CFTYP_RSMIN       (:) = '      '
CFTYP_GAMMA       (:) = '      '
CFTYP_WRMAX_CF    (:) = '      '
CFTYP_RGL         (:) = '      '
CFTYP_CV          (:) = '      '
CFTYP_Z0_O_Z0H    (:) = '      '
CFTYP_ALBNIR_VEG  (:,:) = '      '
CFTYP_ALBVIS_VEG  (:,:) = '      '
CFTYP_ALBUV_VEG   (:,:) = '      '
CFTYP_ALBNIR_SOIL (:,:) = '      '
CFTYP_ALBVIS_SOIL (:,:) = '      '
CFTYP_ALBUV_SOIL  (:,:) = '      '
CFTYP_GMES        (:) = '      '
CFTYP_BSLAI       (:) = '      '
CFTYP_LAIMIN      (:) = '      '
CFTYP_SEFOLD      (:) = '      '
CFTYP_GC          (:) = '      '
CFTYP_DMAX        (:) = '      '
CFTYP_F2I         (:) = '      '
CFTYP_H_TREE      (:) = '      '
CFTYP_RE25        (:) = '      '
CFTYP_CE_NITRO    (:) = '      '
CFTYP_CF_NITRO    (:) = '      '
CFTYP_CNA_NITRO   (:) = '      '
!
CFTYP_IRRIG       (:,:) = '      '
CFTYP_WATSUP      (:,:) = '      '
CFTYP_SEED_M      (:) = '      '
CFTYP_SEED_D      (:) = '      '
CFTYP_REAP_M      (:) = '      '
CFTYP_REAP_D      (:) = '      '
!
CFTYP_Z0LITTER    (:,:) = '      '
CFTYP_GNDLITTER   (:,:) = '      '
CFTYP_H_VEG       (:,:) = '      '
!
CFTYP_CONDSAT     (:) = '      '
CFTYP_MPOTSAT     (:) = '      '
CFTYP_BCOEF       (:) = '      '
CFTYP_WWILT       (:) = '      '
CFTYP_WFC         (:) = '      '
CFTYP_WSAT        (:) = '      '
!
isize_lmeb_patch=count(io%LMEB_PATCH(:))
!-------------------------------------------------------------------------------
!
!*    2.      Input file for cover types
!             --------------------------
!
 CALL get_luout(hprogram,iluout)
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_DATA_ISBA',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_data_isba)
!
 CALL close_namelist(hprogram,ilunam)
!
dtv%NTIME = ntime
!
!-------------------------------------------------------------------------------
IF (nvegtype_max < nvegtype) THEN
  WRITE(iluout,*) '------------------------------------'
  WRITE(iluout,*) 'Please update pgd_isba_par.f90 routine : '
  WRITE(iluout,*) 'The maximum number of VEGTYPE  '
  WRITE(iluout,*) 'in the declaration of the namelist variables '
  WRITE(iluout,*) 'must be increased to : ', nvegtype
  WRITE(iluout,*) '------------------------------------'
  CALL abor1_sfx('PGD_ISBA_PAR: MAXIMUM NUMBER OF VEGTYPE MUST BE INCREASED IN NAMELIST DECLARATION')
END IF
!-------------------------------------------------------------------------------
IF (nground_max < io%NGROUND_LAYER) THEN
  WRITE(iluout,*) '------------------------------------'
  WRITE(iluout,*) 'Please update pgd_isba_par.f90 routine : '
  WRITE(iluout,*) 'The maximum number of soil layers  '
  WRITE(iluout,*) 'in the declaration of the namelist variables '
  WRITE(iluout,*) 'must be increased to : ', io%NGROUND_LAYER
  WRITE(iluout,*) '------------------------------------'
  CALL abor1_sfx('PGD_ISBA_PAR: MAXIMUM NUMBER OF SOIL LAYERS MUST BE INCREASED IN NAMELIST DECLARATION')
END IF
!-------------------------------------------------------------------------------
!
IF (ntime/=36 .AND. ntime/=12 .AND. ntime/=2 .AND. ntime/=1) &
   CALL abor1_sfx('PGD_ISBA_PAR: WRONG VALUE FOR NTIME (POSSIBLE VALUES ARE 1, 12 OR 36')
!
!-------------------------------------------------------------------------------
!
!*    3.      Uniform fields are prescribed
!             -----------------------------
!
!-------------------------------------vegtypes-----------------------------------------
!
ALLOCATE(dtv%XPAR_VEGTYPE(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss,  &
                       hprogram,'ARI','VEGTYPE: vegetation type','NAT',cfnam_vegtype,   &
                       cftyp_vegtype,xunif_vegtype,dtv%XPAR_VEGTYPE,dtv%LDATA_VEGTYPE) 
IF (.NOT. dtv%LDATA_VEGTYPE ) THEN
  DO jveg = 1,nvegtype
    CALL av_pgd(dtco, dtv%XPAR_VEGTYPE(:,jveg),s%XCOVER,dtco%XDATA_VEGTYPE(:,jveg),'NAT','ARI',s%LCOVER)
  ENDDO
ENDIF
!
IF (.NOT.io%LECOCLIMAP .AND. .NOT.dtv%LDATA_VEGTYPE) THEN
  !
  WRITE(iluout,*) ' '
  WRITE(iluout,*) '***********************************************************'
  WRITE(iluout,*) '* Error in PGD field preparation of field VEGTYPE         *'
  WRITE(iluout,*) '* There is no prescribed value and no input file          *'
  WRITE(iluout,*) '* Without ECOCLIMAP, this field must be prescribed        *'
  WRITE(iluout,*) '***********************************************************'
  WRITE(iluout,*) ' '
  CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR VEGTYPE')
  !
ELSEIF (dtv%LDATA_VEGTYPE) THEN
  !
  WHERE (dtv%XPAR_VEGTYPE(:,:)==xundef) dtv%XPAR_VEGTYPE(:,:)=0.
  WHERE (dtv%XPAR_VEGTYPE(:,:)/=0.) dtv%XPAR_VEGTYPE(:,:) = dtv%XPAR_VEGTYPE(:,:) / &
                                                    spread(sum(dtv%XPAR_VEGTYPE(:,:),2),2,nvegtype)
  !  
ENDIF
!
!--------------------------------temporal fields-----------------------------------
!
ALLOCATE(dtv%LDATA_VEG(ntime*nvegtype))
ALLOCATE(dtv%XPAR_VEG(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','VEG: vegetation fraction','NAT',cfnam_veg,cftyp_veg,xunif_veg, &
                        dtv%XPAR_VEG,dtv%LDATA_VEG)
IF (all(.NOT.dtv%LDATA_VEG)) DEALLOCATE(dtv%XPAR_VEG)
!
ALLOCATE(dtv%LDATA_LAI(ntime*nvegtype))
ALLOCATE(dtv%XPAR_LAI(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','LAI: leaf area index','NAT',cfnam_lai,cftyp_lai,xunif_lai, &
                        dtv%XPAR_LAI,dtv%LDATA_LAI)
IF (all(.NOT.dtv%LDATA_LAI)) DEALLOCATE(dtv%XPAR_LAI)
!
ALLOCATE(dtv%LDATA_H_VEG(ntime*nvegtype))
ALLOCATE(dtv%XPAR_H_VEG(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','H_VEG: vegetation height','NAT',cfnam_h_veg,cftyp_h_veg,xunif_h_veg, &
                        dtv%XPAR_H_VEG,dtv%LDATA_H_VEG) 
IF (all(.NOT.dtv%LDATA_H_VEG)) DEALLOCATE(dtv%XPAR_H_VEG)
!
ALLOCATE(dtv%LDATA_Z0(ntime*nvegtype))
ALLOCATE(dtv%XPAR_Z0(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'CDN','Z0: roughness length','NAT',cfnam_z0,cftyp_z0,xunif_z0, &
                        dtv%XPAR_Z0,dtv%LDATA_Z0)
IF (all(.NOT.dtv%LDATA_Z0)) DEALLOCATE(dtv%XPAR_Z0)
!
ALLOCATE(dtv%LDATA_EMIS(ntime*nvegtype))
ALLOCATE(dtv%XPAR_EMIS(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','EMIS: emissivity','NAT',cfnam_emis,cftyp_emis,xunif_emis, &
                        dtv%XPAR_EMIS,dtv%LDATA_EMIS)
IF (all(.NOT.dtv%LDATA_EMIS)) DEALLOCATE(dtv%XPAR_EMIS)
!
IF (.NOT.io%LECOCLIMAP .AND. .NOT.(any(dtv%LDATA_VEG) .AND. any(dtv%LDATA_LAI) .AND. &
                                   any(dtv%LDATA_Z0) .AND. any(dtv%LDATA_EMIS))) THEN
  !
  WRITE(iluout,*) ' '
  WRITE(iluout,*) '***********************************************************'
  WRITE(iluout,*) '* Error in PGD field preparation of temporal fields       *'
  WRITE(iluout,*) '* There is no prescribed value and no input file :        *'
  IF (all(.NOT.dtv%LDATA_VEG )) WRITE(iluout,*) '* for VEG                            *'
  IF (all(.NOT.dtv%LDATA_LAI )) WRITE(iluout,*) '* for LAI                            *'
  IF (all(.NOT.dtv%LDATA_Z0  )) WRITE(iluout,*) '* for Z0                             *'
  IF (all(.NOT.dtv%LDATA_EMIS)) WRITE(iluout,*) '* for EMIS                           *'
  WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
  WRITE(iluout,*) '***********************************************************'
  WRITE(iluout,*) ' '
  CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR TEMPORAL PARAMETERS')
  !
ENDIF
!
ALLOCATE(dtv%LDATA_ALBNIR_VEG(ntime*nvegtype))
ALLOCATE(dtv%XPAR_ALBNIR_VEG(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ALBNIR_VEG: NIR albedo of vegetation','NAT',cfnam_albnir_veg,   &
                       cftyp_albnir_veg,xunif_albnir_veg,dtv%XPAR_ALBNIR_VEG,dtv%LDATA_ALBNIR_VEG)
IF (all(.NOT.dtv%LDATA_ALBNIR_VEG)) DEALLOCATE(dtv%XPAR_ALBNIR_VEG)
!
ALLOCATE(dtv%LDATA_ALBVIS_VEG(ntime*nvegtype))
ALLOCATE(dtv%XPAR_ALBVIS_VEG(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ALBVIS_VEG: VIS albedo of vegetation','NAT',cfnam_albvis_veg,   &
                       cftyp_albvis_veg,xunif_albvis_veg,dtv%XPAR_ALBVIS_VEG,dtv%LDATA_ALBVIS_VEG)  
IF (all(.NOT.dtv%LDATA_ALBVIS_VEG)) DEALLOCATE(dtv%XPAR_ALBVIS_VEG)
!
ALLOCATE(dtv%LDATA_ALBUV_VEG(ntime*nvegtype))
ALLOCATE(dtv%XPAR_ALBUV_VEG(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ALBUV_VEG: UV albedo of vegetation','NAT',cfnam_albuv_veg,   &
                       cftyp_albuv_veg,xunif_albuv_veg,dtv%XPAR_ALBUV_VEG,dtv%LDATA_ALBUV_VEG)
IF (all(.NOT.dtv%LDATA_ALBUV_VEG)) DEALLOCATE(dtv%XPAR_ALBUV_VEG)
!
ALLOCATE(dtv%LDATA_ALBNIR_SOIL(ntime*nvegtype))
ALLOCATE(dtv%XPAR_ALBNIR_SOIL(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ALBNIR_SOIL: NIR albedo of SOIL','NAT',cfnam_albnir_soil,   &
                       cftyp_albnir_soil,xunif_albnir_soil,dtv%XPAR_ALBNIR_SOIL,dtv%LDATA_ALBNIR_SOIL)  
IF (all(.NOT.dtv%LDATA_ALBNIR_SOIL)) DEALLOCATE(dtv%XPAR_ALBNIR_SOIL)
!
ALLOCATE(dtv%LDATA_ALBVIS_SOIL(ntime*nvegtype))
ALLOCATE(dtv%XPAR_ALBVIS_SOIL(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ALBVIS_SOIL: VIS albedo of SOIL','NAT',cfnam_albvis_soil,   &
                       cftyp_albvis_soil,xunif_albvis_soil,dtv%XPAR_ALBVIS_SOIL,dtv%LDATA_ALBVIS_SOIL)  
IF (all(.NOT.dtv%LDATA_ALBVIS_SOIL)) DEALLOCATE(dtv%XPAR_ALBVIS_SOIL)
!
ALLOCATE(dtv%LDATA_ALBUV_SOIL(ntime*nvegtype))
ALLOCATE(dtv%XPAR_ALBUV_SOIL(kdim,ntime,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ALBUV_SOIL: UV albedo of SOIL','NAT',cfnam_albuv_soil,   &
                       cftyp_albuv_soil,xunif_albuv_soil,dtv%XPAR_ALBUV_SOIL,dtv%LDATA_ALBUV_SOIL)  
IF (all(.NOT.dtv%LDATA_ALBUV_SOIL)) DEALLOCATE(dtv%XPAR_ALBUV_SOIL)
!
!
! ------------ Begin MEB parameters ---------------------
!
ALLOCATE(dtv%LDATA_GNDLITTER(ntime*nvegtype))
ALLOCATE(dtv%LDATA_Z0LITTER (ntime*nvegtype))
dtv%LDATA_GNDLITTER(:) = .false.
dtv%LDATA_Z0LITTER (:) = .false.
!
IF(isize_lmeb_patch>0) THEN
  !
  ALLOCATE(dtv%XPAR_GNDLITTER(kdim,ntime,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','GNDLITTER: ground litter fraction','NAT',&
                         cfnam_gndlitter,cftyp_gndlitter,xunif_gndlitter,dtv%XPAR_GNDLITTER,dtv%LDATA_GNDLITTER)
  IF (all(.NOT.dtv%LDATA_GNDLITTER)) DEALLOCATE(dtv%XPAR_GNDLITTER)
  !
  ALLOCATE(dtv%XPAR_Z0LITTER(kdim,ntime,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'CDN','Z0LITTER: ground litter roughness length','NAT', &
                         cfnam_z0litter,cftyp_z0litter,xunif_z0litter,dtv%XPAR_Z0LITTER,dtv%LDATA_Z0LITTER)
  IF (all(.NOT.dtv%LDATA_Z0LITTER)) DEALLOCATE(dtv%XPAR_Z0LITTER)
  !
ENDIF
! ------------ End MEB parameters ---------------------
!
!--------------------------------depths fields-----------------------------------
!
ALLOCATE(dtv%XPAR_DG(kdim,io%NGROUND_LAYER,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','DG: ground depth','NAT',cfnam_dg,cftyp_dg,&
                         xunif_dg,dtv%XPAR_DG,dtv%LDATA_DG(1:nvegtype))
IF (all(.NOT.dtv%LDATA_DG)) DEALLOCATE(dtv%XPAR_DG)
!  
ALLOCATE(dtv%XPAR_ROOT_DEPTH(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','ROOT_DEPTH: root depth','NAT',cfnam_root_depth,cftyp_root_depth,&
                       xunif_root_depth,dtv%XPAR_ROOT_DEPTH,dtv%LDATA_ROOT_DEPTH)
IF (all(.NOT.dtv%LDATA_ROOT_DEPTH)) DEALLOCATE(dtv%XPAR_ROOT_DEPTH)
ALLOCATE(dtv%XPAR_GROUND_DEPTH(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','GROUND_DEPTH: ground depth','NAT',cfnam_ground_depth,cftyp_ground_depth,&
      xunif_ground_depth,dtv%XPAR_GROUND_DEPTH,dtv%LDATA_GROUND_DEPTH)
IF (all(.NOT.dtv%LDATA_GROUND_DEPTH)) DEALLOCATE(dtv%XPAR_GROUND_DEPTH)
!
IF(io%CISBA=='DIF')THEN 
  ! 
  ALLOCATE(dtv%XPAR_ROOTFRAC(kdim,io%NGROUND_LAYER,nvegtype))  
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','ROOTFRAC: root fraction','NAT',&
                        cfnam_rootfrac,cftyp_rootfrac,xunif_rootfrac,&
                        dtv%XPAR_ROOTFRAC,dtv%LDATA_ROOTFRAC(1:nvegtype))
  IF (all(.NOT.dtv%LDATA_ROOTFRAC)) DEALLOCATE(dtv%XPAR_ROOTFRAC)
  !
  ALLOCATE(dtv%XPAR_ROOT_EXTINCTION(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','ROOT_EXTINCTION: root extinction','NAT',cfnam_root_extinction,cftyp_root_extinction,&
        xunif_root_extinction,dtv%XPAR_ROOT_EXTINCTION,dtv%LDATA_ROOT_EXTINCTION)
  IF (all(.NOT.dtv%LDATA_ROOT_EXTINCTION)) DEALLOCATE(dtv%XPAR_ROOT_EXTINCTION)
  !      
  ALLOCATE(dtv%XPAR_ROOT_LIN(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','ROOT_LIN: root linear','NAT',cfnam_root_lin,cftyp_root_lin,&
        xunif_root_lin,dtv%XPAR_ROOT_LIN,dtv%LDATA_ROOT_LIN)
  IF (all(.NOT.dtv%LDATA_ROOT_LIN)) DEALLOCATE(dtv%XPAR_ROOT_LIN)
  !
  IF (.NOT.io%LECOCLIMAP) THEN
    IF(all(.NOT.dtv%LDATA_DG) .AND. all(.NOT.dtv%LDATA_ROOTFRAC) .AND. &
       (all(.NOT.dtv%LDATA_ROOT_DEPTH).OR.all(.NOT.dtv%LDATA_ROOT_EXTINCTION).OR.all(.NOT.dtv%LDATA_ROOT_LIN))) THEN
      WRITE(iluout,*) ' '
      WRITE(iluout,*) '****************************************************************************'
      WRITE(iluout,*) '* Without ECOCLIMAP, Error in PGD field preparation for ISBA-DIF           *'
      WRITE(iluout,*) '* There is no prescribed value and no input file                           *'
      WRITE(iluout,*) '*  (1) XUNIF_ROOTFRAC must be given.                                       *'
      WRITE(iluout,*) '*  (2) Other solution, give all these fields:                              *'
      WRITE(iluout,*) '*     - XUNIF_ROOT_DEPTH      (soil root depth)                            *'
      WRITE(iluout,*) '*     - XUNIF_ROOT_EXTINCTION (root extinction parameter [Jackson 1996])   *'
      WRITE(iluout,*) '*     - XUNIF_ROOT_LIN        (0.05 usually; 1=uniform root distribution!!)*'
      WRITE(iluout,*) '****************************************************************************'
      WRITE(iluout,*) ' '
      CALL abor1_sfx("PGD_ISBA_PAR: PROBLEM IN SOIL GRID COMPUTATION") 
    ELSEIF( .NOT.all(io%XSOILGRID(:)==xundef) .AND. &
           (all(.NOT.dtv%LDATA_GROUND_DEPTH)   .OR.all(.NOT.dtv%LDATA_ROOT_DEPTH).OR. &
            all(.NOT.dtv%LDATA_ROOT_EXTINCTION).OR.all(.NOT.dtv%LDATA_ROOT_LIN)  )) THEN
      WRITE(iluout,*) ' '
      WRITE(iluout,*) '****************************************************************************'
      WRITE(iluout,*) '* Without ECOCLIMAP, Error in PGD field preparation for ISBA-DIF           *'
      WRITE(iluout,*) '* There is no prescribed value and no input file.                          *'
      WRITE(iluout,*) '* When XSOILGRID is given, other field are needed :                        *'
      WRITE(iluout,*) '*     - XUNIF_GROUND_DEPTH    (soil ground depth for moisture)             *'
      WRITE(iluout,*) '*     - XUNIF_ROOT_DEPTH      (soil root depth)                            *'
      WRITE(iluout,*) '*     - XUNIF_ROOT_LIN        (0.05 usually; 1=uniform root distribution!!)*'
      WRITE(iluout,*) '*     - XUNIF_ROOT_EXTINCTION (root extinction parameter [Jackson 1996])   *'
      WRITE(iluout,*) '****************************************************************************'
      WRITE(iluout,*) ' '
      CALL abor1_sfx("PGD_ISBA_PAR: PROBLEM IN SOIL GRID COMPUTATION") 
    ENDIF
    IF(all(.NOT.dtv%LDATA_DG) .AND.all(io%XSOILGRID(:)==xundef))THEN
      WRITE(iluout,*) ' '
      WRITE(iluout,*) '****************************************************************************'
      WRITE(iluout,*) '* Without ECOCLIMAP, Error in PGD field preparation for ISBA-DIF           *'
      WRITE(iluout,*) '* There is no prescribed value to compute vertical soil grid.              *'
      WRITE(iluout,*) '* 2 solutions:                                                             *'
      WRITE(iluout,*) '* (1) Give XUNIF_DG in NAM_DATA_ISBA.                                      *'
      WRITE(iluout,*) '*  OR                                                                      *'
      WRITE(iluout,*) '* (2) Give XSOILGRID in NAM_ISBA                                           *'
      WRITE(iluout,*) '****************************************************************************'
      WRITE(iluout,*) ' '       
      CALL abor1_sfx("PGD_ISBA_PAR: PROBLEM IN SOIL GRID COMPUTATION") 
    ENDIF
  ENDIF
  !
ELSE   
  !
  IF ( .NOT.io%LECOCLIMAP .AND. all(.NOT.dtv%LDATA_DG) .AND. &
      (all(.NOT.dtv%LDATA_GROUND_DEPTH).OR.all(.NOT.dtv%LDATA_ROOT_DEPTH)) ) THEN
    WRITE(iluout,*) ' '
    WRITE(iluout,*) '****************************************************************************'
    WRITE(iluout,*) '* Without ECOCLIMAP, Error in PGD field preparation                        *'
    WRITE(iluout,*) '* There is no prescribed value and no input file                           *'
    WRITE(iluout,*) '* XUNIF_DG or both XUNIF_GROUND_DEPTH and XUNIF_ROOT_DEPTH must be given.  *'
    WRITE(iluout,*) '****************************************************************************'
    WRITE(iluout,*) ' '       
    CALL abor1_sfx("PGD_ISBA_PAR: PROBLEM IN SOIL GRID COMPUTATION") 
  ENDIF
  !
ENDIF
!
ALLOCATE(dtv%XPAR_DICE(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','DICE: ice depth for runoff','NAT',cfnam_dice,cftyp_dice,&
                       xunif_dice,dtv%XPAR_DICE,dtv%LDATA_DICE)
!
IF (.NOT.io%LECOCLIMAP.AND.all(.NOT.dtv%LDATA_DICE)) THEN
  dtv%LDATA_DICE(:) = .false.
  IF(io%CISBA/='DIF' .AND. (any(dtv%LDATA_DG).OR.any(dtv%LDATA_ROOT_DEPTH))) THEN
    DO jveg = 1,nvegtype
      IF(dtv%LDATA_DG(jveg))THEN
        dtv%LDATA_DICE(jveg)=.true.
        WHERE(dtv%XPAR_DG(:,2,jveg)/=xundef) dtv%XPAR_DICE(:,jveg) = max(0.2,0.8*dtv%XPAR_DG(:,2,jveg))
      ELSEIF(dtv%LDATA_ROOT_DEPTH(jveg))THEN
        dtv%LDATA_DICE(jveg)=.true.
        WHERE(dtv%XPAR_ROOT_DEPTH(:,jveg)/=xundef) dtv%XPAR_DICE(:,jveg) = max(0.2,0.8*dtv%XPAR_ROOT_DEPTH(:,jveg))
      ENDIF
    ENDDO
  ELSEIF (io%CISBA=='DIF') THEN
    dtv%XPAR_DICE(:,:) = 0.0
    dtv%LDATA_DICE(:)  =.true.
  ELSE
    WRITE(iluout,*) ' '
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) '* Error in PGD field preparation of field DICE            *'
    WRITE(iluout,*) '* There is no prescribed value and no input file          *'
    WRITE(iluout,*) '* Without ECOCLIMAP, this field must be prescribed        *'
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) ' '
    CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR DICE')
  ENDIF
ENDIF  
IF (all(.NOT.dtv%LDATA_DICE)) DEALLOCATE(dtv%XPAR_DICE)
!
!---------------------classical fields---------------------------------------------
!
ALLOCATE(dtv%XPAR_RSMIN(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'INV','RSMIN: minimal stomatal resistance','NAT',cfnam_rsmin,   &
                       cftyp_rsmin,xunif_rsmin,dtv%XPAR_RSMIN,dtv%LDATA_RSMIN)
IF (all(.NOT.dtv%LDATA_RSMIN)) DEALLOCATE(dtv%XPAR_RSMIN)
!
ALLOCATE(dtv%XPAR_GAMMA(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','GAMMA: gamma coefficient','NAT',cfnam_gamma,   &
                       cftyp_gamma,xunif_gamma,dtv%XPAR_GAMMA,dtv%LDATA_GAMMA)
IF (all(.NOT.dtv%LDATA_GAMMA)) DEALLOCATE(dtv%XPAR_GAMMA)
!
ALLOCATE(dtv%XPAR_WRMAX_CF(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','WRMAX_CF: coeff. for max WR','NAT',cfnam_wrmax_cf,   &
                       cftyp_wrmax_cf,xunif_wrmax_cf,dtv%XPAR_WRMAX_CF,dtv%LDATA_WRMAX_CF)
IF (all(.NOT.dtv%LDATA_WRMAX_CF)) DEALLOCATE(dtv%XPAR_WRMAX_CF)
!
ALLOCATE(dtv%XPAR_RGL(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','RGL: max SW rad. for photosynthesis','NAT',cfnam_rgl,   &
                       cftyp_rgl,xunif_rgl,dtv%XPAR_RGL,dtv%LDATA_RGL)  
IF (all(.NOT.dtv%LDATA_RGL)) DEALLOCATE(dtv%XPAR_RGL)
!
ALLOCATE(dtv%XPAR_CV(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'INV','CV: thermal inertia for vegetation','NAT',cfnam_cv,   &
                       cftyp_cv,xunif_cv,dtv%XPAR_CV,dtv%LDATA_CV)  
IF (all(.NOT.dtv%LDATA_CV)) DEALLOCATE(dtv%XPAR_CV)
!
ALLOCATE(dtv%XPAR_Z0_O_Z0H(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','Z0_O_Z0H: ratio of roughness lengths','NAT',cfnam_z0_o_z0h,   &
                       cftyp_z0_o_z0h,xunif_z0_o_z0h,dtv%XPAR_Z0_O_Z0H,dtv%LDATA_Z0_O_Z0H)  
IF (all(.NOT.dtv%LDATA_Z0_O_Z0H)) DEALLOCATE(dtv%XPAR_Z0_O_Z0H)
!
IF (.NOT.io%LECOCLIMAP .AND. .NOT.(any(dtv%LDATA_RSMIN).AND.any(dtv%LDATA_GAMMA).AND.any(dtv%LDATA_WRMAX_CF).AND.&
         any(dtv%LDATA_RGL).AND.any(dtv%LDATA_CV).AND.any(dtv%LDATA_Z0_O_Z0H).AND.any(dtv%LDATA_ALBNIR_VEG).AND. &
         any(dtv%LDATA_ALBVIS_VEG).AND.any(dtv%LDATA_ALBUV_VEG).AND.any(dtv%LDATA_ALBNIR_SOIL).AND.&
         any(dtv%LDATA_ALBVIS_SOIL).AND.any(dtv%LDATA_ALBUV_SOIL))) THEN
  !
  WRITE(iluout,*) ' '
  WRITE(iluout,*) '***********************************************************'
  WRITE(iluout,*) '* Error in PGD field preparation of classical fields      *'
  WRITE(iluout,*) '* There is no prescribed value and no input file :        *'
  IF (all(.NOT.dtv%LDATA_RSMIN       )) WRITE(iluout,*) '* for RSMIN                  *'
  IF (all(.NOT.dtv%LDATA_GAMMA       )) WRITE(iluout,*) '* for GAMMA                  *'
  IF (all(.NOT.dtv%LDATA_WRMAX_CF    )) WRITE(iluout,*) '* for WRMAX_CF               *'
  IF (all(.NOT.dtv%LDATA_RGL         )) WRITE(iluout,*) '* for RGL                    *'
  IF (all(.NOT.dtv%LDATA_CV          )) WRITE(iluout,*) '* for CV                     *'
  IF (all(.NOT.dtv%LDATA_Z0_O_Z0H    )) WRITE(iluout,*) '* for Z0_O_Z0H               *'
  IF (all(.NOT.dtv%LDATA_ALBNIR_VEG  )) WRITE(iluout,*) '* for ALBNIR_VEG             *'
  IF (all(.NOT.dtv%LDATA_ALBVIS_VEG  )) WRITE(iluout,*) '* for ALBVIS_VEG             *'
  IF (all(.NOT.dtv%LDATA_ALBUV_VEG   )) WRITE(iluout,*) '* for ALBUV_VEG              *'
  IF (all(.NOT.dtv%LDATA_ALBNIR_SOIL )) WRITE(iluout,*) '* for ALBNIR_SOIL            *'
  IF (all(.NOT.dtv%LDATA_ALBVIS_SOIL )) WRITE(iluout,*) '* for ALBVIS_SOIL            *'
  IF (all(.NOT.dtv%LDATA_ALBUV_SOIL  )) WRITE(iluout,*) '* for ALBUV_SOIL             *'
  WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
  WRITE(iluout,*) '***********************************************************'
  WRITE(iluout,*) ' '
  CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR CLASSICAL PARAMETERS')
  !
ENDIF
!
!--------------------------------------AGS parameters----------------------------
!
IF (io%CPHOTO/='NON' .OR. &
    (all(.NOT.dtv%LDATA_Z0).AND.(any(dtv%LDATA_LAI).OR.dtv%LDATA_VEGTYPE)) .OR. isize_lmeb_patch>0) THEN
  !
  ALLOCATE(dtv%XPAR_H_TREE(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                        hprogram,'ARI','H_TREE: height of trees','NAT',cfnam_h_tree,   &
                        cftyp_h_tree,xunif_h_tree,dtv%XPAR_H_TREE,dtv%LDATA_H_TREE)  
  IF (all(.NOT.dtv%LDATA_H_TREE)) THEN
    DEALLOCATE(dtv%XPAR_H_TREE)
  ELSE
    IF (u%LECOSG) THEN
      dtv%LDATA_H_TREE(1:3)   = .false.
      dtv%LDATA_H_TREE(13:18) = .false.
      dtv%LDATA_H_TREE(20)    = .false.
    ELSE
      dtv%LDATA_H_TREE (1:3)  = .false.
      dtv%LDATA_H_TREE (7:12) = .false.
      dtv%LDATA_H_TREE(18:19) = .false.
    ENDIF
  ENDIF
  !
ENDIF

IF (io%CPHOTO/='NON' .OR. isize_lmeb_patch>0) THEN
  ALLOCATE(dtv%XPAR_BSLAI(kdim,nvegtype))
 CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                       hprogram,'ARI','BSLAI: biomass over LAI','NAT',cfnam_bslai,   &
                       cftyp_bslai,xunif_bslai,dtv%XPAR_BSLAI,dtv%LDATA_BSLAI)  
  IF (all(.NOT.dtv%LDATA_BSLAI)) DEALLOCATE(dtv%XPAR_BSLAI)
ENDIF
!
IF (.NOT.io%LECOCLIMAP .AND. isize_lmeb_patch>0 ) THEN
  IF (.NOT.(any(dtv%LDATA_H_TREE).AND.any(dtv%LDATA_GNDLITTER).AND.any(dtv%LDATA_Z0LITTER)  &
      .AND. any(dtv%LDATA_BSLAI))) THEN
    WRITE(iluout,*) ' '
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) '* Error in PGD field preparation of MEB fields            *'
    WRITE(iluout,*) '* There is no prescribed value and no input file :        *'
    IF (all(.NOT.dtv%LDATA_GNDLITTER   )) WRITE(iluout,*) '* for GNDLITTER              *'
    IF (all(.NOT.dtv%LDATA_Z0LITTER    )) WRITE(iluout,*) '* for Z0LITTER               *'
    IF (all(.NOT.dtv%LDATA_H_TREE      )) WRITE(iluout,*) '* for H_TREE                 *'
    IF (all(.NOT.dtv%LDATA_BSLAI       )) WRITE(iluout,*) '* for BSLAI                  *'
    WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) ' '
    CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR MEB PARAMETERS')
  ENDIF
  !
ENDIF
!
ALLOCATE(dtv%LDATA_IRRIG (ntime*nvegtype))      
ALLOCATE(dtv%LDATA_WATSUP(ntime*nvegtype))
dtv%LDATA_IRRIG (:) = .false.
dtv%LDATA_WATSUP(:) = .false.
!
IF (io%CPHOTO/='NON') THEN
  !
  ALLOCATE(dtv%XPAR_RE25 (kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','RE25: ecosystem respiration','NAT',cfnam_re25,   &
                         cftyp_re25,xunif_re25,dtv%XPAR_RE25,dtv%LDATA_RE25)  
  IF (all(.NOT.dtv%LDATA_RE25)) DEALLOCATE(dtv%XPAR_RE25)  
  !
  ALLOCATE(dtv%XPAR_LAIMIN(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','LAIMIN: minimum LAI','NAT',cfnam_laimin,   &
                         cftyp_laimin,xunif_laimin,dtv%XPAR_LAIMIN,dtv%LDATA_LAIMIN)  
  IF (all(.NOT.dtv%LDATA_LAIMIN)) DEALLOCATE(dtv%XPAR_LAIMIN)          
  !
  ALLOCATE(dtv%XPAR_SEFOLD(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','SEFOLD: e-folding time for senescence','NAT',cfnam_sefold,   &
                         cftyp_sefold,xunif_sefold,dtv%XPAR_SEFOLD,dtv%LDATA_SEFOLD)  
  IF (all(.NOT.dtv%LDATA_SEFOLD)) DEALLOCATE(dtv%XPAR_SEFOLD)
  !  
  ALLOCATE(dtv%XPAR_GMES(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','GMES: mesophyl conductance','NAT',cfnam_gmes,   &
                         cftyp_gmes,xunif_gmes,dtv%XPAR_GMES,dtv%LDATA_GMES)
  IF (all(.NOT.dtv%LDATA_GMES)) DEALLOCATE(dtv%XPAR_GMES)
  !
  ALLOCATE(dtv%XPAR_GC(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','GC: cuticular conductance','NAT',cfnam_gc,   &
                         cftyp_gc,xunif_gc,dtv%XPAR_GC,dtv%LDATA_GC)  
  IF (all(.NOT.dtv%LDATA_GC)) DEALLOCATE(dtv%XPAR_GC)
  !
  IF (.NOT.io%LECOCLIMAP .AND. .NOT.(any(dtv%LDATA_H_TREE).AND.any(dtv%LDATA_RE25).AND.&
          any(dtv%LDATA_LAIMIN).AND.any(dtv%LDATA_BSLAI).AND.any(dtv%LDATA_SEFOLD).AND.&
          any(dtv%LDATA_GMES).AND.any(dtv%LDATA_GC))) THEN
    !
    WRITE(iluout,*) ' '
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) '* Error in PGD field preparation of AGS fields            *'
    WRITE(iluout,*) '* There is no prescribed value and no input file :        *'
    IF (all(.NOT.dtv%LDATA_H_TREE )) WRITE(iluout,*) '* for H_TREE                      *'
    IF (all(.NOT.dtv%LDATA_RE25   )) WRITE(iluout,*) '* for RE25                        *'    
    IF (all(.NOT.dtv%LDATA_LAIMIN )) WRITE(iluout,*) '* for LAIMIN                      *'    
    IF (all(.NOT.dtv%LDATA_BSLAI  )) WRITE(iluout,*) '* for BSLAI                       *'
    IF (all(.NOT.dtv%LDATA_SEFOLD )) WRITE(iluout,*) '* for SEFOLD                      *'
    IF (all(.NOT.dtv%LDATA_GMES   )) WRITE(iluout,*) '* for GMES                        *'
    IF (all(.NOT.dtv%LDATA_GC     )) WRITE(iluout,*) '* for GC                          *'
    WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) ' '
    CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR AGS PARAMETERS')
    !
  ENDIF
  !  
  !--------------------------------------AGS Stress parameters----------------------------
  !  
  ALLOCATE(dtv%XPAR_F2I(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','F2I: critical normalized soil water content (stress)','NAT',cfnam_f2i,   &
                         cftyp_f2i,xunif_f2i,dtv%XPAR_F2I,dtv%LDATA_F2I)
  IF (all(.NOT.dtv%LDATA_F2I)) DEALLOCATE(dtv%XPAR_F2I)
  !
  ALLOCATE(dtv%XPAR_DMAX(kdim,nvegtype))
  CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                         hprogram,'ARI','DMAX: maximum air saturation deficit','NAT',cfnam_dmax,   &
                         cftyp_dmax,xunif_dmax,dtv%XPAR_DMAX,dtv%LDATA_DMAX)  
  IF (all(.NOT.dtv%LDATA_DMAX)) DEALLOCATE(dtv%XPAR_DMAX)
  !
  ALLOCATE(dtv%LPAR_STRESS(kdim,nvegtype))
  DO jveg=1,nvegtype
    gpar_stress = lunif_stress(jveg)
    IF (xstress(jveg)<1.) gpar_stress = .false.
    IF (xstress(jveg)==1. .AND. .NOT.gpar_stress) dtv%LDATA_STRESS=.true.
    dtv%LPAR_STRESS(:,jveg) = gpar_stress
  ENDDO
  IF (all(.NOT.dtv%LDATA_STRESS)) DEALLOCATE(dtv%LPAR_STRESS)
  !
  IF (.NOT.io%LECOCLIMAP .AND. .NOT.(any(dtv%LDATA_F2I).AND.any(dtv%LDATA_DMAX))) THEN
    !
    WRITE(iluout,*) ' '
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) '* Error in PGD field preparation of AGS Stress fields     *'
    WRITE(iluout,*) '* There is no prescribed value and no input file :        *'
    IF (all(.NOT.dtv%LDATA_F2I  )) WRITE(iluout,*) '* for F2I                           *'
    IF (all(.NOT.dtv%LDATA_DMAX )) WRITE(iluout,*) '* for DMAX                          *'
    WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
    WRITE(iluout,*) '***********************************************************'
    WRITE(iluout,*) ' '
    CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR AGS STRESS PARAMETERS')
    !
  ENDIF
  !    
  !--------------------------------------AGS Nitrogen parameters----------------------------
  !  
  IF (io%CPHOTO=='NIT' .OR. io%CPHOTO=='NCB') THEN
    !
    ALLOCATE(dtv%XPAR_CE_NITRO(kdim,nvegtype))
    CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                           hprogram,'ARI','CE_NITRO: leaf area ratio sensitivity to nitrogen ccion','NAT',&
                           cfnam_ce_nitro, cftyp_ce_nitro,xunif_ce_nitro,dtv%XPAR_CE_NITRO,dtv%LDATA_CE_NITRO)  
    IF (all(.NOT.dtv%LDATA_CE_NITRO)) DEALLOCATE(dtv%XPAR_CE_NITRO)
    !
    ALLOCATE(dtv%XPAR_CF_NITRO(kdim,nvegtype))
    CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                           hprogram,'ARI','CF_NITRO: lethal minimum value of leaf area ratio','NAT',&
                           cfnam_cf_nitro,cftyp_cf_nitro,xunif_cf_nitro,dtv%XPAR_CF_NITRO,dtv%LDATA_CF_NITRO)
    IF (all(.NOT.dtv%LDATA_CF_NITRO)) DEALLOCATE(dtv%XPAR_CF_NITRO)
    !
    ALLOCATE(dtv%XPAR_CNA_NITRO(kdim,nvegtype))
    CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                           hprogram,'ARI','CNA_NITRO: nitrogen ccion of active biomass','NAT',&
                           cfnam_cna_nitro,cftyp_cna_nitro,xunif_cna_nitro,dtv%XPAR_CNA_NITRO,dtv%LDATA_CNA_NITRO)
    IF (all(.NOT.dtv%LDATA_CNA_NITRO)) DEALLOCATE(dtv%XPAR_CNA_NITRO)
    !
    IF (.NOT.io%LECOCLIMAP .AND. &
         .NOT.(any(dtv%LDATA_CE_NITRO).AND.any(dtv%LDATA_CF_NITRO).AND.any(dtv%LDATA_CNA_NITRO))) THEN
      !
      WRITE(iluout,*) ' '
      WRITE(iluout,*) '***********************************************************'
      WRITE(iluout,*) '* Error in PGD field preparation of AGS Nitrogen fields   *'
      WRITE(iluout,*) '* There is no prescribed value and no input file :        *'
      IF (all(.NOT.dtv%LDATA_CE_NITRO  )) WRITE(iluout,*) '* for CE_NITRO                 *'
      IF (all(.NOT.dtv%LDATA_CF_NITRO  )) WRITE(iluout,*) '* for CF_NITRO                 *'
      IF (all(.NOT.dtv%LDATA_CNA_NITRO )) WRITE(iluout,*) '* for CNA_NITRO                *'
      WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
      WRITE(iluout,*) '***********************************************************'
      WRITE(iluout,*) ' '
      CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR AGS NITROGEN PARAMETERS')
      !
    ENDIF
    !    
    IF (lagrip) THEN
      !
      ALLOCATE(dtv%XPAR_IRRIG(kdim,dtv%NTIME,nvegtype))  
      CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                             hprogram,'ARI','IRRIG: flag for irrigation','NAT',&
                             cfnam_irrig,cftyp_irrig,xunif_irrig,dtv%XPAR_IRRIG,dtv%LDATA_IRRIG)
      IF (all(.NOT.dtv%LDATA_IRRIG)) DEALLOCATE(dtv%XPAR_IRRIG)
      !
      ALLOCATE(dtv%XPAR_WATSUP(kdim,dtv%NTIME,nvegtype))
      CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                             hprogram,'ARI','WATSUP: water supply during irr. (mm)','NAT',&
                             cfnam_watsup,cftyp_watsup,xunif_watsup,dtv%XPAR_WATSUP,dtv%LDATA_WATSUP)
      IF (all(.NOT.dtv%LDATA_WATSUP)) DEALLOCATE(dtv%XPAR_WATSUP)
      !
      ALLOCATE(dtv%XPAR_SEED_M(kdim,nvegtype))
      CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                             hprogram,'MAJ','SEED_M: month of seeding','NAT',&
                             cfnam_seed_m,cftyp_seed_m,xunif_seed_m,dtv%XPAR_SEED_M,dtv%LDATA_SEED_M)
      IF (all(.NOT.dtv%LDATA_SEED_M)) DEALLOCATE(dtv%XPAR_SEED_M)
      !
      ALLOCATE(dtv%XPAR_SEED_D(kdim,nvegtype))
      CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                             hprogram,'MAJ','SEED_D: day of seeding','NAT',&
                             cfnam_seed_d,cftyp_seed_d,xunif_seed_d,dtv%XPAR_SEED_D,dtv%LDATA_SEED_D)
      IF (all(.NOT.dtv%LDATA_SEED_D)) DEALLOCATE(dtv%XPAR_SEED_D)
      !
      ALLOCATE(dtv%XPAR_REAP_M(kdim,nvegtype))
      CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                             hprogram,'MAJ','REAP_M: month of reaping','NAT',&
                             cfnam_reap_m,cftyp_reap_m,xunif_reap_m,dtv%XPAR_REAP_M,dtv%LDATA_REAP_M)
      IF (all(.NOT.dtv%LDATA_REAP_M)) DEALLOCATE(dtv%XPAR_REAP_M)
      !
      ALLOCATE(dtv%XPAR_REAP_D(kdim,nvegtype))
      CALL ini_var_from_data(dtco, ug, u, uss, dtv%XPAR_VEGTYPE, &
                             hprogram,'MAJ','REAP_D: day of reaping','NAT',&
                             cfnam_reap_d,cftyp_reap_d,xunif_reap_d,dtv%XPAR_REAP_D,dtv%LDATA_REAP_D)
      IF (all(.NOT.dtv%LDATA_REAP_D)) DEALLOCATE(dtv%XPAR_REAP_D)
      !      
    ENDIF
    !
    IF ((any(dtv%LDATA_IRRIG).OR.any(dtv%LDATA_WATSUP).OR.any(dtv%LDATA_SEED_M).OR.&
         any(dtv%LDATA_SEED_D).OR.any(dtv%LDATA_REAP_M).OR.any(dtv%LDATA_REAP_D)) .AND. &
        (all(.NOT.dtv%LDATA_IRRIG).OR.all(.NOT.dtv%LDATA_WATSUP).OR.all(.NOT.dtv%LDATA_SEED_M).OR.&
         all(.NOT.dtv%LDATA_SEED_D).OR.all(.NOT.dtv%LDATA_REAP_M).OR.all(.NOT.dtv%LDATA_REAP_D))) THEN
      !
      WRITE(iluout,*) ' '
      WRITE(iluout,*) '***********************************************************'
      WRITE(iluout,*) '* Error in PGD field preparation of irrigation parameters   *'
      WRITE(iluout,*) '* There is no prescribed value and no input file :          *'
      WRITE(iluout,*) '* Some are prescribed and some are not. If you prescribe    *'
      WRITE(iluout,*) '* one of IRRIG, WATSUP, SEED_M, SEED_D, REAP_M and REAP_D,  *'
      WRITE(iluout,*) '* you need to prescribe all the others.  *'
      WRITE(iluout,*) '***********************************************************'
      WRITE(iluout,*) ' '
      CALL abor1_sfx('PGD_ISBA_PAR: MISSING PRESCRIBED VALUE OR INPUT FILE FOR IRRIGATION PARAMETERS')
      !
    ENDIF
    !
    IF (.NOT.io%LECOCLIMAP .AND. .NOT.(any(dtv%LDATA_IRRIG).AND.any(dtv%LDATA_WATSUP).AND.&
            any(dtv%LDATA_SEED_M).AND.any(dtv%LDATA_SEED_D).AND.any(dtv%LDATA_REAP_M).AND.&
            any(dtv%LDATA_REAP_D))) THEN
      !
      WRITE(iluout,*) ' '
      WRITE(iluout,*) '***********************************************************'
      WRITE(iluout,*) '* Error in PGD field preparation of irrigation fields      *'
      WRITE(iluout,*) '* There is no prescribed value and no input file :         *'
      IF (all(.NOT.dtv%LDATA_IRRIG   )) WRITE(iluout,*) '* for IRRIG                   *'
      IF (all(.NOT.dtv%LDATA_WATSUP  )) WRITE(iluout,*) '* for WATSUP                  *'
      IF (all(.NOT.dtv%LDATA_SEED_M  )) WRITE(iluout,*) '* for SEED_M                  *'
      IF (all(.NOT.dtv%LDATA_SEED_D  )) WRITE(iluout,*) '* for SEED_D                  *'
      IF (all(.NOT.dtv%LDATA_REAP_M  )) WRITE(iluout,*) '* for REAP_M                  *'
      IF (all(.NOT.dtv%LDATA_REAP_D  )) WRITE(iluout,*) '* for REAP_D                  *'      
      WRITE(iluout,*) '* Without ECOCLIMAP, these fields must be prescribed      *'
      WRITE(iluout,*) '***********************************************************'
      WRITE(iluout,*) ' '
      CALL abor1_sfx('PGD_ISBA_PAR: NO PRESCRIBED VALUE NOR INPUT FILE FOR IRRIGATION PARAMETERS')
      !
    ENDIF
    
    !
  ENDIF
  !
ENDIF
!
!--------------------------------hydrological fields-----------------------------------
!
ALLOCATE(dtv%XPAR_CONDSAT(kdim,io%NGROUND_LAYER))
  CALL ini_var_from_data(dtco, ug, u, uss,  &
                         hprogram,'ARI','CONDSAT: ','NAT',cfnam_condsat,cftyp_condsat,&
                         xunif_condsat,dtv%XPAR_CONDSAT,dtv%LDATA_CONDSAT)
IF (.NOT.dtv%LDATA_CONDSAT) DEALLOCATE(dtv%XPAR_CONDSAT)
!
ALLOCATE(dtv%XPAR_MPOTSAT(kdim,io%NGROUND_LAYER))
  CALL ini_var_from_data(dtco, ug, u, uss,  &
                         hprogram,'ARI','MPOTSAT: ','NAT',cfnam_mpotsat,cftyp_mpotsat,&
                         xunif_mpotsat,dtv%XPAR_MPOTSAT,dtv%LDATA_MPOTSAT)
IF (.NOT.dtv%LDATA_MPOTSAT) DEALLOCATE(dtv%XPAR_MPOTSAT)
!
ALLOCATE(dtv%XPAR_BCOEF(kdim,io%NGROUND_LAYER))
  CALL ini_var_from_data(dtco, ug, u, uss,  &
                         hprogram,'ARI','BCOEF: ','NAT',cfnam_bcoef,cftyp_bcoef,&
                         xunif_bcoef,dtv%XPAR_BCOEF,dtv%LDATA_BCOEF)
IF (.NOT.dtv%LDATA_BCOEF) DEALLOCATE(dtv%XPAR_BCOEF)
!
ALLOCATE(dtv%XPAR_WWILT(kdim,io%NGROUND_LAYER))
  CALL ini_var_from_data(dtco, ug, u, uss,  &
                         hprogram,'ARI','WWILT: ','NAT',cfnam_wwilt,cftyp_wwilt,&
                         xunif_wwilt,dtv%XPAR_WWILT,dtv%LDATA_WWILT)
IF (.NOT.dtv%LDATA_WWILT) DEALLOCATE(dtv%XPAR_WWILT)
!
ALLOCATE(dtv%XPAR_WFC(kdim,io%NGROUND_LAYER))
  CALL ini_var_from_data(dtco, ug, u, uss,  &
                         hprogram,'ARI','WFC: ','NAT',cfnam_wfc,cftyp_wfc,&
                         xunif_wfc,dtv%XPAR_WFC,dtv%LDATA_WFC)
IF (.NOT.dtv%LDATA_WFC) DEALLOCATE(dtv%XPAR_WFC)
!
ALLOCATE(dtv%XPAR_WSAT(kdim,io%NGROUND_LAYER))
  CALL ini_var_from_data(dtco, ug, u, uss,  &
                         hprogram,'ARI','WSAT: ','NAT',cfnam_wsat,cftyp_wsat,&
                         xunif_wsat,dtv%XPAR_WSAT,dtv%LDATA_WSAT)
IF (.NOT.dtv%LDATA_WSAT) DEALLOCATE(dtv%XPAR_WSAT)
!
!----------------------------------------------------------------------------------------
!
IF (io%LECOCLIMAP .AND. dtv%LDATA_VEGTYPE) THEN
  !  
  CALL extrapol_fields(dtco, dtv, kdim, io, s, ug, u, hprogram,iluout)
  !
ENDIF
!

!----------------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PGD_ISBA_PAR',1,zhook_handle)
!
END SUBROUTINE pgd_isba_par