Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
No results found
Show changes
Hide whitespace changes
Inline Side-by-side
Showing
with 395 additions and 1258 deletions
cfgs/ORCA2_ICE_ABL/EXPREF/field_def_nemo-oce.xml deleted 120000 → 0
View file @ e3e3786b
../../SHARED/field_def_nemo-oce.xml
\ No newline at end of file
cfgs/ORCA2_ICE_ABL/EXPREF/file_def_nemo-oce.xml deleted 100644 → 0
View file @ e3e3786b
<?xml version="1.0"?>
<!--
============================================================================================================
= output files definition =
= Define your own filesfor ocean dynamics context =
= put the variables you want... =
============================================================================================================
-->
<file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="5d" min_digits="4">
<file_group id="5d" output_freq="5d" output_level="10" enabled=".TRUE."> <!-- 5d files -->
<file id="file11" name_suffix="_grid_T" description="ocean T grid variables" >
<field field_ref="e3t" />
<field field_ref="toce" name="thetao" />
<field field_ref="soce" name="so" />
<field field_ref="sst" name="tos" />
<field field_ref="sss" name="sos" />
<field field_ref="ssh" name="zos" />
<field field_ref="mldkz5" />
<field field_ref="mldr10_1" />
<field field_ref="sbt" />
<field field_ref="heatc" name="heatc" />
<field field_ref="saltc" name="saltc" />
<field field_ref="empmr" name="wfo" />
<field field_ref="qsr_oce" name="qsr_oce" />
<field field_ref="qns_oce" name="qns_oce" />
<field field_ref="qt_oce" name="qt_oce" />
<field field_ref="saltflx" name="sfx" />
<field field_ref="taum" name="taum" />
<field field_ref="wspd" name="windsp" />
<field field_ref="precip" name="precip" />
<!-- ice and snow -->
<field field_ref="snowpre" />
</file>
<file id="file12" name_suffix="_grid_U" description="ocean U grid variables" >
<field field_ref="e3u" />
<field field_ref="ssu" name="uos" />
<field field_ref="uoce" name="uo" />
<field field_ref="utau" name="tauuo" />
</file>
<file id="file13" name_suffix="_grid_V" description="ocean V grid variables" >
<field field_ref="e3v" />
<field field_ref="ssv" name="vos" />
<field field_ref="voce" name="vo" />
<field field_ref="vtau" name="tauvo" />
</file>
<file id="file14" name_suffix="_grid_ABL" description="ABL grid variables" >
<field field_ref="u_abl" />
<field field_ref="v_abl" />
<field field_ref="t_abl" />
<field field_ref="q_abl" />
<field field_ref="uvz1_abl" />
<field field_ref="tz1_abl" />
<field field_ref="qz1_abl" />
<field field_ref="uvz1_dta" />
<field field_ref="tz1_dta" />
<field field_ref="qz1_dta" />
<field field_ref="pblh" />
</file>
</file_group>
</file_definition>
cfgs/ORCA2_ICE_ABL/EXPREF/grid_def_nemo.xml deleted 120000 → 0
View file @ e3e3786b
../../SHARED/grid_def_nemo.xml
\ No newline at end of file
cfgs/ORCA2_ICE_ABL/EXPREF/namelist_cfg deleted 100644 → 0
View file @ e3e3786b
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE Configuration namelist : overwrite default values defined in SHARED/namelist_ref
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! ORCA2 - ICE - PISCES configuration !!
!!======================================================================
!! *** Domain & Run management namelists *** !!
!! !!
!! namrun parameters of the run
!! namdom space and time domain
!! namcfg parameters of the configuration (default: user defined GYRE)
!! namwad Wetting and drying (default: OFF)
!! namtsd data: temperature & salinity (default: OFF)
!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
!! namc1d 1D configuration options (ln_c1d =T)
!! namc1d_dyndmp 1D newtonian damping applied on currents (ln_c1d =T)
!! namc1d_uvd 1D data (currents) (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
cn_exp = "ORCA2" ! experience name
nn_it000 = 1 ! first time step
nn_itend = 16 ! last time step (std 5475)
nn_date0 = 20130101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
nn_write = 4 ! frequency of write in the output file (modulo referenced to nn_it000)
nn_istate = 0 ! output the initial state (1) or not (0)
ln_mskland = .false. ! mask land points in NetCDF outputs
/
!-----------------------------------------------------------------------
&namdom ! time and space domain
!-----------------------------------------------------------------------
rn_Dt = 5400. ! time step for the dynamics and tracer
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg)
!-----------------------------------------------------------------------
ln_read_cfg = .true. ! (=T) read the domain configuration file
cn_domcfg = "ORCA_R2_zps_domcfg" ! domain configuration filename
!
ln_closea = .false. ! F => suppress closed seas (defined by closea_mask field)
! ! from the bathymetry at runtime.
/
!-----------------------------------------------------------------------
&namtile ! parameters of the tiling
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtsd ! Temperature & Salinity Data (init/dmp) (default: OFF)
!-----------------------------------------------------------------------
! ! =T read T-S fields for:
ln_tsd_init = .true. ! ocean initialisation
ln_tsd_dmp = .true. ! T-S restoring (see namtra_dmp)
cn_dir = './' ! root directory for the T-S data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_tem = 'data_1m_potential_temperature_nomask', -1. ,'votemper', .true. , .true. , 'yearly' , '' , '' , ''
sn_sal = 'data_1m_salinity_nomask' , -1. ,'vosaline', .true. , .true. , 'yearly' , '' , '' , ''
/
!!======================================================================
!! *** Surface Boundary Condition namelists *** !!
!! !!
!! namsbc surface boundary condition manager (default: NO selection)
!! namsbc_flx flux formulation (ln_flx =T)
!! namsbc_blk Bulk formulae formulation (ln_blk =T)
!! namsbc_cpl CouPLed formulation ("key_oasis3" )
!! namsbc_sas Stand-Alone Surface module (SAS_SRC only)
!! namsbc_iif Ice-IF: use observed ice cover (nn_ice = 1 )
!! namtra_qsr penetrative solar radiation (ln_traqsr =T)
!! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T)
!! namsbc_rnf river runoffs (ln_rnf =T)
!! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)
!! namsbc_isf ice shelf melting/freezing (ln_isfcav =T : read (ln_read_cfg=T) or set or usr_def_zgr )
!! namsbc_iscpl coupling option between land ice model and ocean (ln_isfcav =T)
!! namsbc_wave external fields from wave model (ln_wave =T)
!! namberg iceberg floats (ln_icebergs=T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namsbc ! Surface Boundary Condition manager (default: NO selection)
!-----------------------------------------------------------------------
nn_fsbc = 1 ! frequency of SBC module call
! (also = the frequency of sea-ice & iceberg model call)
! Type of air-sea fluxes
ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk )
ln_abl = .true. ! ABL formulation (T => fill namsbc_abl )
! Sea-ice :
nn_ice = 2 ! =0 no ice boundary condition
! ! =1 use observed ice-cover ( => fill namsbc_iif )
! ! =2 or 3 for SI3 and CICE, respectively
! Misc. options of sbc :
ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr)
ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave
ln_rnf = .true. ! runoffs (T => fill namsbc_rnf)
nn_fwb = 2 ! FreshWater Budget:
! ! =2 annual global mean of e-p-r set to zero
ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave)
ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
! ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
! ! = 1 Phillips: v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
! ! = 2 Phillips as (1) but using the wave frequency from a wave model
ln_tauwoc = .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
ln_tauw = .false. ! Activate ocean stress components from wave model
ln_stcor = .false. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
/
!-----------------------------------------------------------------------
&namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T)
!-----------------------------------------------------------------------
! ! bulk algorithm :
ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008)
!
cn_dir = './' ! root directory for the bulk data location
!___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_wndi = 'uwnd_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d', 24., 'uwnd' , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic' , 'U10' , ''
sn_wndj = 'vwnd_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d', 24., 'vwnd' , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic' , 'V10' , ''
sn_tair = 'tair_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d', 24., 'tair' , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bilinear' , '' , ''
sn_humi = 'humi_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d', 24., 'humi' , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bilinear' , '' , ''
sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24., 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24., 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
sn_prec = 'ncar_precip.15JUNE2009_fill' , -1., 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
sn_snow = 'ncar_precip.15JUNE2009_fill' , -1., 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
sn_slp = 'slp.15JUNE2009_fill' , 6., 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
sn_hpgi = 'uhpg_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d', 24., 'uhpg' , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic' , 'UG' , ''
sn_hpgj = 'vhpg_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d', 24., 'vhpg' , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic' , 'VG' , ''
/
!-----------------------------------------------------------------------
&namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T)
!-----------------------------------------------------------------------
cn_dir = './' ! root directory for the location of the ABL grid file
cn_dom = 'dom_cfg_abl_L25Z10.nc'
ln_hpgls_frc = .true.
ln_geos_winds = .false.
nn_dyn_restore = 1
rn_ldyn_min = 7.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour]
rn_ldyn_max = 1.5 ! magnitude of the nudging on ABL dynamics at the top of the ABL [hour]
rn_ltra_min = 7.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour]
rn_ltra_max = 1.5 ! magnitude of the nudging on ABL tracers at the top of the ABL [hour]
/
!-----------------------------------------------------------------------
&namtra_qsr ! penetrative solar radiation (ln_traqsr =T)
!-----------------------------------------------------------------------
! ! type of penetration (default: NO selection)
ln_qsr_rgb = .true. ! RGB light penetration (Red-Green-Blue)
!
nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0)
cn_dir = './' ! root directory for the chlorophyl data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_chl ='chlorophyll' , -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , ''
/
!-----------------------------------------------------------------------
&namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr =T)
!-----------------------------------------------------------------------
nn_sssr = 2 ! add a damping term to the surface freshwater flux
rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day]
ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2)
rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day]
/
!-----------------------------------------------------------------------
&namsbc_rnf ! runoffs (ln_rnf =T)
!-----------------------------------------------------------------------
ln_rnf_mouth = .true. ! specific treatment at rivers mouths
rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T)
rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
rn_rfact = 1.e0 ! multiplicative factor for runoff
cn_dir = './' ! root directory for the location of the runoff files
!___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________!
! ! 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 !
sn_rnf = 'runoff_core_monthly', -1. , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , ''
sn_cnf = 'runoff_core_monthly', 0. , 'socoefr0', .false. , .true. , 'yearly' , '' , '' , ''
sn_s_rnf = 'runoffs' , 24. , 'rosaline', .true. , .true. , 'yearly' , '' , '' , ''
sn_t_rnf = 'runoffs' , 24. , 'rotemper', .true. , .true. , 'yearly' , '' , '' , ''
sn_dep_rnf = 'runoffs' , 0. , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , ''
/
!-----------------------------------------------------------------------
&namsbc_wave ! External fields from wave model (ln_wave=T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namberg ! iceberg parameters (default: OFF)
!-----------------------------------------------------------------------
/
!!======================================================================
!! *** Lateral boundary condition *** !!
!! !!
!! namlbc lateral momentum boundary condition (default: NO selection)
!! namagrif agrif nested grid (read by child model only) ("key_agrif")
!! nam_tide Tidal forcing (default: OFF)
!! nambdy Unstructured open boundaries (default: OFF)
!! nambdy_dta Unstructured open boundaries - external data (see nambdy)
!! nambdy_tide tidal forcing at open boundaries (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namlbc ! lateral momentum boundary condition (default: NO selection)
!-----------------------------------------------------------------------
rn_shlat = 2. ! no slip
/
!-----------------------------------------------------------------------
&namagrif ! AGRIF zoom ("key_agrif")
!-----------------------------------------------------------------------
ln_spc_dyn = .true. ! use 0 as special value for dynamics
rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s]
rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s]
ln_chk_bathy = .false. ! =T check the parent bathymetry
/
!!======================================================================
!! *** Top/Bottom boundary condition *** !!
!! !!
!! namdrg top/bottom drag coefficient (default: NO selection)
!! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T)
!! namdrg_bot bottom friction (ln_drg_OFF=F)
!! nambbc bottom temperature boundary condition (default: OFF)
!! nambbl bottom boundary layer scheme (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namdrg ! top/bottom drag coefficient (default: NO selection)
!-----------------------------------------------------------------------
ln_lin = .true. ! linear drag: Cd = Cd0 Uc0
/
!-----------------------------------------------------------------------
&nambbc ! bottom temperature boundary condition (default: OFF)
!-----------------------------------------------------------------------
ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom
nn_geoflx = 2 ! geothermal heat flux: = 2 read variable flux [mW/m2]
cn_dir = './' ! root directory for the geothermal data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_qgh ='geothermal_heating.nc' , -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , ''
/
!-----------------------------------------------------------------------
&nambbl ! bottom boundary layer scheme (default: OFF)
!-----------------------------------------------------------------------
ln_trabbl = .true. ! Bottom Boundary Layer parameterisation flag
nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0)
nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0)
rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s]
rn_gambbl = 10. ! advective bbl coefficient [s]
/
!!======================================================================
!! Tracer (T-S) namelists !!
!! !!
!! nameos equation of state (default: NO selection)
!! namtra_adv advection scheme (default: NO selection)
!! namtra_ldf lateral diffusion scheme (default: NO selection)
!! namtra_mle mixed layer eddy param. (Fox-Kemper param.) (default: OFF)
!! namtra_eiv eddy induced velocity param. (default: OFF)
!! namtra_dmp T & S newtonian damping (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&nameos ! ocean Equation Of Seawater (default: NO selection)
!-----------------------------------------------------------------------
ln_eos80 = .true. ! = Use EOS80
/
!-----------------------------------------------------------------------
&namtra_adv ! advection scheme for tracer (default: NO selection)
!-----------------------------------------------------------------------
ln_traadv_fct = .true. ! FCT scheme
nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order
nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order
/
!-----------------------------------------------------------------------
&namtra_ldf ! lateral diffusion scheme for tracers (default: NO selection)
!-----------------------------------------------------------------------
ln_traldf_lap = .true. ! laplacian operator
ln_traldf_iso = .true. ! iso-neutral (Standard operator)
ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators)
! ! Coefficients:
nn_aht_ijk_t = 20 ! space/time variation of eddy coef
! ! = 20 aht = 1/2 Ud. max(e1,e2)
rn_Ud = 0.018 ! 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_mle ! mixed layer eddy parametrisation (Fox-Kemper) (default: OFF)
!-----------------------------------------------------------------------
ln_mle = .true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation
/
!-----------------------------------------------------------------------
&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 the eiv coeficient
! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
! ! time invariant coefficients: aei0 = 1/2 Ue*Le
rn_Ue = 0.03 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10)
!
ln_ldfeiv_dia =.true. ! diagnose eiv stream function and velocities
/
!-----------------------------------------------------------------------
&namtra_dmp ! tracer: T & S newtonian damping (default: OFF)
!-----------------------------------------------------------------------
ln_tradmp = .true. ! add a damping term (using resto.nc coef.)
nn_zdmp = 0 ! vertical shape =0 damping throughout the water column
/
!!======================================================================
!! *** Dynamics namelists *** !!
!! !!
!! nam_vvl vertical coordinate options (default: z-star)
!! namdyn_adv formulation of the momentum advection (default: NO selection)
!! namdyn_vor advection scheme (default: NO selection)
!! namdyn_hpg hydrostatic pressure gradient (default: NO selection)
!! namdyn_spg surface pressure gradient (default: NO selection)
!! namdyn_ldf lateral diffusion scheme (default: NO selection)
!! namdta_dyn offline TOP: dynamics read in files (OFF_SRC only)
!!======================================================================
!
!-----------------------------------------------------------------------
&namdyn_adv ! formulation of the momentum advection (default: NO selection)
!-----------------------------------------------------------------------
ln_dynadv_vec = .true. ! vector form - 2nd centered scheme
nn_dynkeg = 0 ! grad(KE) scheme: =0 C2 ; =1 Hollingsworth correction
/
!-----------------------------------------------------------------------
&namdyn_vor ! Vorticity / Coriolis scheme (default: NO selection)
!-----------------------------------------------------------------------
ln_dynvor_een = .true. ! energy & enstrophy scheme
nn_een_e3f = 0 ! =0 e3f = mean masked e3t divided by 4
/
!-----------------------------------------------------------------------
&namdyn_hpg ! Hydrostatic pressure gradient option (default: NO selection)
!-----------------------------------------------------------------------
ln_hpg_sco = .true. ! s-coordinate (standard jacobian formulation)
/
!-----------------------------------------------------------------------
&namdyn_spg ! surface pressure gradient (default: NO selection)
!-----------------------------------------------------------------------
ln_dynspg_ts = .true. ! split-explicit free surface
/
!-----------------------------------------------------------------------
&namdyn_ldf ! lateral diffusion on momentum (default: NO selection)
!-----------------------------------------------------------------------
ln_dynldf_lap = .true. ! laplacian operator
ln_dynldf_lev = .true. ! iso-level
nn_ahm_ijk_t = -30 ! =-30 3D coeff. read in eddy_diffusivity_3D.nc
/
!!======================================================================
!! vertical physics namelists !!
!! !!
!! namzdf vertical physics manager (default: NO selection)
!! namzdf_ric richardson number vertical mixing (ln_zdfric=T)
!! namzdf_tke TKE vertical mixing (ln_zdftke=T)
!! namzdf_gls GLS vertical mixing (ln_zdfgls=T)
!! namzdf_osm OSM vertical diffusion (ln_zdfosm=T)
!! namzdf_iwm tidal mixing parameterization (ln_zdfiwm=T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namzdf ! vertical physics manager (default: NO selection)
!-----------------------------------------------------------------------
ln_zdftke = .true. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke)
ln_zdfevd = .true. ! Enhanced Vertical Diffusion scheme
nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1)
rn_evd = 100. ! evd mixing coefficient [m2/s]
ln_zdfddm = .true. ! double diffusive mixing
rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity)
rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio
ln_zdfiwm = .true. ! internal wave-induced mixing (T => fill namzdf_iwm)
! ! Coefficients
rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F)
rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)
nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)
nn_havtb = 1 ! horizontal shape for avtb (=1) or not (=0)
/
!-----------------------------------------------------------------------
&namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T)
!-----------------------------------------------------------------------
nn_mxl = 3 ! mixing length: = 0 bounded by the distance to surface and bottom
! ! = 1 bounded by the local vertical scale factor
! ! = 2 first vertical derivative of mixing length bounded by 1
! ! = 3 as =2 with distinct dissipative an mixing length scale
nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to NIWs
! = 0 none ; = 1 add a tke source below the ML
! = 2 add a tke source just at the base of the ML
! = 3 as = 1 applied on HF part of the stress (ln_cpl=T)
rn_eice = 0 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4
/
!-----------------------------------------------------------------------
&namzdf_iwm ! internal wave-driven mixing parameterization (ln_zdfiwm =T)
!-----------------------------------------------------------------------
nn_zpyc = 2 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency
ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F)
/
!!======================================================================
!! *** Diagnostics namelists *** !!
!! !!
!! namtrd dynamics and/or tracer trends (default: OFF)
!! namhsb Heat and salt budgets (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namflo float parameters (default: OFF)
!! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)
!! nam_diadct transports through some sections (default: OFF)
!! nam_diatmb Top Middle Bottom Output (default: OFF)
!! nam_dia25h 25h Mean Output (default: OFF)
!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
!!======================================================================
!
!!======================================================================
!! *** Observation & Assimilation *** !!
!! !!
!! namobs observation and model comparison (default: OFF)
!! nam_asminc assimilation increments ('key_asminc')
!!======================================================================
!
!!======================================================================
!! *** Miscellaneous namelists *** !!
!! !!
!! nammpp Massively Parallel Processing
!! namctl Control prints (default: OFF)
!! namsto Stochastic parametrization of EOS (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&nammpp ! Massively Parallel Processing
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namctl ! Control prints (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsto ! Stochastic parametrization of EOS (default: OFF)
!-----------------------------------------------------------------------
/
cfgs/ORCA2_ICE_ABL/EXPREF/namelist_ice_cfg deleted 100644 → 0
View file @ e3e3786b
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! SI3 configuration namelist: Overwrites SHARED/namelist_ice_ref
!! 1 - Generic parameters (nampar)
!! 2 - Ice thickness discretization (namitd)
!! 3 - Ice dynamics (namdyn)
!! 4 - Ice ridging/rafting (namdyn_rdgrft)
!! 5 - Ice rheology (namdyn_rhg)
!! 6 - Ice advection (namdyn_adv)
!! 7 - Ice surface boundary conditions (namsbc)
!! 8 - Ice thermodynamics (namthd)
!! 9 - Ice heat diffusion (namthd_zdf)
!! 10 - Ice lateral melting (namthd_da)
!! 11 - Ice growth in open water (namthd_do)
!! 12 - Ice salinity (namthd_sal)
!! 13 - Ice melt ponds (namthd_pnd)
!! 14 - Ice initialization (namini)
!! 15 - Ice/snow albedos (namalb)
!! 16 - Ice diagnostics (namdia)
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!------------------------------------------------------------------------------
&nampar ! Generic parameters
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namitd ! Ice discretization
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namdyn ! Ice dynamics
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namdyn_rdgrft ! Ice ridging/rafting
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namdyn_rhg ! Ice rheology
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namdyn_adv ! Ice advection
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namsbc ! Ice surface boundary conditions
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd ! Ice thermodynamics
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_zdf ! Ice heat diffusion
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_da ! Ice lateral melting
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_do ! Ice growth in open water
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_sal ! Ice salinity
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namthd_pnd ! Melt ponds
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namini ! Ice initialization
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namalb ! albedo parameters
!------------------------------------------------------------------------------
/
!------------------------------------------------------------------------------
&namdia ! Diagnostics
!------------------------------------------------------------------------------
/
cfgs/ORCA2_ICE_ABL/EXPREF/namelist_ice_ref deleted 120000 → 0
View file @ e3e3786b
../../SHARED/namelist_ice_ref
\ No newline at end of file
cfgs/ORCA2_ICE_ABL/EXPREF/namelist_ref deleted 120000 → 0
View file @ e3e3786b
../../SHARED/namelist_ref
\ No newline at end of file
cfgs/ORCA2_ICE_ABL/cpp_ORCA2_ICE_ABL.fcm deleted 100644 → 0
View file @ e3e3786b
bld::tool::fppkeys key_si3 key_xios key_qco
cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_top_cfg
View file @ 4a9f4b18
......@@ -114,7 +114,7 @@
sn_trcsbc(14) = 'dust.orca.new' , -1 , 'dustfer' , .true. , .true. , 'yearly' , '' , '' , ''
sn_trcsbc(23) = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' , ''
rn_trsfac(5) = 3.774194e-02 ! ( 1E-3 / 31. * 117 )
rn_trsfac(7) = 9.572954e-03 ! ( 8.8 / 28.1 )
rn_trsfac(7) = 9.572954e-03 ! ( 0.269 / 28.1 )
rn_trsfac(14) = 6.2667860e-04 ! ( 0.035 / 55.85 )
rn_trsfac(23) = 5.2232143e-01 ! ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.3125/14 )
rn_sbc_time = 1. ! Time scaling factor for SBC and CBC data (seconds in a day)
......
cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_top_cfg
View file @ 4a9f4b18
......@@ -114,7 +114,7 @@
sn_trcsbc(14) = 'dust.orca.new' , -1 , 'dustfer' , .true. , .true. , 'yearly' , '' , '' , ''
sn_trcsbc(23) = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' , ''
rn_trsfac(5) = 3.774194e-02 ! ( 1E-3 / 31. * 117 )
rn_trsfac(7) = 9.572954e-03 ! ( 8.8 / 28.1 )
rn_trsfac(7) = 9.572954e-03 ! ( 0.269 / 28.1 )
rn_trsfac(14) = 6.2667860e-04 ! ( 0.035 / 55.85 )
rn_trsfac(23) = 5.2232143e-01 ! ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.3125/14 )
rn_sbc_time = 1. ! Time scaling factor for SBC and CBC data (seconds in a day)
......
cfgs/README.rst deleted 100644 → 0
View file @ e3e3786b
********************************
Run the Reference configurations
********************************
.. todo::
Lack of illustrations for ref. cfgs, and more generally in the guide.
NEMO is distributed with a set of reference configurations allowing both
the user to set up his own first applications and
the developer to test/validate his NEMO developments (using SETTE package).
.. contents::
:local:
:depth: 1
.. attention::
Concerning the configurations,
the NEMO System Team is only in charge of the so-called reference configurations described below.
.. hint::
Configurations developed by external research projects or initiatives that
make use of NEMO are welcome to be publicized through the website by
filling up the form :website:`to add an associated project<projects/add>`.
How to compile an experiment from a reference configuration
===========================================================
To compile the ORCA2_ICE_PISCES_ reference configuration using :file:`makenemo`,
one should use the following, by selecting among available architecture file or
providing a user defined one:
.. code-block:: console
$ ./makenemo -r 'ORCA2_ICE_PISCES' -m 'my_arch' -j '4'
A new ``EXP00`` folder will be created within the selected reference configurations,
namely ``./cfgs/ORCA2_ICE_PISCES/EXP00``.
It will be necessary to uncompress the archives listed in the above table for
the given reference configuration that includes input & forcing files.
Then it will be possible to launch the execution of the model through a runscript
(opportunely adapted to the user system).
List of Configurations
======================
All forcing files listed below in the table are available from |DOI data|_
=================== === === === === === ==================================
Configuration Component(s) Archives (input & forcing files)
------------------- ------------------- ----------------------------------
Name O S T P A
=================== === === === === === ==================================
AGRIF_DEMO_ X X X AGRIF_DEMO_v4.0.tar,
ORCA2_ICE_v4.0.tar
AMM12_ X AMM12_v4.0.tar
C1D_PAPA_ X INPUTS_C1D_PAPA_v4.0.tar
GYRE_BFM_ X X *none*
GYRE_PISCES_ X X X *none*
ORCA2_ICE_PISCES_ X X X X ORCA2_ICE_v4.0.tar,
INPUTS_PISCES_v4.0.tar
ORCA2_OFF_PISCES_ X X ORCA2_OFF_v4.0.tar,
INPUTS_PISCES_v4.0.tar
ORCA2_OFF_TRC_ X ORCA2_OFF_v4.0.tar
ORCA2_SAS_ICE_ X ORCA2_ICE_v4.0.tar,
INPUTS_SAS_v4.0.tar
SPITZ12_ X X SPITZ12_v4.0.tar
=================== === === === === === ==================================
.. admonition:: Legend for component combination
O for OCE, S for SI\ :sup:`3`, T for TOP, P for PISCES and A for AGRIF
AGRIF_DEMO
----------
``AGRIF_DEMO`` is based on the ``ORCA2_ICE_PISCES`` global configuration at 2° of resolution with
the inclusion of 3 online nested grids to demonstrate the overall capabilities of AGRIF in
a realistic context (including the nesting of sea ice models).
The configuration includes a 1:1 grid in the Pacific and two successively nested grids with
odd and even refinement ratios over the Arctic ocean,
with the finest grid spanning the whole Svalbard archipelago that is of
particular interest to test sea ice coupling.
.. image:: _static/AGRIF_DEMO_no_cap.jpg
:scale: 66%
:align: center
The 1:1 grid can be used alone as a benchmark to check that
the model solution is not corrupted by grid exchanges.
Note that since grids interact only at the baroclinic time level,
numerically exact results can not be achieved in the 1:1 case.
Perfect reproducibility is obtained only by switching to a fully explicit setup instead of
a split explicit free surface scheme.
AMM12
-----
``AMM12`` stands for *Atlantic Margin Model at 12 km* that is
a regional configuration covering the Northwest European Shelf domain on
a regular horizontal grid of ~12 km of resolution (see :cite:`ODEA2012`).
.. image:: _static/AMM_domain.png
:align: center
This configuration allows to tests several features of NEMO specifically addressed to the shelf seas.
In particular, ``AMM12`` accounts for vertical s-coordinates system, GLS turbulence scheme,
tidal lateral boundary conditions using a flather scheme (see more in ``BDY``).
Boundaries may be completely omitted by setting ``ln_bdy = .false.`` in ``nambdy``.
Sample surface fluxes, river forcing and an initial restart file are included to test a realistic model run
(``AMM12_v4.0.tar``).
Note that, the Baltic boundary is included within the river input file and is specified as a river source,
but unlike ordinary river points the Baltic inputs also include salinity and temperature data.
C1D_PAPA
--------
.. figure:: _static/Papa2015.jpg
:height: 225px
:align: left
``C1D_PAPA`` is a 1D configuration for the `PAPA station`_ located in
the northern-eastern Pacific Ocean at 50.1°N, 144.9°W.
See :gmd:`Reffray et al. (2015) <8/69/2015>` for the description of
its physical and numerical turbulent-mixing behaviour.
| The water column setup, called NEMO1D, is activated by
setting ``ln_c1d = .true.`` in ``namdom`` and
has a horizontal domain of 1x1 grid point.
| This reference configuration uses 75 vertical levels grid (1m at the surface),
GLS turbulence scheme with K-epsilon closure and the NCAR bulk formulae.
Data provided with ``INPUTS_C1D_PAPA_v4.2.tar`` file account for:
- :file:`forcing_PAPASTATION_1h_y201[0-1].nc`:
ECMWF operational analysis atmospheric forcing rescaled to 1h
(with long and short waves flux correction) for years 2010 and 2011
- :file:`init_PAPASTATION_m06d15.nc`: Initial Conditions from
observed data and Levitus 2009 climatology
- :file:`chlorophyll_PAPASTATION.nc`: surface chlorophyll file from Seawifs data
GYRE_BFM
--------
``GYRE_BFM`` shares the same physical setup of GYRE_PISCES_,
but NEMO is coupled with the `BFM`_ biogeochemical model as described in ``./cfgs/GYRE_BFM/README``.
GYRE_PISCES
-----------
``GYRE_PISCES`` is an idealized configuration representing a Northern hemisphere double gyres system,
in the Beta-plane approximation with a regular 1° horizontal resolution and 31 vertical levels,
with PISCES BGC model :cite:`gmd-8-2465-2015`.
Analytical forcing for heat, freshwater and wind-stress fields are applied.
This configuration acts also as demonstrator of the **user defined setup**
(``ln_read_cfg = .false.``) and grid setting are handled through
the ``&namusr_def`` controls in :file:`namelist_cfg`:
.. literalinclude:: ../../../cfgs/GYRE_PISCES/EXPREF/namelist_cfg
:language: fortran
:lines: 35-41
Note that, the default grid size is 30x20 grid points (with ``nn_GYRE = 1``) and
vertical levels are set by ``jpkglo``.
The specific code changes can be inspected in :file:`./src/OCE/USR`.
.. rubric:: Running GYRE as a benchmark
| This simple configuration can be used as a benchmark since it is easy to increase resolution,
with the drawback of getting results that have a very limited physical meaning.
| GYRE grid resolution can be increased at runtime by setting a different value of ``nn_GYRE``
(integer multiplier scaling factor), as described in the following table:
=========== ============ ============ ============ ===============
``nn_GYRE`` ``jpiglo`` ``jpjglo`` ``jpkglo`` Equivalent to
=========== ============ ============ ============ ===============
1 30 20 31 GYRE 1°
25 750 500 101 ORCA 1/2°
50 1500 1000 101 ORCA 1/4°
150 4500 3000 101 ORCA 1/12°
200 6000 4000 101 ORCA 1/16°
=========== ============ ============ ============ ===============
| Note that, it is necessary to set ``ln_bench = .true.`` in ``&namusr_def`` to
avoid problems in the physics computation and that
the model timestep should be adequately rescaled.
| For example if ``nn_GYRE = 150``, equivalent to an ORCA 1/12° grid,
the timestep ``rn_rdt`` should be set to 1200 seconds
Differently from previous versions of NEMO, the code uses by default the time-splitting scheme and
internally computes the number of sub-steps.
ORCA2_ICE_PISCES
----------------
``ORCA2_ICE_PISCES`` is a reference configuration for the global ocean with
a 2°x2° curvilinear horizontal mesh and 31 vertical levels,
distributed using z-coordinate system and with 10 levels in the top 100m.
ORCA is the generic name given to global ocean Mercator mesh,
(i.e. variation of meridian scale factor as cosinus of the latitude),
with two poles in the northern hemisphere so that
the ratio of anisotropy is nearly one everywhere
This configuration uses the three components
- |OCE|, the ocean dynamical core
- |ICE|, the thermodynamic-dynamic sea ice model.
- |MBG|, passive tracer transport module and PISCES BGC model :cite:`gmd-8-2465-2015`
All components share the same grid.
The model is forced with CORE-II normal year atmospheric forcing and
it uses the NCAR bulk formulae.
.. rubric:: Ocean Physics
:horizontal diffusion on momentum:
the eddy viscosity coefficient depends on the geographical position.
It is taken as 40000 m\ :sup:`2`/s, reduced in the equator regions (2000 m\ :sup:`2`/s)
excepted near the western boundaries.
:isopycnal diffusion on tracers:
the diffusion acts along the isopycnal surfaces (neutral surface) with
an eddy diffusivity coefficient of 2000 m\ :sup:`2`/s.
:Eddy induced velocity parametrization:
With a coefficient that depends on the growth rate of baroclinic instabilities
(it usually varies from 15 m\ :sup:`2`/s to 3000 m\ :sup:`2`/s).
:lateral boundary conditions:
Zero fluxes of heat and salt and no-slip conditions are applied through lateral solid boundaries.
:bottom boundary condition:
Zero fluxes of heat and salt are applied through the ocean bottom.
The Beckmann [19XX] simple bottom boundary layer parameterization is applied along
continental slopes.
A linear friction is applied on momentum.
:convection:
The vertical eddy viscosity and diffusivity coefficients are increased to 1 m\ :sup:`2`/s in
case of static instability.
:time step: is 5400sec (1h30') so that there is 16 time steps in one day.
ORCA2_OFF_PISCES
----------------
``ORCA2_OFF_PISCES`` shares the same general offline configuration of ``ORCA2_ICE_TRC``,
but only PISCES model is an active component of TOP.
ORCA2_OFF_TRC
-------------
| ``ORCA2_OFF_TRC`` is based on the ORCA2 global ocean configuration
(see ORCA2_ICE_PISCES_ for general description) along with
the tracer passive transport module (TOP),
but dynamical fields are pre-calculated and read with specific time frequency.
| This enables for an offline coupling of TOP components,
here specifically inorganic carbon compounds (CFC11, CFC12, SF6, C14) and water age module (age).
See :file:`namelist_top_cfg` to inspect the selection of
each component with the dedicated logical keys.
Pre-calculated dynamical fields are provided to NEMO using
the namelist ``&namdta_dyn`` in :file:`namelist_cfg`,
in this case with a 5 days frequency (120 hours):
.. literalinclude:: ../../namelists/namdta_dyn
:language: fortran
Input dynamical fields for this configuration (:file:`ORCA2_OFF_v4.0.tar`) comes from
a 2000 years long climatological simulation of ORCA2_ICE using ERA40 atmospheric forcing.
| Note that,
this configuration default uses linear free surface (``ln_linssh = .true.``) assuming that
model mesh is not varying in time and
it includes the bottom boundary layer parameterization (``ln_trabbl = .true.``) that
requires the provision of BBL coefficients through ``sn_ubl`` and ``sn_vbl`` fields.
| It is also possible to activate PISCES model (see ``ORCA2_OFF_PISCES``) or
a user defined set of tracers and source-sink terms with ``ln_my_trc = .true.``
(and adaptation of ``./src/TOP/MY_TRC`` routines).
In addition, the offline module (OFF) allows for the provision of further fields:
1. **River runoff** can be provided to TOP components by setting ``ln_dynrnf = .true.`` and
by including an input datastream similarly to the following:
.. code-block:: fortran
sn_rnf = 'dyna_grid_T', 120, 'sorunoff' , .true., .true., 'yearly', '', '', ''
2. **VVL dynamical fields**, in the case input data were produced by a dyamical core using
variable volume (``ln_linssh = .false.``)
it is necessary to provide also diverce and E-P at before timestep by
including input datastreams similarly to the following
.. code-block:: fortran
sn_div = 'dyna_grid_T', 120, 'e3t' , .true., .true., 'yearly', '', '', ''
sn_empb = 'dyna_grid_T', 120, 'sowaflupb', .true., .true., 'yearly', '', '', ''
More details can be found by inspecting the offline data manager in
the routine :file:`./src/OFF/dtadyn.F90`.
ORCA2_SAS_ICE
-------------
| ORCA2_SAS_ICE is a demonstrator of the Stand-Alone Surface (SAS) module and
it relies on ORCA2 global ocean configuration (see ORCA2_ICE_PISCES_ for general description).
| The standalone surface module allows surface elements such as sea-ice, iceberg drift, and
surface fluxes to be run using prescribed model state fields.
It can profitably be used to compare different bulk formulae or
adjust the parameters of a given bulk formula.
More informations about SAS can be found in :doc:`NEMO manual <cite>`.
SPITZ12
-------
``SPITZ12`` is a regional configuration around the Svalbard archipelago
at 1/12° of horizontal resolution and 75 vertical levels.
See :gmd:`Rousset et al. (2015) <8/2991/2015>` for more details.
This configuration references to year 2002,
with atmospheric forcing provided every 2 hours using NCAR bulk formulae,
while lateral boundary conditions for dynamical fields have 3 days time frequency.
.. rubric:: References
.. bibliography:: cfgs.bib
:all:
:style: unsrt
:labelprefix: C
cfgs/SHARED/README.namelists deleted 100644 → 0
View file @ e3e3786b
Simple style rules for namelists
--------------------------------
NEMO reference namelists should adhere to the following simple style rules:
1. Comments outside a namelist block start with !! in column 1
2. Each namelist block starts with 3 lines of the form:
!-----------------------------------------------------------------------
&namblockname ! short description of block
!-----------------------------------------------------------------------
with all ! and & 's starting in column 1
3. The closing / for each namelist block is in column 1
4. Comments within namelist blocks never start with !- . Use ! followed
by space or != etc.
These conventions make it possible to construct empty configuration namelists.
For example, a namelist_cfg template can be produced from namelist_ref with
the following grep command; e.g.:
grep -E '^!-|^&|^/' namelist_ref > namelist_cfg.template
head namelist_cfg.template
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdom ! time and space domain
!-----------------------------------------------------------------------
/
.
.
If all configuration namelists are produced and maintained using this
strategy then standard, side-by-side comaparators, such as vimdiff or xxdiff,
can be used to compare and transfer lines from the reference namelist to a
configuration namelist when setting up a new configuration.
Tips and tricks
---------------
1. The following bash function is useful when checking which namelist blocks
are in active use in a configuration namelist:
function list_used_nl(){ grep -n -E '^&|^/' "$1" | sed -e 's/:/ /' \
| awk ' BEGIN { x = 0 } \
{if ( NR % 2 == 0 && $1 - x > 2 ) printf("%3d %s\n", $1 - x , n) ; \
else x = $1; n = $2}' \
| sort -k 2;}
which (assuming the namelist adheres to the conventions) will list the number
of entries in each non-empty namelist block. The list is sorted on the block
name to ease comparisons. For example:
list_used_nl ORCA2_LIM3_PISCES/EXP00/namelist_cfg
1 &nambbc
5 &nambbl
30 &namberg
10 &namcfg
4 &namctl
3 &namdom
1 &namdrg
5 &namdyn_adv
1 &namdyn_hpg
22 &namdyn_ldf
1 &namdyn_spg
5 &namdyn_vor
3 &nameos
1 &namhsb
4 &namrun
1 &namsbc
1 &namsbc_blk
3 &namtra_adv
28 &namtra_ldf
10 &namtra_ldfeiv
25 &namzdf
3 &namzdf_iwm
2. vimdiff can give garish colours in some terminals. Usually this is because
vim assumes, incorrectly, that the terminal only supports 8 colours. Try forcing
256 colours with:
:set t_Co=256
to produce more pastel shades (add this to ~/.vimrc if successful).
3. Switching between vsplit panes in vim is a multi-key sequence. The tool is
much easier to use if the sequence is mapped to a spare key. Here I use the
§ and ± key on my Mac keyboard (add to ~/.vimrc):
map § ^Wl
map ± ^Wh
4. With easy switching between panes, constructing namelists in vimdiff just
requires the following commands in addition to normal editing:
]c - Go to next block of the diff
dp - Push version of the block under cursor into the other pane
do - Pull version of the block under cursor from the other pane
cfgs/SHARED/README.rst deleted 100644 → 0
View file @ e3e3786b
***********
Diagnostics
***********
.. todo::
.. contents::
:local:
Output of diagnostics in NEMO is usually done using XIOS.
This is an efficient way of writing diagnostics because
the time averaging, file writing and even some simple arithmetic or regridding is carried out in
parallel to the NEMO model run.
This page gives a basic introduction to using XIOS with NEMO.
Much more information is available from the :xios:`XIOS homepage<>` above and from the NEMO manual.
Use of XIOS for diagnostics is activated using the pre-compiler key ``key_xios``.
Extracting and installing XIOS
==============================
1. Install the NetCDF4 library.
If you want to use single file output you will need to compile the HDF & NetCDF libraries to
allow parallel IO.
2. Download the version of XIOS that you wish to use.
The recommended version is now XIOS 2.5:
.. code-block:: console
$ svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/branchs/xios-2.5
and follow the instructions in :xios:`XIOS documentation <wiki/documentation>` to compile it.
If you find problems at this stage, support can be found by subscribing to
the :xios:`XIOS mailing list <../mailman/listinfo.cgi/xios-users>` and sending a mail message to it.
XIOS Configuration files
------------------------
XIOS is controlled using XML input files that should be copied to
your model run directory before running the model.
Examples of these files can be found in the reference configurations (:file:`./cfgs`).
The XIOS executable expects to find a file called :file:`iodef.xml` in the model run directory.
In NEMO we have made the decision to use include statements in the :file:`iodef.xml` file to include:
- :file:`field_def_nemo-oce.xml` (for physics),
- :file:`field_def_nemo-ice.xml` (for ice),
- :file:`field_def_nemo-pisces.xml` (for biogeochemistry) and
- :file:`domain_def.xml` from the :file:`./cfgs/SHARED` directory.
Most users will not need to modify :file:`domain_def.xml` or :file:`field_def_nemo-???.xml` unless
they want to add new diagnostics to the NEMO code.
The definition of the output files is organized into separate :file:`file_definition.xml` files which
are included in the :file:`iodef.xml` file.
Modes
=====
Detached Mode
-------------
In detached mode the XIOS executable is executed on separate cores from the NEMO model.
This is the recommended method for using XIOS for realistic model runs.
To use this mode set ``using_server`` to ``true`` at the bottom of the :file:`iodef.xml` file:
.. code-block:: xml
<variable id="using_server" type="boolean">true</variable>
Make sure there is a copy (or link to) your XIOS executable in the working directory and
in your job submission script allocate processors to XIOS.
Attached Mode
-------------
In attached mode XIOS runs on each of the cores used by NEMO.
This method is less efficient than the detached mode but can be more convenient for testing or
with small configurations.
To activate this mode simply set ``using_server`` to false in the :file:`iodef.xml` file
.. code-block:: xml
<variable id="using_server" type="boolean">false</variable>
and don't allocate any cores to XIOS.
.. note::
Due to the different domain decompositions between XIOS and NEMO,
if the total number of cores is larger than the number of grid points in the ``j`` direction then
the model run will fail.
Adding new diagnostics
======================
If you want to add a NEMO diagnostic to the NEMO code you will need to do the following:
1. Add any necessary code to calculate you new diagnostic in NEMO
2. Send the field to XIOS using ``CALL iom_put( 'field_id', variable )`` where
``field_id`` is a unique id for your new diagnostics and
variable is the fortran variable containing the data.
This should be called at every model timestep regardless of how often you want to output the field.
No time averaging should be done in the model code.
3. If it is computationally expensive to calculate your new diagnostic
you should also use "iom_use" to determine if it is requested in the current model run.
For example,
.. code-block:: fortran
IF iom_use('field_id') THEN
!Some expensive computation
!...
!...
iom_put('field_id', variable)
ENDIF
4. Add a variable definition to the :file:`field_def_nemo-???.xml` file.
5. Add the variable to the :file:`iodef.xml` or :file:`file_definition.xml` file.
cfgs/SHARED/domain_def_nemo.xml
View file @ 4a9f4b18
......@@ -189,12 +189,14 @@
<domain id="grid_F_inner" long_name="grid F inner"/>
<!-- zonal mean grid -->
<domain_group id="gznl">
<domain id="gznl" long_name="gznl"/>
<domain id="ptr" domain_ref="gznl" >
<zoom_domain id="ptr" ibegin="0000" jbegin="0" ni="1" nj="0000" />
</domain>
</domain_group>
<domain id="gznl" long_name="gznl"/>
<domain id="ptr" domain_ref="gznl" >
<zoom_domain id="ptr" ibegin="0000" jbegin="0" ni="1" nj="0000" />
</domain>
<domain id="znl_T" domain_ref="gznl" > <zoom_domain id="znl_T"/> </domain>
<domain id="znl_W" domain_ref="gznl" > <zoom_domain id="znl_W"/> </domain>
<!-- other grids -->
......
cfgs/SHARED/field_def_nemo-ice.xml
View file @ 4a9f4b18
......@@ -50,10 +50,11 @@
<field id="icevpnd" long_name="melt pond volume" standard_name="sea_ice_meltpond_volume" unit="m" />
<field id="icehlid" long_name="melt pond lid depth" standard_name="sea_ice_meltpondlid_depth" unit="m" />
<field id="icevlid" long_name="melt pond lid volume" standard_name="sea_ice_meltpondlid_volume" unit="m" />
<field id="dvpn_mlt" long_name="pond volume tendency due to surface melt" standard_name="sea_ice_pondvolume_tendency_melt" unit="kg/m2/s" />
<field id="dvpn_lid" long_name="pond volume tendency due to exchanges with lid" standard_name="sea_ice_pondvolume_tendency_lids" unit="kg/m2/s" />
<field id="dvpn_rnf" long_name="pond volume tendency due to runoff" standard_name="sea_ice_pondvolume_tendency_runoff" unit="kg/m2/s" />
<field id="dvpn_drn" long_name="pond volume tendency due to drainage" standard_name="sea_ice_pondvolume_tendency_drainage" unit="kg/m2/s" />
<field id="iceepnd" long_name="melt pond effective concentration" standard_name="sea_ice_meltpond_effective_concentration" unit="" />
<field id="dvpn_mlt" long_name="pond volume tendency due to surface melt" standard_name="sea_ice_pondvolume_tendency_melt" unit="cm/d" />
<field id="dvpn_lid" long_name="pond volume tendency due to exchanges with lid" standard_name="sea_ice_pondvolume_tendency_lids" unit="cm/d" />
<field id="dvpn_rnf" long_name="pond volume tendency due to runoff" standard_name="sea_ice_pondvolume_tendency_runoff" unit="cm/d" />
<field id="dvpn_drn" long_name="pond volume tendency due to drainage" standard_name="sea_ice_pondvolume_tendency_drainage" unit="cm/d" />
<!-- heat -->
<field id="icetemp" long_name="Mean ice temperature" unit="degC" detect_missing_value="true" />
......@@ -95,6 +96,10 @@
<field id="yield12" long_name="yield surface tensor component 12" standard_name="yield12" unit="N/m" />
<field id="beta_evp" long_name="Relaxation parameter of ice rheology (beta)" standard_name="relaxation_parameter_of_ice_rheology" unit="" />
<field id="isig1" long_name="1st principal stress component for EVP rhg" unit="" />
<field id="isig2" long_name="2nd principal stress component for EVP rhg" unit="" />
<field id="isig3" long_name="convergence measure for EVP rheology (must be around 1)" unit="" />
<!-- surface heat fluxes -->
<field id="qt_ice" long_name="total heat flux at ice surface" standard_name="surface_downward_heat_flux_in_air" unit="W/m2" />
<field id="qsr_ice" long_name="solar heat flux at ice surface" standard_name="surface_downwelling_shortwave_flux_in_air" unit="W/m2" />
......@@ -189,8 +194,15 @@
<field id="icedrift_heat" long_name="Ice heat drift (conservation check)" unit="W/m2" />
<!-- sbcssm variables -->
<field id="sst_m" unit="degC" />
<field id="sss_m" unit="psu" />
<field id="sst_m_pot" unit="degC" />
<!-- EOS-80 -->
<field id="sss_m_pra" unit="psu" />
<!-- TEOS-10 -->
<field id="sss_m_abs" unit="g/kg" />
<!-- SEOS -->
<field id="sss_m_seos" unit="psu" />
<field id="ssu_m" unit="m/s" />
<field id="ssv_m" unit="m/s" />
<field id="ssh_m" unit="m" />
......@@ -407,9 +419,17 @@
<field field_ref="iceapnd" name="siapnd" />
<field field_ref="icehpnd" name="sihpnd" />
<field field_ref="icevpnd" name="sivpnd" />
<field field_ref="iceepnd" name="siepnd" />
<field field_ref="iceage" name="siage" />
<field field_ref="sst_m" name="sst_m" />
<field field_ref="sss_m" name="sss_m" />
<field id="sst_m_pot" unit="degC" />
<!-- EOS-80 -->
<field id="sss_m_pra" unit="psu" />
<!-- TEOS-10 -->
<field id="sss_m_abs" unit="g/kg" />
<!-- SEOS -->
<field id="sss_m_seos" unit="psu" />
<!-- heat -->
<field field_ref="icetemp" name="sitemp" />
......@@ -435,7 +455,7 @@
<field field_ref="sheastr" name="sheastr" />
<field field_ref="sig1_pnorm" name="sig1_pnorm"/>
<field field_ref="sig2_pnorm" name="sig2_pnorm"/>
<field field_ref="icedlt" name="sidelt" />
<field field_ref="icedlt" name="sidelta" />
<!-- heat fluxes -->
<field field_ref="qt_oce_ai" name="qt_oce_ai" />
......
cfgs/SHARED/field_def_nemo-oce.xml
View file @ 4a9f4b18
......@@ -109,19 +109,35 @@ that are available in the tidal-forcing implementation (see
<!-- T grid -->
<field_group id="grid_T" grid_ref="grid_T_2D" >
<field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D" />
<field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D" />
<field id="e3ts" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC"/>
<field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D" />
<field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/>
<field id="e3t_300" field_ref="e3t" grid_ref="grid_T_zoom_300" detect_missing_value="true" />
<field id="e3t_vsum300" field_ref="e3t_300" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="masscello" long_name="Sea Water Mass per unit area" standard_name="sea_water_mass_per_unit_area" unit="kg/m2" grid_ref="grid_T_3D"/>
<field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/>
<field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D" />
<field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field >
<field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/>
<field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field >
<!-- EOS80 -->
<field id="toce_pot" long_name="potential temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="toce_pot_e3t" long_name="potential temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce_pot * e3t </field >
<field id="soce_pra" long_name="practical salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/>
<field id="soce_pra_e3t" long_name="practical salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce_pra * e3t </field >
<!-- TEOS10 -->
<field id="toce_con" long_name="conservative temperature" standard_name="sea_water_conservative_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="toce_con_e3t" long_name="conservative temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce_con * e3t </field >
<field id="soce_abs" long_name="absolute salinity" standard_name="sea_water_absolute_salinity" unit="g/kg" grid_ref="grid_T_3D"/>
<field id="soce_abs_e3t" long_name="absolute salinity (thickness weighted)" unit="g/kg" grid_ref="grid_T_3D" > soce_abs * e3t </field >
<!-- SEOS -->
<field id="toce_seos" long_name="temperature" standard_name="sea_water_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="toce_seos_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce_seos * e3t </field >
<field id="soce_seos" long_name="salinity" standard_name="sea_water_salinity" unit="1e-3" grid_ref="grid_T_3D"/>
<field id="soce_seos_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce_seos * e3t </field >
<field id="toce_e3t_300" field_ref="toce_e3t" unit="degree_C" grid_ref="grid_T_zoom_300" detect_missing_value="true" />
<field id="toce_e3t_vsum300" field_ref="toce_e3t_300" unit="degress_C*m" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="toce_vmean300" field_ref="toce_e3t_vsum300" unit="degree_C" grid_ref="grid_T_vsum" detect_missing_value="true" > toce_e3t_vsum300/e3t_vsum300 </field>
......@@ -159,8 +175,72 @@ that are available in the tidal-forcing implementation (see
<field id="sbs" long_name="sea bottom salinity" unit="0.001" grid_ref="grid_T_2D_inner" />
<field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" />
<field id="taubot" long_name="bottom stress module" unit="N/m2" grid_ref="grid_T_2D_inner" />
<!-- EOS80 -->
<field id="sst_pot" long_name="sea surface potential temperature" standard_name="sea_surface_temperature" unit="degC" />
<field id="sst2_pot" long_name="square of sea surface potential temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst_pot * sst_pot </field >
<field id="sstmax_pot" long_name="max of sea surface potential temperature" field_ref="sst_pot" operation="maximum" />
<field id="sstmin_pot" long_name="min of sea surface potential temperature" field_ref="sst_pot" operation="minimum" />
<field id="sstgrad_pot" long_name="module of potential sst gradient" unit="degC/m" grid_ref="grid_T_2D_inner" />
<field id="sstgrad2_pot" long_name="square of module of potential sst gradient" unit="degC2/m2" grid_ref="grid_T_2D_inner" />
<field id="sbt_pot" long_name="sea bottom potential temperature" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="tosmint_pot" long_name="vertical integral of potential temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner" />
<field id="sst_wl_pot" long_name="Delta potential SST of warm layer" unit="degC" />
<field id="sst_cs_pot" long_name="Delta potential SST of cool skin" unit="degC" />
<field id="temp_3m_pot" long_name="potential temperature at 3m" unit="degC" />
<field id="sss_pra" long_name="sea surface practical salinity" standard_name="sea_surface_practical_salinity" unit="1e-3" />
<field id="sss2_pra" long_name="square of sea surface practical salinity" unit="1e-6" > sss_pra * sss_pra </field >
<field id="sssmax_pra" long_name="max of sea surface practical salinity" field_ref="sss_pra" operation="maximum" />
<field id="sssmin_pra" long_name="min of sea surface practical salinity" field_ref="sss_pra" operation="minimum" />
<field id="sbs_pra" long_name="sea bottom practical salinity" unit="0.001" grid_ref="grid_T_2D_inner" />
<field id="somint_pra" long_name="vertical integral of practical salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_practical_salinity" unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" />
<!-- TEOS10 -->
<field id="sst_con" long_name="sea surface conservative temperature" standard_name="sea_surface_conservative_temperature" unit="degC" />
<field id="sst2_con" long_name="square of sea surface conservative temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst_con * sst_con </field >
<field id="sstmax_con" long_name="max of sea surface conservative temperature" field_ref="sst_con" operation="maximum" />
<field id="sstmin_con" long_name="min of sea surface conservative temperature" field_ref="sst_con" operation="minimum" />
<field id="sstgrad_con" long_name="module of conservative sst gradient" unit="degC/m" grid_ref="grid_T_2D_inner" />
<field id="sstgrad2_con" long_name="square of module of conservative sst gradient" unit="degC2/m2" grid_ref="grid_T_2D_inner" />
<field id="sbt_con" long_name="sea bottom conservative temperature" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="tosmint_con" long_name="vertical integral of conservative temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_conservative_temperature" unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner" />
<field id="sst_wl_con" long_name="Delta conservative SST of warm layer" unit="degC" />
<field id="sst_cs_con" long_name="Delta conservative SST of cool skin" unit="degC" />
<field id="temp_3m_con" long_name="conservative temperature at 3m" unit="degC" />
<field id="sss_abs" long_name="sea surface absolute salinity" standard_name="sea_surface_absolute_salinity" unit="g/kg" />
<field id="sss2_abs" long_name="square of sea surface absolute salinity" unit="1e-6" > sss_abs * sss_abs </field >
<field id="sssmax_abs" long_name="max of sea surface absolute salinity" field_ref="sss_abs" operation="maximum" />
<field id="sssmin_abs" long_name="min of sea surface absolute salinity" field_ref="sss_abs" operation="minimum" />
<field id="sbs_abs" long_name="sea bottom absolute salinity" unit="g/kg" grid_ref="grid_T_2D_inner" />
<field id="somint_abs" long_name="vertical integral of absolute salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_absolute_salinity" unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" />
<!-- SEOS -->
<field id="sst_seos" long_name="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" />
<field id="sst2_seos" long_name="square of sea surface temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst_seos * sst_seos </field >
<field id="sstmax_seos" long_name="max of sea surface temperature" field_ref="sst_seos" operation="maximum" />
<field id="sstmin_seos" long_name="min of sea surface temperature" field_ref="sst_seos" operation="minimum" />
<field id="sstgrad_seos" long_name="module of sst gradient" unit="degC/m" grid_ref="grid_T_2D_inner" />
<field id="sstgrad2_seos" long_name="square of module of sst gradient" unit="degC2/m2" grid_ref="grid_T_2D_inner" />
<field id="sbt_seos" long_name="sea bottom temperature" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="tosmint_seos" long_name="vertical integral of temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_temperature" unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner" />
<field id="sst_wl_seos" long_name="Delta SST of warm layer" unit="degC" />
<field id="sst_cs_seos" long_name="Delta SST of cool skin" unit="degC" />
<field id="temp_3m_seos" long_name="temperature at 3m" unit="degC" />
<field id="sss_seos" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" />
<field id="sss2_seos" long_name="square of sea surface salinity" unit="1e-6" > sss_seos * sss_seos </field >
<field id="sssmax_seos" long_name="max of sea surface salinity" field_ref="sss_seos" operation="maximum" />
<field id="sssmin_seos" long_name="min of sea surface salinity" field_ref="sss_seos" operation="minimum" />
<field id="sbs_seos" long_name="sea bottom salinity" unit="0.001" grid_ref="grid_T_2D_inner" />
<field id="somint_seos" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" />
<field id="taubot" long_name="bottom stress module" unit="N/m2" grid_ref="grid_T_2D_inner" />
<field id="htau" long_name="htau length scale" unit="m"/>
<!-- Case EOS = TEOS-10 : output potential temperature -->
<field id="toce_pot" long_name="Sea Water Potential Temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="sst_pot" long_name="potential sea surface temperature" standard_name="sea_surface_temperature" unit="degC" />
......@@ -176,16 +256,34 @@ that are available in the tidal-forcing implementation (see
<field id="mldr10_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
<field id="mldr10_1max" long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="maximum" />
<field id="mldr10_1min" long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="minimum" />
<field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D_inner" />
<field id="saltc" long_name="Salt content vertically integrated" unit="PSU*kg/m2" grid_ref="grid_T_2D_inner" />
<field id="salt2c" long_name="square of Salt content vertically integrated" unit="PSU2*kg/m2" grid_ref="grid_T_2D_inner" />
<field id="mldzint_1" long_name="Mixed Layer Depth interpolated" standard_name="ocean_mixed_layer_thickness" unit="m" grid_ref="grid_T_2D" />
<field id="mldzint_2" long_name="Mixed Layer Depth interpolated" standard_name="ocean_mixed_layer_thickness" unit="m" grid_ref="grid_T_2D" />
<field id="mldzint_3" long_name="Mixed Layer Depth interpolated" standard_name="ocean_mixed_layer_thickness" unit="m" grid_ref="grid_T_2D" />
<field id="mldzint_4" long_name="Mixed Layer Depth interpolated" standard_name="ocean_mixed_layer_thickness" unit="m" grid_ref="grid_T_2D" />
<field id="mldzint_5" long_name="Mixed Layer Depth interpolated" standard_name="ocean_mixed_layer_thickness" unit="m" grid_ref="grid_T_2D" />
<field id="mldhtc_1" long_name="Mixed Layer Depth integrated heat content" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D" />
<field id="mldhtc_2" long_name="Mixed Layer Depth integrated heat content" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D" />
<field id="mldhtc_3" long_name="Mixed Layer Depth integrated heat content" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D" />
<field id="mldhtc_4" long_name="Mixed Layer Depth integrated heat content" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D" />
<field id="mldhtc_5" long_name="Mixed Layer Depth integrated heat content" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D" />
<field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D_inner" />
<field id="saltc" long_name="Salt content vertically integrated" unit="1e-3.m" grid_ref="grid_T_2D_inner" />
<field id="salt2c" long_name="square of Salt content vertically integrated" unit="1e-6.m2" grid_ref="grid_T_2D_inner" />
<!-- EOS -->
<field id="alpha" long_name="thermal expansion" unit="degC-1" grid_ref="grid_T_3D" />
<field id="beta" long_name="haline contraction" unit="1e3" grid_ref="grid_T_3D" />
<field id="rhop" long_name="potential density (sigma0)" standard_name="sea_water_sigma_theta" unit="kg/m3" grid_ref="grid_T_3D" />
<field id="rhop_e3t" long_name="potential density (sigma0) thickness weighted" standard_name="sea_water_sigma_theta" unit="kg/m3" grid_ref="grid_T_3D"> rhop * e3t </field>
<field id="N_2d" long_name="Depth-mean N" unit="m/s" />
<field id="modslp" long_name="sqrt( slpi^2 + slpj^2 )" unit="1" grid_ref="grid_T_3D" />
<field id="RossRad" long_name="Rossby radius" unit="m" />
<field id="RossRadlim" long_name="Rossby radius (limited)" unit="m" />
<field id="Tclinic_recip" long_name="recip of baroclinic timescale" unit="s-1" />
<field id="RR_GS" long_name="Rossby radius / min(dx,dy)" unit="1" />
<!-- Energy - horizontal divergence -->
<field id="eken" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" />
<field id="sKE" long_name="surface kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_2D_inner" />
<field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D" />
......@@ -263,6 +361,10 @@ that are available in the tidal-forcing implementation (see
<field id="qhcisf3d_par" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qconisf" long_name="Conductive heat flux through the ice shelf ( from isf to oce )" unit="W/m2" />
<field id="qlatisf" long_name="Ice shelf latent heat flux" unit="W/m2" />
<field id="qhcisf" long_name="Ice shelf heat content flux" unit="W/m2" />
<field id="fwfisf" long_name="Ice shelf melting" unit="kg/m2/s" />
<!-- top boundary layer properties -->
<field id="isftfrz_cav" long_name="freezing point temperature at ocean/isf interface" unit="degC" />
<field id="isftfrz_par" long_name="freezing point temperature in the parametrization boundary layer" unit="degC" />
......@@ -330,7 +432,12 @@ that are available in the tidal-forcing implementation (see
<field id="zustar" long_name="friction velocity" unit="m/s" />
<field id="zwstrl" long_name="langmuir velocity scale" unit="m/s" />
<field id="zvstr" long_name="mixed velocity scale" unit="m/s" />
<field id="zla" long_name="langmuir number" unit="m/s" />
<field id="zla" long_name="langmuir number" unit="#" />
<field id="wmp" long_name="wave mean period" unit="s" />
<field id="hsw" long_name="significant wave height" unit="m" />
<field id="wmp_NP" long_name="wave mean period from NP spectrum" unit="s" />
<field id="hsw_NP" long_name="significant wave height from NP spectrum" unit="m" />
<field id="wndm" long_name="U_10" unit="m/s" />
<field id="wind_wave_abs_power" long_name="\rho |U_s| x u*^2" unit="mW" />
<field id="wind_wave_power" long_name="U_s \dot tau" unit="mW" />
<field id="wind_power" long_name="\rho u*^3" unit="mW" />
......@@ -426,7 +533,7 @@ that are available in the tidal-forcing implementation (see
<field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" />
<field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" />
<field id="qsr" long_name="Shortwave Radiation" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" />
<field id="qsr3d" long_name="Shortwave Radiation 3D distribution" standard_name="downwelling_shortwave_flux_in_sea_water" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qsr3d" long_name="Shortwave Radiation 3D distribution" standard_name="downwelling_shortwave_flux_in_sea_water" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qrp" long_name="Surface Heat Flux: Damping" standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation" unit="W/m2" />
<field id="qclosea" long_name="closed sea heat content flux" standard_name="closea_heat_content_downward_flux" unit="W/m2" />
<field id="erp" long_name="Surface Water Flux: Damping" standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation" unit="kg/m2/s" />
......@@ -521,6 +628,15 @@ that are available in the tidal-forcing implementation (see
<!-- sbcssm variables -->
<field id="sst_m" unit="degC" />
<field id="sss_m" unit="psu" />
<field id="sst_m_pot" unit="degC" />
<!-- EOS-80 -->
<field id="sss_m_pra" unit="psu" />
<!-- TEOS-10 -->
<field id="sss_m_abs" unit="g/kg" />
<!-- SEOS -->
<field id="sss_m_seos" unit="psu" />
<field id="ssu_m" unit="m/s" />
<field id="ssv_m" unit="m/s" />
<field id="ssh_m" unit="m" />
......@@ -592,8 +708,8 @@ that are available in the tidal-forcing implementation (see
<field id="uocetr_vsum_op" long_name="ocean current along i-axis * e3u * e2u summed on the vertical" read_access="true" freq_op="1mo" field_ref="e2u" unit="m3/s"> @uocetr_vsum </field>
<field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" />
<field id="msftbarot" long_name="ocean_barotropic_mass_streamfunction" unit="kg s-1" > uocetr_vsum_cumul * $rho0 </field>
<field id="uoce2_e3u" long_name="ocean current along i-axis squared (thickness weighted)" unit="m3/s2" grid_ref="grid_U_3D" > uoce * uoce * e3u </field>
<field id="ssu" long_name="ocean surface current along i-axis" unit="m/s" />
<field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" grid_ref="grid_U_2D_inner" />
<field id="ubar" long_name="ocean barotropic current along i-axis" unit="m/s" />
......@@ -654,6 +770,8 @@ that are available in the tidal-forcing implementation (see
<field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" />
<field id="voce" long_name="ocean current along j-axis" standard_name="sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" />
<field id="voce_e3v" long_name="ocean current along j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > voce * e3v </field>
<field id="voce2_e3v" long_name="ocean current along j-axis squared (thickness weighted)" unit="m3/s2" grid_ref="grid_V_3D" > voce * voce * e3v </field>
<field id="ssv" long_name="ocean surface current along j-axis" unit="m/s" />
<field id="sbv" long_name="ocean bottom current along j-axis" unit="m/s" grid_ref="grid_V_2D_inner" />
<field id="vbar" long_name="ocean barotropic current along j-axis" unit="m/s" />
......@@ -738,6 +856,8 @@ that are available in the tidal-forcing implementation (see
<field id="mf_app" long_name="convective area" standard_name="mf_convective_area" unit="%" grid_ref="grid_W_3D" />
<field id="mf_wp" long_name="convective velocity" standard_name="mf_convective_velo" unit="m/s" grid_ref="grid_W_3D" />
<!-- av_tmx: available with ln_zdftmx -->
<field id="av_tmx" long_name="vertical diffusivity due to tidal mixing" unit="m2/s" />
<!-- avt_tide: available with ln_zdfiwm=T -->
<field id="av_ratio" long_name="S over T diffusivity ratio" standard_name="salinity_over_temperature_diffusivity_ratio" unit="1" />
......@@ -760,7 +880,7 @@ that are available in the tidal-forcing implementation (see
<field id="avt_k" long_name="vertical eddy diffusivity from closure schemes" standard_name="ocean_vertical_eddy_diffusivity" unit="m2/s" />
<field id="avm_k" long_name="vertical eddy viscosity from closure schemes" standard_name="ocean_vertical_eddy_viscosity" unit="m2/s" />
<field id="ediss_k" long_name="Kolmogorov energy dissipation (tke scheme)" standard_name="Kolmogorov_energy_dissipation" unit="W/kg" />
<field id="eshear_k" long_name="energy source from vertical shear" standard_name="energy_source_from_shear" unit="W/kg" />
<field id="eshear_k" long_name="energy source from vertical shear" standard_name="energy_source_from_shear" unit="W/kg" grid_ref="grid_W_3D_inner" />
<field id="estrat_k" long_name="energy sink from stratification" standard_name="energy_sink_from_stratification" unit="W/kg" />
</field_group>
......@@ -783,6 +903,31 @@ that are available in the tidal-forcing implementation (see
<!-- f-eddy viscosity coefficients (ldfdyn) -->
<field id="ahmf_2d" long_name=" surface f-eddy viscosity coefficient" unit="m2/s or m4/s" />
<field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/>
<!-- product fields -->
<field_group id="diaprod">
<field id="ut" long_name="product_of_sea_water_x_velocity_and_potential_temperature" unit="degree_C m/s" grid_ref="grid_U_3D" />
<field id="ut_e3u" long_name="product_of_sea_water_x_velocity_and_potential_temperature * e3u" unit="degree_C m2/s" grid_ref="grid_U_3D" > ut * e3u </field >
<field id="us" long_name="product_of_sea_water_x_velocity_and_salinity" unit="PSU m/s" grid_ref="grid_U_3D" />
<field id="us_e3u" long_name="product_of_sea_water_x_velocity_and_salinity * e3u" unit="PSU m2/s" grid_ref="grid_U_3D" > us * e3u </field >
<field id="urhop" long_name="product_of_sea_water_x_velocity_and_potential_density" unit="(kg/m3).(m/s)" grid_ref="grid_U_3D" />
<field id="urhop_e3u" long_name="product_of_sea_water_x_velocity_and_potential_density * e3u" unit="(kg/m3).(m2/s)" grid_ref="grid_U_3D" > urhop * e3u </field >
<field id="vt" long_name="product_of_sea_water_y_velocity_and_potential_temperature" unit="degree_C m/s" grid_ref="grid_V_3D" />
<field id="vt_e3v" long_name="product_of_sea_water_y_velocity_and_potential_temperature * e3v" unit="degree_C m2/s" grid_ref="grid_V_3D" > vt * e3v </field >
<field id="vs" long_name="product_of_sea_water_y_velocity_and_salinity" unit="PSU m/s" grid_ref="grid_V_3D" />
<field id="vs_e3v" long_name="product_of_sea_water_y_velocity_and_salinity * e3t" unit="PSU m2/s" grid_ref="grid_V_3D" > vs * e3v </field >
<field id="vrhop" long_name="product_of_sea_water_y_velocity_and_potential_density" unit="(kg/m3).(m/s)" grid_ref="grid_V_3D" />
<field id="vrhop_e3v" long_name="product_of_sea_water_y_velocity_and_potential_density * e3t" unit="(kg/m3).(m2/s)" grid_ref="grid_V_3D" > vrhop * e3v </field >
<field id="wt" long_name="product_of_upward_sea_water_velocity_and_potential_temperature" unit="degree_C m/s" grid_ref="grid_W_3D" />
<field id="ws" long_name="product_of_upward_sea_water_velocity_and_salinity" unit="PSU m/s" grid_ref="grid_W_3D" />
<field id="wrhop" long_name="product_of_upward_sea_water_velocity_and_potential_density" unit="(kg/m3).(m/s)" grid_ref="grid_W_3D" />
<field id="uv" long_name="product_of_sea_water_x_velocity_and_sea_water_y_velocity" unit="m2/s2 " grid_ref="grid_T_3D" />
<field id="uw" long_name="product_of_upward_sea_water_velocity_and_sea_water_x_velocity" unit="m2/s2 " grid_ref="grid_W_3D" />
<field id="vw" long_name="product_of_upward_sea_water_velocity_and_sea_water_y_velocity" unit="m2/s2" grid_ref="grid_W_3D" />
</field_group>
<field_group id="scalar" grid_ref="grid_scalar" >
<!-- Need to have a salinity reference climatological file : sali_ref_clim_monthly -->
......@@ -1147,6 +1292,17 @@ that are available in the tidal-forcing implementation (see
<field_group id="mooring" >
<field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" />
<field field_ref="soce" name="so" long_name="sea_water_salinity" />
<!-- EOS80 -->
<field field_ref="toce_pot" name="thetao_pot" long_name="sea_water_potential_temperature" />
<field field_ref="soce_pra" name="so_pra" long_name="sea_water_practical_salinity" />
<!-- TEOS10 -->
<field field_ref="toce_con" name="thetao_con" long_name="sea_water_conservative_temperature" />
<field field_ref="soce_abs" name="so_con" long_name="sea_water_absolute_salinity" />
<!-- SEOS -->
<field field_ref="toce_seos" name="thetao_seos" long_name="sea_water_temperature" />
<field field_ref="soce_seos" name="so_seos" long_name="sea_water_salinity" />
<field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" />
<field field_ref="voce" name="vo" long_name="sea_water_y_velocity" />
<field field_ref="woce" name="wo" long_name="sea_water_z_velocity" />
......@@ -1157,6 +1313,22 @@ that are available in the tidal-forcing implementation (see
<field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" />
<field field_ref="sstgrad" name="tosgrad" long_name="module_of_sea_surface_temperature_gradient" />
<field field_ref="sss" name="sos" long_name="sea_surface_salinity" />
<!-- EOS80 -->
<field field_ref="sst_pot" name="tos_pot" long_name="sea_surface_potential_temperature" />
<field field_ref="sst2_pot" name="tossq_pot" long_name="square_of_sea_surface_potential_temperature" />
<field field_ref="sstgrad_pot" name="tosgrad_pot" long_name="module_of_sea_surface_potential_temperature_gradient" />
<field field_ref="sss_pra" name="sos_pra" long_name="sea_surface_absolute_salinity" />
<!-- TEOS10 -->
<field field_ref="sst_con" name="tos_con" long_name="sea_surface_conservative_temperature" />
<field field_ref="sst2_con" name="tossq_con" long_name="square_of_sea_surface_conservative_temperature" /> <field field_ref="sstgrad_con" name="tosgrad_con" long_name="module_of_sea_surface_conservative_temperature_gradient" />
<field field_ref="sss_abs" name="sos_abs" long_name="sea_surface_absolute_salinity" />
<!-- SEOS -->
<field field_ref="sst_seos" name="tos_seos" long_name="sea_surface_temperature" />
<field field_ref="sst2_seos" name="tossq_seos" long_name="square_of_sea_surface_temperature" />
<field field_ref="sstgrad_seos" name="tosgrad_seos" long_name="module_of_sea_surface_temperature_gradient" />
<field field_ref="sss_seos" name="sos_seos" long_name="sea_surface_salinity" />
<field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" />
<field field_ref="empmr" name="wfo" long_name="water_flux_into_sea_water" />
<field field_ref="qsr" name="rsntds" long_name="surface_net_downward_shortwave_flux" />
......@@ -1184,6 +1356,27 @@ that are available in the tidal-forcing implementation (see
<field field_ref="sst" name="tos" long_name="sea_surface_temperature" />
<field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" />
<field field_ref="sss" name="sos" long_name="sea_surface_salinity" />
<!-- EOS80 -->
<field id="toce_pot" long_name="sea_water_potential_temperature" grid_ref="grid_T_3D" />
<field id="soce_pra" long_name="sea_water_practical_salinity" grid_ref="grid_T_3D" />
<field id="sst_pot" long_name="sea_surface_potential_temperature" grid_ref="grid_T_2D" />
<field id="sst2_pot" long_name="square_of_sea_surface_potential_temperature" grid_ref="grid_T_2D" />
<field id="sss_pra" long_name="sea_surface_practical_salinity" grid_ref="grid_T_2D" />
<!-- TEOS10 -->
<field id="toce_con" long_name="sea_water_conservative_temperature" grid_ref="grid_T_3D" />
<field id="soce_abs" long_name="sea_water_absolute_salinity" grid_ref="grid_T_3D" />
<field id="sst_con" long_name="sea_surface_conservative_temperature" grid_ref="grid_T_2D" unit="degC" />
<field id="sst2_con" long_name="square_of_sea_surface_conservative_temperature" grid_ref="grid_T_2D" />
<field id="sss_abs" long_name="sea_surface_absolute_salinity" grid_ref="grid_T_2D" />
<!-- SEOS -->
<field id="toce_seos" long_name="sea_water_temperature" grid_ref="grid_T_3D" />
<field id="soce_seos" long_name="sea_water_salinity" grid_ref="grid_T_3D" />
<field id="sst_seos" long_name="sea_surface_temperature" grid_ref="grid_T_2D" />
<field id="sst2_seos" long_name="square_of_sea_surface_temperature" grid_ref="grid_T_2D" />
<field id="sss_seos" long_name="sea_surface_salinity" grid_ref="grid_T_2D" />
<field field_ref="rhop" name="sigma0" long_name="potential density" />
<field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" />
<field field_ref="empmr" name="wfo" long_name="water_flux_into_sea_water" />
<field field_ref="qsr" name="rsntds" long_name="surface_net_downward_shortwave_flux" />
......@@ -1217,6 +1410,30 @@ that are available in the tidal-forcing implementation (see
<!-- TMB diagnostic output -->
<field_group id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant">
<!-- EOS80 -->
<field id="top_temp_pot" name="votemper_top_pot" unit="degC" />
<field id="mid_temp_pot" name="votemper_mid_pot" unit="degC" />
<field id="bot_temp_pot" name="votemper_bot_pot" unit="degC" />
<field id="top_sal_pra" name="vosaline_top_pra" unit="psu" />
<field id="mid_sal_pra" name="vosaline_mid_pra" unit="psu" />
<field id="bot_sal_pra" name="vosaline_bot_pra" unit="psu" />
<!-- TEOS10 -->
<field id="top_temp_con" name="votemper_top_con" unit="degC" />
<field id="mid_temp_con" name="votemper_mid_con" unit="degC" />
<field id="bot_temp_con" name="votemper_bot_con" unit="degC" />
<field id="top_sal_abs" name="vosaline_top_abs" unit="g/kg" />
<field id="mid_sal_abs" name="vosaline_mid_abs" unit="g/kg" />
<field id="bot_sal_abs" name="vosaline_bot_abs" unit="g/kg" />
<!-- SEOS -->
<field id="top_temp_seos" name="votemper_top_seos" unit="degC" />
<field id="mid_temp_seos" name="votemper_mid_seos" unit="degC" />
<field id="bot_temp_seos" name="votemper_bot_seos" unit="degC" />
<field id="top_sal_seos" name="vosaline_top_seos" unit="psu" />
<field id="mid_sal_seos" name="vosaline_mid_seos" unit="psu" />
<field id="bot_sal_seos" name="vosaline_bot_seos" unit="psu" />
<field id="top_temp" name="votemper_top" unit="degC" />
<field id="mid_temp" name="votemper_mid" unit="degC" />
<field id="bot_temp" name="votemper_bot" unit="degC" />
......@@ -1246,6 +1463,16 @@ that are available in the tidal-forcing implementation (see
<field id="tempis25h" name="insitu temperature 25h mean" unit="degC" />
<field id="salin25h" name="salinity 25h mean" unit="psu" />
<field id="ssh25h" name="sea surface height 25h mean" grid_ref="grid_T_2D_inner" unit="m" />
<!-- EOS80 -->
<field id="temper25h_pot" name="potential temperature 25h mean" unit="degC" />
<field id="salin25h_pra" name="practical salinity 25h mean" unit="psu" />
<!-- TEOS10 -->
<field id="temper25h_con" name="conservative temperature 25h mean" unit="degC" />
<field id="salin25h_abs" name="absolute salinity 25h mean" unit="g/kg" />
<!-- SEOS -->
<field id="temper25h_seos" name="temperature 25h mean" unit="degC" />
<field id="salin25h_seos" name="salinity 25h mean" unit="psu" />
</field_group>
<field_group id="25h_grid_U" grid_ref="grid_U_3D_inner" operation="instant" >
......
cfgs/SHARED/field_def_nemo-pisces.xml
View file @ 4a9f4b18
......@@ -183,6 +183,9 @@
<field id="PPNEWN" long_name="New Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPNEWP" long_name="New Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPNEWD" long_name="New Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="GPPHYN" long_name="Gross Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="GPPHYP" long_name="Gross Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="GPPHYD" long_name="Gross Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PBSi" long_name="Primary production of Si diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PFeN" long_name="Primary production of nano iron" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PFeP" long_name="Primary production of pico iron" unit="molC/m3/s" grid_ref="grid_T_3D" />
......@@ -276,20 +279,31 @@
<!-- dbio_T on T grid : variables available with diaar5 -->
<field id="TPP" long_name="Total Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D" />
<field id="TPNEW" long_name="New Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D" />
<field id="TPBFE" long_name="Total biogenic iron production" unit="mol/m3/s" grid_ref="grid_T_3D" />
<field id="INTDIC" long_name="DIC content" unit="kg/m2" />
<field id="O2MIN" long_name="Oxygen minimum concentration" unit="mol/m3" />
<field id="ZO2MIN" long_name="Depth of oxygen minimum concentration" unit="m" />
<field id="INTNFIX" long_name="Nitrogen fixation rate : vert. integrated" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > Nfix * e3t </field >
<field id="INTPPPHYN" long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PPPHYN * e3t </field >
<field id="INTPPPHYD" long_name="Vertically integrated primary production by diatom" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PPPHYD * e3t </field >
<field id="INTPPPHYP" long_name="Vertically integrated primary production by picophy" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PPPHYP * e3t </field >
<field id="INTPP" long_name="Vertically integrated primary production by phyto" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > TPP * e3t </field >
<field id="INTPNEW" long_name="Vertically integrated new primary production" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > TPNEW * e3t </field >
<field id="INTPBFE" long_name="Vertically integrated of biogenic iron production" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > TPBFE * e3t </field >
<field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PBSi * e3t </field >
<field id="TPP" long_name="Total Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D" />
<field id="TPNEW" long_name="New Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D" />
<field id="TPBFE" long_name="Total biogenic iron production" unit="mol/m3/s" grid_ref="grid_T_3D" />
<field id="INTDIC" long_name="DIC content" unit="kg/m2" />
<field id="O2MIN" long_name="Oxygen minimum concentration" unit="mol/m3" />
<field id="ZO2MIN" long_name="Depth of oxygen minimum concentration" unit="m" />
<field id="Nfix_e3t" long_name="Nfix * e3t" unit="molN/m2/s" grid_ref="grid_T_3D" > Nfix * e3t </field >
<field id="PPPHYN_e3t" long_name="PPPHYN * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > PPPHYN * e3t </field >
<field id="PPPHYD_e3t" long_name="PPPHYD * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > PPPHYD * e3t </field >
<field id="PPPHYP_e3t" long_name="PPPHYP * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > PPPHYP * e3t </field >
<field id="PBSi_e3t" long_name="PBSi * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > PBSi * e3t </field >
<field id="TPP_e3t" long_name="TPP * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > TPP * e3t </field >
<field id="TPNEW_e3t" long_name="TPNEW * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > TPNEW * e3t </field >
<field id="TPBFE_e3t" long_name="TPBFE * e3t" unit="mol/m2/s" grid_ref="grid_T_3D" > TPBFE * e3t </field >
<field id="INTNFIX" long_name="Nitrogen fixation rate : vert. integrated" field_ref="Nfix_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPPPHYN" long_name="Vertically integrated primary production by nanophy" field_ref="PPPHYN_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPPPHYD" long_name="Vertically integrated primary production by diatom" field_ref="PPPHYD_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPPPHYP" long_name="Vertically integrated primary production by picophy" field_ref="PPPHYP_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPP" long_name="Vertically integrated primary production by phyto" field_ref="TPP_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPNEW" long_name="Vertically integrated new primary production" field_ref="TPNEW_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPBFE" long_name="Vertically integrated of biogenic iron production" field_ref="TPBFE_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" field_ref="PBSi_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<!-- PISCES light : variables available with key_pisces_reduced -->
<field id="FNO3PHY" long_name="FNO3PHY" unit="" grid_ref="grid_T_3D" />
......
cfgs/SHARED/grid_def_nemo.xml
View file @ 4a9f4b18
......@@ -418,5 +418,43 @@
<duplicate_scalar />
</axis>
</grid>
<grid id="grid_EqT" >
<domain id="EqT" />
</grid>
<!-- -->
<grid id="gznl_T_2D">
<domain id="ptr" />
</grid>
<!-- -->
<grid id="gznl_T_3D">
<domain id="ptr" />
<axis axis_ref="deptht" />
</grid>
<!-- -->
<grid id="gznl_W_2D">
<domain id="ptr" />
</grid>
<!-- -->
<grid id="gznl_W_3D">
<domain id="ptr" />
<axis axis_ref="depthw" />
</grid>
<grid id="vert_sum">
<domain id="grid_T"/>
<scalar>
<reduce_axis operation="sum" />
</scalar>
</grid>
<grid id="zoom_300">
<domain id="grid_T" />
<axis axis_ref="deptht300"/>
</grid>
<grid id="zoom_300_sum">
<domain id="grid_T" />
<scalar>
<reduce_axis operation="sum" />
</scalar>
</grid>
</grid_definition>
cfgs/SHARED/namelist_ice_ref
View file @ 4a9f4b18
......@@ -24,10 +24,10 @@
jpl = 5 ! number of ice categories
nlay_i = 2 ! number of ice layers
nlay_s = 2 ! number of snow layers
ln_virtual_itd = .false. ! virtual ITD mono-category parameterization (jpl=1 only)
! i.e. enhanced thermal conductivity & virtual thin ice melting
ln_icedyn = .true. ! ice dynamics (T) or not (F)
ln_icethd = .true. ! ice thermo (T) or not (F)
ln_virtual_itd = .false., ! virtual ITD mono-category parameterization (jpl=1 only)
! i.e. enhan.false.ced thermal conductivity & virtual thin ice melting
ln_icedyn = .true., ! ice dynamics (T) or not (F)
ln_icethd = .true., ! ice thermo (T) or not (F)
rn_amax_n = 0.997 ! maximum tolerated ice concentration NH
rn_amax_s = 0.997 ! maximum tolerated ice concentration SH
cn_icerst_in = "restart_ice" ! suffix of ice restart name (input)
......@@ -38,9 +38,9 @@
!------------------------------------------------------------------------------
&namitd ! Ice discretization
!------------------------------------------------------------------------------
ln_cat_hfn = .true. ! ice categories are defined by a function following rn_himean**(-0.05)
ln_cat_hfn = .true., ! ice categories are defined by a function following rn_himean**(-0.05)
rn_himean = 2.0 ! expected domain-average ice thickness (m)
ln_cat_usr = .false. ! ice categories are defined by rn_catbnd below (m)
ln_cat_usr = .false., ! ice categories are defined by rn_catbnd below (m)
rn_catbnd = 0.,0.45,1.1,2.1,3.7,6.0
rn_himin = 0.1 ! minimum ice thickness (m) allowed
rn_himax = 99.0 ! maximum ice thickness (m) allowed
......@@ -48,14 +48,14 @@
!------------------------------------------------------------------------------
&namdyn ! Ice dynamics
!------------------------------------------------------------------------------
ln_dynALL = .true. ! dyn.: full ice dynamics (rheology + advection + ridging/rafting + correction)
ln_dynRHGADV = .false. ! dyn.: no ridge/raft & no corrections (rheology + advection)
ln_dynADV1D = .false. ! dyn.: only advection 1D (Schar & Smolarkiewicz 1996 test case)
ln_dynADV2D = .false. ! dyn.: only advection 2D w prescribed vel.(rn_uvice + advection)
ln_dynALL = .true., ! dyn.: full ice dynamics (rheology + advection + ridging/rafting + correction)
ln_dynRHGADV = .false., ! dyn.: no ridge/raft & no corrections (rheology + advection)
ln_dynADV1D = .false., ! dyn.: only advection 1D (Schar & Smolarkiewicz 1996 test case)
ln_dynADV2D = .false., ! dyn.: only advection 2D w prescribed vel.(rn_uvice + advection)
rn_uice = 0.5 ! prescribed ice u-velocity
rn_vice = 0.5 ! prescribed ice v-velocity
rn_ishlat = 2. ! lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2)
ln_landfast_L16 = .false. ! landfast: parameterization from Lemieux 2016
ln_landfast_L16 = .false., ! landfast: parameterization from Lemieux 2016
rn_lf_depfra = 0.125 ! fraction of ocean depth that ice must reach to initiate landfast
! recommended range: [0.1 ; 0.25]
rn_lf_bfr = 15. ! maximum bottom stress per unit volume [N/m3]
......@@ -72,30 +72,30 @@
&namdyn_rdgrft ! Ice ridging/rafting
!------------------------------------------------------------------------------
! -- ice_rdgrft_strength -- !
ln_str_H79 = .true. ! ice strength param.: Hibler_79 => P = pstar*<h>*exp(-c_rhg*A)
ln_str_H79 = .true., ! ice strength param.: Hibler_79 => P = pstar*<h>*exp(-c_rhg*A)
rn_pstar = 2.0e+04 ! ice strength thickness parameter [N/m2]
rn_crhg = 20.0 ! ice strength conc. parameter (-)
ln_str_R75 = .false. ! ice strength param.: Rothrock_75 => P = fn of potential energy
ln_str_R75 = .false., ! ice strength param.: Rothrock_75 => P = fn of potential energy
rn_pe_rdg = 17.0 ! coef accouting for frictional dissipation
ln_str_CST = .false. ! ice strength param.: Constant
ln_str_CST = .false., ! ice strength param.: Constant
rn_str = 0.0 ! ice strength value
ln_str_smooth = .true. ! spatial smoothing of the ice strength
ln_str_smooth = .true., ! spatial smoothing of the ice strength
! -- ice_rdgrft -- !
ln_distf_lin = .true. ! redistribution function of ridged ice: linear (Hibler 1980)
ln_distf_exp = .false. ! redistribution function of ridged ice: exponential => not coded yet
ln_distf_lin = .true., ! redistribution function of ridged ice: linear (Hibler, 1980)
ln_distf_exp = .false., ! redistribution function of ridged ice: exponential (Lipscomb et al., 2007)
rn_murdg = 3.0 ! e-folding scale of ridged ice (m**.5)
rn_csrdg = 0.5 ! fraction of shearing energy contributing to ridging
! -- ice_rdgrft_prep -- !
ln_partf_lin = .false. ! Linear ridging participation function (Thorndike et al, 1975)
ln_partf_lin = .false., ! Linear ridging participation function (Thorndike et al., 1975)
rn_gstar = 0.15 ! fractional area of thin ice being ridged
ln_partf_exp = .true. ! Exponential ridging participation function (Lipscomb, 2007)
ln_partf_exp = .true., ! Exponential ridging participation function (Lipscomb et al., 2007)
rn_astar = 0.03 ! exponential measure of ridging ice fraction [set to 0.05 if hstar=100]
ln_ridging = .true. ! ridging activated (T) or not (F)
ln_ridging = .true., ! ridging activated (T) or not (F)
rn_hstar = 25.0 ! determines the maximum thickness of ridged ice [m] (Hibler, 1980)
rn_porordg = 0.3 ! porosity of newly ridged ice (Lepparanta et al., 1995)
rn_fsnwrdg = 0.5 ! snow volume fraction that survives in ridging
rn_fpndrdg = 1.0 ! pond fraction that survives in ridging (small a priori)
ln_rafting = .true. ! rafting activated (T) or not (F)
ln_rafting = .true., ! rafting activated (T) or not (F)
rn_hraft = 0.75 ! threshold thickness for rafting [m]
rn_craft = 5.0 ! squeezing coefficient used in the rafting function
rn_fsnwrft = 0.5 ! snow volume fraction that survives in rafting
......@@ -104,9 +104,9 @@
!------------------------------------------------------------------------------
&namdyn_rhg ! Ice rheology
!------------------------------------------------------------------------------
ln_rhg_EVP = .true. ! EVP rheology
ln_rhg_EAP = .false. ! EAP rheology
ln_aEVP = .true. ! adaptive rheology (Kimmritz et al. 2016 & 2017)
ln_rhg_EVP = .true., ! EVP rheology
ln_rhg_EAP = .false., ! EAP rheology
ln_aEVP = .true., ! adaptive rheology (Kimmritz et al. 2016 & 2017)
rn_creepl = 2.0e-9 ! creep limit [1/s]
rn_ecc = 2.0 ! eccentricity of the elliptical yield curve
nn_nevp = 100 ! number of EVP subcycles
......@@ -117,7 +117,7 @@
! = 1 check at the main time step (output xml: uice_cvg)
! = 2 check at both main and rheology time steps (additional output: ice_cvg.nc)
! this option 2 asks a lot of communications between cpu
ln_rhg_VP = .false. ! VP rheology
ln_rhg_VP = .false., ! VP rheology
nn_vp_nout = 10 ! number of outer iterations
nn_vp_ninn = 1500 ! number of inner iterations
nn_vp_chkcvg = 5 ! iteration step for convergence check
......@@ -125,8 +125,8 @@
!------------------------------------------------------------------------------
&namdyn_adv ! Ice advection
!------------------------------------------------------------------------------
ln_adv_Pra = .true. ! Advection scheme (Prather)
ln_adv_UMx = .false. ! Advection scheme (Ultimate-Macho)
ln_adv_Pra = .true., ! Advection scheme (Prather)
ln_adv_UMx = .false., ! Advection scheme (Ultimate-Macho)
nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order)
/
!------------------------------------------------------------------------------
......@@ -144,8 +144,8 @@
! = 0 Average N(cat) fluxes then apply the average over the N(cat) ice
! = 1 Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity
! = 2 Redistribute a single flux over categories
ln_cndflx = .false. ! Use conduction flux as surface boundary conditions (i.e. for Jules coupling)
ln_cndemulate = .false. ! emulate conduction flux (if not provided in the inputs)
ln_cndflx = .false., ! Use conduction flux as surface boundary conditions (i.e. for Jules coupling)
ln_cndemulate = .false., ! emulate conduction flux (if not provided in the inputs)
nn_qtrice = 0 ! Solar flux transmitted thru the surface scattering layer:
! = 0 Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow)
! = 1 Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities)
......@@ -153,26 +153,26 @@
!------------------------------------------------------------------------------
&namthd ! Ice thermodynamics
!------------------------------------------------------------------------------
ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F)
ln_icedA = .true. ! activate lateral melting param. (T) or not (F)
ln_icedO = .true. ! activate ice growth in open-water (T) or not (F)
ln_icedS = .true. ! activate brine drainage (T) or not (F)
ln_icedH = .true., ! activate ice thickness change from growing/melting (T) or not (F)
ln_icedA = .true., ! activate lateral melting param. (T) or not (F)
ln_icedO = .true., ! activate ice growth in open-water (T) or not (F)
ln_icedS = .true., ! activate brine drainage (T) or not (F)
!
ln_leadhfx = .true. ! heat in the leads is used to melt sea-ice before warming the ocean
ln_leadhfx = .true., ! heat in the leads is used to melt sea-ice before warming the ocean
/
!------------------------------------------------------------------------------
&namthd_zdf ! Ice heat diffusion
!------------------------------------------------------------------------------
ln_zdf_BL99 = .true. ! Heat diffusion follows Bitz and Lipscomb 1999
ln_cndi_U64 = .false. ! sea ice thermal conductivity: k = k0 + beta.S/T (Untersteiner, 1964)
ln_cndi_P07 = .true. ! sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007)
ln_zdf_BL99 = .true., ! Heat diffusion follows Bitz and Lipscomb 1999
ln_cndi_U64 = .false., ! sea ice thermal conductivity: k = k0 + beta.S/T (Untersteiner, 1964)
ln_cndi_P07 = .true., ! sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007)
rn_cnd_s = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971)
! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)
rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m]
rn_kappa_s = 10.0 ! nn_qtrice = 0: radiation attenuation coefficient in snow [1/m]
rn_kappa_smlt = 7.0 ! nn_qtrice = 1: radiation attenuation coefficient in melting snow [1/m]
rn_kappa_sdry = 10.0 ! radiation attenuation coefficient in dry snow [1/m]
ln_zdf_chkcvg = .false. ! check convergence of heat diffusion scheme (outputs: tice_cvgerr, tice_cvgstp)
ln_zdf_chkcvg = .false., ! check convergence of heat diffusion scheme (outputs: tice_cvgerr, tice_cvgstp)
/
!------------------------------------------------------------------------------
&namthd_da ! Ice lateral melting
......@@ -189,7 +189,7 @@
&namthd_do ! Ice growth in open water
!------------------------------------------------------------------------------
rn_hinew = 0.1 ! thickness for new ice formation in open water (m), must be larger than rn_himin
ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind)
ln_frazil = .false., ! Frazil ice parameterization (ice collection as a function of wind)
rn_maxfraz = 1.0 ! maximum fraction of frazil ice collecting at the ice base
rn_vfraz = 0.417 ! thresold drift speed for frazil ice collecting at the ice bottom (m/s)
rn_Cfraz = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom
......@@ -212,22 +212,25 @@
!------------------------------------------------------------------------------
&namthd_pnd ! Melt ponds
!------------------------------------------------------------------------------
ln_pnd = .true. ! activate melt ponds or not
ln_pnd_TOPO = .false. ! topographic melt ponds
ln_pnd_LEV = .true. ! level ice melt ponds
ln_pnd = .true., ! activate melt ponds or not
ln_pnd_TOPO = .false., ! topographic melt ponds
ln_pnd_LEV = .true., ! level ice melt ponds
rn_apnd_min = 0.15 ! minimum meltwater fraction contributing to pond growth (TOPO and LEV)
rn_apnd_max = 0.85 ! maximum meltwater fraction contributing to pond growth (TOPO and LEV)
rn_pnd_flush= 0.1 ! pond flushing efficiency (tuning parameter) (LEV)
ln_pnd_CST = .false. ! constant melt ponds
ln_pnd_CST = .false., ! constant melt ponds
rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC
rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC
ln_pnd_lids = .true. ! frozen lids on top of the ponds (only for ln_pnd_LEV)
ln_pnd_alb = .true. ! effect of melt ponds on ice albedo
ln_pnd_lids = .true., ! frozen lids on top of the ponds (only for ln_pnd_LEV)
ln_pnd_alb = .true., ! effect of melt ponds on ice albedo
nn_pnd_brsal = 0 ! brine salinity formulation 0 = Consistent expression with SI3
! (linear liquidus)
! 1 = used in GOSI9
/
!------------------------------------------------------------------------------
&namini ! Ice initialization
!------------------------------------------------------------------------------
ln_iceini = .true. ! activate ice initialization (T) or not (F)
ln_iceini = .true., ! activate ice initialization (T) or not (F)
nn_iceini_file = 0 ! 0 = Initialise sea ice based on SSTs
! 1 = Initialise sea ice from single category netcdf file
! 2 = Initialise sea ice from multi category restart file
......@@ -280,12 +283,12 @@
!------------------------------------------------------------------------------
&namdia ! Diagnostics
!------------------------------------------------------------------------------
ln_icediachk = .false. ! check online heat, mass & salt budgets
ln_icediachk = .false., ! check online heat, mass & salt budgets
! ! rate of ice spuriously gained/lost at each time step => rn_icechk=1 <=> 1.e-6 m/hour
rn_icechk_cel = 1. ! check at each gridcell (1.e-06m/h)=> stops the code if violated (and writes a file)
rn_icechk_glo = 1.e-04 ! check over the entire ice cover (1.e-10m/h)=> only prints warnings
ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F)
ln_icectl = .false. ! ice points output for debug (T or F)
ln_icediahsb = .false., ! output the heat, mass & salt budgets (T) or not (F)
ln_icectl = .false., ! ice points output for debug (T or F)
iiceprt = 10 ! i-index for debug
jiceprt = 10 ! j-index for debug
/
cfgs/SHARED/namelist_pisces_ref
View file @ 4a9f4b18
......@@ -152,7 +152,6 @@
!-----------------------------------------------------------------------
pislopen = 2. ! P-I slope
pisloped = 2. ! P-I slope for diatoms
xadap = 0. ! Adaptation factor to low light
excretn = 0.05 ! excretion ratio of phytoplankton
excretd = 0.05 ! excretion ratio of diatoms
bresp = 0.033 ! Basal respiration rate
......@@ -172,7 +171,6 @@
excretn = 0.05 ! excretion ratio of phytoplankton
excretp = 0.05 ! excretion ratio of picophytoplankton
excretd = 0.05 ! excretion ratio of diatoms
xadap = 0. ! Adaptation factor to low light
bresp = 0.02 ! Basal respiration rate
thetannm = 0.3 ! Maximum Chl/N in nanophytoplankton
thetanpm = 0.3 ! Maximum Chl/N in picophytoplankton
......