co2_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_init_n (IO, S, PK, PEK, KSIZE, PCO2   )
!     #####################
!
!!****  *CO2_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_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY : isba_p_t, isba_s_t, isba_pe_t
!
USE modd_surfex_mpi, ONLY : nrank,npio
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_options_t), INTENT(INOUT) :: IO
TYPE(isba_s_t), INTENT(INOUT) :: S
TYPE(isba_p_t), INTENT(INOUT) :: PK
TYPE(isba_pe_t), INTENT(INOUT) :: PEK
!
INTEGER, INTENT(IN) :: KSIZE
!
REAL, DIMENSION(:), INTENT(IN) :: PCO2 ! air CO2 concentration (kg/kg)
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(:),   ALLOCATABLE :: ZP_CO2            ! air CO2 concentration (kg/kg)
!
INTEGER :: IPATCH
INTEGER :: INBIOMASS
INTEGER :: JP    ! loop on tiles
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('CO2_INIT_N',0,zhook_handle)
!
inbiomass = SIZE(pk%XINCREASE,2)
!
IF (ksize == 0 ) THEN
  IF (lhook) CALL dr_hook('CO2_INIT_N',1,zhook_handle)
  RETURN
ENDIF
!
IF (maxval(pek%XGMES(:)).NE.xundef .OR. minval(pek%XGMES(:)).NE.xundef) THEN

  CALL pack_co2_init(pk%NR_P(:),ksize)
  !
  CALL cotwoinit_n(io, s, pk, pek, zp_co2 )

  pk%XINCREASE = 0.
  pk%XTURNOVER = 0.
  !

ENDIF
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('CO2_INIT_N',1,zhook_handle)
CONTAINS
!-------------------------------------------------------------------------------
SUBROUTINE pack_co2_init(KMASK,KSIZE)
IMPLICIT NONE
INTEGER, INTENT(IN)               :: KSIZE
INTEGER, DIMENSION(:), INTENT(IN) :: KMASK
!
INTEGER JJ, JI
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('PACK_CO2_INIT',0,zhook_handle)
!
ALLOCATE(zp_co2(ksize))
!
DO jj=1,ksize
  ji                     =    kmask(jj)
  zp_co2(jj)   =    pco2(ji)
END DO
!
IF (lhook) CALL dr_hook('PACK_CO2_INIT',1,zhook_handle)
!-------------------------------------------------------------------------------
END SUBROUTINE pack_co2_init
!
END SUBROUTINE co2_init_n