co2_teb_greenroof_initn.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 co2_teb_greenroof_init_n (I, TGRP, TVG, &
                                           pco2)
!     #####################
!
!!****  *CO2_TEB_GREENROOF_INIT_n* - routine to initialize ISBA-AGS variables
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2003 
!!      J.C. Calvet 01/2004 Externalization
!!      P Le Moigne 11/2004 cotwoinit changed into cotwoinit_n
!!      P Le Moigne 09/2005 AGS modifs of L. Jarlan
!!      S Lafont    09/2008 Add initialisation of POI and ABC (needed for TORI)
!!      A.L. Gibelin 04/2009 : TAU_WOOD for NCB option 
!!      A.L. Gibelin 04/2009 : Add carbon spinup
!!      A.L. Gibelin 07/2009 : Suppress RDK and transform GPP as a diagnostic
!!      A.L. Gibelin 07/2009 : Suppress PPST and PPSTF as outputs
!!
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
USE modd_isba_n, ONLY : isba_t
USE modd_teb_greenroof_pgd_n, ONLY : teb_greenroof_pgd_t
USE modd_teb_veg_n, ONLY : teb_veg_options_t
!
USE modd_surf_par,       ONLY : xundef
USE modd_data_cover_par, ONLY : nvegtype
!
USE modi_cotwoinit_n
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(teb_greenroof_pgd_t), INTENT(INOUT) :: tgrp
TYPE(teb_veg_options_t), INTENT(INOUT) :: tvg
!
REAL, DIMENSION(:), INTENT(IN) :: pco2 ! air CO2 concentration (kg/kg)
!
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(SIZE(TGRP%XVEGTYPE,1)) :: ztau_wood
INTEGER :: ilu   ! size of arrays
INTEGER :: jp    ! loop on tiles
!
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('CO2_TEB_GREENROOF_INIT_N',0,zhook_handle)
ilu = SIZE(tgrp%XVEGTYPE,1)
!
ALLOCATE(tgrp%XANMAX        (ilu))
ALLOCATE(tgrp%XFZERO        (ilu))
ALLOCATE(tgrp%XEPSO         (ilu))
ALLOCATE(tgrp%XGAMM         (ilu))
ALLOCATE(tgrp%XQDGAMM       (ilu))
ALLOCATE(tgrp%XQDGMES       (ilu))
ALLOCATE(tgrp%XT1GMES       (ilu))
ALLOCATE(tgrp%XT2GMES       (ilu))
ALLOCATE(tgrp%XAMAX         (ilu))
ALLOCATE(tgrp%XQDAMAX       (ilu))
ALLOCATE(tgrp%XT1AMAX       (ilu))
ALLOCATE(tgrp%XT2AMAX       (ilu))
ALLOCATE(tgrp%XAH           (ilu))
ALLOCATE(tgrp%XBH           (ilu))
!
     CALL cotwoinit_n(i, &
                      tvg%CPHOTO, tgrp%XVEGTYPE,tgrp%XGMES,pco2,tgrp%XGC,&
            tgrp%XDMAX,tgrp%XABC,tgrp%XPOI,tgrp%XANMAX, tgrp%XFZERO,           &
            tgrp%XEPSO,tgrp%XGAMM,tgrp%XQDGAMM,tgrp%XQDGMES,tgrp%XT1GMES,      &
            tgrp%XT2GMES,tgrp%XAMAX,tgrp%XQDAMAX,tgrp%XT1AMAX,            &
            tgrp%XT2AMAX,tgrp%XAH,tgrp%XBH,ztau_wood                 )  
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('CO2_TEB_GREENROOF_INIT_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE co2_teb_greenroof_init_n