z0v_from_lai.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_z0v_from_lai
!#######################
!
INTERFACE z0v_from_lai
!
    FUNCTION z0v_from_lai_0d(PLAI,PH_TREE,PVEGTYPE,OAGRI_TO_GRASS) RESULT(PZ0)
!
REAL,                 INTENT(IN) :: PLAI         ! Leaf area Index
REAL,                 INTENT(IN) :: PH_TREE      ! height of trees
REAL,   DIMENSION(:), INTENT(IN) :: PVEGTYPE     ! type of vegetation
LOGICAL,              INTENT(IN) :: OAGRI_TO_GRASS
!
REAL                             :: PZ0          ! vegetation roughness
!
END FUNCTION z0v_from_lai_0d
!
!
    FUNCTION z0v_from_lai_1d(PLAI,PH_TREE,PVEGTYPE,OAGRI_TO_GRASS) RESULT(PZ0)
!
REAL,   DIMENSION(:),   INTENT(IN) :: PLAI         ! Leaf area Index
REAL,   DIMENSION(:),   INTENT(IN) :: PH_TREE      ! height of trees
REAL,   DIMENSION(:,:), INTENT(IN) :: PVEGTYPE     ! type of vegetation
LOGICAL,                INTENT(IN) :: OAGRI_TO_GRASS
!
REAL,   DIMENSION(SIZE(PLAI))      :: PZ0          ! vegetation roughness
!
END FUNCTION z0v_from_lai_1d
!
!
    FUNCTION z0v_from_lai_2d(PLAI,PH_TREE,PVEGTYPE,OAGRI_TO_GRASS) RESULT(PZ0)
!
REAL,   DIMENSION(:,:),   INTENT(IN) :: PLAI         ! Leaf area Index
REAL,   DIMENSION(:,:),   INTENT(IN) :: PH_TREE      ! height of trees
REAL,   DIMENSION(:,:,:), INTENT(IN) :: PVEGTYPE     ! type of vegetation
LOGICAL,                  INTENT(IN) :: OAGRI_TO_GRASS
!
REAL,   DIMENSION(SIZE(PLAI,1),SIZE(PLAI,2)) :: PZ0  ! vegetation roughness
!
END FUNCTION z0v_from_lai_2d
!
    FUNCTION z0v_from_lai_vegtype(PLAI,PH_TREE,OAGRI_TO_GRASS) RESULT(PZ0)
!
REAL,   DIMENSION(:),   INTENT(IN) :: PLAI         ! Leaf area Index
REAL,   DIMENSION(:),   INTENT(IN) :: PH_TREE      ! height of trees
LOGICAL,                INTENT(IN) :: OAGRI_TO_GRASS
!
REAL,   DIMENSION(SIZE(PLAI)) :: PZ0  ! vegetation roughness
!
END FUNCTION z0v_from_lai_vegtype
!
END INTERFACE
!
END MODULE modi_z0v_from_lai
!

!   ###########################################################
    FUNCTION z0v_from_lai_0d(PLAI,PH_TREE,PVEGTYPE,OAGRI_TO_GRASS) RESULT(PZ0)
!   ###########################################################
!!
!!    PURPOSE
!!    -------
!
!     Calculates vegetation roughness from leaf
!    area index and type of vegetation
!    (most of types; forest and vineyards; grassland)
!              
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      
!!    none
!!
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!
!!      V. Masson and A. Boone          * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    25/03/99
!!      P. Aumond   10/10/2009     Because drag force applied in atmospheric 
!!                                 model, Z0tree -> z0grass
!!      R. Alkama    05/2012   : Extantion from 12 to 19 vegtypes 
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_surf_par,       ONLY : xundef
USE modi_veg_height_from_lai
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
REAL,                 INTENT(IN) :: PLAI         ! Leaf area Index
REAL,                 INTENT(IN) :: PH_TREE      ! height of trees
REAL,   DIMENSION(:), INTENT(IN) :: PVEGTYPE     ! type of vegetation
LOGICAL,              INTENT(IN) :: OAGRI_TO_GRASS
!
REAL                             :: PZ0          ! vegetation roughness
!
!*      0.2    declarations of local variables
!
REAL                            :: ZALLEN_H    ! Allen formula for height
REAL                            :: ZLAI        ! LAI for vegetated areas
!
REAL, DIMENSION(SIZE(PVEGTYPE)) :: ZH_VEG          ! height for each type
REAL                            :: ZAVG_H      ! averaged height
REAL                            :: ZZREF       ! reference height        
!
INTEGER                         :: JTYPE       ! loop counter
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-----------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_0D',0,zhook_handle)
!
zh_veg(:) = veg_height_from_lai(plai,ph_tree,pvegtype,oagri_to_grass)
!
zzref = 10.
zavg_h = 0.
DO jtype=1,SIZE(pvegtype)
  zavg_h = zavg_h + pvegtype(jtype) / (log(0.13*zh_veg(jtype)/zzref))**2
END DO
zavg_h = max(zavg_h,0.00001)

zavg_h = zzref / 0.13 * exp(-1./sqrt(zavg_h))
!
pz0  = max(0.001, 0.13*zavg_h)
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_0D',1,zhook_handle)
!
!-----------------------------------------------------------------
!
END FUNCTION z0v_from_lai_0d
!
!   ###########################################################
    FUNCTION z0v_from_lai_1d(PLAI,PH_TREE,PVEGTYPE,OAGRI_TO_GRASS) RESULT(PZ0)
!   ###########################################################
!!
!!    PURPOSE
!!    -------
!
!     Calculates vegetation roughness from leaf
!    area index and type of vegetation
!    (most of types; forest and vineyards; grassland)
!              
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      
!!    none
!!
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!
!!      V. Masson and A. Boone          * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    25/03/99
!!
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_surf_par,       ONLY : xundef
USE modi_veg_height_from_lai
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
REAL,   DIMENSION(:),   INTENT(IN) :: PLAI         ! Leaf area Index
REAL,   DIMENSION(:),   INTENT(IN) :: PH_TREE      ! height of trees
REAL,   DIMENSION(:,:), INTENT(IN) :: PVEGTYPE     ! type of vegetation
LOGICAL,                INTENT(IN) :: OAGRI_TO_GRASS
!
REAL,   DIMENSION(SIZE(PLAI))      :: PZ0          ! vegetation roughness
!
!*      0.2    declarations of local variables
!
REAL, DIMENSION(SIZE(PLAI))                  :: ZALLEN_H    ! Allen formula for height
REAL, DIMENSION(SIZE(PLAI))                  :: ZLAI        ! LAI for vegetated areas
!
REAL, DIMENSION(SIZE(PLAI),SIZE(PVEGTYPE,2)) :: ZH_VEG          ! height for each type
REAL, DIMENSION(SIZE(PLAI))                  :: ZAVG_H      ! averaged height
REAL                                         :: ZZREF       ! reference height        
!
INTEGER                                      :: JTYPE       ! loop counter
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-----------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_1D',0,zhook_handle)
zh_veg(:,:) = veg_height_from_lai(plai,ph_tree,pvegtype,oagri_to_grass)
!
zzref = 10.
zavg_h(:) = 0.
DO jtype=1,SIZE(pvegtype,2)
  zavg_h(:) = zavg_h(:) + pvegtype(:,jtype) / (log(0.13*zh_veg(:,jtype)/zzref))**2
END DO

zavg_h = max(zavg_h,0.00001)

zavg_h(:) = zzref / 0.13 * exp(-1./sqrt(zavg_h(:)))
!
pz0(:) = max(0.001, 0.13*zavg_h(:))
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_1D',1,zhook_handle)
!-----------------------------------------------------------------
!
END FUNCTION z0v_from_lai_1d
!
!   ###########################################################
    FUNCTION z0v_from_lai_2d(PLAI,PH_TREE,PVEGTYPE,OAGRI_TO_GRASS) RESULT(PZ0)
!   ###########################################################
!!
!!    PURPOSE
!!    -------
!
!     Calculates vegetation roughness from leaf
!    area index and type of vegetation
!    (most of types; forest and vineyards; grassland)
!              
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      
!!    none
!!
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!
!!      V. Masson and A. Boone          * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    25/03/99
!!
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_surf_par,       ONLY : xundef
USE modi_veg_height_from_lai
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
REAL,   DIMENSION(:,:),   INTENT(IN) :: PLAI         ! Leaf area Index
REAL,   DIMENSION(:,:),   INTENT(IN) :: PH_TREE      ! height of trees
REAL,   DIMENSION(:,:,:), INTENT(IN) :: PVEGTYPE     ! type of vegetation
LOGICAL,                  INTENT(IN) :: OAGRI_TO_GRASS
!
REAL,   DIMENSION(SIZE(PLAI,1),SIZE(PLAI,2)) :: PZ0  ! vegetation roughness
!
!*      0.2    declarations of local variables
!

REAL, DIMENSION(SIZE(PLAI,1),SIZE(PLAI,2))                  :: ZALLEN_H ! Allen formula for height
REAL, DIMENSION(SIZE(PLAI,1),SIZE(PLAI,2))                  :: ZLAI     ! LAI for vegetated areas
!
REAL, DIMENSION(SIZE(PLAI,1),SIZE(PLAI,2),SIZE(PVEGTYPE,3)) :: ZH_VEG       ! height for each type
REAL, DIMENSION(SIZE(PLAI,1),SIZE(PLAI,2))                  :: ZAVG_H   ! averaged height
REAL                                                        :: ZZREF    ! reference height        
!
INTEGER                                                     :: JTYPE    ! loop counter
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-----------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_2D',0,zhook_handle)
zh_veg(:,:,:) = veg_height_from_lai(plai,ph_tree,pvegtype,oagri_to_grass)
!
zzref = 10.
zavg_h(:,:) = 0.
DO jtype=1,SIZE(pvegtype,3)
  zavg_h(:,:) = zavg_h(:,:) + pvegtype(:,:,jtype) / (log(0.13*zh_veg(:,:,jtype)/zzref))**2
END DO
zavg_h(:,:) = max(zavg_h(:,:),0.00001)
zavg_h(:,:) = zzref / 0.13 * exp(-1./sqrt(zavg_h(:,:)))
!
pz0(:,:) = max(0.001, 0.13*zavg_h(:,:))
!
WHERE (plai(:,:) == xundef)
  pz0(:,:) = xundef
END WHERE
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_2D',1,zhook_handle)
!-----------------------------------------------------------------
!
END FUNCTION z0v_from_lai_2d
!
!
!
!   ###########################################################
    FUNCTION z0v_from_lai_vegtype(PLAI,PH_TREE,OAGRI_TO_GRASS) RESULT(PZ0)
!   ###########################################################
!!
!!    PURPOSE
!!    -------
!
!     Calculates vegetation roughness from leaf
!    area index and type of vegetation for each patch
!    (most of types; forest and vineyards; grassland)
!              
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      
!!    none
!!
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!        F.Solmon
!!      V. Masson and A. Boone          * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    25/03/99
!!      
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_surf_par,       ONLY : xundef
USE modi_veg_height_from_lai
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
REAL,   DIMENSION(:),   INTENT(IN) :: PLAI         ! Leaf area Index
REAL,   DIMENSION(:),   INTENT(IN) :: PH_TREE      ! height of trees
LOGICAL,                INTENT(IN) :: OAGRI_TO_GRASS
!
REAL,   DIMENSION(SIZE(PLAI))      :: PZ0          ! vegetation roughness
!
!*      0.2    declarations of local variables
!
REAL, DIMENSION(SIZE(PLAI)) :: ZALLEN_H    ! Allen formula for height
!
REAL, DIMENSION(SIZE(PLAI)) :: ZH_VEG          ! height for each type
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-----------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_VEGTYPE',0,zhook_handle)
zh_veg(:) = veg_height_from_lai(plai,ph_tree,oagri_to_grass)
!
pz0(:) = xundef
!
WHERE(zh_veg(:)/=xundef) pz0(:) = max(0.001, 0.13*zh_veg(:)) ! rugosite pour chaque vegtype
!-----------------------------------------------------------------
!
WHERE (plai(:) == xundef)
  pz0(:) = xundef
END WHERE
IF (lhook) CALL dr_hook('MODI_Z0V_FROM_LAI:Z0V_FROM_LAI_VEGTYPE',1,zhook_handle)
!
!
END FUNCTION z0v_from_lai_vegtype
!