convert_cover_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 convert_cover_isba (DTCO, HALBEDO, &
                                        HISBA,OTR_ML,KDECADE,PCOVER,OCOVER,HPHOTO,   &
                                         HSFTYPE,PVEG,                       &
                                         PLAI,PRSMIN,PGAMMA,PWRMAX_CF,       &
                                         PRGL,PCV,PSOILGRID,PPERM,           &
                                         PDG,KWG_LAYER,PDROOT,PDG2,          &
                                         PD_ICE,PZ0,PZ0_O_Z0H,               &
                                         PALBNIR_VEG,PALBVIS_VEG,PALBUV_VEG, &
                                         PEMIS_ECO,                          &
                                         PVEGTYPE,PROOTFRAC,                 &
                                         PGMES,PBSLAI,PLAIMIN,PSEFOLD,PGC,   &
                                         PDMAX, PF2I, OSTRESS, PH_TREE,PRE25,&
                                         PCE_NITRO, PCF_NITRO, PCNA_NITRO,   &
                                         TPSEED, TPREAP, PWATSUP, PIRRIG,    &
                                         PGNDLITTER, PZ0LITTER, PH_VEG,      &
                                         OWATER_TO_NATURE, OTOWN_TO_ROCK,    &
                                         OGARDEN     )  
!     ##############################################################
!
!!**** *CONVERT_COVER* convert surface cover classes into secondary 
!!                     physiographic variables for ISBA
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original   01/2004
!!    
!!      P Le Moigne 09/2005 AGS modifs of L. Jarlan
!!      B. Decharme 04/2013 Add CDGAVG (average depth)
!!                          Soil depth = Root depth with ISBA-DF
!!                          except for bare soil pft (but limited to 1m)
!!      P Samuelsson 10/2014 MEB
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_data_cover_n, ONLY : data_cover_t, data_cover_init
!
USE modd_data_cover,     ONLY : xdata_lai, xdata_h_tree,                  &
                                  xdata_veg, xdata_z0, xdata_z0_o_z0h,    &
                                  xdata_emis_eco, xdata_gamma, xdata_cv,  &
                                  xdata_rgl, xdata_rsmin,                 &
                                  xdata_albnir_veg, xdata_albvis_veg,     &
                                  xdata_albuv_veg, xdata_dice,            &
                                  xdata_alb_veg_nir, xdata_alb_veg_vis,   &
                                  xdata_alb_soil_nir, xdata_alb_soil_vis, &                                  
                                  xdata_gmes, xdata_bslai, xdata_laimin,  &
                                  xdata_sefold, xdata_gc, xdata_wrmax_cf, &
                                  xdata_stress,                           &
                                  xdata_dmax, xdata_f2i, xdata_re25,      &
                                  xdata_ce_nitro, xdata_cf_nitro,         &
                                  xdata_cna_nitro,                        &
                                  xdata_gmes_st, xdata_bslai_st,          &
                                  xdata_sefold_st, xdata_gc_st,           &
                                  xdata_dmax_st, xdata_watsup,            &
                                  tdata_seed, tdata_reap,xdata_irrig,     &
                                  xdata_root_depth, xdata_ground_depth,   &
                                  xdata_root_extinction, xdata_root_lin,  &
                                  xdata_gndlitter, xdata_z0litter,        &
                                  xdata_h_veg, xdata_lai_all_years, xdata_vegtype
!


USE modd_data_cover_par, ONLY : nvegtype, jpcover, nvt_no, nvt_rock, nvt_snow,   &
                                nvt_tebd, nvt_bone, nvt_trbe, nvt_c3,   &
                                nvt_c4, nvt_irr, nvt_gras, nvt_trog,    &
                                nvt_park, nvt_trbd, nvt_tebe, nvt_tene, &
                                nvt_bobd, nvt_bond, nvt_bogr, nvt_shrb, &
                                nvt_c3w, nvt_c3s, nvt_fltr, nvt_flgr
USE modd_type_date_surf
!
USE modi_av_pgd
USE modi_default_data_cover
USE modi_ini_data_param
!
USE modi_default_lai_eco2_y2002_01
USE modi_default_lai_eco2_y2002_02
USE modi_default_lai_eco2_y2002_03
USE modi_default_lai_eco2_y2002_04
USE modi_default_lai_eco2_y2002_05
USE modi_default_lai_eco2_y2002_06
USE modi_default_lai_eco2_y2002_07
USE modi_default_lai_eco2_y2002_08
USE modi_default_lai_eco2_y2002_09
USE modi_default_lai_eco2_y2002_10
USE modi_default_lai_eco2_y2002_11
USE modi_default_lai_eco2_y2002_12
USE modi_default_lai_eco2_y2002_13
USE modi_default_lai_eco2_y2002_14
USE modi_default_lai_eco2_y2002_15
USE modi_default_lai_eco2_y2002_16
USE modi_default_lai_eco2_y2002_17
USE modi_default_lai_eco2_y2002_18
USE modi_default_lai_eco2_y2002_19
!
USE modi_default_lai_eco2_y2003_01
USE modi_default_lai_eco2_y2003_02
USE modi_default_lai_eco2_y2003_03
USE modi_default_lai_eco2_y2003_04
USE modi_default_lai_eco2_y2003_05
USE modi_default_lai_eco2_y2003_06
USE modi_default_lai_eco2_y2003_07
USE modi_default_lai_eco2_y2003_08
USE modi_default_lai_eco2_y2003_09
USE modi_default_lai_eco2_y2003_10
USE modi_default_lai_eco2_y2003_11
USE modi_default_lai_eco2_y2003_12
USE modi_default_lai_eco2_y2003_13
USE modi_default_lai_eco2_y2003_14
USE modi_default_lai_eco2_y2003_15
USE modi_default_lai_eco2_y2003_16
USE modi_default_lai_eco2_y2003_17
USE modi_default_lai_eco2_y2003_18
USE modi_default_lai_eco2_y2003_19
!
USE modi_default_lai_eco2_y2004_01
USE modi_default_lai_eco2_y2004_02
USE modi_default_lai_eco2_y2004_03
USE modi_default_lai_eco2_y2004_04
USE modi_default_lai_eco2_y2004_05
USE modi_default_lai_eco2_y2004_06
USE modi_default_lai_eco2_y2004_07
USE modi_default_lai_eco2_y2004_08
USE modi_default_lai_eco2_y2004_09
USE modi_default_lai_eco2_y2004_10
USE modi_default_lai_eco2_y2004_11
USE modi_default_lai_eco2_y2004_12
USE modi_default_lai_eco2_y2004_13
USE modi_default_lai_eco2_y2004_14
USE modi_default_lai_eco2_y2004_15
USE modi_default_lai_eco2_y2004_16
USE modi_default_lai_eco2_y2004_17
USE modi_default_lai_eco2_y2004_18
USE modi_default_lai_eco2_y2004_19
!
USE modi_default_lai_eco2_y2005_01
USE modi_default_lai_eco2_y2005_02
USE modi_default_lai_eco2_y2005_03
USE modi_default_lai_eco2_y2005_04
USE modi_default_lai_eco2_y2005_05
USE modi_default_lai_eco2_y2005_06
USE modi_default_lai_eco2_y2005_07
USE modi_default_lai_eco2_y2005_08
USE modi_default_lai_eco2_y2005_09
USE modi_default_lai_eco2_y2005_10
USE modi_default_lai_eco2_y2005_11
USE modi_default_lai_eco2_y2005_12
USE modi_default_lai_eco2_y2005_13
USE modi_default_lai_eco2_y2005_14
USE modi_default_lai_eco2_y2005_15
USE modi_default_lai_eco2_y2005_16
USE modi_default_lai_eco2_y2005_17
USE modi_default_lai_eco2_y2005_18
USE modi_default_lai_eco2_y2005_19
!
USE modi_default_lai_eco2_y2006_01
USE modi_default_lai_eco2_y2006_02
USE modi_default_lai_eco2_y2006_03
USE modi_default_lai_eco2_y2006_04
USE modi_default_lai_eco2_y2006_05
USE modi_default_lai_eco2_y2006_06
USE modi_default_lai_eco2_y2006_07
USE modi_default_lai_eco2_y2006_08
USE modi_default_lai_eco2_y2006_09
USE modi_default_lai_eco2_y2006_10
USE modi_default_lai_eco2_y2006_11
USE modi_default_lai_eco2_y2006_12
USE modi_default_lai_eco2_y2006_13
USE modi_default_lai_eco2_y2006_14
USE modi_default_lai_eco2_y2006_15
USE modi_default_lai_eco2_y2006_16
USE modi_default_lai_eco2_y2006_17
USE modi_default_lai_eco2_y2006_18
USE modi_default_lai_eco2_y2006_19
!
USE modi_default_lai_eco1_01
USE modi_default_lai_eco1_02
USE modi_default_lai_eco1_03
USE modi_default_lai_eco1_04
USE modi_default_lai_eco1_05
USE modi_default_lai_eco1_06
USE modi_default_lai_eco1_07
USE modi_default_lai_eco1_08
USE modi_default_lai_eco1_09
USE modi_default_lai_eco1_10
USE modi_default_lai_eco1_11
USE modi_default_lai_eco1_12
USE modi_default_lai_eco1_13
USE modi_default_lai_eco1_14
USE modi_default_lai_eco1_15
USE modi_default_lai_eco1_16
USE modi_default_lai_eco1_17
USE modi_default_lai_eco1_18
USE modi_default_lai_eco1_19
!
USE modi_ecoclimap2_lai
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
!
 CHARACTER(LEN=*), INTENT(IN) :: HALBEDO
 CHARACTER(LEN=*),       INTENT(IN)    :: HISBA   ! type of soil (Force-Restore OR Diffusion)
 LOGICAL, INTENT(IN) :: OTR_ML
INTEGER,                INTENT(IN)    :: KDECADE
REAL, DIMENSION(:,:),   INTENT(IN)    :: PCOVER
LOGICAL, DIMENSION(:), INTENT(IN)     :: OCOVER
 CHARACTER(LEN=*),       INTENT(IN)    :: HPHOTO  ! type of photosynthesis
 CHARACTER(LEN=*),       INTENT(IN)    :: HSFTYPE ! nature / garden
!
REAL, DIMENSION(:)  ,   OPTIONAL, INTENT(IN)    :: PSOILGRID
REAL, DIMENSION(:)  ,   OPTIONAL, INTENT(IN)    :: PPERM
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PVEG
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PLAI
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PRSMIN
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PGAMMA
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PWRMAX_CF
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PRGL
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PCV
REAL, DIMENSION(:,:,:), OPTIONAL, INTENT(OUT)   :: PDG
INTEGER, DIMENSION(:,:),OPTIONAL, INTENT(OUT)   :: KWG_LAYER
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PDROOT
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PDG2
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PD_ICE
REAL, DIMENSION(:,:,:), OPTIONAL, INTENT(OUT)   :: PROOTFRAC
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PZ0
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PZ0_O_Z0H
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PALBNIR_VEG
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PALBVIS_VEG
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PALBUV_VEG
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PEMIS_ECO
!
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PGNDLITTER
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PZ0LITTER
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PH_VEG
!
REAL, DIMENSION(:,:), OPTIONAL, INTENT(OUT)   :: PVEGTYPE
!
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PGMES
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PRE25
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PBSLAI
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PLAIMIN
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PSEFOLD
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PGC
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PDMAX
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PF2I
LOGICAL, DIMENSION(:,:),OPTIONAL, INTENT(OUT)   :: OSTRESS
!
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PH_TREE
!
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PCE_NITRO
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PCF_NITRO
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PCNA_NITRO
!
TYPE(date_time), DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: TPSEED
TYPE(date_time), DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: TPREAP
!
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PWATSUP
REAL, DIMENSION(:,:),   OPTIONAL, INTENT(OUT)   :: PIRRIG
!
LOGICAL, OPTIONAL, INTENT(IN) :: OWATER_TO_NATURE
LOGICAL, OPTIONAL, INTENT(IN) :: OTOWN_TO_ROCK
LOGICAL, OPTIONAL, INTENT(IN) :: OGARDEN
!
!*    0.2    Declaration of local variables
!            ------------------------------
! calculation of veg from lai in the pixel
!
TYPE(data_cover_t) :: YDTCO
!
REAL, DIMENSION (:,:), ALLOCATABLE           :: ZWORK      ! work array
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZDATA_VEG, ZDATA_LAI, ZDATA_LAI_ALL_YEARS
REAL, DIMENSION(:,:), ALLOCATABLE :: ZDATA_VEGTYPE
!
REAL, DIMENSION (SIZE(OCOVER),NVEGTYPE) :: ZDATA_GROUND_DEPTH, ZDATA_ROOT_DEPTH,  &
                                           ZDATA_ROOT_LIN, ZDATA_ROOT_EXTINCTION
REAL, DIMENSION (SIZE(OCOVER),3,NVEGTYPE) :: ZDATA_DG
!
 CHARACTER(LEN=3)  :: YTREE, YNAT, YLAI, YVEG, YDIF
!
LOGICAL :: GWATER_TO_NATURE, GTOWN_TO_ROCK, GGARDEN
INTEGER :: JLAYER ! loop counter on surface layers
INTEGER :: JVEG, JP   ! loop counter on vegetation types
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
!*    2.      SECONDARY VARIABLES
!             -------------------
!
IF (lhook) CALL dr_hook('CONVERT_COVER_ISBA',0,zhook_handle)
!
gwater_to_nature = .false.
gtown_to_rock = .false.
IF (PRESENT(owater_to_nature)) gwater_to_nature = owater_to_nature
IF (PRESENT(otown_to_rock)) gtown_to_rock = otown_to_rock
IF (PRESENT(ogarden)) ggarden = ogarden
!
IF (ASSOCIATED(dtco%XDATA_WEIGHT)) DEALLOCATE(dtco%XDATA_WEIGHT)
!
IF (hsftype=='NAT') THEN
  ynat='NAT'
  ytree='TRE'
  ylai='LAI'
  yveg='VEG'
  ydif='DVG'
ELSEIF (hsftype=='GRD'.OR.hsftype=='GD'.OR.hsftype=='GR') THEN
  ynat='GRD'
  ytree='GRT'
  ylai='GRL'
  yveg='GRV'
  ydif='GDV'
ENDIF
!
 CALL data_cover_init(ydtco)
ALLOCATE(ydtco%XDATA_VEGTYPE(SIZE(ocover),nvegtype))
ALLOCATE(ydtco%XDATA_NATURE(SIZE(ocover)))
ALLOCATE(ydtco%XDATA_GARDEN(SIZE(ocover)))
ALLOCATE(ydtco%XDATA_TOWN  (SIZE(ocover)))
!
! if we run with ECOSG, and the init file was done with ECOCLIMAP
IF (SIZE(ocover)/=SIZE(xdata_ground_depth,1)) THEN

  ALLOCATE(ydtco%XDATA_WATER (SIZE(ocover)))
  CALL default_data_cover(pdata_vegtype=ydtco%XDATA_VEGTYPE, &
                          pdata_nature=ydtco%XDATA_NATURE,pdata_water=ydtco%XDATA_WATER,&
                          pdata_town=ydtco%XDATA_TOWN,pdata_garden=ydtco%XDATA_GARDEN,&
                          pdata_ground_depth=zdata_ground_depth,&
                          pdata_root_depth=zdata_root_depth)

  IF (gwater_to_nature) THEN
    WHERE(ydtco%XDATA_WATER(:)<1.) ydtco%XDATA_NATURE(:) = ydtco%XDATA_NATURE(:) + ydtco%XDATA_WATER(:)
  ENDIF

  ! save current values for these arrays
  ALLOCATE(zdata_lai_all_years(SIZE(xdata_lai_all_years,1),SIZE(xdata_lai_all_years,2),SIZE(xdata_lai_all_years,3)))
  ALLOCATE(zdata_vegtype(SIZE(xdata_vegtype,1),SIZE(xdata_vegtype,2)))
  ALLOCATE(zdata_veg(SIZE(xdata_veg,1),SIZE(xdata_veg,2),SIZE(xdata_veg,3)))
  ALLOCATE(zdata_lai(SIZE(xdata_lai,1),SIZE(xdata_lai,2),SIZE(xdata_lai,3)))
  zdata_lai_all_years(:,:,:) = xdata_lai_all_years(:,:,:)
  zdata_vegtype(:,:) = xdata_vegtype(:,:)
  zdata_veg(:,:,:) = xdata_veg(:,:,:)
  zdata_lai(:,:,:) = xdata_lai(:,:,:)

  DEALLOCATE(xdata_veg,xdata_lai,xdata_lai_all_years,xdata_vegtype)

  ! re-allocate official arrays
  ALLOCATE(xdata_vegtype(SIZE(ydtco%XDATA_VEGTYPE,1),SIZE(ydtco%XDATA_VEGTYPE,2)))
  ALLOCATE(xdata_veg(SIZE(ocover),36,nvegtype))
  ALLOCATE(xdata_lai(SIZE(ocover),36,nvegtype))
  ALLOCATE(xdata_lai_all_years(SIZE(ocover),36*5,nvegtype))
  ! xdata_vegtype is used by DEFAULT_LAI....
  xdata_vegtype(:,:) = ydtco%XDATA_VEGTYPE
!
  CALL default_lai_eco1_01
  CALL default_lai_eco1_02
  CALL default_lai_eco1_03
  CALL default_lai_eco1_04
  CALL default_lai_eco1_05
  CALL default_lai_eco1_06
  CALL default_lai_eco1_07
  CALL default_lai_eco1_08
  CALL default_lai_eco1_09
  CALL default_lai_eco1_10
  CALL default_lai_eco1_11
  CALL default_lai_eco1_12
  CALL default_lai_eco1_13
  CALL default_lai_eco1_14
  CALL default_lai_eco1_15
  CALL default_lai_eco1_16
  CALL default_lai_eco1_17
  CALL default_lai_eco1_18
  CALL default_lai_eco1_19

  CALL default_lai_eco2_y2002_01
  CALL default_lai_eco2_y2002_02
  CALL default_lai_eco2_y2002_03
  CALL default_lai_eco2_y2002_04
  CALL default_lai_eco2_y2002_05
  CALL default_lai_eco2_y2002_06
  CALL default_lai_eco2_y2002_07
  CALL default_lai_eco2_y2002_08
  CALL default_lai_eco2_y2002_09
  CALL default_lai_eco2_y2002_10
  CALL default_lai_eco2_y2002_11
  CALL default_lai_eco2_y2002_12
  CALL default_lai_eco2_y2002_13
  CALL default_lai_eco2_y2002_14
  CALL default_lai_eco2_y2002_15
  CALL default_lai_eco2_y2002_16
  CALL default_lai_eco2_y2002_17
  CALL default_lai_eco2_y2002_18
  CALL default_lai_eco2_y2002_19
!
  CALL default_lai_eco2_y2003_01
  CALL default_lai_eco2_y2003_02
  CALL default_lai_eco2_y2003_03
  CALL default_lai_eco2_y2003_04
  CALL default_lai_eco2_y2003_05
  CALL default_lai_eco2_y2003_06
  CALL default_lai_eco2_y2003_07
  CALL default_lai_eco2_y2003_08
  CALL default_lai_eco2_y2003_09
  CALL default_lai_eco2_y2003_10
  CALL default_lai_eco2_y2003_11
  CALL default_lai_eco2_y2003_12
  CALL default_lai_eco2_y2003_13
  CALL default_lai_eco2_y2003_14
  CALL default_lai_eco2_y2003_15
  CALL default_lai_eco2_y2003_16
  CALL default_lai_eco2_y2003_17
  CALL default_lai_eco2_y2003_18
  CALL default_lai_eco2_y2003_19
!
  CALL default_lai_eco2_y2004_01
  CALL default_lai_eco2_y2004_02
  CALL default_lai_eco2_y2004_03
  CALL default_lai_eco2_y2004_04
  CALL default_lai_eco2_y2004_05
  CALL default_lai_eco2_y2004_06
  CALL default_lai_eco2_y2004_07
  CALL default_lai_eco2_y2004_08
  CALL default_lai_eco2_y2004_09
  CALL default_lai_eco2_y2004_10
  CALL default_lai_eco2_y2004_11
  CALL default_lai_eco2_y2004_12
  CALL default_lai_eco2_y2004_13
  CALL default_lai_eco2_y2004_14
  CALL default_lai_eco2_y2004_15
  CALL default_lai_eco2_y2004_16
  CALL default_lai_eco2_y2004_17
  CALL default_lai_eco2_y2004_18
  CALL default_lai_eco2_y2004_19
!
  CALL default_lai_eco2_y2005_01
  CALL default_lai_eco2_y2005_02
  CALL default_lai_eco2_y2005_03
  CALL default_lai_eco2_y2005_04
  CALL default_lai_eco2_y2005_05
  CALL default_lai_eco2_y2005_06
  CALL default_lai_eco2_y2005_07
  CALL default_lai_eco2_y2005_08
  CALL default_lai_eco2_y2005_09
  CALL default_lai_eco2_y2005_10
  CALL default_lai_eco2_y2005_11
  CALL default_lai_eco2_y2005_12
  CALL default_lai_eco2_y2005_13
  CALL default_lai_eco2_y2005_14
  CALL default_lai_eco2_y2005_15
  CALL default_lai_eco2_y2005_16
  CALL default_lai_eco2_y2005_17
  CALL default_lai_eco2_y2005_18
  CALL default_lai_eco2_y2005_19
!
  CALL default_lai_eco2_y2006_01
  CALL default_lai_eco2_y2006_02
  CALL default_lai_eco2_y2006_03
  CALL default_lai_eco2_y2006_04
  CALL default_lai_eco2_y2006_05
  CALL default_lai_eco2_y2006_06
  CALL default_lai_eco2_y2006_07
  CALL default_lai_eco2_y2006_08
  CALL default_lai_eco2_y2006_09
  CALL default_lai_eco2_y2006_10
  CALL default_lai_eco2_y2006_11
  CALL default_lai_eco2_y2006_12
  CALL default_lai_eco2_y2006_13
  CALL default_lai_eco2_y2006_14
  CALL default_lai_eco2_y2006_15
  CALL default_lai_eco2_y2006_16
  CALL default_lai_eco2_y2006_17
  CALL default_lai_eco2_y2006_18
  CALL default_lai_eco2_y2006_19

  nvt_c3w  = 0
  nvt_c3s  = 0
  nvt_fltr = 0
  nvt_flgr = 0
  nvt_no   = 1   ! 1  ! no vegetation (smooth)
  nvt_rock = 2   ! 2  ! no vegetation (rocks)
  nvt_snow = 3   ! 3  ! permanent snow and ice
  nvt_tebd = 4   ! 4  ! temperate broadleaf cold-deciduous summergreen (TREE)
  nvt_bone = 5   ! 5  ! boreal needleleaf evergreen  (CONI)
  nvt_trbe = 6   ! 6  ! tropical broadleaf evergreen (EVER)
  nvt_c3   = 7   ! 7  ! C3 cultures types
  nvt_c4   = 8   ! 8  ! C4 cultures types
  nvt_irr  = 9   ! 9  ! irrigated crops
  nvt_gras =10   !10  ! grassland
  nvt_trog =11   !11  ! tropical grassland
  nvt_park =12   !12  ! peat bogs, parks and gardens (irrigated grass)
  nvt_trbd =13   ! 4  ! tropical broadleaf deciduous (TREE)
  nvt_tebe =14   ! 4  ! temperate broadleaf evergreen (TREE)
  nvt_tene =15   ! 5  ! temperate needleleaf evergreen (CONI)
  nvt_bobd =16   ! 4  ! boreal broadleaf cold-deciduous summergreen (TREE)
  nvt_bond =17   ! 5  ! boreal needleleaf cold-deciduous summergreen (CONI)
  nvt_bogr =18   !10  ! boreal grass (GRAS)
  nvt_shrb =19   ! 4  ! shrub (TREE)

  ! calculate XDATA_LAI
  CALL ecoclimap2_lai(dtco%NYEAR)

  ! no more needed
  DEALLOCATE(xdata_lai_all_years,xdata_vegtype)
  ! we give back their values to data_lai_all_years dans data_vegtype
  ALLOCATE(xdata_lai_all_years(SIZE(zdata_lai_all_years,1),SIZE(zdata_lai_all_years,2),SIZE(zdata_lai_all_years,3)))
  ALLOCATE(xdata_vegtype(SIZE(zdata_vegtype,1),SIZE(zdata_vegtype,2)))
  xdata_lai_all_years(:,:,:) = zdata_lai_all_years(:,:,:)
  xdata_vegtype(:,:) = zdata_vegtype(:,:)
  DEALLOCATE(zdata_lai_all_years,zdata_vegtype)
  ! calculated XDATA_VEG
  CALL ini_data_param(proot_lin=zdata_root_lin,proot_extinction=zdata_root_extinction, &
                      plai=xdata_lai,pveg_out=xdata_veg)

ELSE
  ydtco%XDATA_VEGTYPE(:,:) = dtco%XDATA_VEGTYPE(:,:)
  ydtco%XDATA_NATURE(:) = dtco%XDATA_NATURE(:)
  ydtco%XDATA_TOWN(:) = dtco%XDATA_TOWN(:)
  ydtco%XDATA_GARDEN(:) = dtco%XDATA_GARDEN(:)
  zdata_ground_depth(:,:) = xdata_ground_depth(:,:)
  zdata_root_depth(:,:) = xdata_root_depth(:,:)
  zdata_root_lin(:,:) = xdata_root_lin(:,:)
  zdata_root_extinction(:,:) = xdata_root_extinction(:,:)
ENDIF
!
!*    2.1     fields on natural surfaces only, taking into account patches/ 
!             -------------------------------
!
!
! Leaf Aera Index
! ---------------
!
IF (PRESENT(plai)) THEN
  CALL av_pgd(ydtco, plai ,pcover ,xdata_lai(:,kdecade,:),yveg,'ARI',ocover,kdecade=kdecade)
ENDIF
!
!* 1/Rsmin
!
IF (PRESENT(prsmin)) THEN
  IF (SIZE(prsmin)>0) CALL av_pgd(ydtco, prsmin,pcover ,xdata_rsmin,ylai,'INV',ocover,kdecade=kdecade)  
END IF
!
IF (PRESENT(ph_tree)) CALL av_pgd(ydtco, ph_tree ,pcover ,xdata_h_tree(:,:) ,ytree,'ARI',ocover)  
!
DO jveg=1,nvegtype
  IF (PRESENT(pvegtype)) CALL av_pgd(ydtco, pvegtype(:,jveg),pcover ,ydtco%XDATA_VEGTYPE(:,jveg),ynat,'ARI',ocover)  
END DO
!
!
! vegetation fraction
! -------------------
!
IF (PRESENT(pveg)) CALL av_pgd(ydtco, pveg ,pcover ,xdata_veg(:,kdecade,:),ynat,'ARI',ocover)  
!
!
IF (PRESENT(pgndlitter)) CALL av_pgd(ydtco, pgndlitter ,pcover ,xdata_gndlitter(:,kdecade,:),ynat,'ARI',ocover)
!
! roughness length
! ----------------
!
IF (PRESENT(pz0)) CALL av_pgd(ydtco, pz0 ,pcover ,xdata_z0(:,kdecade,:),ynat,'CDN',ocover)  
!
IF (PRESENT(pz0_o_z0h)) CALL av_pgd(ydtco, pz0_o_z0h ,pcover ,xdata_z0_o_z0h(:,:),ynat,'ARI',ocover)  
!
IF (PRESENT(pz0litter)) CALL av_pgd(ydtco, pz0litter ,pcover ,xdata_z0litter(:,kdecade,:),ynat,'CDN',ocover) 
!
!emis-eco
!--------
!
IF (PRESENT(pemis_eco)) CALL av_pgd(ydtco, pemis_eco ,pcover ,xdata_emis_eco(:,kdecade,:),ynat,'ARI',ocover)  
! 
!---------------------------------------------------------------------------------
!
!* other vegetation parameters
!
IF (PRESENT(pgamma)) CALL av_pgd(ydtco, pgamma     ,pcover ,xdata_gamma(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
!
IF (PRESENT(pwrmax_cf)) CALL av_pgd(ydtco, pwrmax_cf  ,pcover ,xdata_wrmax_cf(:,:),yveg,'ARI',ocover,kdecade=kdecade)
!
IF (PRESENT(prgl)) CALL av_pgd(ydtco, prgl       ,pcover ,xdata_rgl(:,:),yveg,'ARI',ocover,kdecade=kdecade)        
IF (PRESENT(pcv)) CALL av_pgd(ydtco, pcv        ,pcover ,xdata_cv(:,:),yveg,'INV',ocover,kdecade=kdecade)  
!
IF (PRESENT(ph_veg)) THEN
 CALL av_pgd(ydtco, ph_veg,pcover,xdata_h_veg(:,kdecade,:),yveg,'ARI',ocover,kdecade=kdecade) 
ENDIF
!

!
!---------------------------------------------------------------------------------
!
!* soil layers
!  -----------
!
IF (PRESENT(pdg)) THEN
!
!* soil layers (and cumulative root fraction for DIF only)
!
  CALL set_cover_dg(SIZE(pdg,1),SIZE(pdg,2),SIZE(pdg,3),PRESENT(pperm),&
                   PRESENT(pdg2),PRESENT(pdroot),PRESENT(kwg_layer),   &
                   PRESENT(prootfrac)          )
!
END IF
!
!---------------------------------------------------------------------------------
!
!* soil ice for runoff
!  -------------------
!
IF (PRESENT(pd_ice)) CALL av_pgd(ydtco, pd_ice,pcover ,xdata_dice(:,:),ynat,'ARI',ocover)
!
!---------------------------------------------------------------------------------
!
IF (PRESENT(palbnir_veg)) THEN
  IF (halbedo=='CM13') THEN
    CALL av_pgd(ydtco, palbvis_veg,pcover,xdata_alb_veg_nir(:,kdecade,:),yveg,'ARI',ocover,kdecade=kdecade)      
  ELSE   
    CALL av_pgd(ydtco, palbnir_veg,pcover ,xdata_albnir_veg(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
  ENDIF
ENDIF
!
IF (PRESENT(palbvis_veg)) THEN
  IF (halbedo=='CM13') THEN
    CALL av_pgd(ydtco, palbvis_veg,pcover,xdata_alb_veg_vis(:,kdecade,:),yveg,'ARI',ocover,kdecade=kdecade)      
  ELSE     
    CALL av_pgd(ydtco, palbvis_veg,pcover ,xdata_albvis_veg(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
  ENDIF
ENDIF
!
IF (PRESENT(palbuv_veg)) THEN
  IF ((halbedo=='CM13'.OR.otr_ml).AND.PRESENT(palbvis_veg)) THEN
    palbuv_veg(:,:)=palbvis_veg(:,:)
  ELSE
    CALL av_pgd(ydtco, palbuv_veg, pcover ,xdata_albuv_veg(:,:),yveg,'ARI',ocover,kdecade=kdecade)
  ENDIF
ENDIF
!  
! parameters for "stress option"
IF (hphoto /='NON' )  THEN

   IF (PRESENT(pgmes)) THEN
     IF (SIZE(pgmes)>0) CALL av_pgd(ydtco, pgmes  ,pcover ,xdata_gmes_st(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF

   IF (PRESENT(pbslai)) THEN
     IF (SIZE(pbslai)>0) CALL av_pgd(ydtco, pbslai ,pcover ,xdata_bslai_st(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF

   IF (PRESENT(psefold)) THEN
     IF (SIZE(psefold)>0) CALL av_pgd(ydtco, psefold,pcover ,xdata_sefold_st(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF

   IF (PRESENT(pgc)) THEN
     IF (SIZE(pgc)>0) CALL av_pgd(ydtco, pgc    ,pcover ,xdata_gc_st(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF

   IF (PRESENT(pdmax)) THEN
     IF (SIZE(pdmax)>0) CALL av_pgd(ydtco, pdmax  ,pcover ,xdata_dmax_st(:,:),ytree,'ARI',ocover,kdecade=kdecade)  
   END IF

ELSE
   IF (PRESENT(pgmes)) THEN
     IF (SIZE(pgmes)>0) CALL av_pgd(ydtco, pgmes  ,pcover ,xdata_gmes(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF
   IF (PRESENT(pbslai)) THEN
     IF (SIZE(pbslai)>0) CALL av_pgd(ydtco, pbslai ,pcover ,xdata_bslai(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF
   IF (PRESENT(psefold)) THEN
     IF (SIZE(psefold)>0) CALL av_pgd(ydtco, psefold,pcover ,xdata_sefold(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF
   IF (PRESENT(pgc)) THEN
     IF (SIZE(pgc)>0) CALL av_pgd(ydtco, pgc    ,pcover ,xdata_gc(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF
   IF (PRESENT(pdmax)) THEN
     IF (SIZE(pdmax)>0) CALL av_pgd(ydtco, pdmax  ,pcover ,xdata_dmax(:,:),ytree,'ARI',ocover,kdecade=kdecade)  
   END IF

ENDIF

IF (PRESENT(pre25)) THEN
  IF (SIZE(pre25)>0) CALL av_pgd(ydtco, pre25  ,pcover ,xdata_re25(:,:),ynat,'ARI',ocover,kdecade=kdecade)  
END IF

IF (PRESENT(plaimin)) THEN
  IF (SIZE(plaimin)>0) CALL av_pgd(ydtco, plaimin,pcover ,xdata_laimin(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
END IF
IF (PRESENT(pce_nitro)) THEN
  IF (SIZE(pce_nitro)>0) &
          CALL av_pgd(ydtco, pce_nitro  ,pcover ,xdata_ce_nitro(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
END IF
IF (PRESENT(pcf_nitro)) THEN
  IF (SIZE(pcf_nitro)>0) CALL av_pgd(ydtco, pcf_nitro  ,pcover ,xdata_cf_nitro(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
END IF
IF (PRESENT(pcna_nitro)) THEN
  IF (SIZE(pcna_nitro)>0) CALL av_pgd(ydtco, pcna_nitro  ,pcover ,xdata_cna_nitro(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
END IF
IF (PRESENT(pf2i)) THEN
  IF (SIZE(pf2i)>0) CALL av_pgd(ydtco, pf2i   ,pcover ,xdata_f2i(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
END IF
!
IF (PRESENT(ostress)) THEN
  IF (SIZE(ostress)>0) THEN
    ALLOCATE(zwork(SIZE(ostress,1),SIZE(ostress,2)))
    CALL av_pgd(ydtco, zwork,pcover ,xdata_stress(:,:),yveg,'ARI',ocover,kdecade=kdecade)
    WHERE (zwork<0.5) 
      ostress = .false.
    ELSEWHERE
      ostress = .true.
    END WHERE
    DEALLOCATE(zwork)
  END IF
END IF
!
IF (hphoto == 'NIT' .OR. hphoto=='NCB')  THEN
   !
   ! date of seeding
   ! ---------------
   !
   IF (PRESENT(tpseed)) THEN
     IF (SIZE(tpseed)>0) &
       CALL av_pgd (tpseed ,pcover ,tdata_seed(:,:),yveg,'MAJ',ocover,kdecade=kdecade)  
   END IF
   !
   ! date of reaping
   ! ---------------
   !
   IF (PRESENT(tpreap)) THEN
     IF (SIZE(tpreap)>0) &
       CALL av_pgd (tpreap ,pcover ,tdata_reap(:,:),yveg,'MAJ',ocover,kdecade=kdecade)  
   END IF
   !
   ! fraction of irrigated surface
   ! ---------------------------
   !
   IF (PRESENT(pirrig)) THEN
     IF (SIZE(pirrig)>0) &
       CALL av_pgd(ydtco, pirrig ,pcover ,xdata_irrig(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF
   !
   ! water supply for irrigation
   ! ---------------------------
   !
   IF (PRESENT(pwatsup)) THEN
     IF (SIZE(pwatsup)>0) &
       CALL av_pgd(ydtco, pwatsup ,pcover ,xdata_watsup(:,:),yveg,'ARI',ocover,kdecade=kdecade)  
   END IF
!
END IF
!
IF (ASSOCIATED(ydtco%XDATA_WEIGHT)) DEALLOCATE(ydtco%XDATA_WEIGHT)
!
! if we run with ECOSG, and the init file was done with ECOCLIMAP
IF (SIZE(ocover)/=SIZE(xdata_ground_depth,1)) THEN
  nvt_c3   = 0
  nvt_park = 0
  nvt_irr  = 0
  nvt_no   = 1   ! 1  ! no vegetation (smooth)
  nvt_rock = 2   ! 2  ! no vegetation (rocks)
  nvt_snow = 3   ! 3  ! permanent snow and ice
  nvt_bobd = 4   ! 4  ! boreal broadleaf cold-deciduous summergreen (TREE)
  nvt_tebd = 5   ! 4  ! temperate broadleaf cold-deciduous summergreen (TREE)
  nvt_trbd = 6   ! 4  ! tropical broadleaf deciduous (TREE)  
  nvt_tebe = 7   ! 4  ! temperate broadleaf evergreen (TREE)  
  nvt_trbe = 8   ! 6  ! tropical broadleaf evergreen (EVER)
  nvt_bone = 9   ! 5  ! boreal needleleaf evergreen  (CONI)  
  nvt_tene =10   ! 5  ! temperate needleleaf evergreen (CONI)
  nvt_bond =11   ! 5  ! boreal needleleaf cold-deciduous summergreen (CONI)  
  nvt_shrb =12   ! 4  ! shrub (TREE)
  nvt_bogr =13   !10  ! boreal grass (GRAS)
  nvt_gras =14   !10  ! grassland
  nvt_trog =15   !11  ! tropical grassland 
  nvt_c3w  =16   ! 7  ! C3W cultures types
  nvt_c3s  =17   ! 7  ! C3S cultures types
  nvt_c4   =18   ! 8  ! C4 cultures types
  nvt_fltr =19   !12  ! flooded trees
  nvt_flgr =20   !12  ! flooded grassland

  DEALLOCATE(xdata_veg,xdata_lai)
  ALLOCATE(xdata_veg(SIZE(zdata_veg,1),SIZE(zdata_veg,2),SIZE(zdata_veg,3)))
  ALLOCATE(xdata_lai(SIZE(zdata_lai,1),SIZE(zdata_lai,2),SIZE(zdata_lai,3)))
  xdata_veg(:,:,:) = zdata_veg(:,:,:)
  xdata_lai(:,:,:) = zdata_lai(:,:,:)
  DEALLOCATE(zdata_lai,zdata_veg)

ENDIF
!
 CALL data_cover_init(ydtco)
!
IF (lhook) CALL dr_hook('CONVERT_COVER_ISBA',1,zhook_handle)
!
!-------------------------------------------------------------------------------
CONTAINS
!-------------------------------------------------------------------------------
!
SUBROUTINE set_cover_dg(KNI,KGROUND,KPATCH,LPERM,LDG2,LDROOT,LWG_LAYER,LROOTFRAC )
!
USE modd_surf_par,       ONLY : xundef, nundef
!
USE modd_reprod_oper,    ONLY : cdgavg, cdgdif
!
USE modi_ini_data_rootfrac
USE modi_ini_data_soil
USE modi_permafrost_depth
!
IMPLICIT NONE
!
INTEGER, INTENT(IN) :: KNI
INTEGER, INTENT(IN) :: KGROUND
INTEGER, INTENT(IN) :: KPATCH
LOGICAL, INTENT(IN) :: LPERM
LOGICAL, INTENT(IN) :: LDG2
LOGICAL, INTENT(IN) :: LDROOT
LOGICAL, INTENT(IN) :: LWG_LAYER
LOGICAL, INTENT(IN) :: LROOTFRAC
!
INTEGER, DIMENSION (KNI,KPATCH) :: IWG_LAYER
REAL, DIMENSION (KNI,KPATCH) :: ZDTOT, ZDROOT      !  work array
REAL, DIMENSION (KNI,KPATCH) :: ZROOT_EXT  !        "
REAL, DIMENSION (KNI,KPATCH) :: ZROOT_LIN  !        "
!
INTEGER :: JPATCH, JJ, JVEGTYPE, JP
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('CONVERT_COVER_ISBA:SET_COVER_DG',0,zhook_handle)
!
zdtot(:,:) = xundef
zdroot(:,:) = xundef
zroot_ext(:,:) = xundef
zroot_lin(:,:) = xundef
iwg_layer(:,:) = nundef
!
!####################################################################################
!
!CDGAVG : old for reprod = 'ARI' Arithmetic average for all depth 
!         recommended    = 'INV' Harmonic average for all depth (default)
!
!CDGDIF : old for reprod = 'SOIL' d3 soil depth from ecoclimap for isba-df
!         recommended    = 'ROOT' d2 soil depth from ecoclimap for isba-df (default)
!
!####################################################################################
!
IF(hisba/='DIF')THEN
  !  
  DO jp = 1,SIZE(zdata_dg,3)
    CALL ini_data_soil(hisba, zdata_dg(:,:,jp),                   &
                     psurf      = ydtco%XDATA_NATURE,             &
                     psurf2     = ydtco%XDATA_GARDEN,             &
                     prootdepth = zdata_root_depth(:,jp),         &
                     psoildepth = zdata_ground_depth(:,jp)       )
  ENDDO
  !
  DO jlayer=1,kground
     CALL av_pgd(ydtco, pdg(:,jlayer,:),pcover,zdata_dg(:,jlayer,:),ynat,'ARI',ocover,kdecade=kdecade)
  ENDDO
  !
ELSE
!
  IF(cdgdif=='ROOT')THEN
    DO jvegtype=1,nvegtype
       IF(jvegtype==nvt_no)THEN
          WHERE(zdata_ground_depth(:,jvegtype)/=xundef)
               zdata_ground_depth(:,jvegtype) = min(1.0,zdata_ground_depth(:,jvegtype))
          ENDWHERE
       ELSEIF(jvegtype/=nvt_rock.AND.jvegtype/=nvt_snow)THEN
         zdata_ground_depth(:,jvegtype) = max(1.0,zdata_root_depth(:,jvegtype))          
       ELSE
         zdata_ground_depth(:,jvegtype) = zdata_root_depth(:,jvegtype)
       ENDIF
    ENDDO
  ENDIF
!
  CALL av_pgd(ydtco, zdtot(:,:),pcover,zdata_ground_depth,ynat,cdgavg,ocover,kdecade=kdecade)
!  
! CALCULATION OF GROUND_DEPTH over Permafrost area
  IF(lperm)THEN
    CALL permafrost_depth(kni,kpatch,pperm,zdtot)
  ENDIF
!  
  IF (ldg2) THEN
     IF(cdgdif=='ROOT')THEN
       CALL av_pgd(ydtco, pdg2(:,:),pcover,zdata_ground_depth,ynat,cdgavg,ocover)
     ELSE
       CALL av_pgd(ydtco, pdg2(:,:),pcover,zdata_root_depth,ynat,cdgavg,ocover)
     ENDIF
  ENDIF
  IF (ldroot .OR. lrootfrac  .OR. (cdgdif=='ROOT')) THEN
    CALL av_pgd(ydtco, zdroot(:,:),pcover,zdata_root_depth,ydif,cdgavg,ocover,kdecade=kdecade)
    IF (ldroot) pdroot(:,:) = zdroot(:,:)
    IF (cdgdif=='ROOT') WHERE(zdroot(:,:).NE.xundef) zdtot(:,:) = max(zdroot(:,:),zdtot(:,:))
  ENDIF
!
  DO jp = 1,SIZE(pdg,3)
    CALL ini_data_soil(hisba, pdg(:,:,jp), psoildepth=zdtot(:,jp), psoilgrid=psoilgrid, &
                     kwg_layer=iwg_layer(:,jp))
  ENDDO
  IF (lwg_layer) kwg_layer(:,:) = iwg_layer(:,:)
!
  IF (lrootfrac) THEN
!      
    CALL av_pgd(ydtco, zroot_lin(:,:),pcover,zdata_root_lin(:,:),ydif,'ARI',ocover,kdecade=kdecade)
    IF (lrootfrac) THEN
      CALL av_pgd(ydtco, zroot_ext(:,:),pcover,zdata_root_extinction(:,:),ydif,'ARI',ocover,kdecade=kdecade)
      DO jp = 1,SIZE(pdg,3)
        CALL ini_data_rootfrac(pdg(:,:,jp),pdroot(:,jp),zroot_ext(:,jp),zroot_lin(:,jp),prootfrac(:,:,jp))
      ENDDO
    ENDIF
!
  ENDIF
!
ENDIF
!
IF (lhook) CALL dr_hook('CONVERT_COVER_ISBA:SET_COVER_DG',1,zhook_handle)
END SUBROUTINE set_cover_dg
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE convert_cover_isba