modn_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 modn_seaflux_n
!     ##################
!
!!****  *MODN_SEAFLUX_n* - declaration of namelist NAM_SEAFLUXn
!!
!!    PURPOSE
!!    -------
!       The purpose of this module is to specify  the namelist NAM_SEAFLUX_n
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!       
!!    AUTHOR
!!    ------
!!      V. Masson    *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004                    
!!      Modified    01/2006 : sea flux parameterization.
!!      Modified    08/2009 : LSURF_BUDGETC
!!      Modified    01/2014 : S. Senesi : introduce sea-ice model
!!      Modified    03/2014 : S. Belamari - add NZ0 (to choose PZ0SEA formulation)
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!

!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
 CHARACTER(LEN=6)  :: CSEA_FLUX
 CHARACTER(LEN=4)  :: CSEA_ALB
REAL     :: XTSTEP
REAL     :: XOUT_TSTEP
REAL     :: XDIAG_TSTEP
INTEGER  :: N2M
LOGICAL  :: L2M_MIN_ZS         
LOGICAL  :: LSURF_BUDGET
LOGICAL  :: LRAD_BUDGET
LOGICAL  :: LSURF_BUDGETC
LOGICAL  :: LRESET_BUDGETC
LOGICAL  :: LCOEF
LOGICAL  :: LSURF_VARS
LOGICAL  :: LPWG
LOGICAL  :: LPRECIP
LOGICAL  :: LPWEBB
LOGICAL  :: LDIAG_OCEAN
LOGICAL  :: LDIAG_SEAICE
INTEGER  :: NZ0
INTEGER  :: NGRVWAVES
REAL     :: XICHCE
 CHARACTER(LEN=6)  :: CCH_DRY_DEP
LOGICAL  :: LPROGSST
LOGICAL  :: LPERTFLUX              ! True = stochastic flux perturbation (default:False)
INTEGER  :: NTIME_COUPLING
REAL :: XOCEAN_TSTEP
 CHARACTER(LEN=6)  :: CINTERPOL_SST
 CHARACTER(LEN=6)  :: CINTERPOL_SSS
 CHARACTER(LEN=6)  :: CINTERPOL_SIC
 CHARACTER(LEN=6)  :: CINTERPOL_SIT
REAL     :: XFREEZING_SST
REAL     :: XSIC_EFOLDING_TIME
REAL     :: XSIT_EFOLDING_TIME
REAL     :: XSEAICE_TSTEP
REAL     :: XCD_ICE_CST
REAL     :: XSI_FLX_DRV
!
namelist/nam_seafluxn/csea_flux,csea_alb, lpwg, lprecip, lpwebb, ngrvwaves, &
                      nz0, lprogsst, ntime_coupling, xocean_tstep, xichce, &
                      cinterpol_sst, cinterpol_sss, lpertflux
namelist/nam_diag_surfn/n2m,l2m_min_zs,lsurf_budget,lrad_budget, &
                          lsurf_budgetc,lreset_budgetc,lcoef,lsurf_vars  
namelist/nam_ch_seafluxn/cch_dry_dep
namelist/nam_diag_oceann/ldiag_ocean
namelist/nam_seaicen/ldiag_seaice, cinterpol_sic, cinterpol_sit,           &
                     xfreezing_sst, xsic_efolding_time, xsit_efolding_time,&
                     xseaice_tstep, xcd_ice_cst, xsi_flx_drv
!
 CONTAINS
!
SUBROUTINE init_nam_seafluxn (O, S)
!
  USE modd_ocean_n, ONLY : ocean_t
  USE modd_seaflux_n, ONLY : seaflux_t
!
  IMPLICIT NONE

!
  TYPE(ocean_t), INTENT(INOUT) :: o
  TYPE(seaflux_t), INTENT(INOUT) :: s
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_SEAFLUXN',0,zhook_handle)
  xtstep = s%XTSTEP
  xout_tstep = s%XOUT_TSTEP
  csea_flux = s%CSEA_FLUX
  csea_alb = s%CSEA_ALB
  lpwg = s%LPWG
  lprecip = s%LPRECIP
  cinterpol_sst = s%CINTERPOL_SST
  cinterpol_sss = s%CINTERPOL_SSS
  lpwebb = s%LPWEBB
  nz0 = s%NZ0
  ngrvwaves = s%NGRVWAVES
  lprogsst = o%LPROGSST
  ntime_coupling = o%NTIME_COUPLING
  xocean_tstep = o%XOCEAN_TSTEP
  xichce = s%XICHCE
  lpertflux = s%LPERTFLUX
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_SEAFLUXN',1,zhook_handle)
END SUBROUTINE init_nam_seafluxn

SUBROUTINE update_nam_seafluxn (O, S)
!
  USE modd_ocean_n, ONLY : ocean_t
  USE modd_seaflux_n, ONLY : seaflux_t
!
  IMPLICIT NONE

!
  TYPE(ocean_t), INTENT(INOUT) :: o
  TYPE(seaflux_t), INTENT(INOUT) :: s
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_SEAFLUXN',0,zhook_handle)
  s%XTSTEP = xtstep
  s%XOUT_TSTEP = xout_tstep  
  s%CSEA_FLUX = csea_flux
  s%CSEA_ALB = csea_alb
  s%LPWG = lpwg
  s%LPRECIP = lprecip
  s%CINTERPOL_SST = cinterpol_sst
  s%CINTERPOL_SSS = cinterpol_sss
  s%LPWEBB = lpwebb
  s%NZ0 = nz0
  s%NGRVWAVES = ngrvwaves
  o%LPROGSST = lprogsst
  o%NTIME_COUPLING = ntime_coupling
  o%XOCEAN_TSTEP = xocean_tstep
  s%XICHCE = xichce
  s%LPERTFLUX = lpertflux
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_SEAFLUXN',1,zhook_handle)
END SUBROUTINE update_nam_seafluxn
!
SUBROUTINE init_nam_diag_surfn (DGS)
!
  USE modd_diag_seaflux_n, ONLY : diag_seaflux_t
!
  IMPLICIT NONE

!
  TYPE(diag_seaflux_t), INTENT(INOUT) :: dgs
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_DIAG_SURFN',0,zhook_handle)
  xdiag_tstep = dgs%XDIAG_TSTEP  
  n2m = dgs%N2M
  l2m_min_zs = dgs%L2M_MIN_ZS
  lsurf_budget = dgs%LSURF_BUDGET
  lrad_budget = dgs%LRAD_BUDGET
  lsurf_budgetc = dgs%LSURF_BUDGETC
  lreset_budgetc = dgs%LRESET_BUDGETC  
  lcoef = dgs%LCOEF
  lsurf_vars = dgs%LSURF_VARS
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_DIAG_SURFN',1,zhook_handle)
END SUBROUTINE init_nam_diag_surfn

SUBROUTINE update_nam_diag_surfn (DGS)
!
  USE modd_diag_seaflux_n, ONLY : diag_seaflux_t
!
  IMPLICIT NONE

!
  TYPE(diag_seaflux_t), INTENT(INOUT) :: dgs
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_DIAG_SURFN',0,zhook_handle)
  dgs%XDIAG_TSTEP = xdiag_tstep  
  dgs%N2M = n2m
  dgs%L2M_MIN_ZS = l2m_min_zs
  dgs%LSURF_BUDGET = lsurf_budget
  dgs%LRAD_BUDGET = lrad_budget 
  dgs%LSURF_BUDGETC = lsurf_budgetc
  dgs%LRESET_BUDGETC = lreset_budgetc  
  dgs%LCOEF = lcoef
  dgs%LSURF_VARS = lsurf_vars
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_DIAG_SURFN',1,zhook_handle)
END SUBROUTINE update_nam_diag_surfn
!
SUBROUTINE init_nam_ch_seafluxn (CHS)
!
  USE modd_ch_seaflux_n, ONLY : ch_seaflux_t
!
  IMPLICIT NONE

!
  TYPE(ch_seaflux_t), INTENT(INOUT) :: chs
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_CH_SEAFLUXN',0,zhook_handle)
  cch_dry_dep = chs%CCH_DRY_DEP
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_CH_SEAFLUXN',1,zhook_handle)
END SUBROUTINE init_nam_ch_seafluxn

SUBROUTINE update_nam_ch_seafluxn (CHS)
!
  USE modd_ch_seaflux_n, ONLY : ch_seaflux_t
!
  IMPLICIT NONE

!
  TYPE(ch_seaflux_t), INTENT(INOUT) :: chs
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_CH_SEAFLUXN',0,zhook_handle)
  chs%CCH_DRY_DEP = cch_dry_dep
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_CH_SEAFLUXN',1,zhook_handle)
END SUBROUTINE update_nam_ch_seafluxn

SUBROUTINE init_nam_diag_oceann (DGO)
!
  USE modd_diag_ocean_n, ONLY : diag_ocean_t
!
  IMPLICIT NONE

!
  TYPE(diag_ocean_t), INTENT(INOUT) :: dgo
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_DIAG_OCEANN',0,zhook_handle)
  ldiag_ocean = dgo%LDIAG_OCEAN
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_DIAG_OCEANN',1,zhook_handle)
END SUBROUTINE init_nam_diag_oceann

SUBROUTINE update_nam_diag_oceann (DGO)
!
  USE modd_diag_ocean_n, ONLY : diag_ocean_t
!
  IMPLICIT NONE

!
  TYPE(diag_ocean_t), INTENT(INOUT) :: dgo
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_DIAG_OCEANN',0,zhook_handle)
  dgo%LDIAG_OCEAN = ldiag_ocean
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_DIAG_OCEANN',1,zhook_handle)
END SUBROUTINE update_nam_diag_oceann

SUBROUTINE init_nam_seaicen (DGSI, S)
!
  USE modd_diag_seaice_n, ONLY : diag_seaice_t
  USE modd_seaflux_n, ONLY : seaflux_t
!
  IMPLICIT NONE

!
  TYPE(diag_seaice_t), INTENT(INOUT) :: dgsi
  TYPE(seaflux_t), INTENT(INOUT) :: s
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_SEAICEN',0,zhook_handle)
  ldiag_seaice = dgsi%LDIAG_SEAICE
  cinterpol_sic = s%CINTERPOL_SIC
  cinterpol_sit = s%CINTERPOL_SIT
  xsic_efolding_time=s%XSIC_EFOLDING_TIME
  xsit_efolding_time=s%XSIT_EFOLDING_TIME
  xseaice_tstep=s%XSEAICE_TSTEP
  xfreezing_sst = s%XFREEZING_SST
  xcd_ice_cst = s%XCD_ICE_CST
  xsi_flx_drv = s%XSI_FLX_DRV
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_SEAICEN',1,zhook_handle)
END SUBROUTINE init_nam_seaicen

SUBROUTINE update_nam_seaicen (DGSI, S)
!
  USE modd_diag_seaice_n, ONLY : diag_seaice_t
  USE modd_seaflux_n, ONLY : seaflux_t
!
  IMPLICIT NONE

!
  TYPE(diag_seaice_t), INTENT(INOUT) :: dgsi
  TYPE(seaflux_t), INTENT(INOUT) :: s
  REAL(KIND=JPRB) :: zhook_handle

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_SEAICEN',0,zhook_handle)
  dgsi%LDIAG_SEAICE = ldiag_seaice
  s%CINTERPOL_SIC = cinterpol_sic
  s%CINTERPOL_SIT = cinterpol_sit
  s%XSIC_EFOLDING_TIME = xsic_efolding_time
  s%XSIT_EFOLDING_TIME = xsit_efolding_time
  s%XSEAICE_TSTEP = xseaice_tstep
  s%XFREEZING_SST = xfreezing_sst
  s%XCD_ICE_CST = xcd_ice_cst
  s%XSI_FLX_DRV = xsi_flx_drv
IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_SEAICEN',1,zhook_handle)
END SUBROUTINE update_nam_seaicen

END MODULE modn_seaflux_n