CNRS - NERC - Met OFFICE - MERCATOR-ocean - INGV - CMCC NEMO team Ocean General Circulation Model version 3.6 (2015) mynode : mpi initialisation ~~~~~~ Namelist nammpp mpi send type cn_mpi_send = I size in bytes of exported buffer nn_buffer = 0 jpni, jpnj and jpnij will be calculated automatically avoid use of mpi_allgather at the north fold ln_nnogather = F Immediate non-blocking send (isend) AAAAAAAA nemo_ctl: Control prints & Benchmark ~~~~~~~ Namelist namctl run control (for debugging) ln_ctl = F level of print nn_print = 0 Start i indice for SUM control nn_ictls = 0 End i indice for SUM control nn_ictle = 0 Start j indice for SUM control nn_jctls = 0 End j indice for SUM control nn_jctle = 0 number of proc. following i nn_isplt = 1 number of proc. following j nn_jsplt = 1 benchmark parameter (0/1) nn_bench = 0 timing activated (0/1) nn_timing = 0 namcfg : configuration initialization through namelist read ~~~~~~~ Namelist namcfg configuration name cp_cfg = gyre configuration zoom name cp_cfz = no zoom configuration resolution jp_cfg = 1 1st lateral dimension ( >= jpi ) jpidta = 32 2nd " " ( >= jpj ) jpjdta = 22 3nd " " jpkdta = 31 1st dimension of global domain in i jpiglo = 32 2nd - - in j jpjglo = 22 left bottom i index of the zoom (in data domain) jpizoom = 1 left bottom j index of the zoom (in data domain) jpizoom = 1 lateral cond. type (between 0 and 6) jperio = 0 use file attribute if exists as i/p j-start ln_use_jattr = F mpp_init : Message Passing MPI ~~~~~~~~ mpp_init: defines mpp subdomains ~~~~~~ ---------------------- iresti= 2 irestj= 2 jpni= 4 jpnj= 2 sum ilcit(i,1)= 32.0000000000000 jpiglo= 32 sum ilcit(1,j)= 22.0000000000000 jpjglo= 22 **************************************************************** * * * * * 2 * 10 x 12 * 10 x 12 * 9 x 12 * 9 x 12 * * 4 * 5 * 6 * 7 * * * * * * **************************************************************** * * * * * 1 * 10 x 12 * 10 x 12 * 9 x 12 * 9 x 12 * * 0 * 1 * 2 * 3 * * * * * * **************************************************************** 1 2 3 4 nproc = 0 nowe = -1 noea = 1 nono = 4 noso = -4 nbondi = -1 nbondj = -1 npolj = 0 nperio = 0 nlci = 10 nlcj = 12 nimpp = 1 njmpp = 1 nreci = 2 npse = -3 nrecj = 2 npsw = -5 jpreci = 1 npne = 5 jprecj = 1 npnw = 3 mpp_init_ioipsl : iloc = 10 12 ~~~~~~~~~~~~~~~ iabsf = 1 1 ihals = 0 0 ihale = 1 1 phy_cst : initialization of ocean parameters and constants ~~~~~~~ Domain info dimension of model Local domain Global domain Data domain jpi : 10 jpiglo : 32 jpidta : 32 jpj : 12 jpjglo : 22 jpjdta : 22 jpk : 31 jpk : 31 jpkdta : 31 jpij : 120 mpp local domain info (mpp) jpni : 4 jpreci : 1 jpnj : 2 jprecj : 1 jpnij : 8 lateral domain boundary condition type : jperio = 0 Constants mathematical constant rpi = 3.14159265358979 day rday = 86400.0000000000 s sideral year rsiyea = 31558149.0101107 s sideral day rsiday = 86164.0996559118 s omega omega = 7.292115083046062E-005 s^-1 nb of months per year raamo = 12.0000000000000 months nb of hours per day rjjhh = 24.0000000000000 hours nb of minutes per hour rhhmm = 60.0000000000000 mn nb of seconds per minute rmmss = 60.0000000000000 s earth radius ra = 6371229.00000000 m gravity grav = 9.80665000000000 m/s^2 triple point of temperature rtt = 273.160000000000 K freezing point of water rt0 = 273.150000000000 K melting point of snow rt0_snow = 273.150000000000 K melting point of ice rt0_ice = 273.050000000000 K reference density and heat capacity now defined in eosbn2.f90 thermal conductivity of pure ice = 2.03439600000000 J/s/m/K fresh ice specific heat = 2093.00000000000 J/kg/K latent heat of fusion of fresh ice / snow = 333700.000000000 J/kg density times specific heat for snow = 690690.000000000 J/m^3/K density times specific heat for ice = 1883700.00000000 J/m^3/K volumetric latent heat fusion of sea ice = 300330000.000000 J/m latent heat of sublimation of snow = 2800000.00000000 J/kg volumetric latent heat fusion of snow = 110121000.000000 J/m^3 density of sea ice = 900.000000000000 kg/m^3 density of snow = 330.000000000000 kg/m^3 emissivity of snow or ice = 0.970000000000000 salinity of ice = 6.00000000000000 psu salinity of sea = 34.7000000000000 psu latent heat of evaporation (water) = 2500000.00000000 J/m^3 correction factor for solar radiation = 0.900000000000000 von Karman constant = 0.400000000000000 Stefan-Boltzmann constant = 5.670000000000000E-008 J/s/m^2/K^4 conversion: degre ==> radian rad = 1.745329251994330E-002 smallest real computer value rsmall = 1.110223024625157E-016 eos_init : equation of state ~~~~~~~~ Namelist nameos : set eos parameters flag for eq. of state and N^2 nn_eos = 0 model does not use Conservative Temperature use of EOS-80 equation of state (pot. temp. and pract. salinity) volumic mass of reference rau0 = 1026.00000000000 kg/m^3 1. / rau0 r1_rau0 = 9.746588693957114E-004 m^3/kg ocean specific heat rcp = 3991.86795711963 J/Kelvin rau0 * rcp rau0_rcp = 4095656.52400474 1. / ( rau0 * rcp ) r1_rau0_rcp = 2.441610994816036E-007 dom_cfg : set the ocean configuration ~~~~~~~ ocean model configuration used : cp_cfg = gyre jp_cfg = 1 global domain lateral boundaries jperio= 0, closed dom_glo : domain: data / local ~~~~~~~ data input domain : jpidta = 32 jpjdta = 22 jpkdta = 31 global or zoom domain: jpiglo = 32 jpjglo = 22 jpk = 31 local domain : jpi = 10 jpj = 12 jpk = 31 south-west indices jpizoom = 1 jpjzoom = 1 conversion local ==> data i-index domain 1 2 3 4 5 6 7 8 9 10 conversion data ==> local i-index domain starting index 1 2 3 4 5 6 7 8 9 10 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 ending index 1 2 3 4 5 6 7 8 9 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 conversion local ==> data j-index domain 1 2 3 4 5 6 7 8 9 10 11 12 conversion data ==> local j-index domain starting index 1 2 3 4 5 6 7 8 9 10 11 12 13 13 13 13 13 13 13 13 13 13 ending index 1 2 3 4 5 6 7 8 9 10 11 12 12 12 12 12 12 12 12 12 12 12 zoom flags : lzoom = F (T = zoom, F = global ) lzoom_e = F (T = forced closed east boundary) lzoom_w = F (T = forced closed west boundary) lzoom_s = F (T = forced closed South boundary) lzoom_n = F (T = forced closed North boundary) dom_init : domain initialization ~~~~~~~~ dom_nam : domain initialization through namelist read ~~~~~~~ Namelist namrun job number nn_no = 0 experiment name for output cn_exp = GYRE file prefix restart input cn_ocerst_in= restart restart input directory cn_ocerst_indir= . file prefix restart output cn_ocerst_out= restart restart output directory cn_ocerst_outdir= . restart logical ln_rstart = F start with forward time step nn_euler = 1 control of time step nn_rstctl = 0 number of the first time step nn_it000 = 1 number of the last time step nn_itend = 4320 initial calendar date aammjj nn_date0 = 10101 leap year calendar (0/1) nn_leapy = 30 initial state output nn_istate = 0 frequency of restart file nn_stock = 4320 frequency of output file nn_write = 60 multi file dimgout ln_dimgnnn = F mask land points ln_mskland = F additional CF standard metadata ln_cfmeta = F overwrite an existing file ln_clobber = T NetCDF chunksize (bytes) nn_chunksz = 0 ===>>> : W A R N I N G =============== ln_rstart =.FALSE., nn_euler is forced to 0 The IOIPSL calendar is "360d", i.e. 360 days in a year Namelist namdom : space & time domain flag read/compute bathymetry nn_bathy = 0 Depth (if =0 bathy=jpkm1) rn_bathy = 0.000000000000000E+000 min depth of the ocean (>0) or rn_hmin = -3.00000000000000 min number of ocean level (<0) minimum thickness of partial rn_e3zps_min = 20.0000000000000 (m) step level rn_e3zps_rat = 0.100000000000000 create mesh/mask file(s) nn_msh = 1 = 0 no file created = 1 mesh_mask = 2 mesh and mask = 3 mesh_hgr, msh_zgr and mask ocean time step rn_rdt = 7200.00000000000 asselin time filter parameter rn_atfp = 0.100000000000000 acceleration of converge nn_acc = 0 nn_acc=1: surface tracer rdt rn_rdtmin = 7200.00000000000 bottom tracer rdt rdtmax = 7200.00000000000 depth of transition rn_rdth = 800.000000000000 suppression of closed seas (=0) nn_closea = 0 online coarsening of dynamical fields ln_crs = F type of horizontal mesh jphgr_msh = 5 longitude of first raw and column T-point ppglam0 = 0.000000000000000E+000 latitude of first raw and column T-point ppgphi0 = 29.0000000000000 zonal grid-spacing (degrees) ppe1_deg = 999999.000000000 meridional grid-spacing (degrees) ppe2_deg = 999999.000000000 zonal grid-spacing (degrees) ppe1_m = 999999.000000000 meridional grid-spacing (degrees) ppe2_m = 999999.000000000 ORCA r4, r2 and r05 coefficients ppsur = -2033.19429528338 ppa0 = 155.832536966415 ppa1 = 146.361591860189 ppkth = 17.2852037241979 ppacr = 5.00000000000000 Minimum vertical spacing ppdzmin = 999999.000000000 Maximum depth pphmax = 999999.000000000 Use double tanf function for vertical coordinates ldbletanh = F Double tanh function parameters ppa2 = 999999.000000000 ppkth2 = 999999.000000000 ppacr2 = 999999.000000000 Namelist namcla cross land advection nn_cla = 0 dom_clo : closed seas ~~~~~~~ dom_hgr : define the horizontal mesh from ithe following par_oce parameters ~~~~~~~ type of horizontal mesh jphgr_msh = 5 position of the first row and ppglam0 = 0.000000000000000E+000 column grid-point (degrees) ppgphi0 = 29.0000000000000 zonal grid-spacing (degrees) ppe1_deg = 999999.000000000 meridional grid-spacing (degrees) ppe2_deg = 999999.000000000 zonal grid-spacing (meters) ppe1_m = 999999.000000000 meridional grid-spacing (meters) ppe2_m = 999999.000000000 beta-plane with regular grid-spacing and rotated domain (GYRE configu ration) given by ppe1_m and ppe2_m longitude and e1 scale factors ------------------------------ 1 -71.52 -71.18 -71.86 -71.52 106000.0000000000 106000.0000000000 106000.0000000000 106000.0000000000 latitude and e2 scale factors ----------------------------- 1 14.85 15.18 15.18 15.52 106000.0000000000 106000.0000000000 106000.0000000000 106000.0000000000 11 21.59 21.92 21.92 22.26 106000.0000000000 106000.0000000000 106000.0000000000 106000.0000000000 Beta-plane and rotated domain : Coriolis parameter varies in this processor from 3.890233185442712E-005 to 5.390852115542532E-005 Coriolis parameter varies globally from 3.890233185442712E-005 to 1.064301837089190E-004 dom_zgr : vertical coordinate ~~~~~~~ Namelist namzgr : set vertical coordinate z-coordinate - full steps ln_zco = T z-coordinate - partial steps ln_zps = F s- or hybrid z-s-coordinate ln_sco = F ice shelf cavities ln_isfcav = F zgr_z : Reference vertical z-coordinates ~~~~~~~ Value of coefficients for vertical mesh: zsur = -2033.19429528338 za0 = 155.832536966415 za1 = 146.361591860189 zkth = 17.2852037241979 zacr = 5.00000000000000 Reference z-coordinate depth and scale factors: level gdept_1d gdepw_1d e3t_1d e3w_1d 1 4.98 0.00 10.00 9.90 2 15.10 10.00 10.26 10.12 3 25.54 20.27 10.65 10.43 4 36.45 30.92 11.22 10.90 5 48.07 42.15 12.07 11.60 6 60.73 54.24 13.34 12.64 7 74.94 67.60 15.20 14.18 8 91.42 82.84 17.93 16.44 9 111.22 100.82 21.92 19.74 10 135.84 122.82 27.66 24.54 11 167.35 150.58 35.81 31.39 12 208.50 186.52 47.09 41.01 13 262.83 233.77 62.26 54.15 14 334.49 296.22 81.81 71.47 15 427.90 378.20 105.69 93.24 16 547.03 484.04 133.03 119.02 17 694.55 617.15 162.12 147.49 18 871.10 779.24 190.71 176.61 19 1075.10 969.84 216.73 204.15 20 1303.24 1186.41 238.81 228.31 21 1551.23 1425.03 256.44 248.18 22 1814.70 1681.29 269.86 263.64 23 2089.75 1951.00 279.70 275.18 24 2373.16 2230.58 286.72 283.52 25 2662.49 2517.21 291.64 289.40 26 2955.93 2808.79 295.03 293.50 27 3252.20 3103.78 297.36 296.31 28 3550.40 3401.10 298.93 298.22 29 3849.90 3700.01 300.00 299.52 30 4150.28 4000.00 300.72 300.40 31 4451.26 4300.71 301.20 300.99 zgr_bat : defines level and meter bathymetry ~~~~~~~ bathymetry field: flat basin Depth = depthw(jpkm1) Minimum ocean depth: 30.9217699403011 minimum number of ocean levels : 3 zgr_bat_ctl : check the bathymetry ~~~~~~~~~~~ suppress isolated ocean grid points ----------------------------------- no isolated ocean grid points mbathy set to 0 along east and west boundary: nperio = 0 zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ~~~~~~~~~~~~~ zgr_top_level : ocean top k-index of T-, U-, V- and W-levels ~~~~~~~~~~~~~ dommsk : ocean mask ~~~~~~ Namelist namlbc lateral momentum boundary cond. rn_shlat = 0.000000000000000E+000 consistency with analytical form ln_vorlat = F ocean lateral free-slip dom_stp : time stepping setting ~~~~~~~ synchronous time stepping dynamics and tracer time step = 2.00000000000000 hours dom_wri : create NetCDF mesh and mask information file(s) ~~~~~~~ iom_nf90_open ~~~ create new file: mesh_mask_0000.nc in WRI TE mode ---> mesh_mask_0000.nc OK iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask define dimensio n variables done iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask write dimension variables done iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umask defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umask written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmask defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmask written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmask defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmask written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmaskutil defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmaskutil written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umaskutil defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umaskutil written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmaskutil defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmaskutil written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmaskutil defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmaskutil written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamt defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamt written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamu defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamu written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamv defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamv written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamf defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamf written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphit defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphit written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiu defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiu written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiv defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiv written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphif defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphif written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1t defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1t written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1u defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1u written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1v defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1v written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1f defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1f written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2t defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2t written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2u defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2u written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2v defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2v written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2f defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2f written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: ff defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: ff written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: mbathy defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: mbathy written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: misf defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: misf written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: isfdraft defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: isfdraft written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_1d defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_1d written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_1d defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_1d written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_1d defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_1d written ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_1d defined ok iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_1d written ok iom_close ~~~ close file: mesh_mask_0000.nc ok dom_ctl : extrema of the masked scale factors ~~~~~~~ e1t maxi: 106000.00 at i = 2 j= 2 e1t mini: 106000.00 at i = 2 j= 2 e2t maxi: 106000.00 at i = 2 j= 2 e2t mini: 106000.00 at i = 2 j= 2 istate_ini : Initialization of the dynamics and tracers ~~~~~~~~~~ dta_tsd_init : Temperature & Salinity data ~~~~~~~~~~~~ Namelist namtsd Initialisation of ocean T & S with T &S input data ln_tsd_init = F damping of ocean T & S toward T &S input data ln_tsd_tradmp = F T & S data not used *** Info used values : date ndastp : 10100 number of elapsed days since the begining of run : 0.000000000000000E+000 =======>> 1/2 time step before the start of the run DATE Y/M/D = 0/12/30 nsec_day: 82800 nsec_week: -522000 ======>> time-step = 1 New day, DATE Y/M/D = 0001/01/01 nday_year = 001 nsec_year = 3600 nsec_month = 3600 nsec_day = 3600 nsec_week = -514800 istate_gyre : initial analytical T and S profil deduced from LEVITUS ~~~~~~~~~~~ Initial temperature and salinity profiles: level gdept_1d temperature salinity 1 4.98 23.46 36.84 2 15.10 23.24 36.83 3 25.54 23.01 36.82 4 36.45 22.75 36.81 5 48.07 22.46 36.81 6 60.73 22.13 36.80 7 74.94 21.74 36.78 8 91.42 21.28 36.77 9 111.22 20.71 36.74 10 135.84 20.03 36.70 11 167.35 19.23 36.63 12 208.50 18.35 36.54 13 262.83 17.40 36.41 14 334.49 16.34 36.24 15 427.90 15.08 36.02 16 547.03 13.42 35.74 17 694.55 11.23 35.48 18 871.10 8.96 35.30 19 1075.10 7.04 35.20 20 1303.24 5.69 35.15 21 1551.23 4.86 35.13 22 1814.70 4.41 35.12 23 2089.75 4.19 35.12 24 2373.16 4.09 35.12 25 2662.49 4.04 35.12 26 2955.93 4.02 35.12 27 3252.20 4.01 35.12 28 3550.40 4.00 35.12 29 3849.90 4.00 35.12 30 4150.28 4.00 35.12 31 4451.26 0.00 0.00 sbc_init : surface boundary condition setting ~~~~~~~~ Namelist namsbc (partly overwritten with CPP key setting) frequency update of sbc (and ice) nn_fsbc = 1 Type of sbc : analytical formulation ln_ana = T flux formulation ln_flx = F CLIO bulk formulation ln_blk_clio = F CORE bulk formulation ln_blk_core = F MFS bulk formulation ln_blk_mfs = F ocean-atmosphere coupled formulation ln_cpl = F forced-coupled mixed formulation ln_mixcpl = F OASIS coupling (with atm or sas) lk_oasis = F components of your executable nn_components = 0 Multicategory heat flux formulation (LIM3) nn_limflx = -1 Misc. options of sbc : Patm gradient added in ocean & ice Eqs. ln_apr_dyn = F ice management in the sbc (=0/1/2/3) nn_ice = 0 ice-ocean embedded/levitating (=0/1/2) nn_ice_embd = 1 daily mean to diurnal cycle qsr ln_dm2dc = F runoff / runoff mouths ln_rnf = F iceshelf formulation nn_isf = 0 Sea Surface Restoring on SST and/or SSS ln_ssr = F FreshWater Budget control (=0/1/2) nn_fwb = 0 closed sea (=0/1) (set in namdom) nn_closea = 0 n. of iterations if land-sea-mask applied nn_lsm = 0 Use of per-category fluxes (nn_limflx = -1) GYRE analytical formulation sbc_ssm : sea surface mean fields, nn_fsbc=1 : instantaneous values ~~~~~~~ default initialisation of ss?_m arrays dyn_nept_init : Simplified Neptune module ~~~~~~~~~~~~~ --> Reading namelist namdyn_nept parameters: ln_neptsimp = F zdf_init: vertical physics ~~~~~~~~ Namelist namzdf : set vertical mixing mixing parameters vertical eddy viscosity rn_avm0 = 1.200000000000000E-004 vertical eddy diffusivity rn_avt0 = 1.200000000000000E-005 constant background or profile nn_avb = 0 horizontal variation for avtb nn_havtb = 0 time splitting / backward scheme ln_zdfexp = F number of time step nn_zdfexp = 3 enhanced vertical diffusion ln_zdfevd = T applied on momentum (=1/0) nn_evdm = 1 vertical coefficient for evd rn_avevd = 100.000000000000 non-penetrative convection (npc) ln_zdfnpc = F npc call frequency nn_npc = 1 npc print frequency nn_npcp = 365 vertical mixing option : TKE dependent eddy coefficients convection : use enhanced vertical dif. scheme use the 1.5 turbulent closure zdf_bfr_init : momentum bottom friction ~~~~~~~~~~~~~ Namelist nam_bfr : set bottom friction parameters quadratic bottom friction friction coef. rn_bfri2 = 1.000000000000000E-003 Max. coef. (log case) rn_bfri2_max = 0.100000000000000 background tke rn_bfeb2 = 2.500000000000000E-003 log formulation ln_bfr2d = F bottom roughness rn_bfrz0 [m] = 3.000000000000000E-003 implicit bottom friction switch ln_bfrimp = T zdf_tke_init : tke turbulent closure scheme - initialisation ~~~~~~~~~~~~ Namelist namzdf_tke : set tke mixing parameters coef. to compute avt rn_ediff = 0.100000000000000 Kolmogoroff dissipation coef. rn_ediss = 0.700000000000000 tke surface input coef. rn_ebb = 67.8300000000000 minimum value of tke rn_emin = 1.000000000000000E-006 surface minimum value of tke rn_emin0 = 1.000000000000000E-004 background shear (>0) rn_bshear = 9.999999999999999E-021 mixing length type nn_mxl = 2 prandl number flag nn_pdl = 1 surface mixing length = F(stress) or not ln_mxl0 = T surface mixing length minimum value rn_mxl0 = 4.000000000000000E-002 flag to take into acc. Langmuir circ. ln_lc = T coef to compute verticla velocity of LC rn_lc = 0.150000000000000 test param. to add tke induced by wind nn_etau = 0 flag for computation of exp. tke profile nn_htau = 1 fraction of en which pene. the thermocline rn_efr = 5.000000000000000E-002 critical Richardson nb with your parameters ri_cri = 0.222222222222222 use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min ldf_tra_init : lateral tracer physics ~~~~~~~~~~~~ Namelist namtra_ldf : lateral mixing parameters (type, direction, coefficien ts) laplacian operator ln_traldf_lap = T bilaplacian operator ln_traldf_bilap = F iso-level ln_traldf_level = F horizontal (geopotential) ln_traldf_hor = F iso-neutral ln_traldf_iso = T iso-neutral (Griffies) ln_traldf_grif = F Griffies strmfn diagnostics ln_traldf_gdia = F lateral eddy diffusivity rn_aht_0 = 1000.00000000000 background hor. diffusivity rn_ahtb_0 = 0.000000000000000E+000 eddy induced velocity coef. rn_aeiv_0 = 0.000000000000000E+000 maximum isoppycnal slope rn_slpmax = 1.000000000000000E-002 pure lateral mixing in ML ln_triad_iso = F lateral mixing on bottom ln_botmix_grif = F tracer mixing coef. = constant (default option) harmonic tracer diffusion (default) constant eddy diffusivity coef. ahtu = ahtv = ahtw = aht0 = 1000.00000000000 ldf_dyn : lateral momentum physics ~~~~~~~ Namelist namdyn_ldf : set lateral mixing parameters laplacian operator ln_dynldf_lap = T bilaplacian operator ln_dynldf_bilap = F iso-level ln_dynldf_level = F horizontal (geopotential) ln_dynldf_hor = T iso-neutral ln_dynldf_iso = F horizontal laplacian eddy viscosity rn_ahm_0_lap = 100000.000000000 background viscosity rn_ahmb_0 = 0.000000000000000E+000 horizontal bilaplacian eddy viscosity rn_ahm_0_blp = 0.000000000000000E+000 upper limit for laplacian eddy visc rn_ahm_m_lap = 40000.0000000000 upper limit for bilap eddy viscosity rn_ahm_m_blp = -1000000000000.00 momentum mixing coef. = constant (default option) harmonic momentum diff. (default) inildf: constant eddy viscosity coef. ~~~~~~ ahm1 = ahm2 = ahm0 = 100000.000000000 ldf_slp_init : direction of lateral mixing ~~~~~~~~~~~~ tra_qsr_init : penetration of the surface solar radiation ~~~~~~~~~~~~ Namelist namtra_qsr : set the parameter of penetration Light penetration (T) or not (F) ln_traqsr = T RGB (Red-Green-Blue) light penetration ln_qsr_rgb = F 2 band light penetration ln_qsr_2bd = T bio-model light penetration ln_qsr_bio = F light penetration for ice-model LIM3 ln_qsr_ice = T RGB : Chl data (=1/2) or cst value (=0) nn_chldta = 0 RGB & 2 bands: fraction of light (rn_si1) rn_abs = 0.580000000000000 RGB & 2 bands: shortess depth of extinction rn_si0 = 0.350000000000000 2 bands: longest depth of extinction rn_si1 = 23.0000000000000 2 bands light penetration level of light extinction = 18 ref depth = 969.843774414690 m tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating ~~~~~~~ Namelist nambbc : set bbc parameters Apply a geothermal heating at ocean bottom ln_trabbc = F type of geothermal flux nn_geoflx = 0 Constant geothermal flux value rn_geoflx_cst = 8.640000000000000E-002 *** no geothermal heat flux tra_dmp_init : T and S newtonian relaxation ~~~~~~~ Namelist namtra_dmp : set relaxation parameters Apply relaxation or not ln_tradmp = F mixed layer damping option nn_zdmp = 0 Damping file name cn_resto = resto.nc tra_adv_init : choice/control of the tracer advection scheme ~~~~~~~~~~~ Namelist namtra_adv : chose a advection scheme for tracers 2nd order advection scheme ln_traadv_cen2 = F TVD advection scheme ln_traadv_tvd = T MUSCL advection scheme ln_traadv_muscl = F MUSCL2 advection scheme ln_traadv_muscl2 = F UBS advection scheme ln_traadv_ubs = F QUICKEST advection scheme ln_traadv_qck = F upstream scheme within muscl ln_traadv_msc_ups = F TVD advection scheme with zts ln_traadv_tvd_zts = F TVD scheme is used tra_adv_mle_init : mixed layer eddy (MLE) advection acting on tracers ~~~~~~~~~~~~~~~~ Namelist namtra_adv_mle : mixed layer eddy advection on tracers use mixed layer eddy (MLE, i.e. Fox-Kemper param) (T/F) ln_mle = T MLE type: =0 standard Fox-Kemper ; =1 new formulation nn_mle = 1 magnitude of the MLE (typical value: 0.06 to 0.08) rn_ce = 6.000000000000000E-002 scale of ML front (ML radius of deformation) (rn_mle=0) rn_lf = 5000.00000000000 m maximum time scale of MLE (rn_mle=0) rn_time = 172800.000000000 s reference latitude (degrees) of MLE coef. (rn_mle=1) rn_lat = 20.0000000000000 deg space interp. of MLD at u-(v-)pts (0=min,1=averaged,2=max) nn_mld_uv = 0 =1 no MLE in case of convection ; =0 always MLE nn_conv = 0 Density difference used to define ML for FK rn_rho_c_mle = 1.000000000000000E-002 Mixed Layer Eddy induced transport added to tracer advection New formulation ML buoyancy criteria = 9.558138401559454E-005 m/s2 associated ML density criteria defined in zdfmxl = 1.000000000000000E-002 kg/m3 tra_ldf_init : lateral tracer diffusive operator ~~~~~~~~~~~ Namelist namtra_ldf already read in ldftra module see ldf_tra_init report for lateral mixing parameters Rotated laplacian operator tra:ldf_ano : lateral diffusion acting on the full fields ~~~~~~~~~~~ tra_zdf_init : vertical tracer physics scheme ~~~~~~~~~~~ Implicit (euler backward) scheme dyn_adv_init : choice/control of the momentum advection scheme ~~~~~~~~~~~ Namelist namdyn_adv : chose a advection formulation & scheme for momentu m Vector/flux form (T/F) ln_dynadv_vec = T = 0 standard scheme ; =1 Hollingsworth scheme nn_dynkeg = 0 2nd order centred advection scheme ln_dynadv_cen2 = F 3rd order UBS advection scheme ln_dynadv_ubs = F Sub timestepping of vertical advection ln_dynzad_zts = F vector form : keg + zad + vor is used with Centered standard keg scheme dyn_vor_init : vorticity term : read namelist and control the consistency ~~~~~~~~~~~~ Namelist namdyn_vor : choice of the vorticity term scheme energy conserving scheme ln_dynvor_ene = T enstrophy conserving scheme ln_dynvor_ens = F mixed enstrophy/energy conserving scheme ln_dynvor_mix = F enstrophy and energy conserving scheme ln_dynvor_een = F enstrophy and energy conserving scheme (old) ln_dynvor_een_old= F Vector form advection : vorticity = Coriolis + relative vorticity vorticity scheme : energy conserving scheme dyn_ldf_init : Choice of the lateral diffusive operator on dynamics ~~~~~~~~~~~ Namelist nam_dynldf : set lateral mixing parameters (type, direction, co efficients) laplacian operator ln_dynldf_lap = T bilaplacian operator ln_dynldf_bilap = F iso-level ln_dynldf_level = F horizontal (geopotential) ln_dynldf_hor = T iso-neutral ln_dynldf_iso = F laplacian operator dyn_hpg_init : hydrostatic pressure gradient initialisation ~~~~~~~~~~~~ Namelist namdyn_hpg : choice of hpg scheme z-coord. - full steps ln_hpg_zco = T z-coord. - partial steps (interpolation) ln_hpg_zps = F s-coord. (standard jacobian formulation) ln_hpg_sco = F s-coord. (standard jacobian formulation) for isf ln_hpg_isf = F s-coord. (Density Jacobian: Cubic polynomial) ln_hpg_djc = F s-coord. (Pressure Jacobian: Cubic polynomial) ln_hpg_prj = F time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = F dyn_zdf_init : vertical dynamics physics scheme ~~~~~~~~~~~ Implicit (euler backward) scheme dyn_spg_init : choice of the surface pressure gradient scheme ~~~~~~~~~~~ Explicit free surface lk_dynspg_exp = F Free surface with time splitting lk_dynspg_ts = F Filtered free surface cst volume lk_dynspg_flt = T filtered free surface file : solver.stat open ok unit = 20 status = REPLACE form = FORMATTED access = SEQUENTIAL solver_init : solver to compute the surface pressure gradient ~~~~~~~~~~~ Namelist namsol : set solver parameters type of elliptic solver nn_solv = 2 absolute/relative (0/1) precision nn_sol_arp = 0 minimum iterations for solver nn_nmin = 210 maximum iterations for solver nn_nmax = 800 frequency for test nn_nmod = 10 absolute precision of solver rn_eps = 1.000000000000000E-006 absolute precision for SOR solver rn_resmax = 1.000000000000000E-010 optimal coefficient of sor rn_sor = 1.96000000000000 a successive-over-relaxation solver with extra outer halo is used with jpr2di = 0 and jpr2dj = 0 icbini : Namelist namberg ln_icebergs = F , NO icebergs used ~~~~~~~~ ===>>> : W A R N I N G =============== W A R N I N G: end of record or file while reading namelist namsto in configur ation namelist iostat = -1 sto_par_init : stochastic parameterization ~~~~~~~~~~~~ Namelist namsto : stochastic parameterization restart stochastic parameters ln_rststo = F read seed of RNG from restart file ln_rstseed = T suffix of sto restart name (input) cn_storst_in = restart_sto suffix of sto restart name (output) cn_storst_out = restart_sto stochastic equation of state ln_sto_eos = F number of degrees of freedom nn_sto_eos = 1 random walk horz. std (in grid points) rn_eos_stdxy = 1.40000000000000 random walk vert. std (in grid points) rn_eos_stdz = 0.700000000000000 random walk tcor (in timesteps) rn_eos_tcor = 1440.00000000000 order of autoregressive processes nn_eos_ord = 1 passes of Laplacian filter nn_eos_flt = 0 limitation factor rn_eos_lim = 2.00000000000000 stochastic parameterization : dia_ptr_init : poleward transport and msf initialization ~~~~~~~~~~~~ Namelist namptr : set ptr parameters Poleward heat & salt transport (T) or not (F) ln_diaptr = F Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = F dia_hsb_init ~~~~~~~~ check the heat and salt budgets (T) or not (F) ln_diahsb = F trd_init : Momentum/Tracers trends ~~~~~~~~~~ Namelist namtrd : set trends parameters global domain averaged dyn & tra trends ln_glo_trd = F U & V trends: 3D output ln_dyn_trd = F U & V trends: Mixed Layer averaged ln_dyn_mxl = F T & S trends: 3D output ln_tra_trd = F T & S trends: Mixed Layer averaged ln_tra_mxl = F Kinetic Energy trends ln_KE_trd = F Potential Energy trends ln_PE_trd = F Barotropic vorticity trends ln_vor_trd = F frequency of trends diagnostics (glo) nn_trd = 365 Euler time step switch is 0 AAAAAAAA