control_temp_func.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.
!   ################
FUNCTION control_temp_func (PTEMP_IN) RESULT (PTEMPFUNC_RESULT)

!
!   ###############################################################
!!**   CONTROL_TEMP_FUNC
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!     Temperature control factor of decomposition.
!!
!!    EXTERNAL
!!    --------
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      
!!    none
!!
!!    REFERENCE
!!    ---------
!!
!!      Krinner et al., Global Biochemical Cycles, 2005
!!      Gibelin et al. 2008, AFM
!!      
!!    AUTHOR
!!    ------
!!
!!      A.-L. Gibelin           * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    23/06/09
!!      B. Decharme 05/2012 : Optimization and ISBA-DIF coupling
!!      
!-------------------------------------------------------------------------------

!
!*       0.     DECLARATIONS
!               ------------
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
REAL, DIMENSION(:), INTENT(IN)  :: ptemp_in ! temperature (K)
!
!*      0.2    declarations of local variables
!
REAL                            :: zcoef1
REAL, PARAMETER                 :: zcoef2 = 10.0
REAL, PARAMETER                 :: zcoef3 = 30.0
!
REAL, DIMENSION(SIZE(PTEMP_IN)) :: ptempfunc_result ! temperature control factor
!
REAL(KIND=JPRB) :: zhook_handle
!
!
!*       1 Calculates temperature control factor
!
IF (lhook) CALL dr_hook('CONTROL_TEMP_FUNC',0,zhook_handle)
!
zcoef1 = log(2.0)
!
ptempfunc_result(:) = exp( (zcoef1/zcoef2) * (ptemp_in(:)-zcoef3) )
!
ptempfunc_result(:) = min( 1., ptempfunc_result(:) )
!
IF (lhook) CALL dr_hook('CONTROL_TEMP_FUNC',1,zhook_handle)

END FUNCTION control_temp_func