modd_surf_atmn.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_n
!     ######################
!
!!****  *MODD_SURF_ATM - declaration of surface parameters
!!
!!    PURPOSE
!!    -------
!     Declaration of surface parameters
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      V. Masson and A. Boone   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       01/2004
!
!*       0.   DECLARATIONS
!             ------------
!
USE modd_type_date_surf
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE

TYPE surf_atm_t
!
!-----------------------------------------------------------------------------------------------------
!
! Type of each surface scheme
!
  CHARACTER(LEN=6) :: ctown      ! name of the urban surface scheme
  CHARACTER(LEN=6) :: cnature    ! name of the soil&vegetation surface scheme
  CHARACTER(LEN=6) :: cwater     ! name of the scheme for inland water
  CHARACTER(LEN=6) :: csea       ! name for the ocean scheme
!
!-----------------------------------------------------------------------------------------------------
!
! Surface/Tile Fractions:
!
  REAL, POINTER, DIMENSION(:)   :: xtown     ! urban surface fraction of the grid box   (-)
  REAL, POINTER, DIMENSION(:)   :: xnature   ! natural surface fraction of the grid box (-)
  REAL, POINTER, DIMENSION(:)   :: xwater    ! inland water fraction of the grid box    (-)
  REAL, POINTER, DIMENSION(:)   :: xsea      ! sea/ocean fraction of the grid box       (-)
!
!-------------------------------------------------------------------------------
!
! type of initialization of vegetation: from cover types (ecoclimap) or parameters prescribed
!
  LOGICAL                        :: lecoclimap ! T: parameters computed from ecoclimap
!                                              ! F: they are read in the file
!
  LOGICAL                        :: lecosg     ! T: parameters computed from ecosg
!-------------------------------------------------------------------------------
!
! change water (not lake) to nature and/or town to rock : arrange cover properly
!
  LOGICAL                        :: lwater_to_nature ! T: Change Wetland treated as inland water into nature 
  LOGICAL                        :: ltown_to_rock    ! T: Change Town into Rock
!
!-------------------------------------------------------------------------------
!
! include urban green areas for urbanized covers
!
  LOGICAL                        :: lgarden    ! T: define urban green areas
!                                              ! F: no urban green areas
!
!-----------------------------------------------------------------------------------------------------
!
! Masks and number of grid elements for each tile surface
!
! Sea/Ocean:
!
  INTEGER                               :: nsize_sea    ! number of grid points by proc containing a
!                                                     ! sea surface                              (-)
  INTEGER                               :: ndim_sea     ! total number of grid points containing a
!                                                     ! sea surface                             (-)
  INTEGER, POINTER, DIMENSION(:)    :: nr_sea       ! sea/ocean surface mask                  (-)
!
! Inland Water:
!
  INTEGER                               :: nsize_water  ! number of grid points containing an 
!                                                     ! inland water surface                    (-)
  INTEGER                               :: ndim_water   ! total number of grid points by proc containing an
!                                                     ! inland surface
  INTEGER, POINTER, DIMENSION(:)    :: nr_water
!
! Town:
!
  INTEGER                               :: nsize_town   ! number of grid points by proc containing an 
!                                                     ! urban surface                           (-)
  INTEGER                               :: ndim_town    ! total number of grid points containing an
!                                                     ! urban surface
  INTEGER, POINTER, DIMENSION(:)    :: nr_town      ! urban surface mask                      (-)
!
! Natural surface:
!
  INTEGER                               :: nsize_nature ! number of grid points by proc containing a 
!                                                     ! natural surface                         (-)
  INTEGER                               :: ndim_nature  ! total number of grid points containing a
!                                                     ! natural surface                         (-)
  INTEGER, POINTER, DIMENSION(:)    :: nr_nature    ! natural surface mask                    (-)
!
! All surfaces:
!
  INTEGER                               :: nsize_full   ! total number of grid points by proc     (-)
  INTEGER                               :: ndim_full    ! total number of grid points             (-)
!
! global sizes of the physical domain
!
  INTEGER                               :: nimax_surf_ll
  INTEGER                               :: njmax_surf_ll
!
! local sizes of the physical domain
!
  INTEGER                               :: nimax_surf_loc
  INTEGER                               :: njmax_surf_loc
!
!-----------------------------------------------------------------------------------------------------
!
! Surface fields (only 1 horizontal dimension)
!
  REAL, POINTER, DIMENSION(:,:) :: xcover    ! fraction of each ecosystem for each grid box (-)
  LOGICAL, POINTER, DIMENSION(:):: lcover    ! GCOVER(i)=T --> ith cover field is not 0.
  REAL, POINTER, DIMENSION(:)   :: xzs       ! orography                                    (m)
!
!-------------------------------------------------------------------------------
!
  type(date_time)                      :: ttime            ! current date and time
!
  REAL                                  :: xout_tstep       ! output writing time step
!
!-----------------------------------------------------------------------------------------------------
!
! physical fields need into the restart file for ARPEGE/ALADIN run
!
  REAL, POINTER, DIMENSION(:)   :: xrain    ! Rainfall rate at surface               (kg/m2/s)
  REAL, POINTER, DIMENSION(:)   :: xsnow    ! snowfall rate at surface               (kg/m2/s)
  REAL, POINTER, DIMENSION(:)   :: xz0      ! surface roughness length for momentum  (m)
  REAL, POINTER, DIMENSION(:)   :: xz0h     ! surface roughness length for heat      (m)
  REAL, POINTER, DIMENSION(:)   :: xqsurf   ! specific humidity at surface           (kg/kg)
!
!-----------------------------------------------------------------------------------------------------
!
!
END TYPE surf_atm_t
!


CONTAINS


!


!
!
SUBROUTINE surf_atm_init(YSURF_ATM)
TYPE(surf_atm_t), INTENT(INOUT) :: YSURF_ATM
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_SURF_ATM_N:SURF_ATM_INIT",0,zhook_handle)
  NULLIFY(ysurf_atm%XTOWN)
  NULLIFY(ysurf_atm%XNATURE)
  NULLIFY(ysurf_atm%XWATER)
  NULLIFY(ysurf_atm%XSEA)
  NULLIFY(ysurf_atm%NR_SEA)
  NULLIFY(ysurf_atm%NR_WATER)
  NULLIFY(ysurf_atm%NR_TOWN)
  NULLIFY(ysurf_atm%NR_NATURE)
  NULLIFY(ysurf_atm%XCOVER)
  NULLIFY(ysurf_atm%LCOVER)
  NULLIFY(ysurf_atm%XZS)
  NULLIFY(ysurf_atm%XRAIN)
  NULLIFY(ysurf_atm%XSNOW)
  NULLIFY(ysurf_atm%XZ0)
  NULLIFY(ysurf_atm%XZ0H)
  NULLIFY(ysurf_atm%XQSURF)
ysurf_atm%CTOWN=' '
ysurf_atm%CNATURE=' '
ysurf_atm%CWATER=' '
ysurf_atm%CSEA=' '
ysurf_atm%LECOCLIMAP=.false.
ysurf_atm%LECOSG=.false.
ysurf_atm%LWATER_TO_NATURE=.false.
ysurf_atm%LTOWN_TO_ROCK=.false.
ysurf_atm%LGARDEN=.false.
ysurf_atm%NSIZE_SEA=0
ysurf_atm%NDIM_SEA=0
ysurf_atm%NSIZE_WATER=0
ysurf_atm%NDIM_WATER=0
ysurf_atm%NSIZE_TOWN=0
ysurf_atm%NDIM_TOWN=0
ysurf_atm%NSIZE_NATURE=0
ysurf_atm%NDIM_NATURE=0
ysurf_atm%NSIZE_FULL=0
ysurf_atm%NDIM_FULL=0
ysurf_atm%NIMAX_SURF_ll=0
ysurf_atm%NJMAX_SURF_ll=0
ysurf_atm%NIMAX_SURF_LOC=0
ysurf_atm%NJMAX_SURF_LOC=0
ysurf_atm%XOUT_TSTEP=0.
IF (lhook) CALL dr_hook("MODD_SURF_ATM_N:SURF_ATM_INIT",1,zhook_handle)
END SUBROUTINE surf_atm_init



END MODULE modd_surf_atm_n