modd_surf_atm_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_surf_atm_sso_n
!     ######################
!
!!****  *MODD_SURF_ATM_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

TYPE surf_atm_sso_t
!
!-----------------------------------------------------------------------------------------------------
!
! Type of roughness
!
 CHARACTER(LEN=4) :: CROUGH     ! type of orographic roughness
!                              ! 'NONE'
                               ! 'Z01D'
                               ! 'Z04D'
                               ! 'BE04'

!-----------------------------------------------------------------------------------------------------
!
! 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 :: 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)
!
  REAL, DIMENSION(:), POINTER   :: XZ0EFFJPDIR    ! heading of J direction (deg from N clockwise)

  REAL, DIMENSION(:), POINTER   :: XZ0REL         ! relief roughness length                 (m)
!
  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         
!-----------------------------------------------------------------------------------------------------
!
!


END TYPE surf_atm_sso_t



 CONTAINS

!




SUBROUTINE surf_atm_sso_init(YSURF_ATM_SSO)
TYPE(surf_atm_sso_t), INTENT(INOUT) :: ysurf_atm_sso
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_SURF_ATM_SSO_N:SURF_ATM_SSO_INIT",0,zhook_handle)
  nullify(ysurf_atm_sso%XAOSIP)
  nullify(ysurf_atm_sso%XAOSIM)
  nullify(ysurf_atm_sso%XAOSJP)
  nullify(ysurf_atm_sso%XAOSJM)
  nullify(ysurf_atm_sso%XHO2IP)
  nullify(ysurf_atm_sso%XHO2IM)
  nullify(ysurf_atm_sso%XHO2JP)
  nullify(ysurf_atm_sso%XHO2JM)
  nullify(ysurf_atm_sso%XZ0EFFIP)
  nullify(ysurf_atm_sso%XZ0EFFIM)
  nullify(ysurf_atm_sso%XZ0EFFJP)
  nullify(ysurf_atm_sso%XZ0EFFJM)
  nullify(ysurf_atm_sso%XZ0EFFJPDIR)  
  nullify(ysurf_atm_sso%XZ0REL)
  nullify(ysurf_atm_sso%XSSO_SLOPE)
  nullify(ysurf_atm_sso%XSSO_ANIS)
  nullify(ysurf_atm_sso%XSSO_DIR)
  nullify(ysurf_atm_sso%XSSO_STDEV)
  nullify(ysurf_atm_sso%XAVG_ZS)
  nullify(ysurf_atm_sso%XSIL_ZS)
  nullify(ysurf_atm_sso%XMAX_ZS)
  nullify(ysurf_atm_sso%XMIN_ZS)
ysurf_atm_sso%CROUGH=' '
ysurf_atm_sso%XFRACZ0=2.
ysurf_atm_sso%XCOEFBE=2.
IF (lhook) CALL dr_hook("MODD_SURF_ATM_SSO_N:SURF_ATM_SSO_INIT",1,zhook_handle)
END SUBROUTINE surf_atm_sso_init


END MODULE modd_surf_atm_sso_n