isba_flood_properties.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 isba_flood_properties(PLAI,PFFLOOD,PFFROZEN,PZ0FLOOD, &
                                       pffg_nosnow,pffv_nosnow         )  
!     ############################################################################
!
!
!!****  *ISBA_FLOOD_PROPERTIES*  
!!
!!    PURPOSE
!!    -------
!      Calculate the Flood fractions without snow and roughness length.  
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!      
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme           * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original      25/05/08 
!!      Modified      09/2009  B. Decharme: limitation of Ri in surface_ri.F90
!!      Modified      02/2014  B. Decharme: Simplification
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_snow_par,   ONLY : xz0sn
USE modd_flood_par,  ONLY : xz0flood
!
USE mode_surf_flood_frac
!
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)       :: pfflood
REAL, DIMENSION(:), INTENT(IN)       :: pffrozen
!
REAL, DIMENSION(:), INTENT(OUT)      :: pz0flood   ! roughness length over floodplains
REAL, DIMENSION(:), INTENT(OUT)      :: pffg_nosnow
REAL, DIMENSION(:), INTENT(OUT)      :: pffv_nosnow
!
!*      0.2    declarations of local variables
!
REAL, DIMENSION(SIZE(PFFLOOD)) :: zpsng
REAL, DIMENSION(SIZE(PFFLOOD)) :: zpsnv
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!       0.     Initializations
!              ---------------
!
IF (lhook) CALL dr_hook('ISBA_FLOOD_PROPERTIES',0,zhook_handle)
!
!-------------------------------------------------------------------------------
!
!       1.     Flood fractions without snow
!              ----------------------------
!
zpsng(:)=0.0
zpsnv(:)=0.0
!
pffg_nosnow(:) = flood_frac_ground(zpsng,pfflood)
pffv_nosnow(:) = flood_frac_veg(plai,zpsnv,pfflood)
!
!-------------------------------------------------------------------------------
!
!       2.     roughness length
!              ----------------
!
WHERE(pffrozen(:)==0.0)
    pz0flood(:) = xz0flood
ELSEWHERE
    pz0flood(:) = min(0.005,10.*xz0flood)
END WHERE
!
IF (lhook) CALL dr_hook('ISBA_FLOOD_PROPERTIES',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE isba_flood_properties