eORCA025.L121-OPM017

Summary

Compare to eORCA025.L121-OPM016, we switched off the AABW restoring, decreased the diffusion as well as the GM coefficient (get back to UKMO value)

Namelist

Only the change compare to the reference (eORCA025.L121-OPM016) a mentioned here:

namelist_oce

namtsd

  • tradmp input file changed to be equal to initial condition as in OPM006 because AABW restoring turned off.

!-----------------------------------------------------------------------
&namtsd_drk    !    Temperature & Salinity Data  (init/dmp)             (default: OFF)
!              !   if key_drakkar, **only**  namtsd_drk is read
!-----------------------------------------------------------------------
   ln_tsd_init   = .true.   !  Initialisation of ocean T & S with T &S input data (T) or not (F)
   ln_tsd_dmp = .false.     !  damping of ocean T & S toward T &S input data (T) or not (F)

   cn_dir        = './'     !  root directory for the location of the temperature and salinity file
   !___________!_____________________________________!___________________!___________!_____________!________!___________!_____________!__________!_______________!
   !           !  file name                          ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask !
   !           !                                     !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !   filename  ! pairing  !    filename   !
   ! data used for initial condition (istate)
   sn_tem_ini  = 'eORCA025.L121_WOA2018_c3.0_d1.0_v19812010.5.1', -1., 'votemper',  .true.   , .true. , 'yearly'  ,  ''         , ' '      , ' '
   sn_sal_ini  = 'eORCA025.L121_WOA2018_c3.0_d1.0_v19812010.5.1', -1., 'vosaline',  .true.   , .true. , 'yearly'  ,  ''         , ' '      , ' '
   ! data used for damping ( tradmp)
   sn_tem_dmp  = 'eORCA025.L121_WOA2018_c3.0_d1.0_v19812010.5.1', -1., 'votemper',  .true.      , .true., 'yearly'   ,  ''         , ' '      , ' '
   sn_sal_dmp  = 'eORCA025.L121_WOA2018_c3.0_d1.0_v19812010.5.1', -1., 'vosaline',  .true.      , .true., 'yearly'   ,  ''         , ' '      , ' '
   !
/
namtra_ldf
  • We get back to UKMO values (150m2/s)

!-----------------------------------------------------------------------
&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection)
!-----------------------------------------------------------------------
   !                       !  Operator type:
   ln_traldf_OFF   = .false.   !  No explicit diffusion
   ln_traldf_lap   = .true.    !    laplacian operator
   ln_traldf_blp   = .false.   !  bilaplacian operator
   !
   !                       !  Direction of action:
   ln_traldf_lev   = .false.   !  iso-level
   ln_traldf_hor   = .false.   !  horizontal  (geopotential)
   ln_traldf_iso   = .true.    !  iso-neutral (standard operator)
   ln_traldf_triad = .false.   !  iso-neutral (triad    operator)
   !
   !                             !  iso-neutral options:
   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators)
   rn_slpmax       =  0.01     !  slope limit                           (both operators)
   ln_triad_iso    = .false.   !  pure horizontal mixing in ML              (triad only)
   rn_sw_triad     = 1         !  =1 switching triad ; =0 all 4 triads used (triad only)
   ln_botmix_triad = .false.   !  lateral mixing on bottom                  (triad only)
   !
   !                       !  Coefficients:
   nn_aht_ijk_t    = 20        !  space/time variation of eddy coefficient:
      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file
      !                             !   =  0           constant
      !                             !   = 10 F(k)      =ldf_c1d
      !                             !   = 20 F(i,j)    =ldf_c2d
      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)
      !                             !                           or   = 1/12 Ud*Ld^3 (blp case)
      rn_Ud        = 0.011          !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
      rn_Ld        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10)
/
namtra_eiv
  • Based on UKMO test and Julie D. discussion we switch one GM (150m2/s) (300m2/s2 in OPM016)

!-----------------------------------------------------------------------
&namtra_eiv    !   eddy induced velocity param.                         (default: OFF)
!-----------------------------------------------------------------------
   ln_ldfeiv   = .true.    ! use eddy induced velocity parameterization
      !
      !                        !  Coefficients:
      nn_aei_ijk_t    = 21          !  space/time variation of eddy coefficient:
      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
      !                             !   =  0           constant
      !                             !   = 10 F(k)      =ldf_c1d
      !                             !   = 20 F(i,j)    =ldf_c2d
      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le
      rn_Ue        = 0.03           !  lateral diffusive velocity [m/s] (nn_aei_ijk_t= 0, 10, 20, 30)
      rn_Le        = 10.e+3         !  lateral diffusive length   [m]   (nn_aei_ijk_t= 0, 10)
      !
      ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities
      nn_ldfeiv_shape = 2           !  shape of bounding coefficient    (nn_aei_ijk_t= 21 only)
/

namtra_dmp_drk

  • we switch off AABW restoring.

!-----------------------------------------------------------------------
&namtra_dmp_drk !   tracer: T & S newtonian damping                      (default: OFF)
!-----------------------------------------------------------------------
   nn_hdmp     =   -2      !  -2 = Drakkar customisation ( use dtacof in tradmp.F90 )
                           !  any other value : Nemo standard code
   nn_file     =    1      !  create a damping.coeff NetCDF file (=1) or not (=0)
                           !  The output file can then be used as input resto file (cn_resto)
   ln_dmpmask  = .false.    !  Read dmp_mask.nc file  when T (between 0 and 1 )
      rn_timsk    =  730.     !  Time scale used for dmp_mask
      cn_dir      =  './'     ! directtory wher to find damping mask
   !___________!____________!___________________!___________!_____________!________!___________!___________!__________!_______________!
   !           !  file name ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights e ! rotation ! land/sea mask !
   !           !            !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !  filename ! pairing  !    filename   !
   sn_dmp      = 'eORCA025.L121_dmpmask_b0.2_c0.3_d1.0_v0.0' ,  -12. ,  'wdmp' , .false.  , .true. , 'yearly'  ,   ''      ,   ''     ,   ''
/

Input files

Only the file changed between the reference (bug_input_eO025-OPM016) and this simulation are described.

Monitoring

see Global indicators