set_rough.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 set_rough(OCANOPY,HROUGH)
!     ########################################################################
!
!!****  *SET_ROUGH* - routine to set default values for roughness ISBA configuration
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_surf_par,   ONLY : xundef
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
LOGICAL,           INTENT(IN)     :: ocanopy    ! True  = SBL scheme
                                                ! False = no SBL scheme
 CHARACTER(LEN=4),  INTENT(INOUT)  :: hrough     ! type of roughness length
REAL(KIND=JPRB) :: zhook_handle
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
!-------------------------------------------------------------------------------
!
!* default value if needed
!
IF (lhook) CALL dr_hook('SET_ROUGH',0,zhook_handle)
IF (hrough  == "UNDE") THEN
  IF (ocanopy) THEN
    hrough = "BE04"
  ELSE
    hrough = "Z04D"
  END IF
END IF
!
!* coherence check
!
IF (ocanopy .AND. (hrough=="Z01D" .OR. hrough=="Z04D")) THEN
  CALL abor1_sfx('If LISBA_CANOPY=T, orographic roughness cannot be computed by CROUGH="Z01D" or "Z04D"')
END IF
!
IF (.NOT. ocanopy .AND. hrough=="BE04") THEN
  CALL abor1_sfx('If LISBA_CANOPY=F, orographic roughness cannot be computed by CROUGH="BE04"')
END IF
IF (lhook) CALL dr_hook('SET_ROUGH',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE set_rough