modd_seafluxn.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_seaflux_n
!     #################
!
!!****  *MODD_SEAFLUX_n - declaration of surface parameters for an inland water surface
!!
!!    PURPOSE
!!    -------
!     Declaration of surface parameters
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original      01/2004
!!      S. Senesi     01/2014  adapt to fractional seaice, and to seaice scheme
!!      S. Belamari   03/2014  Include NZ0
!
!*       0.   DECLARATIONS
!             ------------
!
USE modd_surf_par,   ONLY : xundef
USE modd_type_date_surf
!
USE modd_types_glt,   ONLY : t_glt
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE

TYPE seaflux_t
!
! General surface: 
!
  REAL, POINTER, DIMENSION(:)   :: XZS     ! orography
  REAL, POINTER, DIMENSION(:,:) :: XCOVER  ! fraction of each ecosystem       (-)
  LOGICAL, POINTER, DIMENSION(:):: LCOVER  ! GCOVER(i)=T --> ith cover field is not 0.
  LOGICAL                       :: LSBL    ! T: SBL scheme between sea and atm. forcing level
!                                          ! F: no atmospheric layers below forcing level      
  LOGICAL                       :: LHANDLE_SIC ! T: we do weight seaice and open sea fluxes
  CHARACTER(LEN=6)              :: CSEAICE_SCHEME! Name of the seaice scheme 
  REAL, POINTER, DIMENSION(:)   :: XSEABATHY   ! bathymetry
!
  LOGICAL                       :: LINTERPOL_SST ! Interpolation of monthly SST
  CHARACTER(LEN=6)              :: CINTERPOL_SST ! Interpolation method of monthly SST
  LOGICAL                       :: LINTERPOL_SSS ! Interpolation of monthly SSS
  CHARACTER(LEN=6)              :: CINTERPOL_SSS ! Interpolation method of monthly SSS
  LOGICAL                       :: LINTERPOL_SIC ! Interpolation of monthly SIC
  CHARACTER(LEN=6)              :: CINTERPOL_SIC ! Interpolation method of monthly SIC
  LOGICAL                       :: LINTERPOL_SIT ! Interpolation of monthly SIT
  CHARACTER(LEN=6)              :: CINTERPOL_SIT ! Interpolation method of monthly SIT
  REAL                          :: XFREEZING_SST ! Value marking frozen sea in SST data
  REAL                          :: XSIC_EFOLDING_TIME ! For damping of SIC (days)
  REAL                          :: XSIT_EFOLDING_TIME ! For damping of SIT (days)
  REAL                          :: XSEAICE_TSTEP ! Sea ice model time step
  REAL                          :: XCD_ICE_CST   ! Turbulent exchange coefficient for seaice
  REAL                          :: XSI_FLX_DRV   ! Derivative of fluxes on seaice w.r.t to the temperature (W m-2 K-1)
  
!
! Type of formulation for the fluxes
!
  CHARACTER(LEN=6)                  :: CSEA_FLUX   ! type of flux computation
  CHARACTER(LEN=4)                  :: CSEA_ALB    ! type of albedo
  LOGICAL                           :: LPWG        ! flag for gust
  LOGICAL                           :: LPRECIP     ! flag for precip correction
  LOGICAL                           :: LPWEBB      ! flag for Webb correction
  INTEGER                           :: NZ0         ! set to 0,1 or 2 according to Z0 formulation
                                                   ! 0= ARPEGE / 1= Smith (1988) / 2= Direct
  INTEGER                           :: NGRVWAVES   ! set to 0,1 or 2 according to the 
                                                   ! gravity waves model used in coare30_flux
  REAL                              :: XICHCE      ! CE coef calculation for ECUME
  LOGICAL                           :: LPERTFLUX   ! flag for stochastic flux perturbation
!
! Sea/Ocean:
!
  REAL, POINTER, DIMENSION(:) :: XSST    ! sea surface temperature
  REAL, POINTER, DIMENSION(:) :: XSSS    ! sea surface salinity
  REAL, POINTER, DIMENSION(:) :: XTICE   ! sea ice temperature
  REAL, POINTER, DIMENSION(:) :: XSIC    ! sea ice concentration ( constraint for seaice scheme )
  REAL, POINTER, DIMENSION(:) :: XSST_INI! initial sea surface temperature
  REAL, POINTER, DIMENSION(:) :: XZ0     ! roughness length
  REAL, POINTER, DIMENSION(:) :: XZ0H    ! roughness length for heat
  REAL, POINTER, DIMENSION(:) :: XEMIS   ! emissivity
  REAL, POINTER, DIMENSION(:) :: XDIR_ALB! direct albedo
  REAL, POINTER, DIMENSION(:) :: XSCA_ALB! diffuse albedo
  REAL, POINTER, DIMENSION(:) :: XICE_ALB! sea-ice albedo from seaice model (ESM or embedded)
  REAL, POINTER, DIMENSION(:) :: XUMER   ! U component of sea current (for ESM coupling)
  REAL, POINTER, DIMENSION(:) :: XVMER   ! V component of sea current (for ESM coupling)
!
  REAL, POINTER, DIMENSION(:,:) :: XSST_MTH! monthly sea surface temperature (precedent, current and next)
  REAL, POINTER, DIMENSION(:,:) :: XSSS_MTH! monthly sea surface salinity    (precedent, current and next)
  REAL, POINTER, DIMENSION(:,:) :: XSIC_MTH! monthly sea ice cover           (precedent, current and next)
  REAL, POINTER, DIMENSION(:,:) :: XSIT_MTH! monthly sea ice thickness       (precedent, current and next)
  REAL, POINTER, DIMENSION(:)   :: XFSIC   ! nudging (or forcing) sea ice cover
  REAL, POINTER, DIMENSION(:)   :: XFSIT   ! nudging sea ice thickness
!
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_WIND ! 10m wind speed for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_FWSU ! zonal wind stress for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_FWSV ! meridian wind stress for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_SNET ! Solar net heat flux
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_HEAT ! Non solar net heat flux
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_EVAP ! Evaporation for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_RAIN ! Rainfall for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_SNOW ! Snowfall for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEA_FWSM ! wind stress for ESM coupling
!  
  REAL, POINTER, DIMENSION(:) :: XCPL_SEAICE_SNET ! Solar net heat flux for ESM coupling
  REAL, POINTER, DIMENSION(:) :: XCPL_SEAICE_HEAT ! Non solar net heat flux
  REAL, POINTER, DIMENSION(:) :: XCPL_SEAICE_EVAP ! Ice sublimation for ESM coupling
!
  REAL, POINTER, DIMENSION(:) :: XPERTFLUX     ! Stochastic flux perturbation pattern
!
! Sea-ice :
!
  TYPE(t_glt)                           :: TGLT ! Sea-ice state , diagnostics and auxilliaries
                                             ! for the case of embedded Gelato Seaice model
!
! Date:
!
  TYPE (date_time)                      :: TTIME   ! current date and time
  TYPE (date_time)                      :: TZTIME  
  LOGICAL                               :: LTZTIME_DONE
  INTEGER                               :: JSX  
!
! Time-step:
!
  REAL                                  :: XTSTEP  ! time step
!
  REAL                                  :: XOUT_TSTEP  ! output writing time step
!
!
!
END TYPE seaflux_t



 CONTAINS

!




SUBROUTINE seaflux_init(YSEAFLUX)
TYPE(seaflux_t), INTENT(INOUT) :: yseaflux
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_SEAFLUX_N:SEAFLUX_INIT",0,zhook_handle)
  nullify(yseaflux%XZS)
  nullify(yseaflux%XCOVER)
  nullify(yseaflux%LCOVER)
  nullify(yseaflux%XSEABATHY)
  nullify(yseaflux%XSST)
  nullify(yseaflux%XSSS)
  nullify(yseaflux%XSIC)
  nullify(yseaflux%XTICE)
  nullify(yseaflux%XSST_INI)
  nullify(yseaflux%XZ0)
  nullify(yseaflux%XZ0H)
  nullify(yseaflux%XEMIS)
  nullify(yseaflux%XDIR_ALB)
  nullify(yseaflux%XSCA_ALB)
  nullify(yseaflux%XICE_ALB)
  nullify(yseaflux%XUMER)
  nullify(yseaflux%XVMER)
  nullify(yseaflux%XSST_MTH)
  nullify(yseaflux%XSSS_MTH)
  nullify(yseaflux%XSIC_MTH)
  nullify(yseaflux%XSIT_MTH)
  nullify(yseaflux%XFSIC)
  nullify(yseaflux%XFSIT)
  nullify(yseaflux%XCPL_SEA_WIND)
  nullify(yseaflux%XCPL_SEA_FWSU)
  nullify(yseaflux%XCPL_SEA_FWSV)
  nullify(yseaflux%XCPL_SEA_SNET)
  nullify(yseaflux%XCPL_SEA_HEAT)
  nullify(yseaflux%XCPL_SEA_EVAP)
  nullify(yseaflux%XCPL_SEA_RAIN)
  nullify(yseaflux%XCPL_SEA_SNOW)
  nullify(yseaflux%XCPL_SEA_FWSM)
  nullify(yseaflux%XCPL_SEAICE_SNET)
  nullify(yseaflux%XCPL_SEAICE_HEAT)
  nullify(yseaflux%XCPL_SEAICE_EVAP)
  nullify(yseaflux%XPERTFLUX)
yseaflux%LSBL=.false.
yseaflux%LHANDLE_SIC=.false.
yseaflux%CSEAICE_SCHEME='NONE  '
yseaflux%LINTERPOL_SST=.false.
yseaflux%CINTERPOL_SST=' '
yseaflux%LINTERPOL_SSS=.false.
yseaflux%CINTERPOL_SSS=' '
yseaflux%LINTERPOL_SIC=.false.
yseaflux%CINTERPOL_SIC=' '
yseaflux%LINTERPOL_SIT=.false.
yseaflux%CINTERPOL_SIT=' '
yseaflux%XFREEZING_SST=-1.8
yseaflux%XSIC_EFOLDING_TIME=0. ! means : no damping
yseaflux%XSIT_EFOLDING_TIME=0. ! means : no damping
yseaflux%XSEAICE_TSTEP=xundef 
yseaflux%XCD_ICE_CST=0.
yseaflux%XSI_FLX_DRV=-20. 
yseaflux%CSEA_FLUX=' '
yseaflux%CSEA_ALB=' '
yseaflux%LPWG=.false.
yseaflux%LPRECIP=.false.
yseaflux%LPWEBB=.false.
yseaflux%NZ0=0
yseaflux%NGRVWAVES=0
yseaflux%XICHCE=0.
yseaflux%LPERTFLUX=.false.
yseaflux%JSX=0
yseaflux%LTZTIME_DONE = .false.
yseaflux%XTSTEP=0.
yseaflux%XOUT_TSTEP=0.
IF (lhook) CALL dr_hook("MODD_SEAFLUX_N:SEAFLUX_INIT",1,zhook_handle)
END SUBROUTINE seaflux_init


END MODULE modd_seaflux_n