ini_data_cover.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.
#ifdef RS6K
@process nooptimize
#endif
!     #########################
      SUBROUTINE ini_data_cover (DTCO, U)
!     #########################
!
!!**** *INI_DATA_COVER* initializes cover-field correspondance arrays
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    15/12/97
!!    F.solmon    01/06/00 adaptation for patch approach
!!    B.Decharme  01/03/09 Arrange cover by user
!!    G.Pigeon      08/12 add ROUGH_WALL/ROUGH_ROOF
!!    V. Masson     04/13 merges Arrange cover & garden use option in arrange_cover routine
!!    R.Alkama      05/15 Add 7 new vegtype (19 rather than 12)
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_type_date_surf
!
USE modd_surfex_mpi,     ONLY : wlog_mpi
!
USE modd_surfex_omp,     ONLY : idc
!
USE modd_surf_par,       ONLY : xundef, nundef
!
!
USE modd_data_cover,     ONLY : xdata_town, xdata_nature, xdata_sea, xdata_water, &
                                  xdata_lai, xdata_vegtype, xdata_h_tree,           &
                                  xdata_alb_veg_vis, xdata_alb_veg_nir,             &
                                  xdata_alb_soil_vis, xdata_alb_soil_nir,           &
                                  xdata_albnir_veg, xdata_albvis_veg,               &
                                  xdata_albuv_veg, xdata_rsmin,                     &
                                  xdata_rgl, xdata_cv, xdata_gamma, xdata_gmes,     &
                                  xdata_gc, xdata_bslai, xdata_sefold, xdata_laimin,&
                                  xdata_dmax, xdata_stress, xdata_f2i,              &
                                  xdata_veg, xdata_green, xdata_z0, xdata_z0_o_z0h, &
                                  xdata_emis_eco, xdata_wrmax_cf,                   &
                                  xdata_ce_nitro,xdata_cf_nitro,xdata_cna_nitro,    &
                                  xdata_ground_depth, xdata_root_depth,             &
                                  xdata_root_extinction, xdata_root_lin,            &
                                  xdata_soilrc_so2, xdata_soilrc_o3,                &
                                  xdata_z0_town, xdata_z0h_town, xdata_alb_roof,    &
                                  xdata_emis_roof, xdata_hc_roof, xdata_tc_roof,    &
                                  xdata_d_roof, xdata_alb_road, xdata_emis_road,    &
                                  xdata_hc_road, xdata_tc_road, xdata_d_road,       &
                                  xdata_alb_wall, xdata_emis_wall, xdata_hc_wall,   &
                                  xdata_tc_wall, xdata_d_wall, xdata_bld_height,    &
                                  xdata_wall_o_hor, xdata_bld, xdata_can_hw_ratio,  &
                                  xdata_garden, xdata_dice,                         &
                                  xdata_h_traffic, xdata_le_traffic,                &
                                  xdata_h_industry, xdata_le_industry, xdata_re25,  &
                                  xdata_gmes_st, xdata_bslai_st, xdata_sefold_st,   &
                                  xdata_gc_st, xdata_dmax_st, tdata_seed,           &
                                  tdata_reap, xdata_watsup, xdata_irrig,            &
                                  xdata_lai_all_years, lread_data_cover,            &
                                  xdata_hc_floor, xdata_tc_floor, xdata_d_floor,    & 
                                  xdata_tcool_target, xdata_theat_target,           &
                                  xdata_f_waste_can, xdata_eff_heat, xdata_qin,     &
                                  xdata_qin_frad, xdata_shgc, xdata_u_win, xdata_gr,&
                                  xdata_shgc_sh, xdata_floor_height, xdata_inf,     &
                                  xdata_f_water_cond, xdata_qin_flat,               &
                                  xdata_hr_target, xdata_v_vent, xdata_cap_sys_heat,&
                                  xdata_cap_sys_rat, xdata_t_adp, xdata_m_sys_rat,  &
                                  xdata_cop_rat, xdata_t_size_max, xdata_t_size_min,&
                                  xdata_shade, xdata_natvent, xdata_rough_roof,     &
                                  xdata_rough_wall, xdata_frac_gr,xdata_residential,&
                                  xdata_emis_panel, xdata_alb_panel,                &
                                  xdata_eff_panel, xdata_frac_panel,                &
                                  xdata_gndlitter, xdata_z0litter, xdata_h_veg
!
USE modd_data_cover_par, ONLY : nvegtype, nvegtype_old, nvegtype_ecosg,   &
                                  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, &
                                  jpcover, ncover, ntype, ndata_road_layer,  &
                                  ndata_wall_layer, ndata_roof_layer,     &
                                  ndata_floor_layer, cnames, nbare_soil,  &
                                  nrock, nsea, nwater, npermsnow, nut_cphr, &
                                  nut_cpmr, nut_cplr, nut_ophr, nut_opmr, &
                                  nut_oplr, nut_lwlr, nut_lalr, nut_spar, &
                                  nut_indu, nvt_c3w, nvt_c3s, nvt_fltr, nvt_flgr 
!
USE modd_write_cover_tex,ONLY : cname, clang
!
!
USE mode_pos_surf
!
USE modi_read_covers_param
USE modi_init_types_param
!
USE modi_abor1_sfx
!
USE modi_default_data_cover
!
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_default_alb_soil_eco1
USE modi_default_alb_soil_eco2
!
USE modi_default_alb_veg_eco1_01
USE modi_default_alb_veg_eco1_02
USE modi_default_alb_veg_eco1_03
USE modi_default_alb_veg_eco1_04
USE modi_default_alb_veg_eco1_05
USE modi_default_alb_veg_eco1_06
USE modi_default_alb_veg_eco1_07
USE modi_default_alb_veg_eco1_08
USE modi_default_alb_veg_eco1_09
USE modi_default_alb_veg_eco1_10
USE modi_default_alb_veg_eco1_11
USE modi_default_alb_veg_eco1_12
USE modi_default_alb_veg_eco1_13
USE modi_default_alb_veg_eco1_14
USE modi_default_alb_veg_eco1_15
USE modi_default_alb_veg_eco1_16
USE modi_default_alb_veg_eco1_17
USE modi_default_alb_veg_eco1_18
USE modi_default_alb_veg_eco1_19
!
USE modi_arrange_cover
USE modi_cover301_573
USE modi_ecoclimap2_lai
USE modi_ini_data_param
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(surf_atm_t), INTENT(INOUT) :: U
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER         :: JCOV, JVEG, JDEC       ! loop counters on covers and decades
INTEGER         :: ICPT_SEA, ICPT_WATER
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!*    0.3    Declaration of namelists
!            ------------------------
!
!-------------------------------------------------------------------------------
                                
IF (lhook) CALL dr_hook('INI_DATA_COVER',0,zhook_handle)
!
IF (u%LECOSG) THEN
  jpcover = sum(ntype)
ELSE
  jpcover = ncover
ENDIF
!
IF (ALLOCATED(xdata_town)) THEN
  IF (SIZE(xdata_town)/=jpcover) THEN
    idc = 0
    DEALLOCATE(xdata_town,xdata_garden,xdata_nature,xdata_water,xdata_sea)
    DEALLOCATE(xdata_lai,xdata_lai_all_years,xdata_vegtype)
    DEALLOCATE(xdata_h_tree,xdata_ground_depth,xdata_root_depth,xdata_dice)
    DEALLOCATE(xdata_root_extinction,xdata_root_lin)
    DEALLOCATE(xdata_albnir_veg,xdata_albvis_veg,xdata_albuv_veg)
    DEALLOCATE(xdata_alb_veg_vis,xdata_alb_veg_nir,xdata_alb_soil_vis,xdata_alb_soil_nir)
    DEALLOCATE(xdata_rsmin,xdata_gamma,xdata_wrmax_cf,xdata_rgl,xdata_cv)
    DEALLOCATE(xdata_gmes,xdata_gmes_st,xdata_re25,xdata_gc,xdata_gc_st)
    DEALLOCATE(xdata_f2i,xdata_bslai,xdata_bslai_st,xdata_dmax,xdata_dmax_st,xdata_stress)
    DEALLOCATE(xdata_sefold,xdata_sefold_st,xdata_laimin,xdata_veg,xdata_green)
    DEALLOCATE(xdata_z0,xdata_z0_o_z0h,xdata_emis_eco,xdata_soilrc_so2,xdata_soilrc_o3)
    DEALLOCATE(xdata_ce_nitro,xdata_cf_nitro,xdata_cna_nitro,tdata_seed,tdata_reap)
    DEALLOCATE(xdata_irrig,xdata_watsup,xdata_gndlitter,xdata_z0litter,xdata_h_veg)
    DEALLOCATE(xdata_z0_town,xdata_alb_roof,xdata_alb_road,xdata_alb_wall)
    DEALLOCATE(xdata_emis_roof,xdata_emis_road,xdata_emis_wall)
    DEALLOCATE(xdata_hc_roof,xdata_hc_road,xdata_hc_wall,xdata_hc_floor)
    DEALLOCATE(xdata_tc_roof,xdata_tc_road,xdata_tc_wall,xdata_tc_floor)
    DEALLOCATE(xdata_d_roof,xdata_d_road,xdata_d_wall,xdata_d_floor)
    DEALLOCATE(xdata_bld_height,xdata_wall_o_hor,xdata_bld,xdata_can_hw_ratio)
    DEALLOCATE(xdata_h_traffic,xdata_le_traffic,xdata_h_industry,xdata_le_industry)
    DEALLOCATE(xdata_tcool_target,xdata_theat_target,xdata_f_waste_can)
    DEALLOCATE(xdata_eff_heat,xdata_qin,xdata_qin_frad,xdata_shgc,xdata_u_win,xdata_gr)
    DEALLOCATE(xdata_shgc_sh,xdata_floor_height,xdata_inf,xdata_f_water_cond)
    DEALLOCATE(xdata_qin_flat,xdata_hr_target,xdata_v_vent,xdata_cap_sys_heat)
    DEALLOCATE(xdata_cap_sys_rat,xdata_t_adp,xdata_m_sys_rat,xdata_cop_rat)
    DEALLOCATE(xdata_t_size_max,xdata_t_size_min,xdata_shade,xdata_natvent)
    DEALLOCATE(xdata_rough_roof,xdata_rough_wall,xdata_residential,xdata_frac_gr)
    DEALLOCATE(xdata_emis_panel,xdata_alb_panel,xdata_eff_panel,xdata_frac_panel)
    DEALLOCATE(nsea,nwater,cnames,cname)
  ENDIF
ENDIF
!
IF (idc==0) THEN
!
!*    1.1    artificial surfaces fraction
!            ----------------------------
!
ALLOCATE(xdata_town(jpcover))
!
xdata_town = 0.
!
ALLOCATE(xdata_garden(jpcover))
!
xdata_garden = 0.
!
!-------------------------------------------------------------------------------
!
!*    1.2    natural and cultivated surfaces fraction
!            ----------------------------------------
!
ALLOCATE(xdata_nature(jpcover))
!
xdata_nature = 0.
!
!-------------------------------------------------------------------------------
!
!*    1.3    inland waters surfaces fraction
!            -------------------------------
!
ALLOCATE(xdata_water(jpcover))
!
xdata_water = 0.
!
!-------------------------------------------------------------------------------
!
!*    1.4    sea surface fraction
!            --------------------
!
ALLOCATE(xdata_sea(jpcover))
!
xdata_sea = 0.
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
!!possible patches correspond to  vegetation types
!
!*    2.0    vegetation type fractions
!
!New name  N   Nold   Comments
!-----------------------------
IF (u%LECOSG) THEN
  nvegtype = nvegtype_ecosg
  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
ELSE
  nvegtype = nvegtype_old 
  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)
ENDIF
!
!*    2.1    leaf area index
!            ---------------
!
ALLOCATE(xdata_lai(jpcover,36,nvegtype))
!
xdata_lai(:,:,:) = xundef
!
!*  LAI from 2002 to 2007
ALLOCATE(xdata_lai_all_years(jpcover,5*36,nvegtype))
!
xdata_lai_all_years(:,:,:) = xundef
!
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
!
ALLOCATE(xdata_vegtype(jpcover,nvegtype))
!
xdata_vegtype(:,:) = 0.
!
!-------------------------------------------------------------------------------
!
!*    2.3   height of trees (m)
!            ---------------
!
ALLOCATE(xdata_h_tree(jpcover,nvegtype))
!
xdata_h_tree(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.4    ground depth
!            ------------
!
ALLOCATE(xdata_ground_depth(jpcover,nvegtype))
!
xdata_ground_depth(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.4    root depth
!            ----------
!
ALLOCATE(xdata_root_depth(jpcover,nvegtype))
!
xdata_root_depth(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.4    ice depth for runoff 
!            --------------------
!
ALLOCATE(xdata_dice(jpcover,nvegtype))
!
xdata_dice(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.4    root extinction
!            ---------------
!
ALLOCATE(xdata_root_extinction(jpcover,nvegtype))
!
xdata_root_extinction(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.4    ponderation coefficient between root formulations
!            -------------------------------------------------
!
ALLOCATE(xdata_root_lin(jpcover,nvegtype))
!
xdata_root_lin(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.5    albnir (veg only)
!            ------
!
ALLOCATE(xdata_albnir_veg(jpcover,nvegtype))
!
xdata_albnir_veg(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.6    albvis (veg only)
!            ------
!
ALLOCATE(xdata_albvis_veg(jpcover,nvegtype))
!
xdata_albvis_veg(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.6    albUV (veg only)
!            -----
!
ALLOCATE(xdata_albuv_veg(jpcover,nvegtype))
!
xdata_albuv_veg(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
ALLOCATE(xdata_alb_veg_vis(jpcover,36,nvegtype))
!
xdata_alb_veg_vis(:,:,:) = xundef
!
ALLOCATE(xdata_alb_veg_nir(jpcover,36,nvegtype))
!
xdata_alb_veg_nir(:,:,:) = xundef
!
ALLOCATE(xdata_alb_soil_vis(jpcover,36,nvegtype))
!
xdata_alb_soil_vis(:,:,:) = xundef
!
ALLOCATE(xdata_alb_soil_nir(jpcover,36,nvegtype))
!
xdata_alb_soil_nir(:,:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.7    Rsmin
!            -----
!
ALLOCATE(xdata_rsmin(jpcover,nvegtype))
!
xdata_rsmin(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.8    Gamma
!            -----
!
ALLOCATE(xdata_gamma(jpcover,nvegtype))
!
xdata_gamma(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.8    Wrmax_cf
!            --------
!
ALLOCATE(xdata_wrmax_cf(jpcover,nvegtype))
!
xdata_wrmax_cf(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.9    Rgl
!            ---
!
ALLOCATE(xdata_rgl(jpcover,nvegtype))
!
xdata_rgl(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.10   Cv
!            --
!
ALLOCATE(xdata_cv(jpcover,nvegtype))
!
xdata_cv(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.11   mesophyll conductance (m s-1)
!            -----------------------------
!
ALLOCATE(xdata_gmes(jpcover,nvegtype))
!
xdata_gmes(:,:) = xundef
!
ALLOCATE(xdata_gmes_st(jpcover,nvegtype))
!
xdata_gmes_st(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.11   Ecosystem respiration (kg/kg.m.s-1)
!            -----------------------------------
!
ALLOCATE(xdata_re25(jpcover,nvegtype))
!
xdata_re25(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.11   cuticular conductance (m s-1)
!            -----------------------------
!
ALLOCATE(xdata_gc(jpcover,nvegtype))
!
xdata_gc(:,:) = xundef
!
ALLOCATE(xdata_gc_st(jpcover,nvegtype))
!
xdata_gc_st(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.11   critical normilized soil water content for stress parameterisation
!            ------------------------------------------------------------------
!
ALLOCATE(xdata_f2i(jpcover,nvegtype))
!
xdata_f2i(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.12   ratio d(biomass)/d(lai) (kg/m2)
!            -----------------------
!
ALLOCATE(xdata_bslai(jpcover,nvegtype))
!
xdata_bslai(:,:) = xundef
!
ALLOCATE(xdata_bslai_st(jpcover,nvegtype))
!
xdata_bslai_st(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.12   maximum air saturation deficit tolerate by vegetation (kg/kg)
!            -------------------------------------------------------------
!
ALLOCATE(xdata_dmax(jpcover,nvegtype))
!
xdata_dmax(:,:) = xundef
!
ALLOCATE(xdata_dmax_st(jpcover,nvegtype))
!
xdata_dmax_st(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.11   Defensive/Offensive strategy
!            ----------------------------
!
ALLOCATE(xdata_stress(jpcover,nvegtype))
!
xdata_stress(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.13   e-folding time for senescence (days)
!            ------------------------------------
!
ALLOCATE(xdata_sefold(jpcover,nvegtype))
!
xdata_sefold(:,:) = xundef
!
ALLOCATE(xdata_sefold_st(jpcover,nvegtype))
!
xdata_sefold_st(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.14   Minimum LAI (m2/m2)
!            -------------------
!
ALLOCATE(xdata_laimin(jpcover,nvegtype))
!
xdata_laimin(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.15   vegetation and greeness fraction fractions
!            ------------------------------------------
!
ALLOCATE(xdata_veg(jpcover,36,nvegtype))
!
xdata_veg(:,:,:) = xundef
!
ALLOCATE(xdata_green(jpcover,36,nvegtype))
!
xdata_green(:,:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.16   z0veg
!            -----
!
ALLOCATE(xdata_z0(jpcover,36,nvegtype))
!
xdata_z0(:,:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.17   z0hveg
!            ------
!
ALLOCATE(xdata_z0_o_z0h(jpcover,nvegtype))
!
xdata_z0_o_z0h(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.18   emissivity
!            ----------
!
ALLOCATE(xdata_emis_eco(jpcover,36,nvegtype))
!
xdata_emis_eco(:,:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.19   for chemistry deposition
!            ------------------------
!
ALLOCATE(xdata_soilrc_so2(jpcover,nvegtype))
ALLOCATE(xdata_soilrc_o3(jpcover,nvegtype))
!
xdata_soilrc_so2(:,:) = xundef
xdata_soilrc_o3(:,:) = xundef
!
!------------------------------------------------------------------------
!
!*    2.20   leaf aera ratio sensitivity to nitrogen concentration
!            -----------------------------------------------------
!
ALLOCATE(xdata_ce_nitro(jpcover,nvegtype))
!
xdata_ce_nitro(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.21   lethal minimum value of leaf area ratio
!            ---------------------------------------
!
ALLOCATE(xdata_cf_nitro(jpcover,nvegtype))
!
xdata_cf_nitro(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.22   nitrogen concentration of active biomass
!            ----------------------------------------
!
ALLOCATE(xdata_cna_nitro(jpcover,nvegtype))
!
xdata_cna_nitro(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    2.24   seeding and reaping dates
!            -------------------------
!
ALLOCATE(tdata_seed(jpcover,nvegtype))
!
tdata_seed(:,:)%TDATE%YEAR  = 9999
tdata_seed(:,:)%TDATE%MONTH = nundef                                
tdata_seed(:,:)%TDATE%DAY   = nundef                                
tdata_seed(:,:)%TIME        = 0.                              
!
ALLOCATE(tdata_reap(jpcover,nvegtype))
!
tdata_reap(:,:)%TDATE%YEAR  = 9999
tdata_reap(:,:)%TDATE%MONTH = nundef                                
tdata_reap(:,:)%TDATE%DAY   = nundef                                
tdata_reap(:,:)%TIME        = 0.                              
!
!-------------------------------------------------------------------------------
!
!*    2.25   irrigated fraction
!            ------------------
!
ALLOCATE(xdata_irrig(jpcover,nvegtype))
!
xdata_irrig(:,:) = 0.                                
!
!-------------------------------------------------------------------------------
!
!*    2.25   water supply
!            ------------
!
ALLOCATE(xdata_watsup(jpcover,nvegtype))
!
xdata_watsup(:,:) = 0.                                
!
!-------------------------------------------------------------------------------
!
!*    2.26   For multi-energy balance (MEB)
!            ------------------------------                     
!
ALLOCATE(xdata_gndlitter(jpcover,36,nvegtype))
xdata_gndlitter(:,:,:) = xundef
!
ALLOCATE(xdata_z0litter(jpcover,36,nvegtype))
xdata_z0litter(:,:,:) = xundef
!
ALLOCATE(xdata_h_veg(jpcover,36,nvegtype))
xdata_h_veg(:,:,:) = xundef
!
!-------------------------------------------------------------------------------
!
nut_cphr = sum(ntype(1:3)) + 1
nut_cpmr = sum(ntype(1:3)) + 2
nut_cplr = sum(ntype(1:3)) + 3
nut_ophr = sum(ntype(1:3)) + 4
nut_opmr = sum(ntype(1:3)) + 5
nut_oplr = sum(ntype(1:3)) + 6
nut_lwlr = sum(ntype(1:3)) + 7
nut_lalr = sum(ntype(1:3)) + 8
nut_spar = sum(ntype(1:3)) + 9
nut_indu = sum(ntype(1:3)) + 10
!
!*    3.1    z0 for artificial surfaces
!            --------------------------
!
ALLOCATE(xdata_z0_town(jpcover))
!
xdata_z0_town(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.2    albedo for artificial surfaces
!            ------------------------------
!
ALLOCATE(xdata_alb_roof(jpcover))
!
xdata_alb_roof(:) = xundef
!
ALLOCATE(xdata_alb_road(jpcover))
!
xdata_alb_road(:) = xundef
!
ALLOCATE(xdata_alb_wall(jpcover))
!
xdata_alb_wall(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.4    emissivity for artificial surfaces
!            ----------------------------------
!
ALLOCATE(xdata_emis_roof(jpcover))
!
xdata_emis_roof(:) = xundef
!
ALLOCATE(xdata_emis_road(jpcover))
!
xdata_emis_road(:) = xundef
!
ALLOCATE(xdata_emis_wall(jpcover))
!
xdata_emis_wall(:) = xundef
!
!-------------------------------------------------------------------------------
ndata_roof_layer=3
ndata_road_layer=3
ndata_wall_layer=3
ndata_floor_layer=3
!-------------------------------------------------------------------------------
!
!*    3.5    heat capacity for artificial surfaces
!            -------------------------------------
!
ALLOCATE(xdata_hc_roof(jpcover,ndata_roof_layer))
!
xdata_hc_roof(:,:) = xundef
!
ALLOCATE(xdata_hc_road(jpcover,ndata_road_layer))
!
xdata_hc_road(:,:) = xundef
!
ALLOCATE(xdata_hc_wall(jpcover,ndata_wall_layer))
!
xdata_hc_wall(:,:) = xundef
!
ALLOCATE(xdata_hc_floor(jpcover,ndata_floor_layer))
!
xdata_hc_floor(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.6    thermal conductivity for artificial surfaces
!            --------------------------------------------
!
ALLOCATE(xdata_tc_roof(jpcover,ndata_roof_layer))
!
xdata_tc_roof(:,:) = xundef
!
ALLOCATE(xdata_tc_road(jpcover,ndata_road_layer))
!
xdata_tc_road(:,:) = xundef
!
ALLOCATE(xdata_tc_wall(jpcover,ndata_wall_layer))
!
xdata_tc_wall(:,:) = xundef
!
ALLOCATE(xdata_tc_floor(jpcover,ndata_floor_layer))
!
xdata_tc_floor(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.7    depth for artificial surfaces layers
!            ------------------------------------
!
ALLOCATE(xdata_d_roof(jpcover,ndata_roof_layer))
!
xdata_d_roof(:,:) = xundef
!
ALLOCATE(xdata_d_road(jpcover,ndata_road_layer))
!
xdata_d_road(:,:) = xundef
!
ALLOCATE(xdata_d_wall(jpcover,ndata_wall_layer))
!
xdata_d_wall(:,:) = xundef
!
ALLOCATE(xdata_d_floor(jpcover,ndata_floor_layer))
!
xdata_d_floor(:,:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.8    building height
!            ---------------
!
ALLOCATE(xdata_bld_height(jpcover))
!
xdata_bld_height(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.9    building shape
!            --------------
!
ALLOCATE(xdata_wall_o_hor(jpcover))
!
xdata_wall_o_hor(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.10   building fraction
!            -----------------
!
ALLOCATE(xdata_bld(jpcover))
!
xdata_bld(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.11   canyon shape
!            ------------
!
ALLOCATE(xdata_can_hw_ratio(jpcover))
!
xdata_can_hw_ratio(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.12   anthropogenic fluxes
!            --------------------
!
ALLOCATE(xdata_h_traffic(jpcover))
ALLOCATE(xdata_le_traffic(jpcover))
ALLOCATE(xdata_h_industry(jpcover))
ALLOCATE(xdata_le_industry(jpcover))
!
xdata_h_traffic(:) = xundef
xdata_le_traffic(:) = xundef
xdata_h_industry(:) = xundef
xdata_le_industry(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.13   For TEB-BEM
!            ------------
!
ALLOCATE(xdata_tcool_target(jpcover))
ALLOCATE(xdata_theat_target(jpcover))
!
xdata_tcool_target(:) = xundef
xdata_theat_target(:) = xundef
!
ALLOCATE(xdata_f_waste_can(jpcover))
ALLOCATE(xdata_eff_heat(jpcover))
ALLOCATE(xdata_qin(jpcover))
ALLOCATE(xdata_qin_frad(jpcover))
ALLOCATE(xdata_shgc(jpcover))
ALLOCATE(xdata_u_win(jpcover))
ALLOCATE(xdata_gr(jpcover))
ALLOCATE(xdata_shgc_sh(jpcover))
ALLOCATE(xdata_floor_height(jpcover))
ALLOCATE(xdata_inf(jpcover))
!
xdata_f_waste_can(:) = xundef
xdata_eff_heat(:) = xundef
xdata_qin(:) = xundef
xdata_qin_frad(:) = xundef
xdata_shgc(:) = xundef
xdata_u_win(:) = xundef
xdata_gr(:) = xundef
xdata_shgc_sh(:) = xundef
xdata_floor_height(:) = xundef
xdata_inf(:) = xundef
!
ALLOCATE(xdata_f_water_cond(jpcover))
ALLOCATE(xdata_qin_flat(jpcover))
ALLOCATE(xdata_hr_target(jpcover))
ALLOCATE(xdata_v_vent(jpcover))
ALLOCATE(xdata_cap_sys_heat(jpcover))
ALLOCATE(xdata_cap_sys_rat(jpcover))
ALLOCATE(xdata_t_adp(jpcover))
ALLOCATE(xdata_m_sys_rat(jpcover))
ALLOCATE(xdata_cop_rat(jpcover))
ALLOCATE(xdata_t_size_max(jpcover))
ALLOCATE(xdata_t_size_min(jpcover))
ALLOCATE(xdata_shade(jpcover))
ALLOCATE(xdata_natvent(jpcover))
!
xdata_f_water_cond(:) = xundef
xdata_qin_flat(:) = xundef
xdata_hr_target(:) = xundef
xdata_v_vent(:) = xundef
xdata_cap_sys_heat(:) = xundef
xdata_cap_sys_rat(:) = xundef
xdata_t_adp(:) = xundef
xdata_m_sys_rat(:) = xundef
xdata_cop_rat(:) = xundef
xdata_t_size_max(:) = xundef
xdata_t_size_min(:) = xundef
xdata_shade(:) = 0.
xdata_natvent(:) = 0.
!
ALLOCATE(xdata_rough_roof(jpcover))
ALLOCATE(xdata_rough_wall(jpcover))
xdata_rough_roof(:) = xundef 
xdata_rough_wall(:) = xundef
!
ALLOCATE(xdata_residential(jpcover))
xdata_residential(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*    3.13   For greenroof
!            -------------
!
ALLOCATE(xdata_frac_gr(jpcover))
!
xdata_frac_gr(:) = 0.
!
!-------------------------------------------------------------------------------
!
!*    3.14   For solar panels
!            ----------------
!
ALLOCATE(xdata_emis_panel(jpcover))
ALLOCATE(xdata_alb_panel(jpcover))
ALLOCATE(xdata_eff_panel(jpcover))
ALLOCATE(xdata_frac_panel(jpcover))
!
xdata_emis_panel(:) = xundef
xdata_alb_panel(:) = xundef
xdata_eff_panel(:) = xundef
xdata_frac_panel(:) = xundef
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
!*    4.     CALL to INITIALIZATION ROUTINES
!            -------------------------------
!
!-------------------------------------------------------------------------------
!
!*   4.1   first version of ecoclimap (global)
!          ----------------------------------- 
!
!* Global data set
!
IF (u%LECOSG) THEN
  CALL init_types_param
ELSEIF (lread_data_cover) THEN
  CALL read_covers_param(1)
ELSE
  CALL default_data_cover(xdata_town,xdata_nature,xdata_water,xdata_sea,        &
         xdata_z0_town,xdata_bld_height,xdata_wall_o_hor,xdata_bld,xdata_garden,&
         xdata_alb_roof,xdata_alb_road,xdata_alb_wall,xdata_emis_roof,          &
         xdata_emis_road,xdata_emis_wall,xdata_hc_roof,xdata_tc_roof,           &
         xdata_d_roof,xdata_hc_road,xdata_tc_road,xdata_d_road,                 &
         xdata_hc_wall,xdata_tc_wall,xdata_d_wall,xdata_h_traffic,              &
         xdata_le_traffic,xdata_h_industry,xdata_le_industry,                   &
         xdata_vegtype,xdata_h_tree,xdata_watsup,xdata_irrig,                   &
         xdata_root_depth,xdata_ground_depth,xdata_dice,tdata_seed,             &
         tdata_reap)
  !
  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_alb_soil_eco1
  CALL default_alb_soil_eco2
  !
  CALL default_alb_veg_eco1_01
  CALL default_alb_veg_eco1_02
  CALL default_alb_veg_eco1_03
  CALL default_alb_veg_eco1_04
  CALL default_alb_veg_eco1_05
  CALL default_alb_veg_eco1_06
  CALL default_alb_veg_eco1_07
  CALL default_alb_veg_eco1_08
  CALL default_alb_veg_eco1_09
  CALL default_alb_veg_eco1_10
  CALL default_alb_veg_eco1_11
  CALL default_alb_veg_eco1_12
  CALL default_alb_veg_eco1_13
  CALL default_alb_veg_eco1_14
  CALL default_alb_veg_eco1_15
  CALL default_alb_veg_eco1_16
  CALL default_alb_veg_eco1_17
  CALL default_alb_veg_eco1_18
  CALL default_alb_veg_eco1_19
ENDIF
!
IF (.NOT.u%LECOSG) CALL cover301_573
!
!-------------------------------------------------------------------------------
!
!default values for ECOSG
IF (u%LECOSG) THEN

  !lai
  DO jdec = 1,SIZE(xdata_lai,2)
    WHERE(xdata_vegtype(:,1)>0.) xdata_lai(:,jdec,1) = 0. 
    WHERE(xdata_vegtype(:,2)>0.) xdata_lai(:,jdec,2) = 0.
    WHERE(xdata_vegtype(:,3)>0.) xdata_lai(:,jdec,3) = 0.
    WHERE(xdata_vegtype(:,4)>0.) xdata_lai(:,jdec,4) = 4.
    WHERE(xdata_vegtype(:,5)>0.) xdata_lai(:,jdec,5) = 4.
    WHERE(xdata_vegtype(:,6)>0.) xdata_lai(:,jdec,6) = 4.
    WHERE(xdata_vegtype(:,7)>0.) xdata_lai(:,jdec,7) = 4.
    WHERE(xdata_vegtype(:,8)>0.) xdata_lai(:,jdec,8) = 4.
    WHERE(xdata_vegtype(:,9)>0.) xdata_lai(:,jdec,9) = 4.
    WHERE(xdata_vegtype(:,10)>0.) xdata_lai(:,jdec,10) = 4.
    WHERE(xdata_vegtype(:,11)>0.) xdata_lai(:,jdec,11) = 4.
    WHERE(xdata_vegtype(:,12)>0.) xdata_lai(:,jdec,12) = 4.
    WHERE(xdata_vegtype(:,13)>0.) xdata_lai(:,jdec,13) = 4.
    WHERE(xdata_vegtype(:,14)>0.) xdata_lai(:,jdec,14) = 3.
    WHERE(xdata_vegtype(:,15)>0.) xdata_lai(:,jdec,15) = 3.
    WHERE(xdata_vegtype(:,16)>0.) xdata_lai(:,jdec,16) = 2.
    WHERE(xdata_vegtype(:,17)>0.) xdata_lai(:,jdec,17) = 2.
    WHERE(xdata_vegtype(:,18)>0.) xdata_lai(:,jdec,18) = 2.
    WHERE(xdata_vegtype(:,19)>0.) xdata_lai(:,jdec,19) = 4.
    WHERE(xdata_vegtype(:,20)>0.) xdata_lai(:,jdec,20) = 3.
  ENDDO

  !root_depth
  WHERE(xdata_vegtype(:,1)>0.) xdata_root_depth(:,1) = 0.5
  WHERE(xdata_vegtype(:,2)>0.) xdata_root_depth(:,2) = 0.2
  WHERE(xdata_vegtype(:,3)>0.) xdata_root_depth(:,3) = 0.2
  WHERE(xdata_vegtype(:,4)>0.) xdata_root_depth(:,4) = 2.0
  WHERE(xdata_vegtype(:,5)>0.) xdata_root_depth(:,5) = 2.0
  WHERE(xdata_vegtype(:,6)>0.) xdata_root_depth(:,6) = 3.0
  WHERE(xdata_vegtype(:,7)>0.) xdata_root_depth(:,7) = 2.0
  WHERE(xdata_vegtype(:,8)>0.) xdata_root_depth(:,8) = 5.0
  WHERE(xdata_vegtype(:,9)>0.) xdata_root_depth(:,9) = 2.0
  WHERE(xdata_vegtype(:,10)>0.) xdata_root_depth(:,10) = 2.0
  WHERE(xdata_vegtype(:,11)>0.) xdata_root_depth(:,11) = 2.0
  WHERE(xdata_vegtype(:,12)>0.) xdata_root_depth(:,12) = 2.0
  WHERE(xdata_vegtype(:,13)>0.) xdata_root_depth(:,13) = 1.0
  WHERE(xdata_vegtype(:,14)>0.) xdata_root_depth(:,14) = 1.0
  WHERE(xdata_vegtype(:,15)>0.) xdata_root_depth(:,15) = 1.5
  WHERE(xdata_vegtype(:,16)>0.) xdata_root_depth(:,16) = 1.5
  WHERE(xdata_vegtype(:,17)>0.) xdata_root_depth(:,17) = 1.5
  WHERE(xdata_vegtype(:,18)>0.) xdata_root_depth(:,18) = 1.5
  WHERE(xdata_vegtype(:,19)>0.) xdata_root_depth(:,19) = 2.0
  WHERE(xdata_vegtype(:,20)>0.) xdata_root_depth(:,20) = 1.0
  
  !soil_depth
  WHERE(xdata_vegtype(:,1)>0.) xdata_ground_depth(:,1) = 1.0
  WHERE(xdata_vegtype(:,2)>0.) xdata_ground_depth(:,2) = 0.2
  WHERE(xdata_vegtype(:,3)>0.) xdata_ground_depth(:,3) = 0.2
  WHERE(xdata_vegtype(:,4)>0.) xdata_ground_depth(:,4) = 3.0
  WHERE(xdata_vegtype(:,5)>0.) xdata_ground_depth(:,5) = 3.0
  WHERE(xdata_vegtype(:,6)>0.) xdata_ground_depth(:,6) = 3.0
  WHERE(xdata_vegtype(:,7)>0.) xdata_ground_depth(:,7) = 3.0
  WHERE(xdata_vegtype(:,8)>0.) xdata_ground_depth(:,8) = 8.0
  WHERE(xdata_vegtype(:,9)>0.) xdata_ground_depth(:,9) = 3.0
  WHERE(xdata_vegtype(:,10)>0.) xdata_ground_depth(:,10) = 3.0
  WHERE(xdata_vegtype(:,11)>0.) xdata_ground_depth(:,11) = 3.0
  WHERE(xdata_vegtype(:,12)>0.) xdata_ground_depth(:,12) = 3.0
  WHERE(xdata_vegtype(:,13)>0.) xdata_ground_depth(:,13) = 1.5
  WHERE(xdata_vegtype(:,14)>0.) xdata_ground_depth(:,14) = 1.5
  WHERE(xdata_vegtype(:,15)>0.) xdata_ground_depth(:,15) = 2.0
  WHERE(xdata_vegtype(:,16)>0.) xdata_ground_depth(:,16) = 2.0
  WHERE(xdata_vegtype(:,17)>0.) xdata_ground_depth(:,17) = 2.0
  WHERE(xdata_vegtype(:,18)>0.) xdata_ground_depth(:,18) = 2.0
  WHERE(xdata_vegtype(:,19)>0.) xdata_ground_depth(:,19) = 2.0
  WHERE(xdata_vegtype(:,20)>0.) xdata_ground_depth(:,20) = 1.5
 
  !height_trees
  WHERE(xdata_vegtype(:,4)>0.) xdata_h_tree(:,4) = 15.0
  WHERE(xdata_vegtype(:,5)>0.) xdata_h_tree(:,5) = 10.0
  WHERE(xdata_vegtype(:,6)>0.) xdata_h_tree(:,6) = 30.0
  WHERE(xdata_vegtype(:,7)>0.) xdata_h_tree(:,7) = 10.0
  WHERE(xdata_vegtype(:,8)>0.) xdata_h_tree(:,8) = 30.0
  WHERE(xdata_vegtype(:,9)>0.) xdata_h_tree(:,9) = 15.0
  WHERE(xdata_vegtype(:,10)>0.) xdata_h_tree(:,10) = 10.0
  WHERE(xdata_vegtype(:,11)>0.) xdata_h_tree(:,11) = 15.0
  WHERE(xdata_vegtype(:,12)>0.) xdata_h_tree(:,12) = 1.0
    
  !ice_depth
  WHERE(xdata_vegtype(:,1)>0.) xdata_dice(:,1) = 0.5
  WHERE(xdata_vegtype(:,2)>0.) xdata_dice(:,2) = 0.2
  WHERE(xdata_vegtype(:,3)>0.) xdata_dice(:,3) = 0.2
  WHERE(xdata_vegtype(:,4)>0.) xdata_dice(:,4) = 1.0
  WHERE(xdata_vegtype(:,5)>0.) xdata_dice(:,5) = 1.5
  WHERE(xdata_vegtype(:,6)>0.) xdata_dice(:,6) = 3.0
  WHERE(xdata_vegtype(:,7)>0.) xdata_dice(:,7) = 2.0
  WHERE(xdata_vegtype(:,8)>0.) xdata_dice(:,8) = 8.0
  WHERE(xdata_vegtype(:,9)>0.) xdata_dice(:,9) = 1.0
  WHERE(xdata_vegtype(:,10)>0.) xdata_dice(:,10) = 1.5
  WHERE(xdata_vegtype(:,11)>0.) xdata_dice(:,11) = 0.3
  WHERE(xdata_vegtype(:,12)>0.) xdata_dice(:,12) = 1.5
  WHERE(xdata_vegtype(:,13)>0.) xdata_dice(:,13) = 0.3
  WHERE(xdata_vegtype(:,14)>0.) xdata_dice(:,14) = 1.0
  WHERE(xdata_vegtype(:,15)>0.) xdata_dice(:,15) = 1.5
  WHERE(xdata_vegtype(:,16)>0.) xdata_dice(:,16) = 1.0
  WHERE(xdata_vegtype(:,17)>0.) xdata_dice(:,17) = 1.2
  WHERE(xdata_vegtype(:,18)>0.) xdata_dice(:,18) = 1.5
  WHERE(xdata_vegtype(:,19)>0.) xdata_dice(:,19) = 1.5
  WHERE(xdata_vegtype(:,20)>0.) xdata_dice(:,20) = 1.0
  
  DO jcov = 1,jpcover
    DO jveg = 1,nvegtype
      IF (xdata_vegtype(jcov,jveg)/=0.) THEN
        !alb_soil_nir
        xdata_alb_soil_nir(jcov,:,jveg) = 0.3
        !alb_soil_vis
        xdata_alb_soil_vis(jcov,:,jveg) = 0.1
        !alb_veg_nir
        xdata_alb_veg_nir(jcov,:,jveg) = 0.3  
        !alb_veg_vis
        xdata_alb_veg_vis(jcov,:,jveg) = 0.1
      ENDIF
    ENDDO
  ENDDO

  WHERE(xdata_town(:)>0.)  
    xdata_alb_roof(:) = 0.15                                               
    xdata_alb_road(:) = 0.08                                                
    xdata_alb_wall(:) = 0.25                                               
    xdata_emis_roof(:) = 0.90                                               
    xdata_emis_road(:) = 0.94                                              
    xdata_emis_wall(:) = 0.85                           
    xdata_hc_roof(:,1) = 2.11 * 1.e6                             
    xdata_hc_roof(:,2) = 0.28 * 1.e6                                 
    xdata_hc_roof(:,3) = 0.29 * 1.e6                                         
    xdata_tc_roof(:,1) = 1.5100                                           
    xdata_tc_roof(:,2) = 0.0800                                          
    xdata_tc_roof(:,3) = 0.0500                                           
    xdata_d_roof(:,1) = 0.050                                         
    xdata_d_roof(:,2) = 0.400                                          
    xdata_d_roof(:,3) = 0.100                                      
    xdata_hc_road(:,1) = 1.94 * 1.e6                                    
    xdata_hc_road(:,2) = 1.28 * 1.e6                                  
    xdata_hc_road(:,3) = 1.28 * 1.e6                                 
    xdata_tc_road(:,1) = 0.7454                                 
    xdata_tc_road(:,2) = 0.2513                                
    xdata_tc_road(:,3) = 0.2513                              
    xdata_d_road(:,1) = 0.050                            
    xdata_d_road(:,2) = 0.100                            
    xdata_d_road(:,3) = 1.000                         
    xdata_hc_wall(:,1) = 1.55 * 1.e6                        
    xdata_hc_wall(:,2) = 1.55 * 1.e6                         
    xdata_hc_wall(:,3) = 0.29 * 1.e6                   
    xdata_tc_wall(:,1) = 0.9338  
    xdata_tc_wall(:,2) = 0.9338
    xdata_tc_wall(:,3) = 0.0500    
    xdata_d_wall(:,1) = 0.020   
    xdata_d_wall(:,2) = 0.125
    xdata_d_wall(:,3) = 0.050
    xdata_le_traffic(:) = 0.
    xdata_le_industry(:)= 0.
  END WHERE
  !
  WHERE(xdata_vegtype(:,8)>0.) 
    xdata_watsup(:,8) = 0.   
    xdata_irrig(:,8) = 0.00
  END WHERE
  !                                               
  WHERE(xdata_vegtype(:,9)>0.) 
    xdata_watsup(:,9) = 30.   
    xdata_irrig(:,9) = 1.00  
    tdata_seed(:,9)%TDATE%DAY   = 10
    tdata_seed(:,9)%TDATE%MONTH = 05
    tdata_reap(:,9)%TDATE%DAY   = 01                                                
    tdata_reap(:,9)%TDATE%MONTH = 08
  END WHERE
  !
ELSE

  !For one cover, the soil albedo from CM13 is the same for each vegtype
  DO jveg=2,nvegtype
    DO jcov = 1, jpcover
      xdata_alb_soil_nir(jcov,:,jveg) = xdata_alb_soil_nir(jcov,:,1)
      xdata_alb_soil_vis(jcov,:,jveg) = xdata_alb_soil_vis(jcov,:,1)
    ENDDO
  ENDDO

ENDIF
!
!-------------------------------------------------------------------------------
!
icpt_sea = 0
icpt_water = 0
!
IF (u%LECOSG) THEN
  ALLOCATE(nsea(1))
  ALLOCATE(nwater(2))
ELSE
  ALLOCATE(nsea(3))
  ALLOCATE(nwater(2))        
ENDIF
!
nsea(:) = 0
nwater(:) = 0
!
DO jcov = 1, jpcover
  !
  IF (xdata_nature(jcov)==1.) THEN
    IF (xdata_vegtype(jcov,nvt_no  )==1.) nbare_soil = jcov
    IF (xdata_vegtype(jcov,nvt_rock)==1.) nrock = jcov
    IF (xdata_vegtype(jcov,nvt_snow)==1.) npermsnow = jcov
  END IF
  !
  IF (xdata_sea(jcov)==1.) THEN
    icpt_sea = icpt_sea + 1
    IF(icpt_sea>SIZE(nsea))THEN
      CALL abor1_sfx('INI_DATA_COVER: problem with ecoclimap param : ICPT_SEA > SIZE(NSEA) ')
    ENDIF    
    nsea(icpt_sea) = jcov
  ENDIF
  !
  IF (xdata_water(jcov)==1.) THEN
    icpt_water = icpt_water + 1
    IF(icpt_water>SIZE(nwater))THEN
      CALL abor1_sfx('INI_DATA_COVER: problem with ecoclimap param : ICPT_WATER > SIZE(NWATER) ')
    ENDIF
    nwater(icpt_water) = jcov
  ENDIF
  !
  IF (xdata_town(jcov)==0.) cycle
  !
  xdata_can_hw_ratio(jcov) = 0.5 * xdata_wall_o_hor(jcov) / (1.-xdata_bld(jcov))
  !
  !* Building Energy Model variables 
  !
  xdata_hc_floor(jcov,:) = 2016000.
  xdata_tc_floor(jcov,:) = 1.95
  xdata_d_floor(jcov,1) = 0.01
  xdata_d_floor(jcov,2) = 0.04
  xdata_d_floor(jcov,3) = 0.10
  !
  xdata_tcool_target(jcov) = 297.16
  xdata_theat_target(jcov) = 292.16
  xdata_f_waste_can(jcov)  = 1.0
  xdata_eff_heat(jcov)     = 0.9
  xdata_qin(jcov)          = 5.8
  xdata_qin_frad(jcov)     = 0.2
  xdata_qin_flat(jcov)     = 0.2
  xdata_shgc(jcov)         = 0.763
  xdata_u_win(jcov)        = 2.716
  xdata_gr(jcov)           = 0.3
  xdata_shgc_sh(jcov)      = 0.763
  xdata_floor_height(jcov) = 3.0
  xdata_inf(jcov)          = 0.5
  xdata_f_water_cond(jcov) = 0.
  xdata_qin_flat(jcov)     = 0.2
  xdata_hr_target(jcov)    = 0.5
  xdata_v_vent(jcov)       = 0.0
  xdata_cap_sys_heat(jcov) = 100.
  xdata_cap_sys_rat(jcov)  = 90.
  xdata_t_adp(jcov)        = 285.66
  xdata_m_sys_rat(jcov)    = 0.0067
  xdata_cop_rat(jcov)      = 2.5
  xdata_t_size_max(jcov)   = 301.95
  xdata_t_size_min(jcov)   = 268.96
  xdata_shade(jcov)        = 0.0
  xdata_natvent(jcov)      = 0.0  
  xdata_rough_roof(jcov)   = 1.52  
  xdata_rough_wall(jcov)   = 1.52  
  xdata_residential(jcov)  = 1.
  !
  xdata_emis_panel(jcov) = 0.9
  xdata_alb_panel(jcov) = 0.1
  xdata_eff_panel(jcov) = 0.14
  xdata_frac_panel(jcov) = 0.
  !  
  IF (xdata_garden(jcov)/=0.) THEN
    DO jveg=1,nvegtype
      IF (xdata_vegtype(jcov,jveg)/=0.) THEN
        xdata_alb_soil_nir(jcov,:,jveg) = 0.3
        xdata_alb_soil_vis(jcov,:,jveg) = 0.1
        xdata_alb_veg_nir(jcov,:,jveg) = 0.3
        xdata_alb_veg_vis(jcov,:,jveg) = 0.1    
      ENDIF
    ENDDO
  ENDIF
  !
END DO
!
IF(icpt_sea<SIZE(nsea))THEN
  CALL abor1_sfx('INI_DATA_COVER: problem with ecoclimap param : ICPT_SEA < SIZE(NSEA) ')
ENDIF    
IF(icpt_water<SIZE(nwater))THEN
  CALL abor1_sfx('INI_DATA_COVER: problem with ecoclimap param : ICPT_WATER < SIZE(NWATER) ')
ENDIF
!
!-------------------------------------------------------------------------------
ALLOCATE(cnames(jpcover,2))
cnames(:,:) = ' '
!-------------------------------------------------------------------------------
!
IF (u%LECOSG) THEN

  cnames(1,1) = 'No vegetation (smooth)'
  cnames(2,1) = 'No vegetation (rocks) '
  cnames(3,1) = 'Permanent snow and ice'
  cnames(4,1) = 'temperate broadleaf cold-deciduous summergreen'
  cnames(5,1) = 'boreal needleleaf evergreen'
  cnames(6,1) = 'tropical broadleaf evergreen'
  cnames(7,1) = 'C3 cultures types'
  cnames(8,1) = 'C4 cultures types'
  cnames(9,1) = 'irrigated crops'
  cnames(10,1) = 'temperate grassland'
  cnames(11,1) = 'tropical grassland'
  cnames(12,1) = 'peat bogs, parks and gardens (irrigated grass)'
  cnames(13,1) = 'tropical broadleaf deciduous'
  cnames(14,1) = 'temperate broadleaf evergreen'
  cnames(15,1) = 'temperate needleleaf evergreen'
  cnames(16,1) = 'boreal broadleaf cold-deciduous summergreen'
  cnames(17,1) = 'boreal needleleaf cold-deciduous summergreen'
  cnames(18,1) = 'boreal grass'
  cnames(19,1) = 'shrub'
!
  cnames(20,1) = 'urban areas'
  cnames(21,1) = 'sea and oceans'
  cnames(22,1) = 'lakes'
  cnames(23,1) = 'rivers'
!
!-------------------------------------------------------------------------------
!
  cnames(1,2) = 'Sol nu'
  cnames(2,2) = 'Roche nue'
  cnames(3,2) = 'Neige permanente et glaciers'
  cnames(4,2) = 'Feuillus decidus en zone temperee'
  cnames(5,2) = 'Coniferes sempervirents en zone boreale'
  cnames(6,2) = 'Feuillus sempervirentsn en zone tropicale'
  cnames(7,2) = 'Cultures C3'
  cnames(8,2) = 'Cultures C4'
  cnames(9,2) = 'Cultures irriguées'
  cnames(10,2) = 'Prairies en zone temperee'
  cnames(11,2) = 'Prairies en zone tropicale'
  cnames(12,2) = 'Zone humides, parcs et jardins'
  cnames(13,2) = 'Feuillus decidus en zone tropicale'
  cnames(14,2) = 'Feuillus sempervirents en zone temperee'
  cnames(15,2) = 'Coniferes sempervirents en zone temperee'
  cnames(16,2) = 'Feuillus boreal broadleaf cold-deciduous summergreen'
  cnames(17,2) = 'Coniferes decidus en zone boreale'
  cnames(18,2) = 'Prairies en zone boreale'
  cnames(19,2) = 'Buissons'
!
  cnames(20,2) = 'Zones urbaines'
  cnames(21,2) = 'Mer et oceans'
  cnames(22,2) = 'Lacs'
  cnames(23,2) = 'Rivieres'        

ELSE

  cnames(1,1) = 'Sea and ocean'
  cnames(2,1) = 'Lakes'
  cnames(3,1) = 'Rivers'
  cnames(4,1) = 'Bare land'
  cnames(5,1) = 'Rocks'
  cnames(6,1) = 'Permanent snow and ice'
  cnames(7,1) = 'Urban and built-up'

  cnames(8,1) = 'Tropical undefined islands'
  cnames(9,1) = 'Subpolar undefined islands'
!
  cnames(10,1) = 'S-America cool ENF'
  cnames(11,1) = 'Boreal ENF'
  cnames(12,1) = 'Asia subtropical ENF' 
  cnames(13,1) = 'American Continental ENF' 
  cnames(14,1) = 'American Subtropical ENF'
  cnames(15,1) = 'American Cool Marine ENF'

  cnames(16,1) = 'Africa Equatorial EBF'
  cnames(17,1) = 'Africa Tr. wind EBF'
  cnames(18,1) = 'Oceanian Equatorial EBF'
  cnames(19,1) = 'Asia tropical EBF'
  cnames(20,1) = 'Oceania tropical EBF'
  cnames(21,1) = 'Amazonian EBF'
  cnames(22,1) = 'SH subtropical EBF'
  cnames(23,1) = 'Cent. America Tr. wind EBF'

  cnames(24,1) = 'Asian boreal DNF'

  cnames(25,1) = 'S-America tropical DBF'
  cnames(26,1) = 'N-America humid continental DBF'
  cnames(27,1) = 'Cent. America Tr. wind DBF'
  cnames(28,1) = 'S-America humid subtropical DBF'

  cnames(29,1) = 'Africa dry tropical MF'
  cnames(30,1) = 'S-America cool MF'
  cnames(31,1) = 'NH Subpolar MF'
  cnames(32,1) = 'NH Humid subtropical MF'
  cnames(33,1) = 'NH Continental MF'

  cnames(34,1) = 'NH Africa WL'
  cnames(35,1) = 'SH Africa WL'
  cnames(36,1) = 'Tr. wind humid and subtrop. WL'
  cnames(37,1) = 'Oceanian Equatorial WL'
  cnames(38,1) = 'Asia wet tropical WL'
  cnames(39,1) = 'S-America tropical WL'
  cnames(40,1) = 'S-America humid subtropical WL'
  cnames(41,1) = 'NH Subpolar WL'
  cnames(42,1) = 'NH Continental WL'
  cnames(43,1) = 'Asia humid subtropical WL'
  cnames(44,1) = 'N-America Semi arid WL'
  cnames(45,1) = 'N-America moderate polar WL'
  cnames(46,1) = 'S-America moderate polar WL'
  cnames(47,1) = 'N-America humid subtropical WL'

  cnames(48,1) = 'NH Africa semiarid WG'
  cnames(49,1) = 'NH Africa dry tropical WG'
  cnames(50,1) = 'Africa dry equatorial WG'
  cnames(51,1) = 'SH Africa dry tropical WG'
  cnames(52,1) = 'Oceania tropical WG'
  cnames(53,1) = 'Oceania semiarid WG'
  cnames(54,1) = 'Oceania subtrop. cool marine WG'
  cnames(55,1) = 'Asia humid and subtropical WG'
  cnames(56,1) = 'S-America trop. and subtrop. WG'
  cnames(57,1) = 'S-America Tr. wind WG'
  cnames(58,1) = 'S-America semiarid WG'
  cnames(59,1) = 'NH Subpolar WG'
  cnames(60,1) = 'NH Continental WG'
  cnames(61,1) = 'Asia wet and dry tropical WG'
  cnames(62,1) = 'N-America semi arid WG'
  cnames(63,1) = 'N-America humid subtropical WG'
  cnames(64,1) = 'S-America moderate polar WG'
  cnames(65,1) = 'Cent. Amer. Tr. wind \& trop. WG'
  cnames(66,1) = 'NH Africa dry summer subtrop. WG'

  cnames(67,1) = 'NH Africa arid CS'
  cnames(68,1) = 'NH Africa semiarid CS'
  cnames(69,1) = 'SH Africa semiarid CS'
  cnames(70,1) = 'Oceania arid CS'
  cnames(71,1) = 'Oceania, S-America semiarid CS'
  cnames(72,1) = 'Oceania Tr. wind CS'
  cnames(73,1) = 'SH dry summer subtropical CS'
  cnames(74,1) = 'Asia polar CS'
  cnames(75,1) = 'Asia continental CS'
  cnames(76,1) = 'Asia tropical CS'
  cnames(77,1) = 'N-America polar CS'
  cnames(78,1) = 'N-America continental CS'
  cnames(79,1) = 'NH Africa dry summer subtrop. CS'

  cnames(80,1) = 'NH arid OS'
  cnames(81,1) = 'NH semiarid tropical OS'
  cnames(82,1) = 'SH Africa and Oceania arid OS'
  cnames(83,1) = 'S-America semiarid tropical OS'
  cnames(84,1) = 'Asia dry tropical OS'
  cnames(85,1) = 'NH Polar OS'
  cnames(86,1) = 'N-America Subpolar OS'
  cnames(87,1) = 'N-America semiarid continental OS'

  cnames(88,1) = 'Africa wet Tropical G'
  cnames(89,1) = 'NH Africa Semiarid G'
  cnames(90,1) = 'SH Africa Semiarid G'
  cnames(91,1) = 'S-America, Oceania equatorial G'
  cnames(92,1) = 'S-America, Oceania Semiarid G'
  cnames(93,1) = 'Oceania cool littoral G'
  cnames(94,1) = 'Asia wet and dry tropical G'
  cnames(95,1) = 'NH S-America wet tropical G'
  cnames(96,1) = 'SH S-America wet tropical G'
  cnames(97,1) = 'S-America semiarid G'
  cnames(98,1) = 'S-America moderate polar G'
  cnames(99,1) = 'NH semiarid Continental G'
  cnames(100,1) = 'Asia Subpolar G'
  cnames(101,1) = 'Asia humid Continental G'
  cnames(102,1) = 'Asia semiarid tropical G'
  cnames(103,1) = 'N-America continental G'
  cnames(104,1) = 'Asia humid subtropical G'

  cnames(105,1) = 'NH Africa arid C'
  cnames(106,1) = 'NH Africa, Asia wet and dry trop. C'
  cnames(107,1) = 'SH Africa wet and dry tropical C'
  cnames(108,1) = 'SH Afr. Tr. wind \& semiarid trop. C'
  cnames(109,1) = 'Oceania dry summer subtropical C'
  cnames(110,1) = 'Cent. \& S-Amer., Oceania Tr. wind C'
  cnames(111,1) = 'S-America humid subtropical C'
  cnames(112,1) = 'SH S-America tropical C'
  cnames(113,1) = 'N-Amer., Asia semiarid continental C'
  cnames(114,1) = 'Asia humid continental C'
  cnames(115,1) = 'Asia humid subtropical C'
  cnames(116,1) = 'Asia subpolar C'
  cnames(117,1) = 'Asia semiarid tropical C'
  cnames(118,1) = 'N-America humid continental C'
  cnames(119,1) = 'N-America humid subtropical C'
  cnames(120,1) = 'NH dry summer subtropical C'
  cnames(121,1) = 'NH Africa dry summer subtropical C'
  cnames(122,1) = 'SH Africa dry summer subtropical C'

  cnames(123,1) = 'Bare soil with sparse polar vegetation'

  cnames(124,1) = 'Warm subtropical wetlands'
  cnames(125,1) = 'Subpolar wetlands'
!
  cnames(151,1) = 'Dense urban'
  cnames(152,1) = 'Mediterranean sub-urban'
  cnames(153,1) = 'Temperate sub-urban'
  cnames(154,1) = 'Cold sub-urban'
  cnames(155,1) = 'Industries and commercial areas'
  cnames(156,1) = 'Road and rail networks'
  cnames(157,1) = 'Port facilities'
  cnames(158,1) = 'Airport'
  cnames(159,1) = 'Mineral extraction, construction sites'
  cnames(160,1) = 'Urban parks'
  cnames(161,1) = 'Sport facilities'
!
  cnames(162,1) = 'Spanish crops'
  cnames(163,1) = 'Estremadura crops'
  cnames(164,1) = 'Mediterranean crops'
  cnames(165,1) = 'Atlantic coast crops'
  cnames(166,1) = 'Temperate crops'
  cnames(167,1) = 'Po plain crops'
  cnames(168,1) = 'Warm temperate crops'
  cnames(169,1) = 'Ukrainian crops'
  cnames(170,1) = 'Subpolar crops'
  cnames(171,1) = 'Mountain crops'
  cnames(172,1) = 'Central Europe crops'
  cnames(173,1) = 'Turkish crops'
!
  cnames(174,1) = 'Mediterranean irrigated crops'
  cnames(175,1) = 'Irrigated crops'
  cnames(176,1) = 'Rice fields'
!
  cnames(177,1) = 'Mediterranean vineyards'
  cnames(178,1) = 'Temperate vineyards'
  cnames(179,1) = 'Mediterranean fruit trees'
  cnames(180,1) = 'Temperate fruit trees'
  cnames(181,1) = 'Olive groves'
!
  cnames(182,1) = 'Temperate pastures'
  cnames(183,1) = 'Atlantic border pastures'
  cnames(184,1) = 'Central and Eastern Europe pastures'
  cnames(185,1) = 'Ukrainian pastures'
  cnames(186,1) = 'Subpolar pastures'
!
  cnames(187,1) = 'Spanish complex cultivation pattern'
  cnames(188,1) = 'Mediter. complex cultivation pat.'
  cnames(189,1) = 'Temperate complex cultivation pat.'
  cnames(190,1) = 'French complex cultivation pat.'
  cnames(191,1) = 'Balkanish complex cultivation pat.'
!
  cnames(192,1) = 'Mediterranean crops and woodland'
  cnames(193,1) = 'Crops and woodland'
  cnames(194,1) = 'French crops and woodland'
  cnames(195,1) = 'Balkanish crops and woodland'
  cnames(196,1) = 'Spanish crops and woodland'
  cnames(197,1) = 'Baltic states crops and woodland'
!
  cnames(198,1) = 'Agro-forestry areas'
!
  cnames(199,1) = 'Spanish broad-leaved forest'
  cnames(200,1) = 'Estremadura broad-leaved forest'
  cnames(201,1) = 'Mediterranean broad-leaved forest'
  cnames(202,1) = 'Atlantic coast broad-leaved forest'
  cnames(203,1) = 'Temperate broad-leaved forest'
  cnames(204,1) = 'Moutain broad-leaved forest'
  cnames(205,1) = 'Balkanish broad-leaved forest'
  cnames(206,1) = 'Subpolar broad-leaved forest'
  cnames(207,1) = 'Black Sea broad-leaved forest'
!
  cnames(208,1) = 'Mediterranean pines'
  cnames(209,1) = 'Landes forest'
  cnames(210,1) = 'Moutain coniferous forest'
  cnames(211,1) = 'Temperate coniferous forest'
  cnames(212,1) = 'Subpolar Taiga'
  cnames(213,1) = 'Russian Taiga'
  cnames(214,1) = 'Turkish coniferous forest'
!
  cnames(215,1) = 'Mediterranean mixed forest'
  cnames(216,1) = 'Atlantic coast \& french mixed forest'
  cnames(217,1) = 'Subpolar mixed forest'
  cnames(218,1) = 'Mountain mixed forest'
  cnames(219,1) = 'Eastern Europe mixed forest'
!
  cnames(220,1) = 'Mediterranean GR'
  cnames(221,1) = 'Atlantic coast GR'
  cnames(222,1) = 'Balkanish GR'
  cnames(223,1) = 'Estremadura GR'
  cnames(224,1) = 'Subpolar GR'
  cnames(225,1) = 'Tundra'
!
  cnames(226,1) = 'Turkish moors'
  cnames(227,1) = 'Mediter. moors \& heath lands'
  cnames(228,1) = 'Moutain moors \& heath lands'
  cnames(229,1) = 'Atlantic coast moors \& heath lands'
!
  cnames(230,1) = 'Turkish shrubland'
  cnames(231,1) = 'Mediterranean maquis'
  cnames(232,1) = 'Moutain maquis'
!
  cnames(233,1) = 'Spanish woodland'
  cnames(234,1) = 'Mediterranean woodland'
  cnames(235,1) = 'Temperate woodland'
!
  cnames(236,1) = 'Sparsely vegetated areas'
  cnames(237,1) = 'Burnt areas'
  cnames(238,1) = 'Temperate wetlands'
  cnames(239,1) = 'Subpolar wetlands'
  cnames(240,1) = 'Peat bogs'
  cnames(241,1) = 'Salines and salt marshes'
!
  cnames(242,1) = 'Intertidal flats'
  cnames(243,1) = 'Coastal lagoons'
!
  cnames(301,1)='N SCANDINAVIA TUNDRA1'
  cnames(302,1)='OURAL BF1'
  cnames(303,1)='CARELIE BF1'
  cnames(304,1)='NORTH RUSSIAN TAIGA1'
  cnames(305,1)='NORTH RUSSIAN TAIGA2'
  cnames(306,1)='CARELIE BF2'
  cnames(307,1)='RUSSIAN TAIGA3'
  cnames(308,1)='RUSSIAN BF1'
  cnames(309,1)='RUSSIAN TAIGA4'
  cnames(310,1)='S SCANDINAVIA TAIGA1'
  cnames(311,1)='SOUTH FINLANDIA MF1'
  cnames(312,1)='SOUTH NORWAY MF1'
  cnames(313,1)='BALTIC BF1'
  cnames(314,1)='BALTIC MF1'
  cnames(315,1)='SOUTH SWEDEN CF1'
  cnames(316,1)='BALTIC MF2'
  cnames(317,1)='SOUTH SWEDEN CF2'
  cnames(318,1)='SOUTH SWEDEN CF3'
  cnames(319,1)='SOUTH SWEDEN MF1'
  cnames(320,1)='MOUNTAIN MF1'
  cnames(321,1)='MOUNTAIN BF1'
  cnames(322,1)='TEMPERATE BF1'
  cnames(323,1)='TEMPERATE COMPLEX1'
  cnames(324,1)='MOUNTAIN CF1'
  cnames(325,1)='TEMP HERBACEOUS CF1'
  cnames(326,1)='ATLANTIC COAST BF1'
  cnames(327,1)='TURKISH CF1'
  cnames(328,1)='BALKAN CF1'
  cnames(329,1)='N SPAIN HERBAC MF1'
  cnames(330,1)='TEMP SW HERBAC CF1'
  cnames(331,1)='ATLANTIC COMPLEX1'
  cnames(332,1)='N SPAIN HERBAC MF2'
  cnames(333,1)='MEDITER COMPLEX1'
  cnames(334,1)='MEDITER COMPLEX2'
  cnames(335,1)='MEDITER COMPLEX3'
  cnames(336,1)='MEDITER COMPLEX4'
  cnames(337,1)='MEDITER COMPLEX5'
  cnames(338,1)='BURNT PORT HERBAC CF1'
  cnames(339,1)='BURNT PORT HERBAC BF1'
  cnames(340,1)='EGEE COAST COMPLEX1'
  cnames(341,1)='W MED COAST COMPLEX1'
  cnames(342,1)='MAGHR HERBACEOUS MF1'
  cnames(343,1)='ESTREM HERBACEOUS MF1'

!herbaceous / shrub covers
  cnames(344,1)='POLAR MOUNT TUNDRA1'
  cnames(345,1)='POLAR MOUNT TUNDRA2'
  cnames(346,1)='S SCANDINAVIA TUNDRA1'
  cnames(347,1)='NORTH TUNDRA1'
  cnames(348,1)='S SCANDINAVIA TUNDRA2'
  cnames(349,1)='NORTH RUSSIA TUNDRA1'
  cnames(350,1)='ARAL CONTINENTAL GR1'
  cnames(351,1)='MOUNTAIN TAIGA MOORS1'
  cnames(352,1)='SCOTTISH SWAMP MOORS1'
  cnames(353,1)='ATLANTIC COMPLEX2'
  cnames(354,1)='ATLANTIC GR1'
  cnames(355,1)='IR SCOT SWAMP MOORS1'
  cnames(356,1)='ASIAN SPARSE GR1'
  cnames(357,1)='AS SPARSE SW COMPLEX1'
  cnames(358,1)='N CASPIAN DES OS1'
  cnames(359,1)='ATLAS AS SPARSE COMP1'
  cnames(360,1)='SPARSE SCO CEN EU GR1'
  cnames(361,1)='TEMPERATE COMPLEX2'
  cnames(362,1)='ATLANTIC COMPLEX3'
  cnames(363,1)='ATLANTIC COMPLEX4'
  cnames(364,1)='N ATLANTIC PASTURES1'
  cnames(365,1)='SPARSE SCO CEN EU GR2'
  cnames(366,1)='SPARSE MOUNT E EU GR1'
  cnames(367,1)='TUR N CASP CONT GR1'
  cnames(368,1)='N CASPIAN CONT GR1'
  cnames(369,1)='IRA N CASP CONT GR1'
  cnames(370,1)='TUR IRA MOUNT CONT GR1'
  cnames(371,1)='E CASPIAN DES OS1'
  cnames(372,1)='N CASPIAN COMPLEX1'
  cnames(373,1)='IRAN MOUNT CONT GR1'
  cnames(374,1)='ASIAN SPARSE DES OS1'
  cnames(375,1)='E CASPIAN DES OS2'
  cnames(376,1)='N MEDITER COMPLEX1'
  cnames(377,1)='N MEDITER COMPLEX2'
  cnames(378,1)='ASIAN MEDIT CONT GR1'
  cnames(379,1)='SOUTH RUSSIA CONT GR1'
  cnames(380,1)='BLSEA SPARSE CONT GR1'
  cnames(381,1)='BLSEA SPARSE CONT GR2'
  cnames(382,1)='TURK MOUNT CONT GR1'
  cnames(383,1)='TURKISH COMPLEX1'
  cnames(384,1)='CAUCASIAN COMPLEX1'
  cnames(385,1)='N CASPIAN CONT GR2'
  cnames(386,1)='VOLGA VALLEY CONT GR1'
  cnames(387,1)='VOLGA VALLEY CONT GR2'
  cnames(388,1)='W CASPIAN CONT GR1'
  cnames(389,1)='CAUCASIAN COMPLEX2'
  cnames(390,1)='CAUCASIAN COMPLEX3'
  cnames(391,1)='BLSEA SPARSE CONT GR3'
  cnames(392,1)='CENT MASSIF COMPLEX1'
  cnames(393,1)='CENT MASSIF COMPLEX2'
  cnames(394,1)='TURK COAST COMPLEX1'
  cnames(395,1)='MESOPOTAMIA GR1'
  cnames(396,1)='TURK CILICIA COMPLEX1'
  cnames(397,1)='ASIAN COMPLEX1'
  cnames(398,1)='N MED SPARSE COMPLEX1'
  cnames(399,1)='MEDITER COMPLEX6'
  cnames(400,1)='MEDIT SPARSE COMPLEX1'
  cnames(401,1)='MEDIT SPARSE COMPLEX2'
  cnames(402,1)='MEDIT SPARSE COMPLEX3'
  cnames(403,1)='MEDIT SPARSE COMPLEX4'
  cnames(404,1)='N MED HERBACEOUS CF1'
  cnames(405,1)='ESTREMADURA GR1'
  cnames(406,1)='TUNISIA COMPLEX1'
  cnames(407,1)='TUNISIA HERBACEOUS1'
  cnames(408,1)='ALGERIA HERBACEOUS1'
  cnames(409,1)='DESERTIC HERBACEOUS1'
  cnames(410,1)='DESERTIC HERBACEOUS2'
  cnames(411,1)='SPAIN DES COMPLEX1'
  cnames(412,1)='MED SPARSE COMPLEX5'
  cnames(413,1)='MED SPARSE COMPLEX6'
  cnames(414,1)='MED SPARSE COMPLEX7'
  cnames(415,1)='ME SPARSE DES COMPL1'
  cnames(416,1)='NORTH ARABIA GR1'
  cnames(417,1)='N ARABIA DES COMPLEX1'
  cnames(418,1)='N ARABIA DESERTIC GR1'
  cnames(419,1)='MOROCCO HERBACEOUS1'
  cnames(420,1)='S MED COAST HERBAC1'
  cnames(421,1)='W MEDITER WOODLAND1'
  cnames(422,1)='S MED COAST HERBAC2'
  cnames(423,1)='MESOP DES HERBACEOUS1'
  cnames(424,1)='MAG COAST DES HERBAC1'
  cnames(425,1)='TU AR SPARSE HERBAC1'
  cnames(426,1)='MEDIT SPARSE COMPLEX8'
  cnames(427,1)='MED SPARSE HERBAC1'
  cnames(428,1)='MEDIT SPARSE COMPLEX9'
  cnames(429,1)='SPAIN SPARSE COMPLEX1'
  cnames(430,1)='N MED SPARSE COMPLEX2'
  cnames(431,1)='N MED SPARSE COMPLEX3'
  cnames(432,1)='MAGHRE DES HERBAC1'
  cnames(433,1)='MAGHRE DES HERBAC2'
  cnames(434,1)='MAGHRE DES HERBAC3'
  cnames(435,1)='N ARAB DES HERBAC1'
  cnames(436,1)='MESOPO DES HERBAC2'
  cnames(437,1)='TOURAN DES HERBAC1'
  cnames(438,1)='MESOPO DES HERBAC2'
  cnames(439,1)='TOURAN DES HERBAC2'
  cnames(440,1)='NEW ZEMBLE HERBAC1'
  cnames(441,1)='NEW ZEMBLE HERBAC2'

!crops
  cnames(442,1)='TRANS SIBERIAN CROPS1'
  cnames(443,1)='PO PLAIN CROPS1'
  cnames(444,1)='PO PLAIN CROPS2'
  cnames(445,1)='SPANISH FRENCH CROPS1'
  cnames(446,1)='SPANISH FR ITAL CROPS1'
  cnames(447,1)='DANUBE PLAIN CROPS1'
  cnames(448,1)='N MED SPARSE COMPLEX4'
  cnames(449,1)='BALKAN CROPS1'
  cnames(450,1)='SPAIN FR ITAL CROPS2'
  cnames(451,1)='ATLANTIC CROPS1'
  cnames(452,1)='FR MED SPARSE CROPS1'
  cnames(453,1)='FR MED SPARSE CROPS2'
  cnames(454,1)='ATL MED SPARSE CROPS1'
  cnames(455,1)='BENE BLACK SEA CROPS1'
  cnames(456,1)='FRENCH ITALIAN CROPS1'
  cnames(457,1)='FR MED SPARSE CROPS3'
  cnames(458,1)='MEDITER SPARSE CROPS1'
  cnames(459,1)='ATLANTIC CROPS2'
  cnames(460,1)='NORTH ATLANTIC CROPS1'
  cnames(461,1)='SOUTH RUSSIA CROPS1'
  cnames(462,1)='S RUSSIA BALTIC CROPS1'
  cnames(463,1)='UKRAINIAN CROPS1'
  cnames(464,1)='EAST CARPATES CROPS1'
  cnames(465,1)='E CENT EUROPE CROPS1'
  cnames(466,1)='W CENT EU SW CROPS1'
  cnames(467,1)='HUNGARIAN CROPS1'
  cnames(468,1)='N BLACK SEA CROPS1'
  cnames(469,1)='HUNG BULG CAUC CROPS1'
  cnames(470,1)='SOUTH SWEDEN CROPS1'
  cnames(471,1)='SW RUSSIA CROPS1'
  cnames(472,1)='SOUTH RUSSIA CROPS1'
  cnames(473,1)='IRAN N CASPIAN CROPS1'
  cnames(474,1)='FR TEMP SPARSE CROPS1'
  cnames(475,1)='BULGARIAN CROPS1'
  cnames(476,1)='BULGARIAN CROPS2'
  cnames(477,1)='SP TURK SPARSE CROPS1'
  cnames(478,1)='FRENCH CENT EU CROPS1'
  cnames(479,1)='N BLACK SEA CROPS2'
  cnames(480,1)='BULGARIAN CROPS3'
  cnames(481,1)='POLE CROPS1'
  cnames(482,1)='POLE CROPS2'
  cnames(483,1)='N BLACK SEA CROPS3'
  cnames(484,1)='CENT EU SPARSE CROPS1'
  cnames(485,1)='GERMAN CROPS1'
  cnames(486,1)='BEAUCE CROPS1'
  cnames(487,1)='DANE CROPS1'
  cnames(488,1)='DANE CROPS2'
  cnames(489,1)='NEU ATL SPARSE CROPS1'
  cnames(490,1)='SYRIAN CROPS1'
  cnames(491,1)='GERMAN CROPS2'
  cnames(492,1)='CHANNEL CROPS1'
  cnames(493,1)='CHANNEL CROPS2'
  cnames(494,1)='ITALIAN CROPS1'
  cnames(495,1)='TURKISH CROPS1'
  cnames(496,1)='N MEDIT SPARSE CROPS1'
  cnames(497,1)='SPAIN TUR ARAB CROPS1'
  cnames(498,1)='NORTH SPAIN CROPS1'
  cnames(499,1)='MOROCCO TUNIS CROPS1'
  cnames(500,1)='MOROCCO CROPS1'
  cnames(501,1)='MOROCCO CROPS2'
  cnames(502,1)='ALGERIAN CROPS1'
  cnames(503,1)='MOROCCO CROPS3'
  cnames(504,1)='WEST SPAIN CROPS1'
  cnames(505,1)='MOROCCO CROPS4'
  cnames(506,1)='NORTH MEDITER CROPS1'
  cnames(507,1)='SOUTH SPANISH CROPS1'
  cnames(508,1)='SICILIAN CROPS1'
  cnames(509,1)='MAGHREB SPARSE CROPS1'
  cnames(510,1)='N MEDIT SPARSE CROPS2'
  cnames(511,1)='N MEDIT SPARSE CROPS3'
  cnames(512,1)='SP IT WCOAST CROPS1'
  cnames(513,1)='ESTREMADURA CROPS1'
  cnames(514,1)='ESTREMADURA CROPS2'
  cnames(515,1)='SP IT WCOAST CROPS2'
  cnames(516,1)='ESTREMADURA CROPS3'
  cnames(517,1)='MEDIT ISLANDS CROPS1'
  cnames(518,1)='SPAIN W COAST CROPS1'
  cnames(519,1)='ESTREMADURA CROPS4'
  cnames(520,1)='MECOAST SPARSE CROPS1'
  cnames(521,1)='BRITTANY CROPS1'
  cnames(522,1)='SYRIAN CROPS2'

!irrigated crops
  cnames(523,1)='NIL VALLEY CROPS1'
  cnames(524,1)='NIL VALLEY CROPS2'
  cnames(525,1)='NIL VALLEY CROPS3'
  cnames(526,1)='NIL VALLEY CROPS4'
  cnames(527,1)='SPANISH IRR CROPS1'
  cnames(528,1)='NIL VALLEY CROPS5'
  cnames(529,1)='EGEE IRR CROPS1'
  cnames(530,1)='MEDITER IRR CROPS1'
  cnames(531,1)='S SPAIN IRR CROPS1'
  cnames(532,1)='NIL VALLEY CROPS6'

!bare land      
  cnames(533,1)='BARE ROCK1'
  cnames(534,1)='BARE ROCK2'
  cnames(535,1)='SANDY DESERT1'
  cnames(536,1)='BARE LAND1'
  cnames(537,1)='BARE LAND2'
  cnames(538,1)='BARE LAND3'
  cnames(539,1)='BARE LAND4'
  cnames(540,1)='BARE LAND5'
  cnames(541,1)='BARE LAND6'
  cnames(542,1)='BARE LAND7'
  cnames(543,1)='BARE LAND8'
  cnames(544,1)='BARE LAND9'
  cnames(545,1)='BARE LAND10'
  cnames(546,1)='BARE LAND11'
  cnames(547,1)='BARE LAND12'
  cnames(548,1)='PERMANENT SNOW1'

!Estuaries and swamp areas      
  cnames(549,1)='WADDEN SEA'
  cnames(550,1)='ESTUARY1'
  cnames(551,1)='ESTUARY2'
  cnames(552,1)='POLAR WETLANDS1'
  cnames(553,1)='ESTUARY3'
  cnames(554,1)='ESTUARY4'
  cnames(555,1)='ESTUARY5'
  cnames(556,1)='ESTUARY6'
  cnames(557,1)='POLAR WETLANDS2'
  cnames(558,1)='SUBPOLAR WETLANDS1'
  cnames(559,1)='SUBPOLAR WETLANDS2'
  cnames(560,1)='SUBPOLAR WETLANDS3'

!urban
  cnames(561,1)='TEMPERATE SUBURBAN1'
  cnames(562,1)='TEMPERATE SUBURBAN2'
  cnames(563,1)='TEMPERATE SUBURBAN3'
  cnames(564,1)='TEMPERATE SUBURBAN4'
  cnames(565,1)='TEMPERATE SUBURBAN5'
  cnames(566,1)='COLD SUBURBAN1'
  cnames(567,1)='WARM SUBURBAN1'
  cnames(568,1)='WARM SUBURBAN2'
  cnames(569,1)='TEMPERATE SUBURBAN6'
  cnames(570,1)='TEMPERATE SUBURBAN7'
  cnames(571,1)='WARM SUBURBAN3'

!added classes of permanent crops      
  cnames(572,1)='SPANISH VINEYARDS1'
  cnames(573,1)='LANGUEDOC VINEYARDS1'

!-------------------------------------------------------------------------------
!

  cnames(1,2) = "Mers et oc\'eans"
  cnames(2,2) = "Lacs"
  cnames(3,2) = "rivi\`eres"
  cnames(4,2) = "Sol nu"
  cnames(5,2) = "Rochers"
  cnames(6,2) = "Neiges perp\'etuelles"
  cnames(7,2) = "Ville"

  cnames(8,2) = "\^iles non-d\'efinies (tropicales)"
  cnames(9,2) = "\^iles non-d\'efinies (subpolaires)"
!
  cnames(10,2) = "For\^et persistante de Patagonie "
  cnames(11,2) = 'Taiga'
  cnames(12,2) = "For\^et de conif\`eres d'Extr\^eme-Orient "
  cnames(13,2) = "For\^et de conif\`eres continentale am\'ericaine"
  cnames(14,2) = "For\^et de conif\`eres subtropicale am\'ericaine"
  cnames(15,2) = '"Rain forest"'

  cnames(16,2) = "For\^et \'equatoriale africaine"
  cnames(17,2) = "For\^et de Madagascar (aliz\'es)"
  cnames(18,2) = "For\^et \'equatoriale d'Oc\'eanie"
  cnames(19,2) = "For\^et persistante de Mousson"
  cnames(20,2) = "For\^et tropicale humide d'Oc\'eanie"
  cnames(21,2) = "For\^et amazonienne"
  cnames(22,2) = "For\^et primaire de Tasmanie"
  cnames(23,2) = "For\^et pan-am\'ericaine"

  cnames(24,2) = "Taiga sib\'erienne caduque"

  cnames(25,2) = "For\^et sub-amazonienne"
  cnames(26,2) = "For\^et caduque am\'ericaine"
  cnames(27,2) = "For\^et du Yucatan"
  cnames(28,2) = "For\^et sub-tropicale br\'esilienne"

  cnames(29,2) = "For\^et m\'elang\'ee tropicale s\`eche d'Afrique"
  cnames(30,2) = "For\^et m\'elang\'ee caduque de Patagonie"
  cnames(31,2) = "For\^et m\'elang\'ee sub-polaire"
  cnames(32,2) = "For\^et m\'elang\'ee humide sub-tropicale"
  cnames(33,2) = "For\^et m\'elang\'ee continentale"

  cnames(34,2) = "Prairie bois\'ee sah\'elienne"
  cnames(35,2) = "Prairie bois\'ee d'Afrique m\'eridionale"
  cnames(36,2) = "Prairie bois\'ee \`a feuilles persistante"
  cnames(37,2) = "Prairie bois\'ee \'equatoriale d'Oc\'eanie"
  cnames(38,2) = "Prairie bois\'ee de Mousson"
  cnames(39,2) = "Cerradao"
  cnames(40,2) = "Cerradao d'Argentine"
  cnames(41,2) = "Prairie bois\'ee polaire"
  cnames(42,2) = "Prairie bois\'ee continentale"
  cnames(43,2) = "Prairie bois\'ee d'Extr\^eme-Orient"
  cnames(44,2) = "Prairie bois\'ee semi-aride"
  cnames(45,2) = "Prairie bois\'ee des Rocheuses"
  cnames(46,2) = "Prairie bois\'ee de Patagonie"
  cnames(47,2) = "Prairie bois\'ee des Etats du Sud"

  cnames(48,2) = "Savane arbor\'ee sub-saharienne"
  cnames(49,2) = "Savane arbor\'ee sah\'elienne"
  cnames(50,2) = "Prairie arbor\'ee \'equatoriale africaine"
  cnames(51,2) = "Prairie arbor\'ee subtropicale africaine"
  cnames(52,2) = "Prairie arbor\'ee nord-australienne"
  cnames(53,2) = "Savane arbor\'ee australienne aride"
  cnames(54,2) = "Prairie arbor\'ee n\'eo-z\'elandaise"
  cnames(55,2) = "Prairie arbor\'ee d'Extr\^eme-Orient"
  cnames(56,2) = "Cerrado"
  cnames(57,2) = "Cerrado c\^otier"
  cnames(58,2) = "Sertao"
  cnames(59,2) = "Prairie arbor\'ee polaire"
  cnames(60,2) = "Prairie arbor\'ee continentale"
  cnames(61,2) = "Prairie arbor\'ee de Mousson"
  cnames(62,2) = "Prairie arbor\'ee pan-am\'ericaine"
  cnames(63,2) = "Prairie arbor\'ee des Etats du sud"
  cnames(64,2) = "Prairie arbor\'ee de la cordill\`ere des Andes"
  cnames(65,2) = "Prairie arbor\'ee Amazonienne"
  cnames(66,2) = "Prairie arbor\'ee du Maghreb"

  cnames(67,2) = "Savane arbustive d'Ethiopie"
  cnames(68,2) = "Savane arbustive sah\'elienne"
  cnames(69,2) = '"Bush" africain'
  cnames(70,2) = '"Bush" australien'
  cnames(71,2) = "Savane arbustive semi-d\'esertique"
  cnames(72,2) = "Savane arbustive dense c\^oti\`ere d'Oc\'eanie"
  cnames(73,2) = '"Bush" dense sud-australien'
  cnames(74,2) = "Tundra arbustive dense"
  cnames(75,2) = "Steppe asiatique arbustive dense"
  cnames(76,2) = "Prairie arbustive de Mousson"
  cnames(77,2) = "Tundra dense am\'ericaine"
  cnames(78,2) = "Prairie arbustive d'altitude (Rocheuses)"
  cnames(79,2) = "Prairie arbustive du Maghreb"

  cnames(80,2) = "Savane \'eparse semi-d\'esertique"
  cnames(81,2) = "Savane Sah\'elienne"
  cnames(82,2) = "Bush \'epars"
  cnames(83,2) = "Pampa"
  cnames(84,2) = "Prairie \'eparse de Mousson"
  cnames(85,2) = "Tundra polaire"
  cnames(86,2) = "Tundra sub-polaire"
  cnames(87,2) = "Prairie d'altitude (Rocheuses)"

  cnames(88,2) = "Prairie tropicale africaine"
  cnames(89,2) = "Prairie africaine semi-aride"
  cnames(90,2) = "Prairie d'Afrique m\'eridionale"
  cnames(91,2) = "Prairie \'equatoriale c\^oti\`ere"
  cnames(92,2) = "Prairie semi-d\'esertique"
  cnames(93,2) = "Prairie n\'eo-z\'elandaise"
  cnames(94,2) = "Prairie de Mousson"
  cnames(95,2) = "Prairie v\'en\'ezu\'elienne"
  cnames(96,2) = "Prairie sub-amazonienne"
  cnames(97,2) = "Pampa et prairie"
  cnames(98,2) = "Prairie de Patagonie"
  cnames(99,2) = "Steppes"
  cnames(100,2) = "Prairie polaire"
  cnames(101,2) = "Prairie d'Extr\^eme-Orient"
  cnames(102,2) = "Prairie indienne"
  cnames(103,2) = "Prairie d'altitude (Rocheuses)"
  cnames(104,2) = "Prairie chinoise"

  cnames(105,2) = "Cultures d\'esertique et le long du Nil"
  cnames(106,2) = "Cultures sah\'elienne"
  cnames(107,2) = "Cultures tropicales"
  cnames(108,2) = "Cultures c\^oti\`eres africaines"
  cnames(109,2) = "Cultures sud-australiennes"
  cnames(110,2) = "Cultures pan-am\'ericaines"
  cnames(111,2) = "Cultures de la plaine du Parana"
  cnames(112,2) = "Zones cultiv\'ees sub-amazoniennes"
  cnames(113,2) = "Zones cultiv\'ees trans-siberiennes"
  cnames(114,2) = "Zones cultiv\'ees de Mandchourie"
  cnames(115,2) = "Zones cultiv\'ees de Chine"
  cnames(116,2) = "Zones cultiv\'ees sub-polaires"
  cnames(117,2) = "Zones cultiv\'ees sub-himalayennes"
  cnames(118,2) = "Cultures des Grandes Plaines"
  cnames(119,2) = "Cultures des Etats du Sud"
  cnames(120,2) = "Cultures californiennes"
  cnames(121,2) = "Zones cultiv\'ees du Maghreb"
  cnames(122,2) = "Zones cultiv\'ees d'Afrique du Sud"

  cnames(123,2) = "Tundra \'eparse"

  cnames(124,2) = "Marais tropicaux et sub-tropicaux"
  cnames(125,2) = "Marais polaires et sub-polaires"
!
  cnames(151,2) = "Urbain dense"
  cnames(152,2) = "Suburbain m\'editerran\'een"
  cnames(153,2) = "Suburbain temp\'er\'e"
  cnames(154,2) = "Suburbain froid"
  cnames(155,2) = "Zones industrielles et commerciales"
  cnames(156,2) = "R\'eseaux routiers et ferroviaires"
  cnames(157,2) = "Zones portuaires"
  cnames(158,2) = "A\'eroports"
  cnames(159,2) = "Extraction de mat\'eriaux, sites de construction"
  cnames(160,2) = "Parcs urbains"
  cnames(161,2) = "Equipements sportifs et de loisirs"
!
  cnames(162,2) = "Cultures espagnoles"
  cnames(163,2) = "Cultures d\'Estremadure"
  cnames(164,2) = "Cultures m\'editerran\'eennes"
  cnames(165,2) = "Cultures de la c\^ote atlantique"
  cnames(166,2) = "Cultures temp\'er\'ees"
  cnames(167,2) = "Cultures de la plaine du P\^o"
  cnames(168,2) = "Cultures temp\'er\'ees chaudes"
  cnames(169,2) = "Cultures ukrainiennes"
  cnames(170,2) = "Cultures subpolaires"
  cnames(171,2) = "Cultures montagneuses"
  cnames(172,2) = "Cultures d'Europe Centrale"
  cnames(173,2) = "Cultures turques"
!
  cnames(174,2) = "Cultures m\'editerran\'eennes irrigu\'ees"
  cnames(175,2) = "Cultures irrigu\'ees"
  cnames(176,2) = "Rizi\`eres"
!
  cnames(177,2) = "Vignobles m\'editerran\'eens"
  cnames(178,2) = "Vignobles temp\'er\'es"
  cnames(179,2) = "Arbres fruitiers m\'editerran\'eens"
  cnames(180,2) = "Arbres fruitiers temp\'er\'es"
  cnames(181,2) = "Oliveraies"
!
  cnames(182,2) = "P\^atures temp\'er\'ees"
  cnames(183,2) = "P\^atures c\^oti\`eres atlantiques"
  cnames(184,2) = "P\^atures d'Europe Centrale et orientale"
  cnames(185,2) = "P\^atures ukrainiennes"
  cnames(186,2) = "P\^atures subpolaires"
!
  cnames(187,2) = "Cultures complexes espagnoles"
  cnames(188,2) = "Cultures complexes m\'editerran\'eennes"
  cnames(189,2) = "Cultures complexes temp\'er\'ees"
  cnames(190,2) = "Cultures complexes françaises"
  cnames(191,2) = "Cultures complexes balkaniques"
!
  cnames(192,2) = "Cultures bois\'ees m\'editerran\'eennes"
  cnames(193,2) = "Cultures bois\'ees"
  cnames(194,2) = "Cultures bois\'ees françaises"
  cnames(195,2) = "Cultures bois\'ees m\'editerran\'eennes"
  cnames(196,2) = "Cultures bois\'ees espagnoles"
  cnames(197,2) = "Cultures bois\'ees des pays baltes"
!
  cnames(198,2) = "Zones agro-foresti\`eres"
!
  cnames(199,2) = "For\^et de feuillus espagnole"
  cnames(200,2) = "For\^et de feuillus d\'Estremadure"
  cnames(201,2) = "For\^et de feuillus m\'editerran\'eennes"
  cnames(202,2) = "For\^et de feuillus de la c\^ote atlantique"
  cnames(203,2) = "For\^et de feuillus temp\'er\'ee"
  cnames(204,2) = "For\^et de feuillus montagneuses"
  cnames(205,2) = "For\^et de feuillus balkanique"
  cnames(206,2) = "For\^et de feuillus balkanique"
  cnames(207,2) = "For\^et de feuillus de la mer Noire"
!
  cnames(208,2) = "Pins m\'editerran\'eens"
  cnames(209,2) = "For\^et des Landes"
  cnames(210,2) = "For\^et de conif\`eres montagneuse"
  cnames(211,2) = "For\^et de conif\`eres temp\'er\'ee"
  cnames(212,2) = "Taiga subpolaire"
  cnames(213,2) = "Taiga russe"
  cnames(214,2) = "For\^et de conif\`eres turque"
!
  cnames(215,2) = "For\^et mixte m\'editerran\'eenne"
  cnames(216,2) = "For\^et mixte française \& c\^oti\`ere atlantique"
  cnames(217,2) = "For\^et mixte subpolaire"
  cnames(218,2) = "For\^et mixte montagneuses"
  cnames(219,2) = "For\^et mixte d'Europe orientale"
!
  cnames(220,2) = "Prairie m\'editerran\'eenne"
  cnames(221,2) = "Prairie c\^oti\`ere atlantique"
  cnames(222,2) = "Prairie balkanique"
  cnames(223,2) = "Prairie d\'Estremadure"
  cnames(224,2) = "Prairie subpolaire"
  cnames(225,2) = "Tundra"
!
  cnames(226,2) = "Landes turques"
  cnames(227,2) = "Landes et broussailles m\'editerran\'eennes"
  cnames(228,2) = "Landes et broussailles montagneuses"
  cnames(229,2) = "Landes et broussailles c\^oti\`eres atlantiques"
!
  cnames(230,2) = "V\'eg\'etation arbustive turque"
  cnames(231,2) = "Maquis m\'editerran\'een"
  cnames(232,2) = "Maquis montagneux"
!
  cnames(233,2) = "Zones bois\'ees espagnoles"
  cnames(234,2) = "Zones bois\'ees m\'editerran\'eennes"
  cnames(235,2) = "Zones bois\'ees temp\'er\'ees"
!
  cnames(236,2) = "V\'eg\'etation clairsem\'ee"
  cnames(237,2) = "Zones br\^ul\'ees"
  cnames(238,2) = "Zones humides temp\'er\'ees"
  cnames(239,2) = "Zones humides subpolaires"
  cnames(240,2) = "Tourbi\`eres"
  cnames(241,2) = "Marais salants"

!
  cnames(242,2) = "Zones intertidales"
  cnames(243,2) = "Lagunes littorales"
!
  cnames(301,2)='N SCANDINAVIA TUNDRA1'
  cnames(302,2)='OURAL BF1'
  cnames(303,2)='CARELIE BF1'
  cnames(304,2)='NORTH RUSSIAN TAIGA1'
  cnames(305,2)='NORTH RUSSIAN TAIGA2'
  cnames(306,2)='CARELIE BF2'
  cnames(307,2)='RUSSIAN TAIGA3'
  cnames(308,2)='RUSSIAN BF1'
  cnames(309,2)='RUSSIAN TAIGA4'
  cnames(310,2)='S SCANDINAVIA TAIGA1'
  cnames(311,2)='SOUTH FINLANDIA MF1'
  cnames(312,2)='SOUTH NORWAY MF1'
  cnames(313,2)='BALTIC BF1'
  cnames(314,2)='BALTIC MF1'
  cnames(315,2)='SOUTH SWEDEN CF1'
  cnames(316,2)='BALTIC MF2'
  cnames(317,2)='SOUTH SWEDEN CF2'
  cnames(318,2)='SOUTH SWEDEN CF3'
  cnames(319,2)='SOUTH SWEDEN MF1'
  cnames(320,2)='MOUNTAIN MF1'
  cnames(321,2)='MOUNTAIN BF1'
  cnames(322,2)='TEMPERATE BF1'
  cnames(323,2)='TEMPERATE COMPLEX1'
  cnames(324,2)='MOUNTAIN CF1'
  cnames(325,2)='TEMP HERBACEOUS CF1'
  cnames(326,2)='ATLANTIC COAST BF1'
  cnames(327,2)='TURKISH CF1'
  cnames(328,2)='BALKAN CF1'
  cnames(329,2)='N SPAIN HERBAC MF1'
  cnames(330,2)='TEMP SW HERBAC CF1'
  cnames(331,2)='ATLANTIC COMPLEX1'
  cnames(332,2)='N SPAIN HERBAC MF2'
  cnames(333,2)='MEDITER COMPLEX1'
  cnames(334,2)='MEDITER COMPLEX2'
  cnames(335,2)='MEDITER COMPLEX3'
  cnames(336,2)='MEDITER COMPLEX4'
  cnames(337,2)='MEDITER COMPLEX5'
  cnames(338,2)='BURNT PORT HERBAC CF1'
  cnames(339,2)='BURNT PORT HERBAC BF1'
  cnames(340,2)='EGEE COAST COMPLEX1'
  cnames(341,2)='W MED COAST COMPLEX1'
  cnames(342,2)='MAGHR HERBACEOUS MF1'
  cnames(343,2)='ESTREM HERBACEOUS MF1'

!herbaceous / shrub covers
  cnames(344,2)='POLAR MOUNT TUNDRA1'
  cnames(345,2)='POLAR MOUNT TUNDRA2'
  cnames(346,2)='S SCANDINAVIA TUNDRA1'
  cnames(347,2)='NORTH TUNDRA1'
  cnames(348,2)='S SCANDINAVIA TUNDRA2'
  cnames(349,2)='NORTH RUSSIA TUNDRA1'
  cnames(350,2)='ARAL CONTINENTAL GR1'
  cnames(351,2)='MOUNTAIN TAIGA MOORS1'
  cnames(352,2)='SCOTTISH SWAMP MOORS1'
  cnames(353,2)='ATLANTIC COMPLEX2'
  cnames(354,2)='ATLANTIC GR1'
  cnames(355,2)='IR SCOT SWAMP MOORS1'
  cnames(356,2)='ASIAN SPARSE GR1'
  cnames(357,2)='AS SPARSE SW COMPLEX1'
  cnames(358,2)='N CASPIAN DES OS1'
  cnames(359,2)='ATLAS AS SPARSE COMP1'
  cnames(360,2)='SPARSE SCO CEN EU GR1'
  cnames(361,2)='TEMPERATE COMPLEX2'
  cnames(362,2)='ATLANTIC COMPLEX3'
  cnames(363,2)='ATLANTIC COMPLEX4'
  cnames(364,2)='N ATLANTIC PASTURES1'
  cnames(365,2)='SPARSE SCO CEN EU GR2'
  cnames(366,2)='SPARSE MOUNT E EU GR1'
  cnames(367,2)='TUR N CASP CONT GR1'
  cnames(368,2)='N CASPIAN CONT GR1'
  cnames(369,2)='IRA N CASP CONT GR1'
  cnames(370,2)='TUR IRA MOUNT CONT GR1'
  cnames(371,2)='E CASPIAN DES OS1'
  cnames(372,2)='N CASPIAN COMPLEX1'
  cnames(373,2)='IRAN MOUNT CONT GR1'
  cnames(374,2)='ASIAN SPARSE DES OS1'
  cnames(375,2)='E CASPIAN DES OS2'
  cnames(376,2)='N MEDITER COMPLEX1'
  cnames(377,2)='N MEDITER COMPLEX2'
  cnames(378,2)='ASIAN MEDIT CONT GR1'
  cnames(379,2)='SOUTH RUSSIA CONT GR1'
  cnames(380,2)='BLSEA SPARSE CONT GR1'
  cnames(381,2)='BLSEA SPARSE CONT GR2'
  cnames(382,2)='TURK MOUNT CONT GR1'
  cnames(383,2)='TURKISH COMPLEX1'
  cnames(384,2)='CAUCASIAN COMPLEX1'
  cnames(385,2)='N CASPIAN CONT GR2'
  cnames(386,2)='VOLGA VALLEY CONT GR1'
  cnames(387,2)='VOLGA VALLEY CONT GR2'
  cnames(388,2)='W CASPIAN CONT GR1'
  cnames(389,2)='CAUCASIAN COMPLEX2'
  cnames(390,2)='CAUCASIAN COMPLEX3'
  cnames(391,2)='BLSEA SPARSE CONT GR3'
  cnames(392,2)='CENT MASSIF COMPLEX1'
  cnames(393,2)='CENT MASSIF COMPLEX2'
  cnames(394,2)='TURK COAST COMPLEX1'
  cnames(395,2)='MESOPOTAMIA GR1'
  cnames(396,2)='TURK CILICIA COMPLEX1'
  cnames(397,2)='ASIAN COMPLEX1'
  cnames(398,2)='N MED SPARSE COMPLEX1'
  cnames(399,2)='MEDITER COMPLEX6'
  cnames(400,2)='MEDIT SPARSE COMPLEX1'
  cnames(401,2)='MEDIT SPARSE COMPLEX2'
  cnames(402,2)='MEDIT SPARSE COMPLEX3'
  cnames(403,2)='MEDIT SPARSE COMPLEX4'
  cnames(404,2)='N MED HERBACEOUS CF1'
  cnames(405,2)='ESTREMADURA GR1'
  cnames(406,2)='TUNISIA COMPLEX1'
  cnames(407,2)='TUNISIA HERBACEOUS1'
  cnames(408,2)='ALGERIA HERBACEOUS1'
  cnames(409,2)='DESERTIC HERBACEOUS1'
  cnames(410,2)='DESERTIC HERBACEOUS2'
  cnames(411,2)='SPAIN DES COMPLEX1'
  cnames(412,2)='MED SPARSE COMPLEX5'
  cnames(413,2)='MED SPARSE COMPLEX6'
  cnames(414,2)='MED SPARSE COMPLEX7'
  cnames(415,2)='ME SPARSE DES COMPL1'
  cnames(416,2)='NORTH ARABIA GR1'
  cnames(417,2)='N ARABIA DES COMPLEX1'
  cnames(418,2)='N ARABIA DESERTIC GR1'
  cnames(419,2)='MOROCCO HERBACEOUS1'
  cnames(420,2)='S MED COAST HERBAC1'
  cnames(421,2)='W MEDITER WOODLAND1'
  cnames(422,2)='S MED COAST HERBAC2'
  cnames(423,2)='MESOP DES HERBACEOUS1'
  cnames(424,2)='MAG COAST DES HERBAC1'
  cnames(425,2)='TU AR SPARSE HERBAC1'
  cnames(426,2)='MEDIT SPARSE COMPLEX8'
  cnames(427,2)='MED SPARSE HERBAC1'
  cnames(428,2)='MEDIT SPARSE COMPLEX9'
  cnames(429,2)='SPAIN SPARSE COMPLEX1'
  cnames(430,2)='N MED SPARSE COMPLEX2'
  cnames(431,2)='N MED SPARSE COMPLEX3'
  cnames(432,2)='MAGHRE DES HERBAC1'
  cnames(433,2)='MAGHRE DES HERBAC2'
  cnames(434,2)='MAGHRE DES HERBAC3'
  cnames(435,2)='N ARAB DES HERBAC1'
  cnames(436,2)='MESOPO DES HERBAC2'
  cnames(437,2)='TOURAN DES HERBAC1'
  cnames(438,2)='MESOPO DES HERBAC2'
  cnames(439,2)='TOURAN DES HERBAC2'
  cnames(440,2)='NEW ZEMBLE HERBAC1'
  cnames(441,2)='NEW ZEMBLE HERBAC2'

!crops
  cnames(442,2)='TRANS SIBERIAN CROPS1'
  cnames(443,2)='PO PLAIN CROPS1'
  cnames(444,2)='PO PLAIN CROPS2'
  cnames(445,2)='SPANISH FRENCH CROPS1'
  cnames(446,2)='SPANISH FR ITAL CROPS1'
  cnames(447,2)='DANUBE PLAIN CROPS1'
  cnames(448,2)='N MED SPARSE COMPLEX4'
  cnames(449,2)='BALKAN CROPS1'
  cnames(450,2)='SPAIN FR ITAL CROPS2'
  cnames(451,2)='ATLANTIC CROPS1'
  cnames(452,2)='FR MED SPARSE CROPS1'
  cnames(453,2)='FR MED SPARSE CROPS2'
  cnames(454,2)='ATL MED SPARSE CROPS1'
  cnames(455,2)='BENE BLACK SEA CROPS1'
  cnames(456,2)='FRENCH ITALIAN CROPS1'
  cnames(457,2)='FR MED SPARSE CROPS3'
  cnames(458,2)='MEDITER SPARSE CROPS1'
  cnames(459,2)='ATLANTIC CROPS2'
  cnames(460,2)='NORTH ATLANTIC CROPS1'
  cnames(461,2)='SOUTH RUSSIA CROPS1'
  cnames(462,2)='S RUSSIA BALTIC CROPS1'
  cnames(463,2)='UKRAINIAN CROPS1'
  cnames(464,2)='EAST CARPATES CROPS1'
  cnames(465,2)='E CENT EUROPE CROPS1'
  cnames(466,2)='W CENT EU SW CROPS1'
  cnames(467,2)='HUNGARIAN CROPS1'
  cnames(468,2)='N BLACK SEA CROPS1'
  cnames(469,2)='HUNG BULG CAUC CROPS1'
  cnames(470,2)='SOUTH SWEDEN CROPS1'
  cnames(471,2)='SW RUSSIA CROPS1'
  cnames(472,2)='SOUTH RUSSIA CROPS1'
  cnames(473,2)='IRAN N CASPIAN CROPS1'
  cnames(474,2)='FR TEMP SPARSE CROPS1'
  cnames(475,2)='BULGARIAN CROPS1'
  cnames(476,2)='BULGARIAN CROPS2'
  cnames(477,2)='SP TURK SPARSE CROPS1'
  cnames(478,2)='FRENCH CENT EU CROPS1'
  cnames(479,2)='N BLACK SEA CROPS2'
  cnames(480,2)='BULGARIAN CROPS3'
  cnames(481,2)='POLE CROPS1'
  cnames(482,2)='POLE CROPS2'
  cnames(483,2)='N BLACK SEA CROPS3'
  cnames(484,2)='CENT EU SPARSE CROPS1'
  cnames(485,2)='GERMAN CROPS1'
  cnames(486,2)='BEAUCE CROPS1'
  cnames(487,2)='DANE CROPS1'
  cnames(488,2)='DANE CROPS2'
  cnames(489,2)='NEU ATL SPARSE CROPS1'
  cnames(490,2)='SYRIAN CROPS1'
  cnames(491,2)='GERMAN CROPS2'
  cnames(492,2)='CHANNEL CROPS1'
  cnames(493,2)='CHANNEL CROPS2'
  cnames(494,2)='ITALIAN CROPS1'
  cnames(495,2)='TURKISH CROPS1'
  cnames(496,2)='N MEDIT SPARSE CROPS1'
  cnames(497,2)='SPAIN TUR ARAB CROPS1'
  cnames(498,2)='NORTH SPAIN CROPS1'
  cnames(499,2)='MOROCCO TUNIS CROPS1'
  cnames(500,2)='MOROCCO CROPS1'
  cnames(501,2)='MOROCCO CROPS2'
  cnames(502,2)='ALGERIAN CROPS1'
  cnames(503,2)='MOROCCO CROPS3'
  cnames(504,2)='WEST SPAIN CROPS1'
  cnames(505,2)='MOROCCO CROPS4'
  cnames(506,2)='NORTH MEDITER CROPS1'
  cnames(507,2)='SOUTH SPANISH CROPS1'
  cnames(508,2)='SICILIAN CROPS1'
  cnames(509,2)='MAGHREB SPARSE CROPS1'
  cnames(510,2)='N MEDIT SPARSE CROPS2'
  cnames(511,2)='N MEDIT SPARSE CROPS3'
  cnames(512,2)='SP IT WCOAST CROPS1'
  cnames(513,2)='ESTREMADURA CROPS1'
  cnames(514,2)='ESTREMADURA CROPS2'
  cnames(515,2)='SP IT WCOAST CROPS2'
  cnames(516,2)='ESTREMADURA CROPS3'
  cnames(517,2)='MEDIT ISLANDS CROPS1'
  cnames(518,2)='SPAIN W COAST CROPS1'
  cnames(519,2)='ESTREMADURA CROPS4'
  cnames(520,2)='MECOAST SPARSE CROPS1'
  cnames(521,2)='BRITTANY CROPS1'
  cnames(522,2)='SYRIAN CROPS2'

!irrigated crops
  cnames(523,2)='NIL VALLEY CROPS1'
  cnames(524,2)='NIL VALLEY CROPS2'
  cnames(525,2)='NIL VALLEY CROPS3'
  cnames(526,2)='NIL VALLEY CROPS4'
  cnames(527,2)='SPANISH IRR CROPS1'
  cnames(528,2)='NIL VALLEY CROPS5'
  cnames(529,2)='EGEE IRR CROPS1'
  cnames(530,2)='MEDITER IRR CROPS1'
  cnames(531,2)='S SPAIN IRR CROPS1'
  cnames(532,2)='NIL VALLEY CROPS6'

!bare land      
  cnames(533,2)='BARE ROCK1'
  cnames(534,2)='BARE ROCK2'
  cnames(535,2)='SANDY DESERT1'
  cnames(536,2)='BARE LAND1'
  cnames(537,2)='BARE LAND2'
  cnames(538,2)='BARE LAND3'
  cnames(539,2)='BARE LAND4'
  cnames(540,2)='BARE LAND5'
  cnames(541,2)='BARE LAND6'
  cnames(542,2)='BARE LAND7'
  cnames(543,2)='BARE LAND8'
  cnames(544,2)='BARE LAND9'
  cnames(545,2)='BARE LAND10'
  cnames(546,2)='BARE LAND11'
  cnames(547,2)='BARE LAND12'
  cnames(548,2)='PERMANENT SNOW1'

!Estuaries and swamp areas      
  cnames(549,2)='WADDEN SEA'
  cnames(550,2)='ESTUARY1'
  cnames(551,2)='ESTUARY2'
  cnames(552,2)='POLAR WETLANDS1'
  cnames(553,2)='ESTUARY3'
  cnames(554,2)='ESTUARY4'
  cnames(555,2)='ESTUARY5'
  cnames(556,2)='ESTUARY6'
  cnames(557,2)='POLAR WETLANDS2'
  cnames(558,2)='SUBPOLAR WETLANDS1'
  cnames(559,2)='SUBPOLAR WETLANDS2'
  cnames(560,2)='SUBPOLAR WETLANDS3'

!urban
  cnames(561,2)='TEMPERATE SUBURBAN1'
  cnames(562,2)='TEMPERATE SUBURBAN2'
  cnames(563,2)='TEMPERATE SUBURBAN3'
  cnames(564,2)='TEMPERATE SUBURBAN4'
  cnames(565,2)='TEMPERATE SUBURBAN5'
  cnames(566,2)='COLD SUBURBAN1'
  cnames(567,2)='WARM SUBURBAN1'
  cnames(568,2)='WARM SUBURBAN2'
  cnames(569,2)='TEMPERATE SUBURBAN6'
  cnames(570,2)='TEMPERATE SUBURBAN7'
  cnames(571,2)='WARM SUBURBAN3'

!added classes of permanent crops      
  cnames(572,2)='SPANISH VINEYARDS1'
  cnames(573,2)='LANGUEDOC VINEYARDS1'
ENDIF

IF (.NOT. ALLOCATED(cname)) ALLOCATE(cname(jpcover))
IF (clang=='EN') THEN
  DO jcov=1,jpcover
    cname(jcov) =   cnames(jcov,1)
  ENDDO
ELSE
  DO jcov=1,jpcover
    cname(jcov) =   cnames(jcov,2)
  ENDDO
END IF
!
ENDIF
!-------------------------------------------------------------------------------
!
!*    9.     Arrange cover (optional nam_pgd_arrange_cover & option to use !gardens or not)
!            ------------------------------------------------------------------------------
!
 CALL arrange_cover(dtco, u%LWATER_TO_NATURE, u%LTOWN_TO_ROCK,                   &
                    xdata_nature,xdata_town,xdata_sea,xdata_water,xdata_vegtype, &
                    xdata_garden,u%LGARDEN, xdata_bld, xdata_wall_o_hor            )
!
!-------------------------------------------------------------------------------
!
!*   10.     LAI for ecoclimap2: climatological or not
!            -----------------------------------------
!
IF (.NOT.u%LECOSG) CALL ecoclimap2_lai(dtco%NYEAR)
!
!-------------------------------------------------------------------------------
!
!*    11.    Secondary variables on natural covers
!            -------------------------------------
!
 CALL ini_data_param(ph_tree=xdata_h_tree,plai=xdata_lai, &
                     palbnir_veg=xdata_albnir_veg, palbvis_veg=xdata_albvis_veg,                    &
                     palbuv_veg=xdata_albuv_veg, prsmin=xdata_rsmin,                                &
                     prgl=xdata_rgl, pcv=xdata_cv, pgamma=xdata_gamma,                              &
                     pgmes=xdata_gmes, pgc=xdata_gc, pbslai=xdata_bslai,                            &
                     psefold=xdata_sefold, plaimin_out=xdata_laimin, pdmax=xdata_dmax,              &
                     pstress=xdata_stress, pf2i=xdata_f2i, pveg_out=xdata_veg,                      &
                     pgreen=xdata_green, pz0=xdata_z0, pz0_o_z0h=xdata_z0_o_z0h,                    &
                     pemis_eco=xdata_emis_eco, pwrmax_cf=xdata_wrmax_cf,                            &
                     proot_lin=xdata_root_lin, proot_extinction=xdata_root_extinction,              &
                     psoilrc_so2=xdata_soilrc_so2, psoilrc_o3=xdata_soilrc_o3, pre25=xdata_re25,    &
                     pce_nitro=xdata_ce_nitro,pcf_nitro=xdata_cf_nitro,pcna_nitro=xdata_cna_nitro,  &
                     pgmes_st=xdata_gmes_st, pgc_st=xdata_gc_st, pbslai_st=xdata_bslai_st,          &
                     psefold_st=xdata_sefold_st, pdmax_st=xdata_dmax_st,                            &
                     pgndlitter=xdata_gndlitter, ph_veg=xdata_h_veg, pz0litter=xdata_z0litter      )
!
idc = 1
!
!
! type dominant par cover
! veg roughness length par cover
! sand & clay par cover
! veg par cover
! albedo par cover (min max si varie)
! LAI par cover
! rsmin par cover
! root depth par cover
! orography par covers
!
IF (lhook) CALL dr_hook('INI_DATA_COVER',1,zhook_handle)
!
!------------------------------------------------------------------------------
!
END SUBROUTINE ini_data_cover