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_misc_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_misc_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 (DGO)
!
  USE modd_diag_n, ONLY : diag_options_t
!
  IMPLICIT NONE

!
  TYPE(diag_options_t), INTENT(INOUT) :: DGO
  REAL(KIND=JPRB) :: ZHOOK_HANDLE

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:INIT_NAM_DIAG_SURFN',0,zhook_handle)
  xdiag_tstep = dgo%XDIAG_TSTEP  
  n2m = dgo%N2M
  l2m_min_zs = dgo%L2M_MIN_ZS
  lsurf_budget = dgo%LSURF_BUDGET
  lrad_budget = dgo%LRAD_BUDGET
  lsurf_budgetc = dgo%LSURF_BUDGETC
  lreset_budgetc = dgo%LRESET_BUDGETC  
  lcoef = dgo%LCOEF
  lsurf_vars = dgo%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 (DGO)
!
  USE modd_diag_n, ONLY : diag_options_t
!
  IMPLICIT NONE

!
  TYPE(diag_options_t), INTENT(INOUT) :: DGO
  REAL(KIND=JPRB) :: ZHOOK_HANDLE

  IF (lhook) CALL dr_hook('MODN_SEAFLUX_N:UPDATE_NAM_DIAG_SURFN',0,zhook_handle)
  dgo%XDIAG_TSTEP = xdiag_tstep  
  dgo%N2M = n2m
  dgo%L2M_MIN_ZS = l2m_min_zs
  dgo%LSURF_BUDGET = lsurf_budget
  dgo%LRAD_BUDGET = lrad_budget 
  dgo%LSURF_BUDGETC = lsurf_budgetc
  dgo%LRESET_BUDGETC = lreset_budgetc  
  dgo%LCOEF = lcoef
  dgo%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 (DGMSI, S)
!
  USE modd_diag_misc_seaice_n, ONLY : diag_misc_seaice_t
  USE modd_seaflux_n, ONLY : seaflux_t
!
  IMPLICIT NONE

!
  TYPE(diag_misc_seaice_t), INTENT(INOUT) :: DGMSI
  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_misc_seaice = dgmsi%LDIAG_MISC_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 (DGMSI, S)
!
  USE modd_diag_misc_seaice_n, ONLY : diag_misc_seaice_t
  USE modd_seaflux_n, ONLY : seaflux_t
!
  IMPLICIT NONE
!
  TYPE(diag_misc_seaice_t), INTENT(INOUT) :: DGMSI
  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)
  dgmsi%LDIAG_MISC_SEAICE = ldiag_misc_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