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