canopy_evol_temp.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 canopy_evol_temp(KI,KLVL,PTSTEP,KIMPL,PTHA,PK,PDKDDVDZ,           &
                                   psflux_t,pforc_t,pdforc_tdt,pdz,pdzf,pth,pwth,   &
                                   palfa,pbeta                                      ) 
!     #########################################
!
!!****  *CANOPY_EVOL_TEMP* - evolution of wind in canopy 
!!                        
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    07/2006 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modi_canopy_evol_field
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
INTEGER,                  INTENT(IN)    :: ki        ! number of horizontal points
INTEGER,                  INTENT(IN)    :: klvl      ! number of levels in canopy
REAL,                     INTENT(IN)    :: ptstep    ! time-step                             (s)
INTEGER,                  INTENT(IN)    :: kimpl     ! implicitation code: 
!                                                    ! 1 : computes only alfa and beta coupling
!                                                    !     coefficients for all variables
!                                                    ! 2 : computes temporal evolution of the
!                                                    !     variables
REAL, DIMENSION(KI),      INTENT(IN)    :: ptha      ! potential temp.  at forcing level     (K)
REAL, DIMENSION(KI,KLVL), INTENT(IN)    :: pk        ! mixing exchange coefficient           (m2/s)
REAL, DIMENSION(KI,KLVL), INTENT(IN)    :: pdkddvdz  ! derivative of mixing coefficient as a
!                                                    ! function of vertical gradient of wind
!                                                    ! (at mid levels)                       (m2)
REAL, DIMENSION(KI),      INTENT(IN)    :: psflux_t  ! surface flux w'Th'                    (mK/s)
REAL, DIMENSION(KI,KLVL), INTENT(IN)    :: pforc_t   ! tendency of wind due to canopy drag   (K/s)
REAL, DIMENSION(KI,KLVL), INTENT(IN)    :: pdforc_tdt! formal derivative of the tendency of
!                                                    ! wind due to canopy drag               (1/s)
REAL, DIMENSION(KI,KLVL), INTENT(IN)    :: pdz       ! deltaZ between canopy half levels,
!                                                    ! located at full levels                (m)
REAL, DIMENSION(KI,KLVL), INTENT(IN)    :: pdzf      ! deltaZ between canopy (full) levels,
!                                                    ! located at half levels                (m)
REAL, DIMENSION(KI,KLVL), INTENT(INOUT) :: pth       ! pot. temp. at canopy levels           (K)
REAL, DIMENSION(KI,KLVL), INTENT(OUT)   :: pwth      ! turbulent flux (at half levels)       (mK/s)
REAL, DIMENSION(KI),      INTENT(OUT)   :: palfa     !  V+(1) = alfa F(1) + beta
REAL, DIMENSION(KI),      INTENT(OUT)   :: pbeta     !  V+(1) = alfa F(1) + beta
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER                     :: jlayer   ! loop counter on layers
!
REAL, DIMENSION(KI,KLVL)   :: zdthdz   ! dTh/dz at mid levels
REAL, DIMENSION(KI,KLVL)   :: zf       ! turbulent flux at mid levels
REAL, DIMENSION(KI,KLVL)   :: zdfddvdz ! derivative of turbulent flux as a
!                                       ! function of vertical gradient of wind variable
!                                       ! (at mid levels)
REAL, DIMENSION(KI,KLVL)   :: zext     ! external forcing at full levels
REAL, DIMENSION(KI,KLVL)   :: zdextdv  ! derivative of external forcing as a
!                                       ! function of vertical variable
!                                       ! (at full levels)
REAL, DIMENSION(KI,KLVL)   :: zth      ! work variable : pot. temp at futur instant 
!                                       ! (or past at the end of the routine) 
REAL, DIMENSION(KI)         :: zdthadt  ! dTHa/dt   at forcing level
REAL, DIMENSION(KI)         :: zdwthdz  ! dw'Th'/dz at forcing level
LOGICAL                     :: limpl
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!
!*    1. initializations
!        ---------------
!
!* external forces
!
IF (lhook) CALL dr_hook('CANOPY_EVOL_TEMP',0,zhook_handle)
zext = 0.
zdextdv = 0.
!
!-------------------------------------------------------------------------------
!
!*    5. Forcing due to drag (at full levels)
!        -------------------
!
zext = zext + pforc_t
zdextdv = zdextdv + pdforc_tdt
!
!-------------------------------------------------------------------------------
!
!*    6. External forcing due to large-scale forces (at full levels)
!        ------------------------------------------
!
!
!* forces due to large-scale forcing
!
! These are computed from wind evolution equation at forcing level :
!
! dUa/dt = Large_Scale_Forcing - d(u'w')/dz|z=forcing_level
!
! because vertical derivative of turbulent flux is not available at forcing
! level, one must make the assumption that the turbulent flux is uniform
! between the forcing level and the level just below.
! This means that one assume that the forcing layer is in the inertial sublayer
! (where turbulent fluxes are constant).
!
zdthadt(:) = ( ptha(:) - pth(:,klvl) ) /ptstep
!
zdwthdz(:) = 0.
!
DO jlayer=1,klvl
  zext(:,jlayer) = zext(:,jlayer) + zdwthdz(:) + zdthadt(:)
END DO
!
 CALL  canopy_evol_field(ki, klvl, ptstep, kimpl, pk, pdkddvdz,   &
                        psflux_t, pforc_t, pdforc_tdt, pdz, pdzf,      &
                        zext, zdextdv, pth, pwth, palfa, pbeta         ) 
!
IF (lhook) CALL dr_hook('CANOPY_EVOL_TEMP',1,zhook_handle)
!----------------------------------------------------------------
!
END SUBROUTINE canopy_evol_temp