carbon_spinup.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 carbon_spinup(KMONTH,KDAY,PTIME,                      &
                   ospinupcarbs, ospinupcarbw, pspinmaxs, pspinmaxw,   &
                   knbyearspins, knbyearspinw, knbyearsold, hphoto,    &
                   hrespsl, kspins, kspinw                             )
  
!     #######################################################################
!
!
!!****  *CARBON_SPINUP*  
!!
!!    PURPOSE
!!    -------
!     Number of times the accelerated subroutine is called  
!     for each time step  
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!      
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      R. Alkama           * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original      03/26/2012
!!
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!                              
USE modi_spinup_max
!                              
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
 CHARACTER(LEN=3), INTENT(IN) :: hrespsl                ! Soil Respiration
!                                                      ! 'DEF' = Norman 1992
!                                                      ! 'PRM' = Rivalland PhD Thesis (2003)
!                                                      ! 'CNT' = CENTURY model (Gibelin 2008)
 CHARACTER(LEN=3), INTENT(IN) :: hphoto                 ! type of photosynthesis
!
INTEGER,          INTENT(IN)  :: kmonth                ! Current month
INTEGER,          INTENT(IN)  :: kday                  ! Current day
REAL,             INTENT(IN)  :: ptime                 ! Current time
LOGICAL,          INTENT(IN)  :: ospinupcarbs          ! T: do the soil carb spinup, F: no
LOGICAL,          INTENT(IN)  :: ospinupcarbw          ! T: do the wood carb spinup, F: no
REAL,             INTENT(IN)  :: pspinmaxs             ! max number of times CARBON_SOIL subroutine
REAL,             INTENT(IN)  :: pspinmaxw             ! max number of times the wood is accelerated
INTEGER,          INTENT(IN)  :: knbyearspins          ! nbr years needed to reaches soil equilibrium
INTEGER,          INTENT(IN)  :: knbyearspinw          ! nbr years needed to reaches wood equilibrium
!
INTEGER,          INTENT(INOUT) :: knbyearsold         ! nbr years executed at curent time step
INTEGER,          INTENT(OUT)   :: kspins              ! number of times the soil is accelerated
INTEGER,          INTENT(OUT)   :: kspinw              ! number of times the wood is accelerated
!
!*      0.    declarations of local variables
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('CARBON_SPINUP',0,zhook_handle)
!
!       1.     Initializations
!              ---------------
!
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! number of times CARBON_SOIL subroutine is called for each time step
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
kspins =1
IF ( ospinupcarbs .AND. hphoto/='NON' .AND. hrespsl=='CNT' ) THEN
   CALL spinup_max(pspinmaxs,knbyearspins,knbyearsold,kspins)
ENDIF
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! number of times  WOOD carbon subroutine is called for each time step
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
kspinw=1
IF ( ospinupcarbw .AND. hphoto=='NCB' ) THEN
   CALL spinup_max(pspinmaxw,knbyearspinw,knbyearsold,kspinw)
ENDIF
!
IF (kmonth == 1 .AND. kday==1 .AND. ptime == 0.0 )THEN
   knbyearsold = knbyearsold + 1
ENDIF
!
IF (lhook) CALL dr_hook('CARBON_SPINUP',1,zhook_handle)
!
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE