modd_flaken.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_flake_n
!     ####################
!
!!****  *MODD_FLAKE_n - declaration of surface parameters for the FLake model 
!!                      for inland water surfaces
!!
!!    PURPOSE
!!    -------
!     Declaration of surface parameters
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      V. Masson  *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 flake_t 
!
!-------------------------------------------------------------------------------------
! General surface: 
!-------------------------------------------------------------------------------------
!
  REAL, POINTER, DIMENSION(:) :: XZS       ! orography                     (m)
  REAL, POINTER, DIMENSION(:) :: XZ0       ! roughness length              (m)
  REAL, POINTER, DIMENSION(:) :: XUSTAR    ! air friction velocity         (m/s)
  REAL, POINTER, DIMENSION(:) :: XEMIS     ! water surface emissivity (NOT USED BY FLAKE)
!
  REAL, POINTER, DIMENSION(:,:) :: XCOVER  ! fraction of each ecosystem    (-)
!                                          ! F: no atmospheric layers below forcing level  
!
  LOGICAL, POINTER, DIMENSION(:) :: LCOVER ! GCOVER(i)=T --> ith cover field is not 0.
  LOGICAL                        :: LSBL   ! T: SBL scheme within the Surface Boundary Layer
!
!-------------------------------------------------------------------------------------
! Date and time:
!-------------------------------------------------------------------------------------
!
  TYPE (date_time)                  :: TTIME         ! current date and time
!
  REAL                              :: XTSTEP        ! time step
!
  REAL                              :: XOUT_TSTEP    ! output writing time step
!
!-------------------------------------------------------------------------------------
! FLake switches
!-------------------------------------------------------------------------------------
!
  LOGICAL            :: LSEDIMENTS  ! flag to use or not the bottom sediments
  LOGICAL            :: LSKINTEMP   ! flag to use or not the skin temperature computation
  CHARACTER(LEN=3)   :: CSNOW_FLK   ! FLake snow scheme
  CHARACTER(LEN=5)   :: CFLK_FLUX   ! Type of flux computation
  CHARACTER(LEN=4)   :: CFLK_ALB    ! Type of albedo
!
!-------------------------------------------------------------------------------------
! FLake parameters and variables
!-------------------------------------------------------------------------------------
!
  REAL, POINTER, DIMENSION(:) :: XWATER_DEPTH  ! Lake depth (m)
  REAL, POINTER, DIMENSION(:) :: XWATER_FETCH  ! Lake fetch (m)
  REAL, POINTER, DIMENSION(:) :: XT_BS         ! Temperature at the outer edge of the thermally 
                                               !       active layer of the bottom sediments [K]
  REAL, POINTER, DIMENSION(:) :: XDEPTH_BS     ! Depth of the thermally active layer of the
                                               !       bottom sediments [m]
  REAL, POINTER, DIMENSION(:) :: XCORIO        ! The Coriolis parameter [s^{-1}]
  REAL, POINTER, DIMENSION(:) :: XDIR_ALB      ! Water surface direct albedo
  REAL, POINTER, DIMENSION(:) :: XSCA_ALB      ! Water surface diffuse albedo
  REAL, POINTER, DIMENSION(:) :: XICE_ALB      ! Ice surface albedo (for ESM coupling)
  REAL, POINTER, DIMENSION(:) :: XSNOW_ALB     ! Snow surface albedo
  REAL, POINTER, DIMENSION(:) :: XEXTCOEF_WATER ! Extinction coefficient for the water [m^{-1}]
  REAL, POINTER, DIMENSION(:) :: XEXTCOEF_ICE   ! Extinction coefficient for the ice [m^{-1}]
  REAL, POINTER, DIMENSION(:) :: XEXTCOEF_SNOW  ! Extinction coefficient for the snow [m^{-1}] 
  REAL, POINTER, DIMENSION(:) :: XT_SNOW       ! Temperature at the air-snow interface [K]    
  REAL, POINTER, DIMENSION(:) :: XT_ICE        ! Temperature at the snow-ice or air-ice 
                                               !        interface [K]
  REAL, POINTER, DIMENSION(:) :: XT_MNW        ! Mean temperature of the water column [K]
  REAL, POINTER, DIMENSION(:) :: XT_WML        ! Mixed-layer temperature [K]
  REAL, POINTER, DIMENSION(:) :: XT_BOT        ! Temperature at the water-bottom sediment 
                                               !        interface [K]
  REAL, POINTER, DIMENSION(:) :: XT_B1         ! Temperature at the bottom of the upper 
                                               !        layer of the sediments [K]
  REAL, POINTER, DIMENSION(:) :: XCT           ! Shape factor (thermocline)
  REAL, POINTER, DIMENSION(:) :: XH_SNOW       ! Snow thickness [m]
  REAL, POINTER, DIMENSION(:) :: XH_ICE        ! Ice thickness [m]
  REAL, POINTER, DIMENSION(:) :: XH_ML         ! Thickness of the mixed-layer [m]
  REAL, POINTER, DIMENSION(:) :: XH_B1         ! Thickness of the upper layer of bottom sediments [m]                                    
!
  REAL, POINTER, DIMENSION(:) :: XTS  ! surface temperature  (K)
                                      ! (water or ice or snow)
!
!-------------------------------------------------------------------------------------
! Coupling field for Earth system model
!-------------------------------------------------------------------------------------
!
  REAL, POINTER, DIMENSION(:) :: XCPL_FLAKE_EVAP ! Evaporation for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_FLAKE_RAIN ! Rainfall for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_FLAKE_SNOW ! Snowfall for ESM coupling
!
END TYPE flake_t
!
!-------------------------------------------------------------------------------------
!


 CONTAINS
!
!-------------------------------------------------------------------------------------
!
!
!


!
!
!-------------------------------------------------------------------------------------
!
SUBROUTINE flake_init(YFLAKE)
TYPE(flake_t), INTENT(INOUT) :: yflake
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_FLAKE_N:FLAKE_INIT",0,zhook_handle)
  nullify(yflake%XZS)
  nullify(yflake%XZ0)
  nullify(yflake%XUSTAR)
  nullify(yflake%XCOVER)
  nullify(yflake%LCOVER)
  nullify(yflake%XEMIS)
  nullify(yflake%XWATER_DEPTH)
  nullify(yflake%XWATER_FETCH)
  nullify(yflake%XT_BS)
  nullify(yflake%XDEPTH_BS)
  nullify(yflake%XCORIO)
  nullify(yflake%XDIR_ALB)
  nullify(yflake%XSCA_ALB)
  nullify(yflake%XICE_ALB)
  nullify(yflake%XSNOW_ALB)
  nullify(yflake%XEXTCOEF_WATER)
  nullify(yflake%XEXTCOEF_ICE)
  nullify(yflake%XEXTCOEF_SNOW)
  nullify(yflake%XT_SNOW)
  nullify(yflake%XT_ICE)
  nullify(yflake%XT_MNW)
  nullify(yflake%XT_WML)
  nullify(yflake%XT_BOT)
  nullify(yflake%XT_B1)
  nullify(yflake%XCT)
  nullify(yflake%XH_SNOW)
  nullify(yflake%XH_ICE)
  nullify(yflake%XH_ML)
  nullify(yflake%XH_B1)
  nullify(yflake%XTS)
  nullify(yflake%XCPL_FLAKE_EVAP)
  nullify(yflake%XCPL_FLAKE_RAIN)
  nullify(yflake%XCPL_FLAKE_SNOW)
yflake%LSBL=.false.
yflake%XTSTEP=0.
yflake%XOUT_TSTEP=0.
yflake%LSEDIMENTS=.false.
yflake%LSKINTEMP=.false.
yflake%CSNOW_FLK='   '
yflake%CFLK_ALB='    '
yflake%CFLK_FLUX='     '
IF (lhook) CALL dr_hook("MODD_FLAKE_N:FLAKE_INIT",1,zhook_handle)
END SUBROUTINE flake_init
!
!-------------------------------------------------------------------------------------
!
!
!-------------------------------------------------------------------------------------
!
END MODULE modd_flake_n