Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/domain_def.xml
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/domain_def.xml (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/domain_def.xml (revision 5267)
@@ -0,0 +1,176 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/field_def.xml
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/field_def.xml (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/field_def.xml (revision 5267)
@@ -0,0 +1,659 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/iodef.xml
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/iodef.xml (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/iodef.xml (revision 5267)
@@ -0,0 +1,227 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 50000000
+ 2
+ 0
+ false
+ false
+ oceanx
+
+
+
+
+
+
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/kRGB61.txt
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/kRGB61.txt (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/kRGB61.txt (revision 5267)
@@ -0,0 +1,61 @@
+ 0.010 0.01618 0.07464 0.37807
+ 0.011 0.01654 0.07480 0.37823
+ 0.013 0.01693 0.07499 0.37840
+ 0.014 0.01736 0.07518 0.37859
+ 0.016 0.01782 0.07539 0.37879
+ 0.018 0.01831 0.07562 0.37900
+ 0.020 0.01885 0.07586 0.37923
+ 0.022 0.01943 0.07613 0.37948
+ 0.025 0.02005 0.07641 0.37976
+ 0.028 0.02073 0.07672 0.38005
+ 0.032 0.02146 0.07705 0.38036
+ 0.035 0.02224 0.07741 0.38070
+ 0.040 0.02310 0.07780 0.38107
+ 0.045 0.02402 0.07821 0.38146
+ 0.050 0.02501 0.07866 0.38189
+ 0.056 0.02608 0.07914 0.38235
+ 0.063 0.02724 0.07967 0.38285
+ 0.071 0.02849 0.08023 0.38338
+ 0.079 0.02984 0.08083 0.38396
+ 0.089 0.03131 0.08149 0.38458
+ 0.100 0.03288 0.08219 0.38526
+ 0.112 0.03459 0.08295 0.38598
+ 0.126 0.03643 0.08377 0.38676
+ 0.141 0.03842 0.08466 0.38761
+ 0.158 0.04057 0.08561 0.38852
+ 0.178 0.04289 0.08664 0.38950
+ 0.200 0.04540 0.08775 0.39056
+ 0.224 0.04811 0.08894 0.39171
+ 0.251 0.05103 0.09023 0.39294
+ 0.282 0.05420 0.09162 0.39428
+ 0.316 0.05761 0.09312 0.39572
+ 0.355 0.06130 0.09474 0.39727
+ 0.398 0.06529 0.09649 0.39894
+ 0.447 0.06959 0.09837 0.40075
+ 0.501 0.07424 0.10040 0.40270
+ 0.562 0.07927 0.10259 0.40480
+ 0.631 0.08470 0.10495 0.40707
+ 0.708 0.09056 0.10749 0.40952
+ 0.794 0.09690 0.11024 0.41216
+ 0.891 0.10374 0.11320 0.41502
+ 1.000 0.11114 0.11639 0.41809
+ 1.122 0.11912 0.11984 0.42142
+ 1.259 0.12775 0.12356 0.42500
+ 1.413 0.13707 0.12757 0.42887
+ 1.585 0.14715 0.13189 0.43304
+ 1.778 0.15803 0.13655 0.43754
+ 1.995 0.16978 0.14158 0.44240
+ 2.239 0.18248 0.14701 0.44765
+ 2.512 0.19620 0.15286 0.45331
+ 2.818 0.21102 0.15918 0.45942
+ 3.162 0.22703 0.16599 0.46601
+ 3.548 0.24433 0.17334 0.47313
+ 3.981 0.26301 0.18126 0.48080
+ 4.467 0.28320 0.18981 0.48909
+ 5.012 0.30502 0.19903 0.49803
+ 5.623 0.32858 0.20898 0.50768
+ 6.310 0.35404 0.21971 0.51810
+ 7.079 0.38154 0.23129 0.52934
+ 7.943 0.41125 0.24378 0.54147
+ 8.912 0.44336 0.25725 0.55457
+ 10.000 0.47804 0.27178 0.56870
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist (revision 5267)
@@ -0,0 +1,1210 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OPA : 1 - run manager (namrun, namcfg)
+!! namelists 2 - Domain (namzgr, namzgr_sco, namdom, namtsd)
+!! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas
+!! namsbc_cpl, namtra_qsr, namsbc_rnf,
+!! namsbc_apr, namsbc_ssr, namsbc_alb)
+!! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide)
+!! 5 - bottom boundary (nambfr, nambbc, nambbl)
+!! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_dmp)
+!! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
+!! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx)
+!! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namptr, namhsb)
+!! 10 - miscellaneous (namsol, nammpp, namctl)
+!! 11 - Obs & Assim (namobs, nam_asminc)
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
+!!======================================================================
+!! *** Run management namelists ***
+!!======================================================================
+!! namrun parameters of the run
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namrun ! parameters of the run
+!-----------------------------------------------------------------------
+ nn_no = 0 ! job number (no more used...)
+ cn_exp = "PISCES" ! experience name
+ nn_it000 = 1 ! first time step
+ nn_itend = 3650 ! last time step (std 5475)
+ nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+ nn_leapy = 0 ! Leap year calendar (1) or not (0)
+ ln_rstart = .false. ! start from rest (F) or from a restart file (T)
+ nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=.true.
+ nn_rstctl = 0 ! restart control => activated only if ln_rstart = T
+ ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
+ ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
+ ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart
+ cn_ocerst_in = "restart" ! suffix of ocean restart name (input)
+ cn_ocerst_out = "restart" ! suffix of ocean restart name (output)
+ nn_istate = 0 ! output the initial state (1) or not (0)
+ nn_stock = 3650 ! frequency of creation of a restart file (modulo referenced to 1)
+ nn_write = 3650 ! frequency of write in the output file (modulo referenced to nn_it000)
+ ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)
+ ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%)
+ ln_clobber = .false. ! clobber (overwrite) an existing file
+ nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
+/
+!
+!-----------------------------------------------------------------------
+&namcfg ! default parameters of the configuration
+!-----------------------------------------------------------------------
+ cp_cfg = "orca" ! name of the configuration
+ cp_cfz = '' ! name of the zoom of configuration
+ jp_cfg = 2 ! resolution of the configuration
+ jpidta = 182 ! 1st lateral dimension ( >= jpi )
+ jpjdta = 149 ! 2nd " " ( >= jpj )
+ jpkdta = 31 ! number of levels ( >= jpk )
+ jpiglo = 182 ! 1st dimension of global domain --> i =jpidta
+ jpjglo = 149 ! 2nd - - --> j =jpjdta
+ jpizoom = 1 ! left bottom (i,j) indices of the zoom
+ jpjzoom = 1 ! in data domain indices
+ jperio = 4 ! lateral cond. type (between 0 and 6)
+ ! = 0 closed ; = 1 cyclic East-West
+ ! = 2 equatorial symmetric ; = 3 North fold T-point pivot
+ ! = 4 cyclic East-West AND North fold T-point pivot
+ ! = 5 North fold F-point pivot
+ ! = 6 cyclic East-West AND North fold F-point pivot
+/
+!!======================================================================
+!! *** Domain namelists ***
+!!======================================================================
+!! namzgr vertical coordinate
+!! namzgr_sco s-coordinate or hybrid z-s-coordinate
+!! namdom space and time domain (bathymetry, mesh, timestep)
+!! namtsd data: temperature & salinity
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namzgr ! vertical coordinate
+!-----------------------------------------------------------------------
+ ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined)
+ ln_zps = .true. ! z-coordinate - partial steps (T/F)
+ ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F)
+/
+!-----------------------------------------------------------------------
+&namzgr_sco ! s-coordinate or hybrid z-s-coordinate
+!-----------------------------------------------------------------------
+ ln_s_sh94 = .true. ! Song & Haidvogel 1994 hybrid S-sigma (T)|
+ ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied
+ ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch
+ ! stretching coefficients for all functions
+ rn_sbot_min = 10.0 ! minimum depth of s-bottom surface (>0) (m)
+ rn_sbot_max = 7000.0 ! maximum depth of s-bottom surface (= ocean depth) (>0) (m)
+ rn_hc = 150.0 ! critical depth for transition to stretched coordinates
+ !!!!!!! Envelop bathymetry
+ rn_rmax = 0.3 ! maximum cut-off r-value allowed (00) or min number of ocean level (<0)
+ rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of
+ rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0 fill namsbc_ana )
+ ln_flx = .false. ! flux formulation (T => fill namsbc_flx )
+ ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio)
+ ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core)
+ ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs )
+ ln_cpl = .false. ! Coupled formulation (T => fill namsbc_cpl )
+ ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )
+ nn_ice = 2 ! =0 no ice boundary condition ,
+ ! =1 use observed ice-cover ,
+ ! =2 ice-model used ("key_lim3" or "key_lim2)
+ nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect)
+ ! =1 levitating ice with mass and salt exchange but no presure effect
+ ! =2 embedded sea-ice (full salt and mass exchanges and pressure)
+ ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave
+ ln_rnf = .true. ! runoffs (T => fill namsbc_rnf)
+ ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
+ nn_fwb = 2 ! FreshWater Budget: =0 unchecked
+ ! =1 global mean of e-p-r set to zero at each time step
+ ! =2 annual global mean of e-p-r set to zero
+ ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave)
+ ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave)
+ ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave)
+ cn_iceflx = 'linear' ! redistribution of solar input into ice categories during coupling ice/atm.
+/
+!-----------------------------------------------------------------------
+&namsbc_ana ! analytical surface boundary condition
+!-----------------------------------------------------------------------
+ nn_tau000 = 0 ! gently increase the stress over the first ntau_rst time-steps
+ rn_utau0 = 0.5 ! uniform value for the i-stress
+ rn_vtau0 = 0.e0 ! uniform value for the j-stress
+ rn_qns0 = 0.e0 ! uniform value for the total heat flux
+ rn_qsr0 = 0.e0 ! uniform value for the solar radiation
+ rn_emp0 = 0.e0 ! uniform value for the freswater budget (E-P)
+/
+!-----------------------------------------------------------------------
+&namsbc_flx ! surface boundary condition : flux formulation
+!-----------------------------------------------------------------------
+! ! 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_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the flux files
+/
+!-----------------------------------------------------------------------
+&namsbc_clio ! namsbc_clio CLIO bulk formulae
+!-----------------------------------------------------------------------
+! ! 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_utau = 'taux_1m' , -1 , 'sozotaux', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_vtau = 'tauy_1m' , -1 , 'sometauy', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_wndm = 'flx' , -1 , 'socliowi', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_tair = 'flx' , -1 , 'socliot2', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_humi = 'flx' , -1 , 'socliohu', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ccov = 'flx' , -1 , 'socliocl', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_prec = 'flx' , -1 , 'socliopl', .false. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the bulk files are
+/
+!-----------------------------------------------------------------------
+&namsbc_core ! namsbc_core CORE bulk formulae
+!-----------------------------------------------------------------------
+! ! 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_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , ''
+ sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , ''
+ 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_tair = 't_10.15JUNE2009_fill' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
+ sn_humi = 'q_10.15JUNE2009_fill' , 6 , 'Q_10_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_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the bulk files
+ ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F)
+ ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data
+ ln_bulk2z = .false. ! Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu
+ rn_zqt = 3. ! Air temperature and humidity reference height (m) (ln_bulk2z)
+ rn_zu = 4. ! Wind vector reference height (m) (ln_bulk2z)
+ rn_pfac = 1. ! multiplicative factor for precipitation (total & snow)
+ rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.)
+ rn_vfac = 0. ! multiplicative factor for ocean/ice velocity
+ ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
+/
+!-----------------------------------------------------------------------
+&namsbc_mfs ! namsbc_mfs MFS bulk formulae
+!-----------------------------------------------------------------------
+! ! 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_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' , ''
+ sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' , ''
+ sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false. , 'daily' ,'bilinear.nc', '' , ''
+ sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false. , 'daily' ,'bicubic.nc' , '' , ''
+ sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' , ''
+ sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' , ''
+ sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' , ''
+
+ cn_dir = './ECMWF/' ! root directory for the location of the bulk files
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled")
+!-----------------------------------------------------------------------
+! ! description ! multiple ! vector ! vector ! vector !
+! ! ! categories ! reference ! orientation ! grids !
+! send
+sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , ''
+sn_snd_alb = 'weighted ice' , 'no' , '' , '' , ''
+sn_snd_thick = 'none' , 'no' , '' , '' , ''
+sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T'
+sn_snd_co2 = 'coupled' , 'no' , '' , '' , ''
+! receive
+sn_rcv_w10m = 'none' , 'no' , '' , '' , ''
+sn_rcv_taumod = 'coupled' , 'no' , '' , '' , ''
+sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V'
+sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , ''
+sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , ''
+sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , ''
+sn_rcv_emp = 'conservative' , 'no' , '' , '' , ''
+sn_rcv_rnf = 'coupled' , 'no' , '' , '' , ''
+sn_rcv_cal = 'coupled' , 'no' , '' , '' , ''
+sn_rcv_co2 = 'coupled' , 'no' , '' , '' , ''
+/
+!-----------------------------------------------------------------------
+&namsbc_sas ! analytical surface boundary condition
+!-----------------------------------------------------------------------
+! ! 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_usp = 'sas_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_vsp = 'sas_grid_V' , 120 , 'vomecrty' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_tem = 'sas_grid_T' , 120 , 'sosstsst' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_sal = 'sas_grid_T' , 120 , 'sosaline' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ssh = 'sas_grid_T' , 120 , 'sossheig' , .true. , .true. , 'yearly' , '' , '' , ''
+
+ ln_3d_uv = .true. ! specify whether we are supplying a 3D u,v field
+ cn_dir = './' ! root directory for the location of the bulk files are
+/
+!-----------------------------------------------------------------------
+&namtra_qsr ! penetrative solar radiation
+!-----------------------------------------------------------------------
+! ! 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_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the runoff files
+ ln_traqsr = .false. ! Light penetration (T) or not (F)
+ ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration
+ ln_qsr_2bd = .false. ! 2 bands light penetration
+ ln_qsr_bio = .false. ! bio-model light penetration
+ nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0)
+ rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1)
+ rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction
+ rn_si1 = 23.0 ! 2 bands: longest depth of extinction
+ ln_qsr_ice = .true. ! light penetration for ice-model LIM3
+/
+!-----------------------------------------------------------------------
+&namsbc_rnf ! runoffs namelist surface boundary condition
+!-----------------------------------------------------------------------
+! ! 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' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the runoff files
+ ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F)
+ ln_rnf_mouth = .true. ! specific treatment at rivers mouths
+ rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used
+ rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s]
+ rn_rfact = 1.e0 ! multiplicative factor for runoff
+ ln_rnf_depth = .false. ! read in depth information for runoff
+ ln_rnf_tem = .false. ! read in temperature information for runoff
+ ln_rnf_sal = .false. ! read in salinity information for runoff
+/
+!-----------------------------------------------------------------------
+&namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk
+!-----------------------------------------------------------------------
+! ! 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_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the bulk files
+ rn_pref = 101000._wp ! reference atmospheric pressure [N/m2]/
+ ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F)
+ ln_apr_obc = .false. ! inverse barometer added to OBC ssh data
+/
+!-----------------------------------------------------------------------
+&namsbc_ssr ! surface boundary condition : sea surface restoring
+!-----------------------------------------------------------------------
+! ! 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_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the runoff files
+ nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0)
+ nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2)
+ ! or to SSS only (=1) or no damping term (=0)
+ rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K]
+ 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_alb ! albedo parameters
+!-----------------------------------------------------------------------
+ rn_cloud = 0.06 ! cloud correction to snow and ice albedo
+ rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic
+ rn_alphd = 0.80 ! coefficients for linear interpolation used to
+ rn_alphc = 0.65 ! compute albedo between two extremes values
+ rn_alphdi = 0.72 ! (Pyane, 1972)
+/
+!-----------------------------------------------------------------------
+&namberg ! iceberg parameters
+!-----------------------------------------------------------------------
+ ln_icebergs = .false.
+ ln_bergdia = .true. ! Calculate budgets
+ nn_verbose_level = 1 ! Turn on more verbose output if level > 0
+ nn_verbose_write = 15 ! Timesteps between verbose messages
+ nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage
+ ! Initial mass required for an iceberg of each class
+ rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
+ ! Proportion of calving mass to apportion to each class
+ rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
+ ! Ratio between effective and real iceberg mass (non-dim)
+ ! i.e. number of icebergs represented at a point
+ rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1
+ ! thickness of newly calved bergs (m)
+ rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
+ rn_rho_bergs = 850. ! Density of icebergs
+ rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs
+ ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics
+ rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits
+ rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (00)
+ ! + F(distance-to-coast) + Red and Med Seas
+ nn_zdmp = 0 ! vertical shape =0 damping throughout the water column
+ ! =1 no damping in the mixing layer (kz criteria)
+ ! =2 no damping in the mixed layer (rho crieria)
+ rn_surf = 50. ! surface time scale of damping [days]
+ rn_bot = 360. ! bottom time scale of damping [days]
+ rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters]
+ nn_file = 0 ! create a damping.coeff NetCDF file (=1) or not (=0)
+/
+
+!!======================================================================
+!! *** Dynamics namelists ***
+!!======================================================================
+!! namdyn_adv formulation of the momentum advection
+!! namdyn_vor advection scheme
+!! namdyn_hpg hydrostatic pressure gradient
+!! namdyn_spg surface pressure gradient (CPP key only)
+!! namdyn_ldf lateral diffusion scheme
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namdyn_adv ! formulation of the momentum advection
+!-----------------------------------------------------------------------
+ ln_dynadv_vec = .true. ! vector form (T) or flux form (F)
+ ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme
+ ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme
+/
+!-----------------------------------------------------------------------
+&nam_vvl ! vertical coordinate options
+!-----------------------------------------------------------------------
+ ln_vvl_zstar = .true. ! zstar vertical coordinate
+ ln_vvl_ztilde = .false. ! ztilde vertical coordinate: only high frequency variations
+ ln_vvl_layer = .false. ! full layer vertical coordinate
+ ln_vvl_ztilde_as_zstar = .false. ! ztilde vertical coordinate emulating zstar
+ ln_vvl_zstar_at_eqtor = .false. ! ztilde near the equator
+ rn_ahe3 = 0.0e0 ! thickness diffusion coefficient
+ rn_rst_e3t = 30.e0 ! ztilde to zstar restoration timescale [days]
+ rn_lf_cutoff = 5.0e0 ! cutoff frequency for low-pass filter [days]
+ rn_zdef_max = 0.9e0 ! maximum fractional e3t deformation
+ ln_vvl_dbg = .true. ! debug prints (T/F)
+/
+!-----------------------------------------------------------------------
+&namdyn_vor ! option of physics/algorithm (not control by CPP keys)
+!-----------------------------------------------------------------------
+ ln_dynvor_ene = .false. ! enstrophy conserving scheme
+ ln_dynvor_ens = .false. ! energy conserving scheme
+ ln_dynvor_mix = .false. ! mixed scheme
+ ln_dynvor_een = .true. ! energy & enstrophy scheme
+/
+!-----------------------------------------------------------------------
+&namdyn_hpg ! Hydrostatic pressure gradient option
+!-----------------------------------------------------------------------
+ ln_hpg_zco = .false. ! z-coordinate - full steps
+ ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation)
+ ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation)
+ ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial)
+ ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme)
+ ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T)
+ ! centered time scheme (F)
+/
+!-----------------------------------------------------------------------
+!namdyn_spg ! surface pressure gradient (CPP key only)
+!-----------------------------------------------------------------------
+! ! explicit free surface ("key_dynspg_exp")
+! ! filtered free surface ("key_dynspg_flt")
+! ! split-explicit free surface ("key_dynspg_ts")
+
+!-----------------------------------------------------------------------
+&namdyn_ldf ! lateral diffusion on momentum
+!-----------------------------------------------------------------------
+ ! ! Type of the operator :
+ ln_dynldf_lap = .true. ! laplacian operator
+ ln_dynldf_bilap = .false. ! bilaplacian operator
+ ! ! Direction of action :
+ ln_dynldf_level = .false. ! iso-level
+ ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)
+ ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp")
+ ! ! Coefficient
+ rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s]
+ rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s]
+ rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s]
+ rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity
+ rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity
+ rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity
+ rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp)
+ rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap)
+/
+!-----------------------------------------------------------------------
+&namdta_dyn ! offline dynamics read in files ("key_offline")
+!-----------------------------------------------------------------------
+! ! 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_tem = 'dyna_grid_T' , 120 , 'votemper' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_sal = 'dyna_grid_T' , 120 , 'vosaline' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_mld = 'dyna_grid_T' , 120 , 'somixhgt' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_emp = 'dyna_grid_T' , 120 , 'sowaflup' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_fmf = 'dyna_grid_T' , 120 , 'iowaflup' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ice = 'dyna_grid_T' , 120 , 'soicecov' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_qsr = 'dyna_grid_T' , 120 , 'soshfldo' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_wnd = 'dyna_grid_T' , 120 , 'sowindsp' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_uwd = 'dyna_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_vwd = 'dyna_grid_V' , 120 , 'vomecrty' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_wwd = 'dyna_grid_W' , 120 , 'vovecrtz' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_avt = 'dyna_grid_W' , 120 , 'voddmavs' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ubl = 'dyna_grid_U' , 120 , 'sobblcox' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_vbl = 'dyna_grid_V' , 120 , 'sobblcoy' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ahu = 'dyna_grid_U' , 120 , 'vozoahtu' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ahv = 'dyna_grid_V' , 120 , 'vomeahtv' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ahw = 'dyna_grid_W' , 120 , 'voveahtz' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_eiu = 'dyna_grid_U' , 120 , 'vozoaeiu' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_eiv = 'dyna_grid_V' , 120 , 'vomeaeiv' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_eiw = 'dyna_grid_W' , 120 , 'soleaeiw' , .true. , .true. , 'yearly' , '' , '' , ''
+!
+ cn_dir = './' ! root directory for the location of the dynamical files
+ ln_degrad = .false. ! flag for degradation - requires ("key_degrad")
+ ln_dynwzv = .true. ! computation of vertical velocity instead of using the one read in file
+ ln_dynbbl = .true. ! bbl coef are in files, so read them - requires ("key_trabbl")
+/
+!!======================================================================
+!! Tracers & Dynamics vertical physics namelists
+!!======================================================================
+!! namzdf vertical physics
+!! namzdf_ric richardson number dependent vertical mixing ("key_zdfric")
+!! namzdf_tke TKE dependent vertical mixing ("key_zdftke")
+!! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")
+!! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm")
+!! namzdf_tmx tidal mixing parameterization ("key_zdftmx")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namzdf ! vertical physics
+!-----------------------------------------------------------------------
+ rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst")
+ rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst")
+ nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)
+ nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0)
+ ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F)
+ nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1)
+ rn_avevd = 100. ! evd mixing coefficient [m2/s]
+ ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F)
+ nn_npc = 1 ! frequency of application of npc
+ nn_npcp = 365 ! npc control print frequency
+ ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping
+ nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T
+/
+!-----------------------------------------------------------------------
+&namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" )
+!-----------------------------------------------------------------------
+ rn_avmri = 100.e-4 ! maximum value of the vertical viscosity
+ rn_alp = 5. ! coefficient of the parameterization
+ nn_ric = 2 ! coefficient of the parameterization
+ rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation
+ rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m)
+ rn_mldmax =1000.0 ! maximum allowable mixed-layer depth estimate (m)
+ rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer
+ rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer
+ ln_mldw = .true. ! Flag to use or not the mized layer depth param.
+/
+!-----------------------------------------------------------------------
+&namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke")
+!-----------------------------------------------------------------------
+ rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
+ rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation
+ rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
+ rn_emin = 1.e-6 ! minimum value of tke [m2/s2]
+ rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2]
+ rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it)
+ nn_mxl = 2 ! 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 disspipative an mixing length scale
+ nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
+ ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F)
+ rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value
+ ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002)
+ rn_lc = 0.15 ! coef. associated to Langmuir cells
+ nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves
+ ! = 0 no penetration
+ ! = 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 ("key_coupled")
+ rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
+ nn_htau = 1 ! type of exponential decrease of tke penetration below the ML
+ ! = 0 constant 10 m length scale
+ ! = 1 0.5m at the equator to 30m poleward of 40 degrees
+/
+!------------------------------------------------------------------------
+&namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:
+!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
+ ln_kpprimix = .true. ! shear instability mixing
+ rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]
+ rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]
+ rn_riinfty = 0.8 ! local Richardson Number limit for shear instability
+ rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]
+ rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]
+ rn_difcon = 1. ! maximum mixing in interior convection [m2/s]
+ nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv
+ nn_ave = 1 ! constant (=0) or profile (=1) background on avt
+/
+!-----------------------------------------------------------------------
+&namzdf_gls ! GLS vertical diffusion ("key_zdfgls")
+!-----------------------------------------------------------------------
+ rn_emin = 1.e-6 ! minimum value of e [m2/s2]
+ rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3]
+ ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988)
+ rn_clim_galp = 0.53 ! galperin limit
+ ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation
+ ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case
+ rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux
+ rn_charn = 70000. ! Charnock constant for wb induced roughness length
+ nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg)
+ nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum)
+ nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg)
+ nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum)
+ nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
+ nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
+/
+!-----------------------------------------------------------------------
+&namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm")
+!-----------------------------------------------------------------------
+ rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity)
+ rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio
+/
+!-----------------------------------------------------------------------
+&namzdf_tmx ! tidal mixing parameterization ("key_zdftmx")
+!-----------------------------------------------------------------------
+ rn_htmx = 500. ! vertical decay scale for turbulence (meters)
+ rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1)
+ rn_tfe = 0.333 ! tidal dissipation efficiency
+ rn_me = 0.2 ! mixing efficiency
+ ln_tmx_itf = .true. ! ITF specific parameterisation
+ rn_tfe_itf = 1. ! ITF tidal dissipation efficiency
+/
+
+!!======================================================================
+!! *** Miscellaneous namelists ***
+!!======================================================================
+!! nammpp Massively Parallel Processing ("key_mpp_mpi)
+!! namctl Control prints & Benchmark
+!! namsol elliptic solver / island / free surface
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsol ! elliptic solver / island / free surface
+!-----------------------------------------------------------------------
+ nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)
+ ! =2 successive-over-relaxation (sor)
+ nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test
+ rn_eps = 1.e-6 ! absolute precision of the solver
+ nn_nmin = 300 ! minimum of iterations for the SOR solver
+ nn_nmax = 800 ! maximum of iterations for the SOR solver
+ nn_nmod = 10 ! frequency of test for the SOR solver
+ rn_resmax = 1.e-10 ! absolute precision for the SOR solver
+ rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain)
+/
+!-----------------------------------------------------------------------
+&nammpp ! Massively Parallel Processing ("key_mpp_mpi)
+!-----------------------------------------------------------------------
+ cn_mpi_send = 'I' ! mpi send/recieve type ='S', 'B', or 'I' for standard send,
+ ! buffer blocking send or immediate non-blocking sends, resp.
+ nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation
+ ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold
+ jpni = 0 ! jpni number of processors following i (set automatically if < 1)
+ jpnj = 0 ! jpnj number of processors following j (set automatically if < 1)
+ jpnij = 0 ! jpnij number of local domains (set automatically if < 1)
+/
+!-----------------------------------------------------------------------
+&namctl ! Control prints & Benchmark
+!-----------------------------------------------------------------------
+ ln_ctl = .false. ! trends control print (expensive!)
+ nn_print = 0 ! level of print (0 no extra print)
+ nn_ictls = 0 ! start i indice of control sum (use to compare mono versus
+ nn_ictle = 0 ! end i indice of control sum multi processor runs
+ nn_jctls = 0 ! start j indice of control over a subdomain)
+ nn_jctle = 0 ! end j indice of control
+ nn_isplt = 1 ! number of processors in i-direction
+ nn_jsplt = 1 ! number of processors in j-direction
+ nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench
+ ! (no physical validity of the results)
+ nn_timing = 1 ! timing by routine activated (=1) creates timing.output file, or not (=0)
+/
+!-----------------------------------------------------------------------
+&namc1d ! 1D configuration options ("key_c1d")
+!-----------------------------------------------------------------------
+ rn_lat1d = 30 ! Column latitude
+ rn_lon1d = -64 ! Column longitude
+ ln_c1d_locpt = .true.
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd ! data: U & V currents ("key_c1d")
+!-----------------------------------------------------------------------
+! ! 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_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , ''
+ sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , ''
+!
+ cn_dir = './' ! root directory for the location of the files
+ ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F)
+ ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp ! U & V newtonian damping ("key_c1d")
+!-----------------------------------------------------------------------
+ ln_dyndmp = .false. ! add a damping term (T) or not (F)
+/
+!!======================================================================
+!! *** Diagnostics namelists ***
+!!======================================================================
+!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
+!! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld")
+!! namflo float parameters ("key_float")
+!! namptr Poleward Transport Diagnostics
+!! namhsb Heat and salt budgets
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4")
+!-----------------------------------------------------------------------
+ nn_nchunks_i= 4 ! number of chunks in i-dimension
+ nn_nchunks_j= 4 ! number of chunks in j-dimension
+ nn_nchunks_k= 31 ! number of chunks in k-dimension
+ ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
+ ! is optimal for postprocessing which works exclusively with horizontal slabs
+ ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression
+ ! (F) ignore chunking information and produce netcdf3-compatible files
+/
+!-----------------------------------------------------------------------
+&namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra")
+! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor")
+!-----------------------------------------------------------------------
+ nn_trd = 365 ! time step frequency dynamics and tracers trends
+ nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n /seconds ; =86400. -> /day)
+ cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input)
+ cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output)
+ ln_trdmld_restart = .false. ! restart for ML diagnostics
+ ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S
+/
+!-----------------------------------------------------------------------
+&namflo ! float parameters ("key_float")
+!-----------------------------------------------------------------------
+ jpnfl = 1 ! total number of floats during the run
+ jpnnewflo = 0 ! number of floats for the restart
+ ln_rstflo = .false. ! float restart (T) or not (F)
+ nn_writefl = 75 ! frequency of writing in float output file
+ nn_stockfl = 5475 ! frequency of creation of the float restart file
+ ln_argo = .false. ! Argo type floats (stay at the surface each 10 days)
+ ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T)
+ ! or computed with Blanke' scheme (F)
+ ln_ariane = .true. ! Input with Ariane tool convention(T)
+ ln_flo_ascii = .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T)
+/
+!-----------------------------------------------------------------------
+&namptr ! Poleward Transport Diagnostic
+!-----------------------------------------------------------------------
+ ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F)
+ ln_diaznl = .true. ! Add zonal means and meridional stream functions
+ ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not
+ ! (orca configuration only, need input basins mask file named "subbasins.nc"
+ ln_ptrcomp = .true. ! Add decomposition : overturning
+ nn_fptr = 1 ! Frequency of ptr computation [time step]
+ nn_fwri = 15 ! Frequency of ptr outputs [time step]
+/
+!-----------------------------------------------------------------------
+&namhsb ! Heat and salt budgets
+!-----------------------------------------------------------------------
+ ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm')
+!-----------------------------------------------------------------------
+ nit000_han = 1 ! First time step used for harmonic analysis
+ nitend_han = 75 ! Last time step used for harmonic analysis
+ nstep_han = 15 ! Time step frequency for harmonic analysis
+ tname(1) = 'M2' ! Name of tidal constituents
+ tname(2) = 'K1'
+/
+!-----------------------------------------------------------------------
+&namdct ! transports through sections
+!-----------------------------------------------------------------------
+ nn_dct = 15 ! time step frequency for transports computing
+ nn_dctwri = 15 ! time step frequency for transports writing
+ nn_secdebug = 112 ! 0 : no section to debug
+ ! -1 : debug all section
+ ! 0 < n : debug section number n
+/
+
+!!======================================================================
+!! *** Observation & Assimilation namelists ***
+!!======================================================================
+!! namobs observation and model comparison ('key_diaobs')
+!! nam_asminc assimilation increments ('key_asminc')
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namobs ! observation usage switch ('key_diaobs')
+!-----------------------------------------------------------------------
+ ln_t3d = .false. ! Logical switch for T profile observations
+ ln_s3d = .false. ! Logical switch for S profile observations
+ ln_ena = .false. ! Logical switch for ENACT insitu data set
+ ! ! ln_cor Logical switch for Coriolis insitu data set
+ ln_profb = .false. ! Logical switch for feedback insitu data set
+ ln_sla = .false. ! Logical switch for SLA observations
+
+ ln_sladt = .false. ! Logical switch for AVISO SLA data
+
+ ln_slafb = .false. ! Logical switch for feedback SLA data
+ ! ln_ssh Logical switch for SSH observations
+
+ ln_sst = .true. ! Logical switch for SST observations
+ ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations
+ ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations
+
+ ln_sstfb = .false. ! Logical switch for feedback SST data
+ ! ln_sss Logical switch for SSS observations
+ ! ln_seaice Logical switch for Sea Ice observations
+ ! ln_vel3d Logical switch for velocity observations
+ ! ln_velavcur Logical switch for velocity daily av. cur.
+ ! ln_velhrcur Logical switch for velocity high freq. cur.
+ ! ln_velavadcp Logical switch for velocity daily av. ADCP
+ ! ln_velhradcp Logical switch for velocity high freq. ADCP
+ ! ln_velfb Logical switch for feedback velocity data
+ ! ln_grid_global Global distribtion of observations
+ ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table
+ ! grid_search_file Grid search lookup file header
+ ! enactfiles ENACT input observation file names
+ ! coriofiles Coriolis input observation file name
+ ! ! profbfiles: Profile feedback input observation file name
+ profbfiles = 'profiles_01.nc'
+ ! ln_profb_enatim Enact feedback input time setting switch
+ ! slafilesact Active SLA input observation file name
+ ! slafilespas Passive SLA input observation file name
+ ! ! slafbfiles: Feedback SLA input observation file name
+ slafbfiles = 'sla_01.nc'
+ ! sstfiles GHRSST input observation file name
+ ! ! sstfbfiles: Feedback SST input observation file name
+ sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc'
+ ! seaicefiles Sea Ice input observation file name
+ ! velavcurfiles Vel. cur. daily av. input file name
+ ! velhvcurfiles Vel. cur. high freq. input file name
+ ! velavadcpfiles Vel. ADCP daily av. input file name
+ ! velhvadcpfiles Vel. ADCP high freq. input file name
+ ! velfbfiles Vel. feedback input observation file name
+ ! dobsini Initial date in window YYYYMMDD.HHMMSS
+ ! dobsend Final date in window YYYYMMDD.HHMMSS
+ ! n1dint Type of vertical interpolation method
+ ! n2dint Type of horizontal interpolation method
+ ! ln_nea Rejection of observations near land switch
+ nmsshc = 0 ! MSSH correction scheme
+ ! mdtcorr MDT correction
+ ! mdtcutoff MDT cutoff for computed correction
+ ln_altbias = .false. ! Logical switch for alt bias
+ ln_ignmis = .true. ! Logical switch for ignoring missing files
+ ! endailyavtypes ENACT daily average types
+ ln_grid_global = .true.
+ ln_grid_search_lookup = .false.
+/
+!-----------------------------------------------------------------------
+&nam_asminc ! assimilation increments ('key_asminc')
+!-----------------------------------------------------------------------
+ ln_bkgwri = .false. ! Logical switch for writing out background state
+ ln_trainc = .false. ! Logical switch for applying tracer increments
+ ln_dyninc = .false. ! Logical switch for applying velocity increments
+ ln_sshinc = .false. ! Logical switch for applying SSH increments
+ ln_asmdin = .false. ! Logical switch for Direct Initialization (DI)
+ ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU)
+ nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1]
+ nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1]
+ nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1]
+ nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1]
+ niaufn = 0 ! Type of IAU weighting function
+ ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin
+ salfixmin = -9999 ! Minimum salinity after applying the increments
+ nn_divdmp = 0 ! Number of iterations of divergence damping operator
+/
+!-----------------------------------------------------------------------
+&namsbc_wave ! External fields from wave model
+!-----------------------------------------------------------------------
+! ! 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_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , ''
+!
+ cn_dir_cdg = './' ! root directory for the location of drag coefficient files
+/
+!-----------------------------------------------------------------------
+&namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed)
+!-----------------------------------------------------------------------
+ ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model
+ ln_neptsimp = .false. ! yes/no use simplified neptune
+
+ ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep
+ rn_tslse = 1.2e4 ! value of lengthscale L at the equator
+ rn_tslsp = 3.0e3 ! value of lengthscale L at the pole
+ ! Specify whether to ramp down the Neptune velocity in shallow
+ ! water, and if so the depth range controlling such ramping down
+ ln_neptramp = .true. ! ramp down Neptune velocity in shallow water
+ rn_htrmin = 100.0 ! min. depth of transition range
+ rn_htrmax = 200.0 ! max. depth of transition range
+/
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist_pisces
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist_pisces (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist_pisces (revision 5267)
@@ -0,0 +1,413 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! PISCES (key_pisces) reference namelist (see below for key_pisces_reduced)
+!! 1 - air-sea exchange (nampisext)
+!! 2 - biological parameters (nampisbio)
+!! 3 - parameters for nutrient limitations (nampislim)
+!! 4 - parameters for phytoplankton (nampisprod,nampismort)
+!! 5 - parameters for zooplankton (nampismes,nampiszoo)
+!! 6 - parameters for remineralization (nampisrem)
+!! 7 - parameters for calcite chemistry (nampiscal)
+!! 8 - parameters for inputs deposition (nampissed)
+!! 9 - parameters for Kriest parameterization (nampiskrp, nampiskrs)
+!! 10 - additional 2D/3D diagnostics (nampisdia)
+!! 11 - Damping (nampisdmp)
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisext ! air-sea exchange
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ ln_co2int = .false. ! read atm pco2 from a file (T) or constant (F)
+ atcco2 = 280. ! Constant value atmospheric pCO2 - ln_co2int = F
+ clname = 'atcco2.txt' ! Name of atm pCO2 file - ln_co2int = T
+ nn_offset = 0 ! Offset model-data start year - ln_co2int = T
+! ! If your model year is iyy, nn_offset=(years(1)-iyy)
+! ! then the first atmospheric CO2 record read is at years(1)
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisatm ! Atmospheric prrssure
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+! ! 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_patm = 'presatm' , -1 , 'patm' , .true. , .true. , 'yearly' , '' , '' , ''
+ cn_dir = './' ! root directory for the location of the dynamical files
+!
+ ln_presatm = .false. ! constant atmopsheric pressure (F) or from a file (T)
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisbio ! biological parameters
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ nrdttrc = 2 ! time step frequency for biology
+ wsbio = 2. ! POC sinking speed
+ xkmort = 1.E-8 ! half saturation constant for mortality
+ ferat3 = 10.E-6 ! Fe/C in zooplankton
+ no3rat3 = 0.167 ! N/C ratio in zooplankton
+ po4rat3 = 0.0094 ! P/C ratio in zooplankton
+ wsbio2 = 50. ! Big particles sinking speed
+ niter1max = 1 ! Maximum number of iterations for POC
+ niter2max = 1 ! Maximum number of iterations for GOC
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampislim ! parameters for nutrient limitations
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ concnno3 = 3e-6 ! Nitrate half saturation of nanophytoplankton
+ concpno3 = 1e-6
+ concdno3 = 4E-6 ! Phosphate half saturation for diatoms
+ concnnh4 = 1.5E-6 ! NH4 half saturation for phyto
+ concpnh4 = 3E-7
+ concdnh4 = 2E-6 ! NH4 half saturation for diatoms
+ concnpo4 = 3E-6 ! PO4 half saturation for phyto
+ concppo4 = 1E-6
+ concdpo4 = 4E-6 ! PO4 half saturation for diatoms
+ concnfer = 3E-9 ! Iron half saturation for phyto
+ concpfer = 1E-9
+ concdfer = 4E-9 ! Iron half saturation for diatoms
+ concbfe = 1.E-11 ! Half-saturation for Fe limitation of Bacteria
+ concbnh4 = 1.E-7 ! NH4 half saturation for phyto
+ concbno3 = 5.E-7 ! Phosphate half saturation for diatoms
+ concbpo4 = 1E-7 ! Phosphate half saturation for bacteria
+ xsizedia = 1.E-6 ! Minimum size criteria for diatoms
+ xsizephy = 1.E-6 ! Minimum size criteria for phyto
+ xsizepic = 1.E-6
+ xsizern = 1.0 ! Size ratio for nanophytoplankton
+ xsizerp = 1.0
+ xsizerd = 4.0 ! Size ratio for diatoms
+ xksi1 = 2.E-6 ! half saturation constant for Si uptake
+ xksi2 = 20E-6 ! half saturation constant for Si/C
+ xkdoc = 417.E-6 ! half-saturation constant of DOC remineralization
+ qfnopt = 7.E-6 ! Optimal quota of phyto
+ qfpopt = 7.E-6
+ qfdopt = 7.E-6 ! Optimal quota of diatoms
+ caco3r = 0.4 ! mean rain ratio
+ qnnmin = 0.29
+ qnnmax = 1.525
+ qpnmin = 0.28
+ qpnmax = 1.06
+ qnpmin = 0.42
+ qnpmax = 1.525
+ qppmin = 0.25
+ qppmax = 0.7
+ qndmin = 0.25
+ qndmax = 1.525
+ qpdmin = 0.29
+ qpdmax = 1.32
+ qfnmax = 40E-6
+ qfpmax = 40E-6
+ qfdmax = 40E-6
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisopt ! parameters for optics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+! ! 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_par = 'par.orca' , 24 , 'fr_par' , .true. , .true. , 'yearly' , '' , '' , ''
+ cn_dir = './' ! root directory for the location of the dynamical files
+ ln_varpar = .true. ! boolean for PAR variable
+ parlux = 0.43 ! Fraction of shortwave as PAR
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisprod ! parameters for phytoplankton growth
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ pislope = 3. ! P-I slope
+ pislopep = 4.
+ pislope2 = 3. ! P-I slope for diatoms
+ xadap = 0. ! Adaptation factor to low light
+ excret = 0.05 ! excretion ratio of phytoplankton
+ excretp = 0.05
+ excret2 = 0.05 ! excretion ratio of diatoms
+ bresp = 0.02 ! Basal respiration rate
+ chlcnm = 0.033 ! Minimum Chl/C in nanophytoplankton
+ chlcpm = 0.033
+ chlcdm = 0.05 ! Minimum Chl/C in diatoms
+ chlcmin = 0.004 ! Maximum Chl/c in phytoplankton
+ grosip = 0.131 ! mean Si/C ratio
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampismort ! parameters for phytoplankton sinks
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ wchl = 0.01 ! quadratic mortality of phytoplankton
+ wchlp = 0.001
+ wchld = 0.01 ! maximum quadratic mortality of diatoms
+ wchldm = 0.02 ! maximum quadratic mortality of diatoms
+ mprat = 0.01 ! phytoplankton mortality rate
+ mpratp = 0.01
+ mprat2 = 0.01 ! Diatoms mortality rate
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampismes ! parameters for mesozooplankton
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ part2 = 0.75 ! part of calcite not dissolved in mesozoo guts
+ grazrat2 = 0.85 ! maximal mesozoo grazing rate
+ resrat2 = 0.005 ! exsudation rate of mesozooplankton
+ mzrat2 = 0.02 ! mesozooplankton mortality rate
+ xpref2d = 1. ! zoo preference for phyto
+ xpref2p = 1. ! zoo preference for POC
+ xpref2z = 1. ! zoo preference for zoo
+ xpref2m = 0.2 ! meso preference for zoo
+ xpref2c = 0.3 ! zoo preference for poc
+ xthresh2zoo = 1E-8 ! zoo feeding threshold for mesozooplankton
+ xthresh2dia = 1E-8 ! diatoms feeding threshold for mesozooplankton
+ xthresh2phy = 1E-8 ! nanophyto feeding threshold for mesozooplankton
+ xthresh2mes = 1E-8 ! meso feeding threshold for mesozooplankton
+ xthresh2poc = 1E-8 ! poc feeding threshold for mesozooplankton
+ xthresh2 = 3E-7 ! Food threshold for grazing
+ xkgraz2 = 20.E-6 ! half sturation constant for meso grazing
+ epsher2 = 0.5 ! Efficicency of Mesozoo growth
+ unass2c = 0.3 ! non assimilated fraction of P by mesozoo
+ unass2n = 0.3 ! non assimilated fraction of N by mesozoo
+ unass2p = 0.3 ! non assimilated fraction of P by mesozoo
+ grazflux = 3.e3 ! flux-feeding rate
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampiszoo ! parameters for microzooplankton
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ part = 0.5 ! part of calcite not dissolved in microzoo gutsa
+ grazrat = 2.75 ! maximal zoo grazing rate
+ resrat = 0.03 ! exsudation rate of zooplankton
+ mzrat = 0.005 ! zooplankton mortality rate
+ xprefc = 0.2 ! Microzoo preference for POM
+ xprefn = 1. ! Microzoo preference for Nanophyto
+ xprefp = 1.5
+ xprefd = 1.0 ! Microzoo preference for Diatoms
+ xprefz = 0.3 ! Microzoo preference for microzooplankton
+ xthreshdia = 1.E-8 ! Diatoms feeding threshold for microzooplankton
+ xthreshphy = 1.E-8 ! Nanophyto feeding threshold for microzooplankton
+ xthreshpic = 1.E-8
+ xthreshzoo = 1.E-8 ! Nanophyto feeding threshold for microzooplankton
+ xthreshpoc = 1.E-8 ! POC feeding threshold for microzooplankton
+ xthresh = 3.E-7 ! Food threshold for feeding
+ xkgraz = 20.E-6 ! half sturation constant for grazing
+ epsher = 0.5 ! Efficiency of microzoo growth
+ unassc = 0.3 ! non assimilated fraction of C by zoo
+ unassn = 0.3 ! non assimilated fraction of C by zoo
+ unassp = 0.3 ! non assimilated fraction of C by zoo
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisfer ! parameters for iron chemistry
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ ln_fechem = .false. ! complex iron chemistry ( T/F )
+ ln_ligvar = .false. ! variable ligand concentration
+ xlam1 = 0.005 ! scavenging rate of Iron
+ xlamdust = 150.0 ! Scavenging rate of dust
+ ligand = 0.6E-9 ! Ligands concentration
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisrem ! parameters for remineralization
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ xremikc = 0.35 ! remineralization rate of DOC
+ xremikn = 0.35 ! remineralization rate of DON
+ xremikp = 0.35 ! remineralization rate of DOP
+ xremipc = 0.025 ! remineralisation rate of POC
+ xremipn = 0.03 ! remineralisation rate of PON
+ xremipp = 0.035 ! remineralisation rate of POP
+ nitrif = 0.05 ! NH4 nitrification rate
+ xsirem = 0.003 ! remineralization rate of Si
+ xsiremlab = 0.03 ! fast remineralization rate of Si
+ xsilab = 0.5 ! Fraction of labile biogenic silica
+ oxymin = 1.E-6 ! Half-saturation constant for anoxia
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampiscal ! parameters for Calcite chemistry
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ kdca = 6. ! calcite dissolution rate constant (1/time)
+ nca = 1. ! order of dissolution reaction (dimensionless)
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampissbc ! parameters for inputs deposition
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+! ! 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_dust = 'dust.orca' , -1 , 'dust' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_solub = 'solubility.orca' , -12 , 'solubility1' , .false. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdic = 'river.orca' , 120 , 'riverdic' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdoc = 'river.orca' , 120 , 'riverdoc' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdin = 'river.orca' , 120 , 'riverdin' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdon = 'river.orca' , 120 , 'riverdon' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdip = 'river.orca' , 120 , 'riverdip' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdop = 'river.orca' , 120 , 'riverdop' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_riverdsi = 'river.orca' , 120 , 'riverdsi' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ndepo = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' , ''
+ sn_ironsed = 'bathy.orca' , -12 , 'bathy' , .false. , .true. , 'yearly' , '' , '' , ''
+ sn_hydrofe = 'hydrofe.orca' , -12 , 'epsdb' , .false. , .true. , 'yearly' , '' , '' , ''
+!
+ cn_dir = './' ! root directory for the location of the dynamical files
+ ln_dust = .true. ! boolean for dust input from the atmosphere
+ ln_solub = .true. ! boolean for variable solubility of atm. Iron
+ ln_river = .true. ! boolean for river input of nutrients
+ ln_ndepo = .true. ! boolean for atmospheric deposition of N
+ ln_ironsed = .true. ! boolean for Fe input from sediments
+ ln_ironice = .true. ! boolean for Fe input from sea ice
+ ln_hydrofe = .true. ! boolean for from hydrothermal vents
+ sedfeinput = 1.5e-9 ! Coastal release of Iron
+ dustsolub = 0.02 ! Solubility of the dusta
+ mfrac = 0.035 ! Fe mineral fraction of dust
+ wdust = 2.0 ! Dust sinking speed
+ icefeinput = 15.e-9 ! Iron concentration in sea ice
+ nitrfix = 1.e-7 ! Nitrogen fixation rate
+ diazolight = 30. ! Diazotrophs sensitivity to light (W/m2)
+ concfediaz = 1.e-10 ! Diazotrophs half-saturation Cste for Iron
+ hratio = 1.e+7 ! Fe to 3He ratio assumed for vent iron supply
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampiskrp ! Kriest parameterization : parameters "key_kriest"
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ xkr_eta = 1.17 ! Sinking exponent
+ xkr_zeta = 2.28 ! N content exponent
+ xkr_ncontent = 5.7E-6 ! N content factor
+ xkr_mass_min = 0.0002 ! Minimum mass for Aggregates
+ xkr_mass_max = 1. ! Maximum mass for Aggregates
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampiskrs ! Kriest parameterization : size classes "key_kriest"
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ xkr_sfact = 942. ! Sinking factor
+ xkr_stick = 0.5 ! Stickiness
+ xkr_nnano = 2.337 ! Nbr of cell in nano size class
+ xkr_ndiat = 3.718 ! Nbr of cell in diatoms size class
+ xkr_nmeso = 7.147 ! Nbr of cell in mesozoo size class
+ xkr_naggr = 9.877 ! Nbr of cell in aggregates size class
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisdia ! additional 2D/3D tracers diagnostics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+! ! name ! title of the field ! units !
+! ! ! ! !
+ pisdia2d(1) = 'Cflx ' , 'DIC flux ', 'molC/m2/s '
+ pisdia2d(2) = 'Oflx ' , 'Oxygen flux ', 'molC/m2/s '
+ pisdia2d(3) = 'Kg ' , 'Gas transfer ', 'mol/m2/s/uatm'
+ pisdia2d(4) = 'Delc ' , 'Delta CO2 ', 'uatm '
+ pisdia2d(5) = 'PMO ' , 'POC export ', 'molC/m2/s '
+ pisdia2d(6) = 'PMO2 ' , 'GOC export ', 'molC/m2/s '
+ pisdia2d(7) = 'ExpFe1 ' , 'Nano iron export ', 'molFe/m2/s '
+ pisdia2d(8) = 'ExpFe2 ' , 'Diatoms iron export ', 'molFe/m2/s '
+ pisdia2d(9) = 'ExpSi ' , 'Silicate export ', 'molSi/m2/s '
+ pisdia2d(10) = 'ExpCaCO3 ' , 'Calcite export ', 'molC/m2/s '
+ pisdia2d(11) = 'heup ' , 'euphotic layer depth ', 'm '
+ pisdia2d(12) = 'Fedep ' , 'Iron dep ', 'molFe/m2/s '
+ pisdia2d(13) = 'Nfix ' , 'Nitrogen Fixation ', 'molN/m2/s '
+ pisdia3d(1) = 'PH ' , 'PH ', '- '
+ pisdia3d(2) = 'CO3 ' , 'Bicarbonates ', 'mol/l '
+ pisdia3d(3) = 'CO3sat ' , 'CO3 saturation ', 'mol/l '
+ pisdia3d(4) = 'PAR ' , 'light penetration ', 'W/m2 '
+ pisdia3d(5) = 'PPPHY ' , 'Primary production of nanophyto ', 'molC/m3/s '
+ pisdia3d(6) = 'PPPHY2 ' , 'Primary production of diatoms ', 'molC/m3/s '
+ pisdia3d(7) = 'PPNEWN ' , 'New Primary production of nano ', 'molC/m3/s '
+ pisdia3d(8) = 'PPNEWD ' , 'New Primary production of diat ', 'molC/m3/s '
+ pisdia3d(9) = 'PBSi ' , 'Primary production of Si diatoms ', 'molSi/m3/s '
+ pisdia3d(10) = 'PFeN ' , 'Primary production of nano iron ', 'molFe/m3/s '
+ pisdia3d(11) = 'PFeD ' , 'Primary production of diatoms iron', 'molFe/m3/s '
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisdmp ! Damping
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ ln_pisdmp = .true. ! Relaxation fo some tracers to a mean value
+ nn_pisdmp = 3650 ! Frequency of Relaxation
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampismass ! Mass conservation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ ln_check_mass = .false. ! Check mass conservation
+/
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! PISCES reduced (key_pisces_reduced, ex LOBSTER) : namelists
+!! 1 - biological parameters for phytoplankton (namlobphy)
+!! 2 - biological parameters for nutrients (namlobnut)
+!! 3 - biological parameters for zooplankton (namlobzoo)
+!! 4 - biological parameters for detritus (namlobdet)
+!! 5 - biological parameters for DOM (namlobdom)
+!! 6 - parameters from aphotic layers to sediment (namlobsed)
+!! 7 - general coefficients (namlobrat)
+!! 8 - optical parameters (namlobopt)
+
+!! 10 - biological diagnostics trends (namlobdbi)
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobphy ! biological parameters for phytoplankton
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ tmumax = 1.21e-5 ! maximal phytoplankton growth rate [s-1]
+ rgamma = 0.05 ! phytoplankton exudation fraction [%]
+ fphylab = 0.75 ! NH4 fraction of phytoplankton exsudation
+ tmminp = 5.8e-7 ! minimal phytoplancton mortality rate [0.05/86400 s-1=20 days]
+ aki = 33. ! light photosynthesis half saturation constant[W/m2]
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobnut ! biological parameters for nutrients
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ akno3 = 0.7 ! nitrate limitation half-saturation value [mmol/m3]
+ aknh4 = 0.001 ! ammonium limitation half-saturation value [mmol/m3]
+ taunn = 5.80e-7 ! nitrification rate [s-1]
+ psinut = 3. ! inhibition of nitrate uptake by ammonium
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobzoo ! biological parameters for zooplankton
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ rppz = 0.8 ! zooplankton nominal preference for phytoplancton food [%]
+ taus = 9.26E-6 ! specific zooplankton maximal grazing rate [s-1]
+! ! 0.75/86400 s-1=8.680555E-6 1/86400 = 1.15e-5
+ aks = 1. ! half-saturation constant for total zooplankton grazing [mmolN.m-3]
+ rpnaz = 0.3 ! non-assimilated phytoplankton by zooplancton [%]
+ rdnaz = 0.3 ! non-assimilated detritus by zooplankton [%]
+ tauzn = 8.1e-7 ! zooplancton specific excretion rate [0.1/86400 s-1=10 days]
+ fzoolab = 0.5 ! NH4 fraction of zooplankton excretion
+ fdbod = 0.5 ! zooplankton mortality fraction that goes to detritus
+ tmminz = 2.31e-6 ! minimal zooplankton mortality rate [(mmolN/m3)-1 d-1]
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobdet ! biological parameters for detritus
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ taudn = 5.80e-7 ! detritus breakdown rate [0.1/86400 s-1=10 days]
+ fdetlab = 0. ! NH4 fraction of detritus dissolution
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobdom ! biological parameters for DOM
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ taudomn = 6.43e-8 ! DOM breakdown rate [s-1]
+! ! slow remineralization rate of semi-labile dom to nh4 (1 month)
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobsed ! parameters from aphotic layers to sediment
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ sedlam = 3.86e-7 ! time coefficient of POC remineralization in sediments [s-1]
+ sedlostpoc = 0. ! mass of POC lost in sediments
+ vsed = 3.47e-5 ! detritus sedimentation speed [m/s]
+ xhr = -0.858 ! coeff for martin''s remineralisation profile
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobrat ! general coefficients
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ rcchl = 60. ! Carbone/Chlorophyl ratio [mgC.mgChla-1]
+ redf = 6.56 ! redfield ratio (C:N) for phyto
+ reddom = 6.56 ! redfield ratio (C:N) for DOM
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&namlobopt ! optical parameters
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+ xkg0 = 0.0232 ! green absorption coefficient of water
+ xkr0 = 0.225 ! red absorption coefficent of water
+ xkgp = 0.074 ! green absorption coefficient of chl
+ xkrp = 0.037 ! red absorption coefficient of chl
+ xlg = 0.674 ! green chl exposant for absorption
+ xlr = 0.629 ! red chl exposant for absorption
+ rpig = 0.7 ! chla/chla+pheo ratio
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisdbi ! biological diagnostics trends
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+! ! 2D bio diagnostics units : mmole/m2/s ("key_trdmld_trc")
+! ! name ! title of the field ! units !
+ pisdiabio(1) = 'NO3PHY' , 'Flux from NO3 to PHY ', 'mmole/m3/s'
+ pisdiabio(2) = 'NH4PHY' , 'Flux from NH4 to PHY ', 'mmole/m3/s'
+ pisdiabio(3) = 'PHYNH4' , 'Flux from PHY to NH4 ', 'mmole/m3/s'
+ pisdiabio(4) = 'PHYDOM' , 'Flux from PHY to DOM ', 'mmole/m3/s'
+ pisdiabio(5) = 'PHYZOO' , 'Flux from PHY to ZOO ', 'mmole/m3/s'
+ pisdiabio(6) = 'PHYDET' , 'Flux from PHY to DET ', 'mmole/m3/s'
+ pisdiabio(7) = 'DETZOO' , 'Flux from DET to ZOO ', 'mmole/m3/s'
+ pisdiabio(8) = 'DETSED' , 'Flux from DET to SED ', 'mmole/m3/s'
+ pisdiabio(9) = 'ZOODET' , 'Flux from ZOO to DET ', 'mmole/m3/s'
+ pisdiabio(10) = 'ZOOBOD' , 'Zooplankton closure ', 'mmole/m3/s'
+ pisdiabio(11) = 'ZOONH4' , 'Flux from ZOO to NH4 ', 'mmole/m3/s'
+ pisdiabio(12) = 'ZOODOM' , 'Flux from ZOO to DOM ', 'mmole/m3/s'
+ pisdiabio(13) = 'NH4NO3' , 'Flux from NH4 to NO3 ', 'mmole/m3/s'
+ pisdiabio(14) = 'DOMNH4' , 'Flux from DOM to NH4 ', 'mmole/m3/s'
+ pisdiabio(15) = 'DETNH4' , 'Flux from DET to NH4 ', 'mmole/m3/s'
+ pisdiabio(16) = 'DETDOM' , 'Flux from DET to DOM ', 'mmole/m3/s'
+ pisdiabio(17) = 'SEDNO3' , 'NO3 remineralization from SED', 'mmole/m3/s'
+/
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist_top
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist_top (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/namelist_top (revision 5267)
@@ -0,0 +1,176 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/TOP1 : - tracer run information (namtrc_run)
+!! - tracer definition (namtrc )
+!! - tracer data initialisation (namtrc_dta)
+!! - tracer advection (namtrc_adv)
+!! - tracer lateral diffusion (namtrc_ldf)
+!! - tracer vertical physics (namtrc_zdf)
+!! - tracer newtonian damping (namtrc_dmp)
+!! - dynamical tracer trends (namtrc_trd)
+!! - tracer output diagonstics (namtrc_dia)
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!-----------------------------------------------------------------------
+&namtrc_run ! run information
+!-----------------------------------------------------------------------
+ nn_dttrc = 1 ! time step frequency for passive sn_tracers
+ nn_writetrc = 5475 ! time step frequency for sn_tracer outputs
+ ln_top_euler = .true. ! use Euler time-stepping for TOP
+ ln_rsttr = .true. ! start from a restart file (T) or not (F)
+ nn_rsttr = 0 ! restart control = 0 initial time step is not compared to the restart file value
+ ! = 1 do not use the value in the restart file
+ ! = 2 calendar parameters read in the restart file
+ cn_trcrst_in = "restart_trc" ! suffix of pass. sn_tracer restart name (input)
+ cn_trcrst_out = "restart_trc" ! suffix of pass. sn_tracer restart name (output)
+/
+!-----------------------------------------------------------------------
+&namtrc ! tracers definition
+!-----------------------------------------------------------------------
+ ln_trcdta = .true. ! Initialisation from data input file (T) or not (F)
+ ln_trcdmp = .false. ! add a damping termn (T) or not (F)
+ ln_trcdmp_clo = .false. ! damping term (T) or not (F) on closed seas
+ sn_tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', 'eq/L ' , .true. , .true.
+ sn_tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(5) = 'PO4 ' , 'Phosphate Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(6) = 'POC ' , 'Small organic carbon Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(7) = 'Si ' , 'Silicate Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(13) = 'DSi ' , 'Diatoms Silicate Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(15) = 'BFe ' , 'Big iron particles Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(17) = 'SFe ' , 'Small iron particles Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(19) = 'GSi ' , 'Sinking biogenic Silicate Concentration', 'mol-C/L' , .false. , .true.
+ sn_tracer(20) = 'NFe ' , 'Nano iron Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(23) = 'NO3 ' , 'Nitrates Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(24) = 'NH4 ' , 'Ammonium Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(25) = 'DON ' , 'Dissolved Organic N Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(26) = 'DOP ' , 'Dissolved organic P Concentration ', 'mol-C/L' , .true. , .true.
+ sn_tracer(27) = 'PON ' , 'Small PON Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(28) = 'POP ' , 'Small POP Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(29) = 'GON ' , 'Big PON Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(30) = 'GOP ' , 'Big POP Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(31) = 'PHYN ' , 'PHYN Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(32) = 'PHYP ' , 'PHYN Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(33) = 'DIAN ' , 'PHYN Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(34) = 'DIAP ' , 'PHYN Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(35) = 'PIC ' , 'PICO Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(36) = 'PICN ' , 'PICO N Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(37) = 'PICP ' , 'PICO P Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(38) = 'PFe ' , 'PICO Fe Concentration ', 'mol-C/L' , .false. , .true.
+ sn_tracer(39) = 'PCHL ' , 'PICO Chl Concentration ', 'mol-C/L' , .false. , .true.
+
+/
+!-----------------------------------------------------------------------
+&namtrc_dta ! Initialisation from data input file
+!-----------------------------------------------------------------------
+!
+ sn_trcdta(1) = 'data_1m_DIC_nomask' , -1 , 'DIC' , .false. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(2) = 'data_1m_Alkalini_nomask' , -1 , 'Alkalini', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(3) = 'data_1m_O2_nomask' , -1 , 'O2' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(5) = 'data_1m_PO4_nomask' , -1 , 'PO4' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(7) = 'data_1m_Si_nomask' , -1 , 'Si' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(10) = 'data_1m_DOC_nomask' , -1 , 'DOC' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(14) = 'data_1m_Fer_nomask' , -1 , 'Fer' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(23) = 'data_1m_NO3_nomask' , -1 , 'NO3' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(25) = 'data_1m_DOC_nomask' , -1 , 'DOC' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_trcdta(26) = 'data_1m_DOC_nomask' , -1 , 'DOC' , .true. , .true. , 'yearly' , '' , '' , ''
+ cn_dir = './' ! root directory for the location of the data files
+ rn_trfac(1) = 1.0e-06 ! multiplicative factor
+ rn_trfac(2) = 1.0e-06 ! - - - -
+ rn_trfac(3) = 44.6e-06 ! - - - -
+ rn_trfac(5) = 122.0e-06 ! - - - -
+ rn_trfac(7) = 1.0e-06 ! - - - -
+ rn_trfac(10) = 1.0 ! - - - -
+ rn_trfac(14) = 1.0 ! - - - -
+ rn_trfac(23) = 7.6e-06 ! - - - -
+ rn_trfac(25) = 1.0 ! - - - -
+ rn_trfac(26) = 1.0 ! - - - -
+/
+!-----------------------------------------------------------------------
+&namtrc_adv ! advection scheme for passive tracer
+!-----------------------------------------------------------------------
+ ln_trcadv_cen2 = .false. ! 2nd order centered scheme
+ ln_trcadv_tvd = .false. ! TVD scheme
+ ln_trcadv_muscl = .true. ! MUSCL scheme
+ ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries
+ ln_trcadv_ubs = .false. ! UBS scheme
+ ln_trcadv_qck = .false. ! QUICKEST scheme
+ ln_trcadv_msc_ups = .false. ! use upstream scheme within muscl
+/
+!-----------------------------------------------------------------------
+&namtrc_ldf ! lateral diffusion scheme for passive tracer
+!-----------------------------------------------------------------------
+ ln_trcldf_diff = .true. ! performs lateral diffusion (T) or not (F)
+! ! Type of the operator :
+ ln_trcldf_lap = .true. ! laplacian operator
+ ln_trcldf_bilap = .false. ! bilaplacian operator
+ ! Direction of action :
+ ln_trcldf_level = .false. ! iso-level
+ ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T)
+ ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp")
+! ! Coefficient
+ rn_ahtrc_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s]
+ rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s]
+/
+!-----------------------------------------------------------------------
+&namtrc_zdf ! vertical physics
+!-----------------------------------------------------------------------
+ ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping
+ nn_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T
+/
+!-----------------------------------------------------------------------
+&namtrc_rad ! treatment of negative concentrations
+!-----------------------------------------------------------------------
+ ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namtrc_dmp ! passive tracer newtonian damping
+!-----------------------------------------------------------------------
+ nn_hdmp_tr = -1 ! horizontal shape =-1, damping in Med and Red Seas only
+ ! =XX, damping poleward of XX degrees (XX>0)
+ ! + F(distance-to-coast) + Red and Med Seas
+ nn_zdmp_tr = 1 ! vertical shape =0 damping throughout the water column
+ ! =1 no damping in the mixing layer (kz criteria)
+ ! =2 no damping in the mixed layer (rho crieria)
+ rn_surf_tr = 50. ! surface time scale of damping [days]
+ rn_bot_tr = 360. ! bottom time scale of damping [days]
+ rn_dep_tr = 800. ! depth of transition between rn_surf and rn_bot [meters]
+ nn_file_tr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0)
+/
+!-----------------------------------------------------------------------
+&namtrc_trd ! diagnostics on tracer trends ('key_trdtrc')
+! or mixed-layer trends ('key_trdmld_trc')
+!----------------------------------------------------------------------
+ nn_trd_trc = 5475 ! time step frequency and tracers trends
+ nn_ctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n /seconds ; =86400. -> /day)
+ ln_trdmld_trc_restart = .false. ! restart for ML diagnostics
+ ln_trdmld_trc_instant = .true. ! flag to diagnose trends of instantantaneous or mean ML T/S
+ ln_trdtrc(1) = .true.
+ ln_trdtrc(2) = .true.
+ ln_trdtrc(23) = .true.
+/
+!-----------------------------------------------------------------------
+&namtrc_dia ! parameters for passive tracer additional diagnostics
+!----------------------------------------------------------------------
+ ln_diatrc = .true. ! save additional diag. (T) or not (F)
+ ln_diabio = .true. ! output biological trends
+ nn_writedia = 3650 ! time step frequency for diagnostics
+ nn_writebio = 10 !: frequency of biological outputs
+/
+!----------------------------------------------------------------------
+! namtrc_bc ! data for boundary conditions
+!-----------------------------------------------------------------------
+&namtrc_bc
+!
+ cn_dir = './' ! root directory for the location of the data files
+/
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/xmlio_server.def
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/xmlio_server.def (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/EXP00/xmlio_server.def (revision 5267)
@@ -0,0 +1,35 @@
+! xmlio_server namelist
+!
+! using_server: .TRUE.(.FALSE.) to use(bypass) the io_server
+! using_oasis : .TRUE.(.FALSE.) if nemo is coupled through OASIS
+! client_id : used only for OASIS, NEMO id in the namecouple
+! server_id : used only for OASIS, io_server id in the namecouple
+!
+&coupling_param
+ using_server = .FALSE.
+ using_oasis = .FALSE.
+ client_id = 'nemo.x'
+ server_id = 'ionemo'
+/
+
+! global_mpi_buffer_size: size in Mo of the MPI buffer used by the io_server
+&mpi_param
+ global_mpi_buffer_size = 512
+/
+
+!!======================================================================
+!! namnc4 netcdf4 chunking and compression settings
+!!======================================================================
+!-----------------------------------------------------------------------
+&namnc4 ! netcdf4 chunking and compression settings
+ ! (benign if "key_netcdf4" is not used)
+!-----------------------------------------------------------------------
+ nn_nchunks_i = 4 ! number of chunks in i-dimension
+ nn_nchunks_j = 4 ! number of chunks in j-dimension
+ nn_nchunks_k = 31 ! number of chunks in k-dimension
+ ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
+ ! is optimal for postprocessing which works exclusively with horizontal slabs
+ ln_nc4zip = .TRUE. ! (T) use netcdf4 chunking and compression
+ ! (F) ignore chunking information and produce netcdf3-compatible files
+/
+
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/COMP/pisces.card
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/COMP/pisces.card (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/COMP/pisces.card (revision 5267)
@@ -0,0 +1,65 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+OPA_NDT_DAY=4
+#==============================================================
+#-- Name of directory related to dynamical
+#-- available directories are: v3, cpl
+#-- NOTE: if you use your own directory you have to comment it
+OPA_dyn=v3
+#==============================================================
+
+[InitialStateFiles]
+List= ()
+
+[BoundaryFiles]
+List= ()
+ListNonDel= (${R_BC}/${config_UserChoices_TagName}/DYN/${pisces_UserChoices_OPA_dyn}/NEMOV3_mesh_mask.nc, mesh_mask.nc), \
+ (${R_BC}/${config_UserChoices_TagName}/DYN/${pisces_UserChoices_OPA_dyn}/NEMOV3_5d_21210101_21211231_grid_T.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/DYN/${pisces_UserChoices_OPA_dyn}/NEMOV3_5d_21210101_21211231_grid_U.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/DYN/${pisces_UserChoices_OPA_dyn}/NEMOV3_5d_21210101_21211231_grid_V.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/DYN/${pisces_UserChoices_OPA_dyn}/NEMOV3_5d_21210101_21211231_grid_W.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_DIC_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_Alkalini_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_O2_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_NO3_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_PO4_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_Si_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_DOC_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/INI/data_1m_Fer_nomask.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/ANC/kRGB61.txt, .), \
+ (${R_BC}/${config_UserChoices_TagName}/ANC/dust.orca.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/ANC/bathy.orca.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/ANC/river.orca.nc, .), \
+ (${R_BC}/${config_UserChoices_TagName}/ANC/ndeposition.orca.nc, .)
+
+[ParametersFiles]
+List= (${SUBMIT_DIR}/PARAM/namelist, .), \
+ (${SUBMIT_DIR}/PARAM/namelist_top, .), \
+ (${SUBMIT_DIR}/PARAM/namelist_pisces, .), \
+ (${SUBMIT_DIR}/PARAM/iodef.xml, .), \
+ (${SUBMIT_DIR}/PARAM/xmlio_server.def, .)
+
+[RestartFiles]
+List= (${config_UserChoices_JobName}_${NEMO_END}_restart_trc.nc, restart_trc.nc, restart_trc.nc)
+
+[OutputText]
+List= (ocean.output)
+
+[OutputFiles]
+List= (${PREFIX_WFT1}_ptrc_T.nc, ${R_OUT_MBG_WFT1}/${PREFIX}_${WFT1}_ptrc_T.nc, NONE), \
+ (${PREFIX_WFT1}_diad_T.nc, ${R_OUT_MBG_WFT1}/${PREFIX}_${WFT1}_diad_T.nc, NONE), \
+ (${PREFIX_WFT2}_ptrc_T.nc, ${R_OUT_MBG_WFT2}/${PREFIX}_${WFT2}_ptrc_T.nc, Post_1M_ptrc_T), \
+ (${PREFIX_WFT3}_diad_T.nc, ${R_OUT_MBG_WFT3}/${PREFIX}_${WFT3}_diad_T.nc, Post_5D_diad_T)
+
+[Post_1M_ptrc_T]
+Patches= ()
+GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter)
+TimeSeriesVars = (DCHL, NCHL, NO3)
+
+[Post_5D_diad_T]
+Patches = ()
+GatherWithInternal = (nav_lon, nav_lat, time_counter)
+TimeSeriesVars = (Cflx)
+
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/COMP/pisces.driver
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/COMP/pisces.driver (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/COMP/pisces.driver (revision 5267)
@@ -0,0 +1,192 @@
+#!/bin/ksh
+
+#-----------------------------------
+function MBG_Initialize
+{
+ IGCM_debug_PushStack "MBG_Initialize"
+
+ JOB_NAME=${config_UserChoices_JobName}
+ OPA_NDT_DAY=${pisces_UserChoices_OPA_NDT_DAY:=4}
+
+# Local function to find namelists parameters
+supergrep () {
+ grep "^ *$1 *=" $2 | sed -e "s% *\!.*%%"
+}
+
+ ##--Variables used by OPA --
+
+ # cexper experience name for vairmer format
+ # nit000 number of the first time step
+ # nitend number of the last time step
+ # nleapy leap year calendar (0/1) (30 for 360d)
+ # nwrite frequency of OUTPUT file
+ # ln_rstart boolean term for restart (true or false)
+ # nstock frequency of restart file
+ # ndate0 initial calendar date aammjj
+ # rdt time step in seconds (coming from namelist)
+
+
+ PAT_CEXPER=$( supergrep cn_exp ${SUBMIT_DIR}/PARAM/namelist )
+ PAT_NIT000=$( supergrep nn_it000 ${SUBMIT_DIR}/PARAM/namelist )
+ PAT_NITEND=$( supergrep nn_itend ${SUBMIT_DIR}/PARAM/namelist )
+ PAT_NLEAPY=$( supergrep nn_leapy ${SUBMIT_DIR}/PARAM/namelist )
+ PAT_NSTOCK=$( supergrep nn_stock ${SUBMIT_DIR}/PARAM/namelist )
+ PAT_NDATE0=$( supergrep nn_date0 ${SUBMIT_DIR}/PARAM/namelist )
+
+
+ OPA_RDT=$( supergrep rn_rdt ${SUBMIT_DIR}/PARAM/namelist | sed 's/ *rn_rdt *=//' | sed 's/\. *,//' )
+
+ # Period Length In Days between DateBegin and first day of calendar 0001 01 01
+ # /!\ Needed by OPA namelist to compute file names /!\
+ (( DaysSinceJC = $( IGCM_date_DaysSinceJC ${DateBegin} ) + 1 ))
+
+ ##--Variables used by MBG --
+
+ # lrsttr boolean term for tracer model restart (true or false)
+ # nrsttr control of the time step (0, 1 or 2)
+
+ ##--Variables used by PISCES --
+
+ PAT_TOP_LRSTTR=$( supergrep ln_rsttr ${SUBMIT_DIR}/PARAM/namelist_top )
+ PAT_TOP_NRSTTR=$( supergrep nn_rsttr ${SUBMIT_DIR}/PARAM/namelist_top )
+ PAT_PIS_LRIVER=$( supergrep ln_river ${SUBMIT_DIR}/PARAM/namelist_pisces )
+
+ ##--Variables used for XMLIO_SERVER configuration file
+
+ # using_server: .TRUE.(.FALSE.) to use(bypass) the io_server
+
+ PAT_SERVER=$( supergrep using_server ${SUBMIT_DIR}/PARAM/xmlio_server.def )
+
+
+ IGCM_debug_PopStack "MBG_Initialize"
+}
+
+#-----------------------------------
+function MBG_Update
+{
+ IGCM_debug_PushStack "MBG_Update"
+
+ # Period Length In Days between DateBegin and first day of calendar 0001 01 01
+ # Needed by OPA namelist to compute correct file names
+ (( PeriodDaysSinceJC = $( IGCM_date_DaysSinceJC ${PeriodDateBegin} ) + 1 ))
+
+ ## Count of number of time steps from begin of OPA job
+ (( OPA_NIT000 = ( PeriodDaysSinceJC - DaysSinceJC ) * OPA_NDT_DAY + 1 ))
+ (( OPA_NDT = PeriodLengthInDays * OPA_NDT_DAY ))
+ (( OPA_NITEND = OPA_NIT000 + OPA_NDT - 1))
+
+ ## Forcing to write one restart file at the end of period
+ OPA_NSTOCK="${OPA_NITEND}"
+
+ ## Verification of number of time steps per day
+ (( NB_SEC_DAY_MODEL = OPA_NDT_DAY * OPA_RDT ))
+ (( NB_SEC_DAY = 60 * 60 * 24 ))
+
+ if [ ${NB_SEC_DAY_MODEL} -ne ${NB_SEC_DAY} ]
+ then
+ echo " VERIFY OPA_NDT_DAY in pisces.card "
+ exit
+ fi
+
+ ## nleapy configuration
+ case ${config_UserChoices_CalendarType} in
+ leap)
+ OPA_NLEAPY=1;;
+ noleap)
+ OPA_NLEAPY=0;;
+ 360d)
+ OPA_NLEAPY=30;;
+ esac
+
+
+ typeset -r PRECIS=8
+ NEMO_END=$( echo $( awk "BEGIN { printf \"%0${PRECIS}d\",${OPA_NITEND} }" ) )
+
+ sed -e "s%${PAT_CEXPER}% cn_exp=\"${config_UserChoices_JobName}\"%" \
+ -e "s%${PAT_NIT000}% nn_it000=${OPA_NIT000}%" \
+ -e "s%${PAT_NITEND}% nn_itend=${OPA_NITEND}%" \
+ -e "s%${PAT_NSTOCK}% nn_stock=${OPA_NSTOCK}%" \
+ -e "s%${PAT_NDATE0}% nn_date0=${PeriodDateBegin}%" \
+ -e "s%${PAT_NLEAPY}% nn_leapy=${OPA_NLEAPY}%" \
+ namelist > namelist.tmp
+
+ IGCM_sys_Mv namelist.tmp namelist
+
+ ##--Write Frequency Purpose
+ Nbfreq=$( echo ${config_MBG_WriteFrequency} | wc -w )
+ typeset -Z1 i=1
+ while [ $i -le ${Nbfreq} ]; do
+ topfrequency=$( echo ${config_MBG_WriteFrequency} | awk -v var="$i" "-F " '{print $var}' )
+ topfactor=$( echo ${topfrequency} | sed -e "s/[yYmMdD]//" )
+ case ${topfrequency} in
+ *Y|*y)
+ eval PREFIX_WFT${i}=${config_UserChoices_JobName}_${topfactor}y_${DatesPeriod} ;
+ eval R_OUT_MBG_WFT${i}=${R_OUT_MBG_O_Y} ;
+ eval WFT${i}=${topfactor}Y ;;
+ *M|*m)
+ eval PREFIX_WFT${i}=${config_UserChoices_JobName}_${topfactor}m_${DatesPeriod} ;
+ eval R_OUT_MBG_WFT${i}=${R_OUT_MBG_O_M} ;
+ eval WFT${i}=${topfactor}M ;;
+ *D|*d)
+ eval PREFIX_WFT${i}=${config_UserChoices_JobName}_${topfactor}d_${DatesPeriod} ;
+ eval R_OUT_MBG_WFT${i}=${R_OUT_MBG_O_D} ;
+ eval WFT${i}=${topfactor}D ;;
+ *)
+ esac
+ let i=$i+1
+ done
+
+ ##-- Restart configuration
+ if ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_MBG_Restart}" = "n" ] ) ; then
+
+ #echo "NO RESTART FOR TOP"
+ TOP_LRSTTR=.FALSE.
+ TOP_NRSTTR=0
+
+ PISCES_LRIVER=.FALSE.
+
+ elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_MBG_Restart}" = "y" ] ) ; then
+
+ #echo "RESTART TOP"
+ TOP_LRSTTR=.TRUE.
+ TOP_NRSTTR=0
+
+ PISCES_LRIVER=.TRUE.
+
+ else
+
+ #echo "RESTART TOP"
+ TOP_LRSTTR=.TRUE.
+ TOP_NRSTTR=2
+
+ PISCES_LRIVER=.TRUE.
+
+ fi
+
+# on force les sorties PISCES comme les sorties OPA
+
+ sed -e "s%${PAT_TOP_LRSTTR}% ln_rsttr=${TOP_LRSTTR}%" \
+ -e "s%${PAT_TOP_NRSTTR}% nn_rsttr=${TOP_NRSTTR}%" \
+ namelist_top > namelist_top.tmp
+
+ IGCM_sys_Mv namelist_top.tmp namelist_top
+
+ sed -e "s%${PAT_PIS_LRIVER}% ln_river=${PISCES_LRIVER}%" \
+ namelist_pisces > namelist_pisces.tmp
+
+ IGCM_sys_Mv namelist_pisces.tmp namelist_pisces
+
+
+ IGCM_debug_PopStack "MBG_Update"
+}
+
+#-----------------------------------
+function MBG_Finalize
+{
+ IGCM_debug_PushStack "MBG_Finalize"
+
+ echo FINALIZE MBG !!!
+
+ IGCM_debug_PopStack "MBG_Finalize"
+}
+
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/README
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/README (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/README (revision 5267)
@@ -0,0 +1,58 @@
+ config.card
+ ORCA2_OFF_PISCES
+ 5 years run - trunk of 1 year per job with diagnostics ( key_trc_diaadd & key_trc_dia3d activated )
+
+ - NEMO trunk revision 1656
+ - mpi on 4 processors ( 1x4 domain decomposition )
+ - libIGCM trunk revision 183
+ - dynamical files from IPSLCM5-VA66 runs from Olivier Marti
+
+without monitoring and atlas
+
+Tested on NEC SX9 mercure at CCRT.
+
+Ouputs on : /dmnfs10/cont003/p48ethe/IGCM_OUT/ORCA2_OFF_PISCES/VA66P
+
+========================================================================
+run.card
+
+# $Date: $
+# $Name: $
+# $Revision: $
+# IPSL (2006)
+# This software is governed by the CeCILL licence see libIGCM/libIGCM_CeCILL.LIC
+
+#========================================================================
+[Configuration]
+#last PREFIX
+OldPrefix= VA66P_20051231
+#Compute date of loop
+PeriodDateBegin= 2006-01-01
+PeriodDateEnd= 2006-12-31
+CumulPeriod= 6
+# State of Job "Start", "Running", "OnQueue", "Completed"
+PeriodState= Completed
+
+#========================================================================
+[PostProcessing]
+
+#RebuildRunning=y/n
+#RebuildToComplete=number of dir to rebuild
+#RebuildDateID=PeriodDateEnd correspondant au dernier flag de RebuildFrequency=true
+
+TimeSeriesRunning=n
+TimeSeriesCompleted=
+
+#========================================================================
+[Log]
+# Executables Size
+LastExeSize= ( 11596207 )
+
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# CumulPeriod | PeriodDateBegin | PeriodDateEnd | RunDateBegin | RunDateEnd | RealCpuTime | UserCpuTime | SysCpuTime | ExeDate
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# 1 | 20010101 | 20011231 | 2009-10-14T15:05:34 | 2009-10-14T15:10:39 | 304.19207 | 1191.95484 | 4.41881 | MBG_Oct_14_15:02
+# 2 | 20020101 | 20021231 | 2009-10-14T15:14:18 | 2009-10-14T15:20:13 | 354.59008 | 1412.65806 | 2.16025 | MBG_Oct_14_15:02
+# 3 | 20030101 | 20031231 | 2009-10-14T15:24:17 | 2009-10-14T15:30:10 | 353.09587 | 1407.03031 | 2.22349 | MBG_Oct_14_15:02
+# 4 | 20040101 | 20041231 | 2009-10-14T15:34:16 | 2009-10-14T15:40:18 | 362.49957 | 1430.77230 | 3.31767 | MBG_Oct_14_15:02
+# 5 | 20050101 | 20051231 | 2009-10-14T15:44:16 | 2009-10-14T15:50:17 | 360.82695 | 717.62953 | 1.79535 | MBG_Oct_14_15:02
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/config.card
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/config.card (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/IGCM00/config.card (revision 5267)
@@ -0,0 +1,96 @@
+#
+# This is config.card file for NEMO ORCA2_LIM_PISCES configuration
+#
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#D-- UserChoices -
+[UserChoices]
+#============================
+#-- (<8 chars MAX for JobName)
+JobName=PISCES
+LongName="SCRIPT_V1"
+TagName=ORCA2_OFF_PISCES
+#============================
+#-- leap, noleap, 360d
+CalendarType=noleap
+#-- Begin and end of Job
+#-- "YYYY-MM-DD"
+DateBegin=2001-01-01
+DateEnd=2005-12-31
+#============================
+#-- 1Y, 1M, 5D, 1D
+PeriodLength=1Y
+#============================
+#-- PBS Class (required on brodie mono/multi)
+JobClass=mono
+#JobClass=multi
+#============================
+#-- Total Number of Processors
+JobNumProcTot=1
+#JobNumProcTot=8
+#============================
+#-- Run Options :
+#JobRunOptions=''
+JobRunOptions='"-np "${BATCH_NUM_PROC_TOT}"'
+#============================
+#============================
+# config.card configuration options :
+# config.card configuration options :
+###R_BC=/u/rech/eee/reee451/IGCM/BC
+R_BC=/dmnfs/cont003/p48ethe/IGCM/BC
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+MBG= (pisces, PISCES)
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=nemo
+#D- For each component, Real name of executable
+MBG= (nemo.exe, nemo)
+#========================================================================
+#D-- Restarts -
+[Restarts]
+#D- If you want a GENERAL RULE FOR RESTARTS, put this flag to 'y'
+OverRule=n
+#D- Last day of the experience used as restart
+RestartDate=2000-12-31
+#D- Define restart simulation name
+RestartJobName=OROFP
+#D- Path Server Group Login
+RestartPath=${ARCHIVE}/IGCM_OUT/ORCA2_OFF_PISCES
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission
+#D- BETA FEATURE
+RebuildFrequency=NONE
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+RebuildFromArchive=false
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=NONE
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=NONE
+#========================================================================
+#D-- MBG -
+[MBG]
+##-- frequency
+WriteFrequency="1Y 1M 5D"
+##-- restart
+Restart= n
+##-- Last day of the experience used as restart
+RestartDate=2000-12-31
+# Define restart simulation name
+RestartJobName=OROFP1
+RestartPath=${ARCHIVE}/IGCM_OUT/ORCA2_OFF_PISCES
+# Old component name for restart (if empty, use new name)
+OldName=
+#========================================================================
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/cpp_ORCA2_OFF_PISCES_QUOTA.fcm
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/cpp_ORCA2_OFF_PISCES_QUOTA.fcm (revision 5267)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_QUOTA/cpp_ORCA2_OFF_PISCES_QUOTA.fcm (revision 5267)
@@ -0,0 +1,1 @@
+bld::tool::fppkeys key_trabbl key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces_quota key_iomput key_mpp_mpi
Index: branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/cfg.txt
===================================================================
--- branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/cfg.txt (revision 5266)
+++ branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/CONFIG/cfg.txt (revision 5267)
@@ -4,4 +4,5 @@
GYRE_XIOS OPA_SRC
ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC
+ORCA2_OFF_PISCES_QUOTA OPA_SRC OFF_SRC TOP_SRC
ORCA2_SAS_LIM OPA_SRC SAS_SRC LIM_SRC_2 NST_SRC
C1D_PAPA OPA_SRC