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