permafrost_depth.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 permafrost_depth (KNI,KPATCH,PPERM,PSOILDEPTH)
!     ###################################################
!
!!****  *PERMAFROST_DEPTH*  
!!
!!    PURPOSE
!!    -------
! 
!     Extended  ground depth to 12m over permafrost area
!
!!**  METHOD
!!    ------
!
!     Direct calculation
!
!!    EXTERNAL
!!    --------
!
!     None
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!      
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme     
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    30/08/12 
!-------------------------------------------------------------------------------
!
USE modd_surf_par, ONLY : xundef
USE modd_isba_par, ONLY : xpermfrac, xpermdepth
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
!*      0.1    declarations of arguments
!
IMPLICIT NONE
!
INTEGER,             INTENT(IN   ) :: kni        ! number of point
!
INTEGER,             INTENT(IN   ) :: kpatch     ! patch number
!
REAL, DIMENSION(:),  INTENT(IN   ) :: pperm      ! permafrost area (fraction)
!
REAL, DIMENSION(:,:),INTENT(INOUT) :: psoildepth ! output soil depth distribution (m)
!
!*      0.2    declarations of local variables
!
REAL, DIMENSION(KNI) :: zperm
!
INTEGER :: jj, jpatch
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PERMAFROST_DEPTH',0,zhook_handle)
!
zperm(:)=0.0
WHERE(pperm(:)/=xundef)zperm(:)=pperm(:)
!
DO jpatch=1,kpatch
   DO jj=1,kni
      IF(zperm(jj)>=xpermfrac.AND.psoildepth(jj,jpatch)/=xundef)THEN
         psoildepth(jj,jpatch)=max(psoildepth(jj,jpatch),xpermdepth)
       ENDIF
   ENDDO
ENDDO

!
IF (lhook) CALL dr_hook('PERMAFROST_DEPTH',1,zhook_handle)
!
END SUBROUTINE permafrost_depth