modd_sson.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 modd_sso_n
!     ######################
!
!!****  *MODD_SSO - declaration of surface parameters related to orography
!!
!!    PURPOSE
!!    -------
!     Declaration of surface parameters
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      V. Masson    *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       01/2004
!
!*       0.   DECLARATIONS
!             ------------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
! utilisé uniquement tel quel. 
TYPE sso_t
!
!-----------------------------------------------------------------------------------------------------
!
! Type of roughness
!
 CHARACTER(LEN=4) :: crough     ! type of orographic roughness
!                              ! 'NONE'
                               ! 'Z01D'
                               ! 'Z04D'
                               ! 'BE04'
!
  REAL, DIMENSION(:), POINTER   :: xz0effjpdir    ! heading of J direction (deg from N clockwise)
!
  REAL, DIMENSION(:), POINTER   :: xsso_slope         ! slope of S.S.O.
  REAL, DIMENSION(:), POINTER   :: xsso_anis          ! anisotropy of S.S.O.
  REAL, DIMENSION(:), POINTER   :: xsso_dir           ! direction of S.S.O. (deg from N clockwise) 
  REAL, DIMENSION(:), POINTER   :: xsso_stdev         ! S.S.O. standard deviation           (m)
!
  REAL, DIMENSION(:), POINTER   :: xavg_zs        ! averaged orography                      (m)
  REAL, DIMENSION(:), POINTER   :: xsil_zs        ! silhouette orography                    (m)
  REAL, DIMENSION(:), POINTER   :: xmax_zs        ! maximum subgrid orography               (m)
  REAL, DIMENSION(:), POINTER   :: xmin_zs        ! minimum subgrid orography               (m)
! Zo threshold
  REAL   :: xfracz0                                ! Z0=Min(Z0, Href/XFRACZ0)
  REAL   :: xcoefbe                                ! Beljaars coefficient         
!
!-----------------------------------------------------------------------------------------------------
!
! Subgrid orography parameters
!
  REAL, DIMENSION(:), POINTER :: xaosip, xaosim, xaosjp, xaosjm
! directional A/S quantities in 4 coordinate directions
! (IP: i index up;  IM: i index down;  JP: j index up;  JM: j index down)
! They are used in soil routines to compute effective roughness length
!
  REAL, DIMENSION(:), POINTER :: xho2ip, xho2im, xho2jp, xho2jm
! directional h/2 quantities in 4 coordinate directions
! (IP: i index up;  IM: i index down;  JP: j index up;  JM: j index down)
! They are used in soil routines to compute effective roughness length
!
  REAL, DIMENSION(:), POINTER   :: xz0rel         ! relief roughness length     (m)
!
! utilisé par l'atmosphère 
!
  REAL, DIMENSION(:), POINTER :: xz0effip, xz0effim, xz0effjp, xz0effjm
! directional total roughness lenghts in 4 coordinate directions
! (IP: i index up;  IM: i index down;  JP: j index up;  JM: j index down)
!
END TYPE sso_t
!
! utilisé par isba
!
TYPE sso_np_t
!
TYPE(sso_t), DIMENSION(:), POINTER :: al=>null()
!
END TYPE sso_np_t
!
CONTAINS
!
SUBROUTINE sso_init(YSSO)
TYPE(sso_t), INTENT(INOUT) :: YSSO
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_SSO_N:SSO_INIT",0,zhook_handle)
  NULLIFY(ysso%XZ0EFFJPDIR)  
  NULLIFY(ysso%XSSO_SLOPE)
  NULLIFY(ysso%XSSO_ANIS)
  NULLIFY(ysso%XSSO_DIR)
  NULLIFY(ysso%XSSO_STDEV)
  NULLIFY(ysso%XAVG_ZS)
  NULLIFY(ysso%XSIL_ZS)
  NULLIFY(ysso%XMAX_ZS)
  NULLIFY(ysso%XMIN_ZS)
  NULLIFY(ysso%XAOSIP)
  NULLIFY(ysso%XAOSIM)
  NULLIFY(ysso%XAOSJP)
  NULLIFY(ysso%XAOSJM)
  NULLIFY(ysso%XHO2IP)
  NULLIFY(ysso%XHO2IM)
  NULLIFY(ysso%XHO2JP)
  NULLIFY(ysso%XHO2JM)
  NULLIFY(ysso%XZ0REL)

  NULLIFY(ysso%XZ0EFFIP)
  NULLIFY(ysso%XZ0EFFIM)
  NULLIFY(ysso%XZ0EFFJP)
  NULLIFY(ysso%XZ0EFFJM)
  
ysso%CROUGH=' '
ysso%XFRACZ0=2.
ysso%XCOEFBE=2.
IF (lhook) CALL dr_hook("MODD_SSO_N:SSO_INIT",1,zhook_handle)
END SUBROUTINE sso_init
!
SUBROUTINE sso_np_init(YSSO_NP,KPATCH)
TYPE(sso_np_t), INTENT(INOUT) :: YSSO_NP
INTEGER, INTENT(IN) :: KPATCH
INTEGER :: JP
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_SSO_N:SSO_NP_INIT",0,zhook_handle)
!
IF (ASSOCIATED(ysso_np%AL)) THEN
  DO jp = 1,kpatch
    CALL sso_init(ysso_np%AL(jp))
  ENDDO        
  DEALLOCATE(ysso_np%AL)
ELSE
  ALLOCATE(ysso_np%AL(kpatch))
  DO jp = 1,kpatch
    CALL sso_init(ysso_np%AL(jp))
  ENDDO
ENDIF
!
IF (lhook) CALL dr_hook("MODD_SSO_N:SSO_NP_INIT",1,zhook_handle)
END SUBROUTINE sso_np_init
!
END MODULE modd_sso_n