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.
MODULE modi_permafrost_depth
!
INTERFACE permafrost_depth
!
SUBROUTINE permafrost_depth_1d (KNI,KPATCH,PPERM,PSOILDEPTH)
!
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)
!
END SUBROUTINE permafrost_depth_1d
!
SUBROUTINE permafrost_depth_2d (KNI,KPATCH,PPERM,PSOILDEPTH)
!
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)
!
END SUBROUTINE permafrost_depth_2d
!
END INTERFACE permafrost_depth
!
END MODULE modi_permafrost_depth
!
!     #########################
SUBROUTINE permafrost_depth_1d (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_1D',0,zhook_handle)
!
zperm(:)=0.0
WHERE(pperm(:)/=xundef)zperm(:)=pperm(:)
!
DO jj=1,kni
  IF(zperm(jj)>=xpermfrac.AND.psoildepth(jj)/=xundef)THEN
    psoildepth(jj)=max(psoildepth(jj),xpermdepth)
  ENDIF
ENDDO

!
IF (lhook) CALL dr_hook('PERMAFROST_DEPTH_1D',1,zhook_handle)
!
END SUBROUTINE permafrost_depth_1d
!
!     #########################
SUBROUTINE permafrost_depth_2d (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_2D',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_2D',1,zhook_handle)
!
END SUBROUTINE permafrost_depth_2d