convert_patch_garden.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_patch_garden (DTCO, DTI, I, TGDO, TGDPE, TGDP, TOP, TVG, &
                                 klu,kdecade)
!#############################################################
!
!!****  *CONVERT_PATCH_GARDEN* - routine to initialize GARDEN parameters 
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      A. Lemonsu  *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    09/2009
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_data_isba_n, ONLY : data_isba_t
USE modd_isba_n, ONLY : isba_t
USE modd_teb_garden_option_n, ONLY : teb_garden_options_t
USE modd_teb_garden_pgd_evol_n, ONLY : teb_garden_pgd_evol_t
USE modd_teb_garden_pgd_n, ONLY : teb_garden_pgd_t
USE modd_teb_option_n, ONLY : teb_options_t
USE modd_teb_veg_n, ONLY : teb_veg_options_t
!
USE modd_type_date_surf,    ONLY: date_time
!
USE modi_convert_patch_isba
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: dtco
TYPE(data_isba_t), INTENT(INOUT) :: dti
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(teb_garden_options_t), INTENT(INOUT) :: tgdo
TYPE(teb_garden_pgd_evol_t), INTENT(INOUT) :: tgdpe
TYPE(teb_garden_pgd_t), INTENT(INOUT) :: tgdp
TYPE(teb_options_t), INTENT(INOUT) :: top
TYPE(teb_veg_options_t), INTENT(INOUT) :: tvg
!
INTEGER, INTENT(IN) :: klu     ! number of points
INTEGER, INTENT(IN) :: kdecade ! number of decades
!
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL,             DIMENSION(KLU,1)               :: zlai
REAL,             DIMENSION(KLU,1)               :: zveg
REAL,             DIMENSION(KLU,1)               :: zz0
REAL,             DIMENSION(KLU,1)               :: zemis
REAL,             DIMENSION(KLU,1)               :: zrsmin
REAL,             DIMENSION(KLU,1)               :: zgamma
REAL,             DIMENSION(KLU,1)               :: zwrmax_cf
REAL,             DIMENSION(KLU,1)               :: zrgl
REAL,             DIMENSION(KLU,1)               :: zcv
REAL,             DIMENSION(KLU,1)               :: zgmes
REAL,             DIMENSION(KLU,1)               :: zbslai
REAL,             DIMENSION(KLU,1)               :: zlaimin
REAL,             DIMENSION(KLU,1)               :: zsefold
REAL,             DIMENSION(KLU,1)               :: zgc
REAL,             DIMENSION(KLU,1)               :: zdmax
REAL,             DIMENSION(KLU,1)               :: zf2i
REAL,             DIMENSION(KLU,1)               :: zalbnir_veg
REAL,             DIMENSION(KLU,1)               :: zalbvis_veg
REAL,             DIMENSION(KLU,1)               :: zalbuv_veg
REAL,             DIMENSION(KLU,1)               :: zce_nitro
REAL,             DIMENSION(KLU,1)               :: zcf_nitro
REAL,             DIMENSION(KLU,1)               :: zcna_nitro
REAL,             DIMENSION(KLU,1)               :: zre25
REAL,             DIMENSION(KLU,1)               :: zh_tree
REAL,             DIMENSION(KLU,1)               :: zz0_o_z0h
REAL,             DIMENSION(KLU,1)               :: zd_ice
REAL,             DIMENSION(KLU,TGDO%NGROUND_LAYER,1) :: zrootfrac
REAL,             DIMENSION(KLU,TGDO%NGROUND_LAYER,1) :: zdg
REAL,             DIMENSION(KLU,1)               :: zdroot
REAL,             DIMENSION(KLU,1)               :: zdg2
INTEGER,          DIMENSION(KLU,1)               :: iwg_layer
LOGICAL,          DIMENSION(KLU,1)               :: gstress
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!               Initialisation for IO
!
IF (lhook) CALL dr_hook('CONVERT_PATCH_GARDEN',0,zhook_handle)
!
!-------------------------------------------------------------------------------
!
!
  CALL convert_patch_isba(dtco, dti, i, &
                          tvg%CISBA,kdecade,kdecade,top%XCOVER,top%LCOVER,tvg%CPHOTO,.false.,  &
                        .false.,tvg%LTR_ML,'GRD',pveg=zveg,plai=zlai,              &
                        prsmin=zrsmin,pgamma=zgamma,pwrmax_cf=zwrmax_cf,       &
                        prgl=zrgl,pcv=zcv,psoilgrid=tgdo%XSOILGRID,                 &
                        pdg=zdg,kwg_layer=iwg_layer,pdroot=zdroot,pdg2=zdg2,   &
                        pz0=zz0,pz0_o_z0h=zz0_o_z0h,                           &
                        palbnir_veg=zalbnir_veg,palbvis_veg=zalbvis_veg,       &
                        palbuv_veg=zalbuv_veg,pemis_eco=zemis,                 &
                        pvegtype=tgdp%XVEGTYPE,prootfrac=zrootfrac,                 &
                        pgmes=zgmes,pbslai=zbslai,plaimin=zlaimin,             &
                        psefold=zsefold,pgc=zgc,                               &
                        pdmax=zdmax,pf2i=zf2i,ostress=gstress,ph_tree=zh_tree, &
                        pre25=zre25,pce_nitro=zce_nitro,pcf_nitro=zcf_nitro,   &
                        pcna_nitro=zcna_nitro,pd_ice=zd_ice                    )
!
tgdpe%CUR%XVEG         = zveg(:,1)
tgdpe%CUR%XLAI         = zlai(:,1)
tgdpe%CUR%XZ0          = zz0(:,1)
tgdpe%CUR%XEMIS        = zemis(:,1)
tgdp%XRSMIN       = zrsmin(:,1)
tgdp%XGAMMA       = zgamma(:,1)
tgdp%XWRMAX_CF    = zwrmax_cf(:,1)
tgdp%XRGL         = zrgl(:,1)
tgdp%XCV          = zcv(:,1)
tgdp%XGMES        = zgmes(:,1)
tgdp%XBSLAI       = zbslai(:,1)
tgdp%XLAIMIN      = zlaimin(:,1)
tgdp%XSEFOLD      = zsefold(:,1)
tgdp%XGC          = zgc(:,1)
tgdp%XDMAX        = zdmax(:,1)
tgdp%XF2I         = zf2i(:,1)
tgdp%LSTRESS      = gstress(:,1)
tgdp%XALBNIR_VEG  = zalbnir_veg(:,1)
tgdp%XALBVIS_VEG  = zalbvis_veg(:,1)
tgdp%XALBUV_VEG   = zalbuv_veg(:,1)
tgdp%XCE_NITRO    = zce_nitro(:,1)
tgdp%XCF_NITRO    = zcf_nitro(:,1)
tgdp%XCNA_NITRO   = zcna_nitro(:,1)
tgdp%XH_TREE      = zh_tree(:,1)
tgdp%XRE25        = zre25(:,1)
tgdp%XROOTFRAC    = zrootfrac(:,:,1)
tgdp%XDG          = zdg(:,:,1)
tgdp%XDROOT       = zdroot(:,1)
tgdp%XDG2         = zdg2(:,1)
tgdp%NWG_LAYER    = iwg_layer(:,1)
tgdp%XZ0_O_Z0H    = zz0_o_z0h(:,1)
tgdp%XD_ICE       = zd_ice(:,1)
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('CONVERT_PATCH_GARDEN',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
!
END SUBROUTINE convert_patch_garden