CNRS - NERC - Met OFFICE - MERCATOR-ocean - INGV - CMCC NEMO team Ocean General Circulation Model version 3.4 (2011) 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 = 1 namcfg : configuration initialization through namelist read ~~~~~~~ Namelist namcfg configuration name cp_cfg = orca configuration zoom name cp_cfz = configuration resolution jp_cfg = 1 1st lateral dimension ( >= jpi ) jpidta = 362 2nd " " ( >= jpj ) jpjdta = 292 3nd " " jpkdta = 46 1st dimension of global domain in i jpiglo = 362 2nd - - in j jpjglo = 292 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 = 6 mpp_init : Message Passing MPI ~~~~~~~~ defines mpp subdomains ---------------------- iresti= 3 irestj= 2 jpni = 21 jpnj = 2 **************************************************************** * * * * * 2 * 20 x147 * 20 x147 * 20 x147 * 19 x147 * * * * * * **************************************************************** * * * * * 1 * 20 x147 * 20 x147 * 20 x147 * 19 x147 * * * * * * **************************************************************** 1 2 3 4 **************************************************************** * * * * * 2 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** * * * * * 1 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** 5 6 7 8 **************************************************************** * * * * * 2 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** * * * * * 1 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** 9 10 11 12 **************************************************************** * * * * * 2 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** * * * * * 1 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** 13 14 15 16 **************************************************************** * * * * * 2 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** * * * * * 1 * 19 x147 * 19 x147 * 19 x147 * 19 x147 * * * * * * **************************************************************** 17 18 19 20 **************** * * 2 * 19 x147 * * * **************** * * 1 * 19 x147 * * * **************** 21 sum ilcit(i,1) = 362.000000000000 jpiglo = 362 sum ilcit(1,j) = 292.000000000000 jpjglo = 292 nproc = 0 nowe = 20 noea = 1 nono = 21 noso = -21 nbondi = 0 nbondj = -1 npolj = 0 nperio = 0 nlci = 20 nlcj = 147 nimpp = 1 njmpp = 1 nbse = 0 npse = -20 nbsw = 0 npsw = -1 nbne = 1 npne = 22 nbnw = 1 npnw = 41 mpp_init_ioipsl : iloc = 20 147 ~~~~~~~~~~~~~~~ iabsf = 1 1 ihals = 0 0 ihale = 1 1 file : timing.output open ok unit = 20 status = REPLACE form = FORMATTED access = SEQUENTIAL phy_cst : initialization of ocean parameters and constants ~~~~~~~ Domain info dimension of model Local domain Global domain Data domain jpi : 20 jpiglo : 362 jpidta : 362 jpj : 147 jpjglo : 292 jpjdta : 292 jpk : 46 jpk : 46 jpkdta : 46 jpij : 2940 mpp local domain info (mpp) jpni : 21 jpreci : 1 jpnj : 2 jprecj : 1 jpnij : 42 lateral domain boundary condition type : jperio = 6 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.160000000000 K melting point of ice rt0_ice = 273.160000000000 K volumic mass of pure water rauw = 1000.00000000000 kg/m^3 volumic mass of reference rau0 = 1035.00000000000 kg/m^3 1. / rau0 r1_rau0 = 9.661835748792270E-004 m^3/kg ocean specific heat rcp = 4000.00000000000 J/Kelvin 1. / ( rau0 * rcp ) r1_rau0_rcp = 2.415458937198068E-007 thermal conductivity of the snow = 0.310000000000000 J/s/m/K thermal conductivity of the ice = 2.03439600000000 J/s/m/K fresh ice specific heat = 2067.00000000000 J/kg/K latent heat of fusion of fresh ice / snow = 334000.000000000 J/kg latent heat of subl. of fresh ice / snow = 2834000.00000000 J/kg volumetric latent heat fusion of snow = 110220000.000000 J/m^3 density of sea ice = 917.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 thermal exp. coef. (linear) rn_alpha = 2.000000000000000E-004 saline exp. coef. (linear) rn_beta = 7.700000000000000E-004 use of Jackett & McDougall (1994) equation of state and McDougall (1987) Brunt-Vaisala frequency dom_cfg : set the ocean configuration ~~~~~~~ ocean model configuration used : cp_cfg = orca jp_cfg = 1 global domain lateral boundaries jperio= 6, cyclic east-west and north fold with F-point pivot dom_glo : domain: data / local ~~~~~~~ data input domain : jpidta = 362 jpjdta = 292 jpkdta = 46 global or zoom domain: jpiglo = 362 jpjglo = 292 jpk = 46 local domain : jpi = 20 jpj = 147 jpk = 46 south-west indices jpizoom = 1 jpjzoom = 1 conversion local ==> data i-index domain 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 conversion data ==> local i-index domain starting index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 ending index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 conversion local ==> data j-index domain 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 conversion data ==> local j-index domain starting index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 ending index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 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 = C36B file prefix restart input cn_ocerst_in= restart_oce file prefix restart output cn_ocerst_out= restart_oce restart logical ln_rstart = F start with forward time step nn_euler = 1 control of time step nn_rstctl = 1 number of the first time step nn_it000 = 1 number of the last time step nn_itend = 744 initial calendar date aammjj nn_date0 = 19900101 leap year calendar (0/1) nn_leapy = 1 initial state output nn_istate = 0 frequency of restart file nn_stock = 0 frequency of output file nn_write = -1 multi file dimgout ln_dimgnnn = F mask land points ln_mskland = F overwrite an existing file ln_clobber = F NetCDF chunksize (bytes) nn_chunksz = 0 ===>>> : W A R N I N G =============== ln_rstart =.FALSE., nn_euler is forced to 0 ===>>> : W A R N I N G =============== nstock = 0 it is forced to 744 The IOIPSL calendar is "gregorian", i.e. leap year Namelist namdom : space & time domain flag read/compute bathymetry nn_bathy = 1 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 = 0 = 0 no file created = 1 mesh_mask = 2 mesh and mask = 3 mesh_hgr, msh_zgr and mask ocean time step rn_rdt = 3600.00000000000 asselin time filter parameter rn_atfp = 0.100000000000000 acceleration of converge nn_acc = 0 nn_acc=1: surface tracer rdt rn_rdtmin = 28800.0000000000 bottom tracer rdt rdtmax = 28800.0000000000 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 = 0 longitude of first raw and column T-point ppglam0 = 999999.000000000 latitude of first raw and column T-point ppgphi0 = 999999.000000000 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 = -4762.96143546300 ppa0 = 255.580490704400 ppa1 = 245.581322324900 ppkth = 21.4333619793800 ppacr = 3.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 = 0 position of the first row and ppglam0 = 999999.000000000 column grid-point (degrees) ppgphi0 = 999999.000000000 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 curvilinear coordinate on the sphere read in "coordinate" file hgr_read : read the horizontal coordinates ~~~~~~~~ jpiglo = 362 jpjglo = 292 jpk = 46 iom_nf90_open ~~~ open existing file: coordinates.nc in REA D mode ---> coordinates.nc OK read glamt (rec: 1) in coordinates.nc ok read glamu (rec: 1) in coordinates.nc ok read glamv (rec: 1) in coordinates.nc ok read glamf (rec: 1) in coordinates.nc ok read gphit (rec: 1) in coordinates.nc ok read gphiu (rec: 1) in coordinates.nc ok read gphiv (rec: 1) in coordinates.nc ok read gphif (rec: 1) in coordinates.nc ok read e1t (rec: 1) in coordinates.nc ok read e1u (rec: 1) in coordinates.nc ok read e1v (rec: 1) in coordinates.nc ok read e1f (rec: 1) in coordinates.nc ok read e2t (rec: 1) in coordinates.nc ok read e2u (rec: 1) in coordinates.nc ok read e2v (rec: 1) in coordinates.nc ok read e2f (rec: 1) in coordinates.nc ok iom_close ~~~ close file: coordinates.nc ok orca_r1: Gibraltar : e2u reduced to 20 km orca_r1: Bhosporus : e2u reduced to 10 km orca_r1: Lombok : e1v reduced to 10 km orca_r1: Sumba : e1v reduced to 8 km orca_r1: Ombai : e1v reduced to 13 km orca_r1: Timor Passage : e1v reduced to 20 km orca_r1: W Halmahera : e1v reduced to 30 km orca_r1: E Halmahera : e1v reduced to 50 km longitude and e1 scale factors ------------------------------ 1 72.50 73.00 72.50 73.00 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808 11 82.50 83.00 82.50 83.00 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808 latitude and e2 scale factors ----------------------------- 1 -78.39 -78.39 -78.29 -78.29 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808 11 -76.20 -76.20 -76.08 -76.08 26524.8299354813 26524.8299354813 26750.5179637307 26750.5179637307 21 -73.60 -73.60 -73.46 -73.46 31393.5671869958 31393.5671869958 31657.3897510458 31657.3897510458 31 -70.53 -70.53 -70.36 -70.36 37065.0717616678 37065.0717616678 37371.1271782116 37371.1271782116 41 -66.91 -66.91 -66.71 -66.71 43612.0739372158 43612.0739372158 43963.3187131693 43963.3187131693 51 -62.66 -62.66 -62.43 -62.43 51073.7561269995 51073.7561269995 51470.7880537348 51470.7880537348 61 -57.70 -57.70 -57.43 -57.43 59426.2991489708 59426.2991489708 59865.5913742379 59865.5913742379 71 -51.95 -51.95 -51.64 -51.64 68542.4202550136 68542.4202550136 69014.0413289406 69014.0413289406 81 -45.35 -45.35 -45.00 -45.00 78143.3450408712 78143.3450408712 78628.5295414200 78628.5295414200 91 -37.89 -37.89 -37.50 -37.50 87755.5230030943 87755.5230030943 88225.0221273672 88225.0221273672 101 -29.59 -29.59 -29.15 -29.15 96696.8748118213 96696.8748118213 97111.6439805149 97111.6439805149 111 -20.55 -20.55 -20.08 -20.08 104126.3062508076 104126.3062508076 104442.1936993819 104442.1936993819 121 -11.39 -11.39 -10.98 -10.98 90176.3801300774 90176.3801300774 87910.2293468842 87910.2293468842 131 -5.54 -5.54 -5.34 -5.34 44969.8813203525 44969.8813203525 43942.4027802381 43942.4027802381 141 -2.00 -2.00 -1.83 -1.83 37154.1706497673 37154.1706497673 37119.2489256383 37119.2489256383 dom_zgr : vertical coordinate ~~~~~~~ Namelist namzgr : set vertical coordinate z-coordinate - full steps ln_zco = F z-coordinate - partial steps ln_zps = T s- or hybrid z-s-coordinate ln_sco = F zgr_z : Reference vertical z-coordinates ~~~~~~~ Value of coefficients for vertical mesh: zsur = -4762.96143546300 za0 = 255.580490704400 za1 = 245.581322324900 zkth = 21.4333619793800 zacr = 3.00000000000000 Reference z-coordinate depth and scale factors: level gdept_1d gdepw_1d e3t_1d e3w_1d 1 5.00 0.00 10.00 10.00 2 15.00 10.00 10.00 10.00 3 25.00 20.00 10.00 10.00 4 35.01 30.00 10.01 10.00 5 45.01 40.01 10.01 10.01 6 55.03 50.02 10.02 10.02 7 65.06 60.04 10.04 10.03 8 75.13 70.09 10.09 10.06 9 85.25 80.18 10.17 10.12 10 95.49 90.35 10.33 10.24 11 105.97 100.69 10.65 10.47 12 116.90 111.36 11.27 10.91 13 128.70 122.65 12.47 11.77 14 142.20 135.16 14.78 13.43 15 158.96 150.03 19.23 16.65 16 181.96 169.42 27.66 22.78 17 216.65 197.37 43.26 34.30 18 272.48 241.13 70.88 55.21 19 364.30 312.74 116.11 90.99 20 511.53 429.72 181.55 146.43 21 732.20 611.89 261.03 220.35 22 1033.22 872.87 339.39 301.42 23 1405.70 1211.59 402.26 373.31 24 1830.89 1612.98 444.87 426.00 25 2289.77 2057.13 470.55 459.47 26 2768.24 2527.22 484.95 478.83 27 3257.48 3011.90 492.70 489.44 28 3752.44 3504.46 496.78 495.07 29 4250.40 4001.16 498.90 498.02 30 4749.91 4500.02 500.00 499.54 31 5250.23 5000.00 500.56 500.33 32 5750.95 5500.55 500.86 500.73 33 6251.89 6001.40 501.00 500.94 34 6752.94 6502.40 501.08 501.05 35 7254.04 7003.48 501.12 501.10 36 7755.17 7504.60 501.14 501.13 37 8256.32 8005.74 501.15 501.15 38 8757.47 8506.89 501.16 501.15 39 9258.63 9008.05 501.16 501.16 40 9759.79 9509.21 501.16 501.16 41 10260.95 10010.37 501.16 501.16 42 10762.11 10511.53 501.16 501.16 43 11263.27 11012.69 501.16 501.16 44 11764.43 11513.85 501.16 501.16 45 12265.60 12015.01 501.16 501.16 46 12766.76 12516.18 501.16 501.16 zgr_bat : defines level and meter bathymetry ~~~~~~~ iom_nf90_open ~~~ open existing file: bathy_meter.nc in REA D mode ---> bathy_meter.nc OK read Bathymetry (rec: 1) in bathy_meter.nc ok iom_close ~~~ close file: bathy_meter.nc ok Minimum ocean depth: 30.0032158075635 minimum number of ocean levels : 3 zgr_zps : z-coordinate with partial steps ~~~~~~~ mbathy is recomputed : bathy_level file is NOT used zgr_bat_ctl : check the bathymetry ~~~~~~~~~~~ suppress isolated ocean grid points ----------------------------------- the number of ocean level at grid-point (i,j) = 10 41 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 11 42 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 19 43 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 5 52 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 8 68 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 13 99 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 10 100 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 16 122 is changed from 31 to 30 672 ocean grid points suppressed mbathy set to 0 along east and west boundary: nperio = 0 maximum number of ocean level = 31 < jpk-1 you can decrease jpk to 32 zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ~~~~~~~~~~~~~ dommsk : ocean mask ~~~~~~ Namelist namlbc lateral momentum boundary cond. rn_shlat = 2.00000000000000 consistency with analytical form ln_vorlat = F ocean lateral no-slip orca_r1: increase friction near the following straits : Gibraltar Bhosporus Makassar (Top) Lombok Ombai Timor Passage West Halmahera East Halmahera dom_stp : time stepping setting ~~~~~~~ synchronous time stepping dynamics and tracer time step = 1.00000000000000 hours dom_ctl : extrema of the masked scale factors ~~~~~~~ e1t maxi: 111198.92 at i = 1 j= 147 e1t mini: 25207.90 at i = 92 j= 8 e2t maxi: 105405.52 at i = 1 j= 114 e2t mini: 24988.38 at i = 208 j= 283 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 = T damping of ocean T & S toward T &S input data ln_tsd_tradmp = F dta_tsd : Temperature & Salinity data ~~~~~~~ namtsd Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./temperature variable name: votemper frequency: -1.00000000000000 time interp: F climatology: T weights : pairing : data type: monthly land/sea mask: root filename: ./salinity variable name: vosaline frequency: -1.00000000000000 time interp: F climatology: T weights : pairing : data type: monthly land/sea mask: *** Info used values : date ndastp : 19900100 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 = 1989/12/31 nsec_day: 84600 nsec_week: -1800 ======>> time-step = 1 New day, DATE Y/M/D = 1990/01/01 nday_year = 001 nsec_year = 1800 nsec_month = 1800 nsec_day = 1800 nsec_week = 1800 iom_nf90_open ~~~ open existing file: ./temperature_m01.nc in READ mode ---> ./temperature_m01.nc OK read votemper (rec: 1) in ./temperature_m01.nc ok iom_nf90_open ~~~ open existing file: ./salinity_m01.nc in READ mode ---> ./salinity_m01.nc OK read vosaline (rec: 1) in ./salinity_m01.nc ok fld_read: var votemper kt = 1 ( 0.0208 days), Y/M/D = 1990/01/01, record: 0001 (days 0.0000 <-> 31.0000) fld_read: var vosaline kt = 1 ( 0.0208 days), Y/M/D = 1990/01/01, record: 0001 (days 0.0000 <-> 31.0000) temperature Levitus level = 1 1 147 28.48 127 28.16 107 25.88 87 14.67 67 2.58 47 0.83 27 0.00 7 0.00 level = 23 1 147 7.51 127 6.54 107 7.44 87 6.20 67 2.12 47 1.03 27 0.00 7 0.00 level = 45 1 147 0.00 127 0.00 107 0.00 87 0.00 67 0.00 47 0.00 27 0.00 7 0.00 salinity Levitus level = 1 1 147 35.15 127 34.36 107 35.42 87 34.89 67 33.88 47 33.67 27 0.00 7 0.00 level = 23 1 147 34.97 127 34.78 107 34.54 87 34.45 67 34.67 47 34.71 27 0.00 7 0.00 level = 45 1 147 0.00 127 0.00 107 0.00 87 0.00 67 0.00 47 0.00 27 0.00 7 0.00 dta_tsd: deallocte T & S arrays as they are only use to initialize the run dyn_nept_init : Simplified Neptune module ~~~~~~~~~~~~~ --> Reading namelist namdyn_nept parameters: ln_neptsimp = F 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 = F 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 coupled formulation (T if key_sbc_cpl) ln_cpl = T Flux handling over ice categories cn_iceflx = linear 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 = 3 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 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 Fluxes averaged over all ice categories ln_iceflx_ave = F Fluxes distributed linearly over ice categories ln_iceflx_linear = T coupled formulation sbc_ssm : sea surface mean fields, nn_fsbc=1 : instantaneous values ~~~~~~~ sbc_cpl_init : namsbc_cpl namelist ~~~~~~~~~~~~ received fields (mutiple ice categogies) 10m wind module = none (no) stress module = none (no) surface stress = oce and ice (no) - referential = spherical - orientation = eastward-northward - mesh = T non-solar heat flux sensitivity = coupled (no) solar heat flux = oce and ice (yes) non-solar heat flux = oce and ice (yes) freshwater budget = conservative (no) runoffs = coupled (no) calving = coupled (no) sea ice heat fluxes = none (no) atm co2 = none (no) sent fields (multiple ice categories) surface temperature = oce and ice (yes) albedo = ice (no) ice/snow thickness = ice and snow (yes) surface current = none (no) - referential = spherical - orientation = eastward-northward - mesh = T oce co2 flux = none (no) cpl_prism_define : initialization in coupled ocean/atmosphere case ~~~~~~~~~~~~~~~~~ Define 1 1 OIceFrc.C001 for 20 Define 1 2 OIceFrc.C002 for 20 Define 1 3 OIceFrc.C003 for 20 Define 1 4 OIceFrc.C004 for 20 Define 1 5 OIceFrc.C005 for 20 Define 2 1 O_SSTSST for 20 Define 3 1 OTepIce.C001 for 20 Define 3 2 OTepIce.C002 for 20 Define 3 3 OTepIce.C003 for 20 Define 3 4 OTepIce.C004 for 20 Define 3 5 OTepIce.C005 for 20 Define 5 1 OAlbIce for 20 Define 7 1 OIceTck.C001 for 20 Define 7 2 OIceTck.C002 for 20 Define 7 3 OIceTck.C003 for 20 Define 7 4 OIceTck.C004 for 20 Define 7 5 OIceTck.C005 for 20 Define 8 1 OSnwTck.C001 for 20 Define 8 2 OSnwTck.C002 for 20 Define 8 3 OSnwTck.C003 for 20 Define 8 4 OSnwTck.C004 for 20 Define 8 5 OSnwTck.C005 for 20 Define 1 1 O_OTaux1 for 21 Define 2 1 O_OTauy1 for 21 Define 7 1 O_ITaux1 for 21 Define 8 1 O_ITauy1 for 21 Define 13 1 O_QsrOce for 21 Define 14 1 O_QsrIce.C001 for 21 Define 14 2 O_QsrIce.C002 for 21 Define 14 3 O_QsrIce.C003 for 21 Define 14 4 O_QsrIce.C004 for 21 Define 14 5 O_QsrIce.C005 for 21 Define 16 1 O_QnsOce for 21 Define 17 1 O_QnsIce.C001 for 21 Define 17 2 O_QnsIce.C002 for 21 Define 17 3 O_QnsIce.C003 for 21 Define 17 4 O_QnsIce.C004 for 21 Define 17 5 O_QnsIce.C005 for 21 Define 19 1 OTotRain for 21 Define 20 1 OTotSnow for 21 Define 21 1 OTotEvap for 21 Define 22 1 OIceEvap for 21 Define 27 1 O_dQnsdT for 21 Define 28 1 O_Runoff for 21 Define 29 1 OCalving for 21 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 = 0 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 linear botton friction friction coef. rn_bfri1 = 4.000000000000000E-004 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 = 1 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 zdf_ddm : double diffusive mixing ~~~~~~~ Namelist namzdf_ddm : set dd mixing parameter maximum avs for dd mixing rn_avts = 1.000000000000000E-004 heat/salt buoyancy flux ratio rn_hsbfr = 1.60000000000000 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 = 1000.00000000000 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. = F( latitude, longitude) harmonic tracer diffusion (default) ldf_tra_c2d : 2D eddy diffusivity and eddy ~~~~~~~~~~~ -- induced velocity coefficients 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 = 20000.0000000000 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 = 20000.0000000000 upper limit for bilap eddy viscosity rn_ahm_m_blp = -1000000000000.00 momentum mixing coef. = F( latitude, longitude, depth) harmonic momentum diff. (default) ldf_dyn_c3d : 3D lateral eddy viscosity coefficient ~~~~~~~~~~~ laplacian operator: ahm proportional to e1 maximum grid-spacing = 111198.923448546 maximum value for ahm = 20000.0000000000 ahm profile : jk ahm depth t-level 1 1.0000 4.9999 2 0.9994 15.0003 3 0.9989 25.0018 4 0.9982 35.0054 5 0.9976 45.0133 6 0.9970 55.0295 7 0.9963 65.0618 8 0.9956 75.1255 9 0.9948 85.2504 10 0.9940 95.4943 11 0.9932 105.9699 12 0.9923 116.8962 13 0.9913 128.6979 14 0.9901 142.1952 15 0.9885 158.9606 16 0.9861 181.9628 17 0.9821 216.6479 18 0.9745 272.4767 19 0.9583 364.3030 20 0.9192 511.5348 21 0.8191 732.2009 22 0.6067 1033.2173 23 0.3759 1405.6975 24 0.2766 1830.8850 25 0.2544 2289.7679 26 0.2506 2768.2423 27 0.2501 3257.4789 28 0.2500 3752.4422 29 0.2500 4250.4012 30 0.2500 4749.9133 31 0.2500 5250.2266 32 0.2500 5750.9524 33 0.2500 6251.8903 34 0.2500 6752.9371 35 0.2500 7254.0399 36 0.2500 7755.1714 37 0.2500 8256.3176 38 0.2500 8757.4714 39 0.2500 9258.6291 40 0.2500 9759.7888 41 0.2500 10260.9496 42 0.2500 10762.1108 43 0.2500 11263.2724 44 0.2500 11764.4340 45 0.2500 12265.5958 46 0.2500 12766.7575 ahm profile : jk ahm depth t-level 1 1.0000 4.9999 2 0.9994 15.0003 3 0.9989 25.0018 4 0.9982 35.0054 5 0.9976 45.0133 6 0.9970 55.0295 7 0.9963 65.0618 8 0.9956 75.1255 9 0.9948 85.2504 10 0.9940 95.4943 11 0.9932 105.9699 12 0.9923 116.8962 13 0.9913 128.6979 14 0.9901 142.1952 15 0.9885 158.9606 16 0.9861 181.9628 17 0.9821 216.6479 18 0.9745 272.4767 19 0.9583 364.3030 20 0.9192 511.5348 21 0.8191 732.2009 22 0.6067 1033.2173 23 0.3759 1405.6975 24 0.2766 1830.8850 25 0.2544 2289.7679 26 0.2506 2768.2423 27 0.2501 3257.4789 28 0.2500 3752.4422 29 0.2500 4250.4012 30 0.2500 4749.9133 31 0.2500 5250.2266 32 0.2500 5750.9524 33 0.2500 6251.8903 34 0.2500 6752.9371 35 0.2500 7254.0399 36 0.2500 7755.1714 37 0.2500 8256.3176 38 0.2500 8757.4714 39 0.2500 9258.6291 40 0.2500 9759.7888 41 0.2500 10260.9496 42 0.2500 10762.1108 43 0.2500 11263.2724 44 0.2500 11764.4340 45 0.2500 12265.5958 46 0.2500 12766.7575 ORCA R1, R2 or R4: overwrite the previous definition of ahm ================= ldfdyn_c3d_orca : 3D eddy viscosity coefficient ~~~~~~~~~~~~~~~ orca R1, R2 or R4 configuration: reduced in the surface Eq. strip file : ahmcoef open ok unit = 21 status = OLD form = FORMATTED access = SEQUENTIAL 1D zcoef array ~~~~~~~~~~~~~~ jk zcoef 1 1.000 2 1.000 3 1.000 4 1.000 5 1.000 6 1.000 7 1.000 8 1.000 9 1.000 10 1.000 11 1.000 12 1.000 13 1.000 14 1.000 15 1.000 16 1.000 17 1.000 18 1.000 19 1.000 20 1.000 21 1.000 22 2.000 23 3.000 24 5.000 25 7.000 26 9.000 27 10.000 28 10.000 29 10.000 30 10.000 31 10.000 32 10.000 33 10.000 34 10.000 35 10.000 36 10.000 37 10.000 38 10.000 39 10.000 40 10.000 41 10.000 42 10.000 43 10.000 44 10.000 45 10.000 46 10.000 3D ahm1 array (k=1) 1 147 1.00 127 4.41 107 20.00 87 20.00 67 20.00 47 20.00 27 20.00 7 20.00 3D ahm2 array (k=1) 1 147 1.00 127 4.08 107 20.00 87 20.00 67 20.00 47 20.00 27 20.00 7 20.00 3D ahm2 array (k=jpk) 1 147 10.00 127 20.00 107 20.00 87 20.00 67 20.00 47 20.00 27 20.00 7 20.00 ldf_slp_init : direction of lateral mixing ~~~~~~~~~~~~ Horizontal mixing in s-coordinate: slope = slope of s-surfaces 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 = T 2 band light penetration ln_qsr_2bd = F bio-model light penetration ln_qsr_bio = F light penetration for ice-model LIM3 ln_qsr_ice = T RGB : Chl data (=1) 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 light penetration for ice-model LIM3 ln_qsr_ice = T R-G-B light penetration - Constant Chlorophyll trc_oce_rgb : Initialisation of the optical look-up table ~~~~~~~~~~~ RGB longest depth of extinction r_si2 = 61.8046971569839 jc = 1 Chl = 1.000000000000000E-002 irgb = 1 jc = 2 Chl = 1.100000000000000E-002 irgb = 2 jc = 3 Chl = 1.300000000000000E-002 irgb = 3 jc = 4 Chl = 1.400000000000000E-002 irgb = 4 jc = 5 Chl = 1.600000000000000E-002 irgb = 5 jc = 6 Chl = 1.800000000000000E-002 irgb = 6 jc = 7 Chl = 2.000000000000000E-002 irgb = 7 jc = 8 Chl = 2.200000000000000E-002 irgb = 8 jc = 9 Chl = 2.500000000000000E-002 irgb = 9 jc = 10 Chl = 2.800000000000000E-002 irgb = 10 jc = 11 Chl = 3.200000000000000E-002 irgb = 11 jc = 12 Chl = 3.500000000000000E-002 irgb = 12 jc = 13 Chl = 4.000000000000000E-002 irgb = 13 jc = 14 Chl = 4.500000000000000E-002 irgb = 14 jc = 15 Chl = 5.000000000000000E-002 irgb = 15 jc = 16 Chl = 5.600000000000000E-002 irgb = 16 jc = 17 Chl = 6.300000000000000E-002 irgb = 17 jc = 18 Chl = 7.099999999999999E-002 irgb = 18 jc = 19 Chl = 7.900000000000000E-002 irgb = 19 jc = 20 Chl = 8.900000000000000E-002 irgb = 20 jc = 21 Chl = 0.100000000000000 irgb = 21 jc = 22 Chl = 0.112000000000000 irgb = 22 jc = 23 Chl = 0.126000000000000 irgb = 23 jc = 24 Chl = 0.141000000000000 irgb = 24 jc = 25 Chl = 0.158000000000000 irgb = 25 jc = 26 Chl = 0.178000000000000 irgb = 26 jc = 27 Chl = 0.200000000000000 irgb = 27 jc = 28 Chl = 0.224000000000000 irgb = 28 jc = 29 Chl = 0.251000000000000 irgb = 29 jc = 30 Chl = 0.282000000000000 irgb = 30 jc = 31 Chl = 0.316000000000000 irgb = 31 jc = 32 Chl = 0.355000000000000 irgb = 32 jc = 33 Chl = 0.398000000000000 irgb = 33 jc = 34 Chl = 0.447000000000000 irgb = 34 jc = 35 Chl = 0.501000000000000 irgb = 35 jc = 36 Chl = 0.562000000000000 irgb = 36 jc = 37 Chl = 0.631000000000000 irgb = 37 jc = 38 Chl = 0.708000000000000 irgb = 38 jc = 39 Chl = 0.794000000000000 irgb = 39 jc = 40 Chl = 0.891000000000000 irgb = 40 jc = 41 Chl = 1.00000000000000 irgb = 41 jc = 42 Chl = 1.12200000000000 irgb = 42 jc = 43 Chl = 1.25900000000000 irgb = 43 jc = 44 Chl = 1.41300000000000 irgb = 44 jc = 45 Chl = 1.58500000000000 irgb = 45 jc = 46 Chl = 1.77800000000000 irgb = 46 jc = 47 Chl = 1.99500000000000 irgb = 47 jc = 48 Chl = 2.23900000000000 irgb = 48 jc = 49 Chl = 2.51200000000000 irgb = 49 jc = 50 Chl = 2.81800000000000 irgb = 50 jc = 51 Chl = 3.16200000000000 irgb = 51 jc = 52 Chl = 3.54800000000000 irgb = 52 jc = 53 Chl = 3.98100000000000 irgb = 53 jc = 54 Chl = 4.46700000000000 irgb = 54 jc = 55 Chl = 5.01200000000000 irgb = 55 jc = 56 Chl = 5.62300000000000 irgb = 56 jc = 57 Chl = 6.31000000000000 irgb = 57 jc = 58 Chl = 7.07900000000000 irgb = 58 jc = 59 Chl = 7.94300000000000 irgb = 59 jc = 60 Chl = 8.91200000000000 irgb = 60 jc = 61 Chl = 10.0000000000000 irgb = 61 level of light extinction = 25 ref depth = 2527.21691545123 m Constant Chlorophyll concentration = 0.05 fixed volume: light distribution computed one for all tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating ~~~~~~~ Namelist nambbc : set bbc parameters Apply a geothermal heating at ocean bottom ln_trabbc = T type of geothermal flux nn_geoflx = 1 Constant geothermal flux value rn_geoflx_cst = 8.640000000000000E-002 *** constant heat flux = 8.640000000000000E-002 tra_bbl_init : bottom boundary layer initialisation ~~~~~~~~~~~~ Namelist nambbl : set bbl parameters diffusive bbl (=1) or not (=0) nn_bbl_ldf = 1 advective bbl (=1/2) or not (=0) nn_bbl_adv = 0 diffusive bbl coefficient rn_ahtbbl = 1000.00000000000 m2/s advective bbl coefficient rn_gambbl = 10.0000000000000 s tra_dmp_init : T and S newtonian damping ~~~~~~~ Namelist namtra_dmp : set damping parameter add a damping term or not ln_tradmp = F T and S damping option nn_hdmp = -1 mixed layer damping option nn_zdmp = 0 (non-C1D zoom: forced to 0) surface time scale (days) rn_surf = 50.0000000000000 bottom time scale (days) rn_bot = 360.000000000000 depth of transition (meters) rn_dep = 800.000000000000 create a damping.coeff file nn_file = 0 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 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 = F 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 parametrisation NOT used 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 2nd order centred advection scheme ln_dynadv_cen2 = F 3rd order UBS advection scheme ln_dynadv_ubs = F vector form : keg + zad + vor is used 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 = F enstrophy conserving scheme ln_dynvor_ens = F mixed enstrophy/energy conserving scheme ln_dynvor_mix = F enstrophy and energy conserving scheme ln_dynvor_een = T Vector form advection : vorticity = Coriolis + relative vorticity vorticity scheme : energy and enstrophy 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 = F z-coord. - partial steps (interpolation) ln_hpg_zps = T s-coord. (standard jacobian formulation) ln_hpg_sco = 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 = 22 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 = 1 absolute/relative (0/1) precision nn_sol_arp = 0 minimum iterations for solver nn_nmin = 300 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.92000000000000 a preconditioned conjugate gradient solver is used icbini : Namelist namberg ln_icebergs = F , NO icebergs used ~~~~~~~~ dia_ptr_init : poleward transport and msf initialization ~~~~~~~~~~~~ Namelist namptr : set ptr parameters Poleward heat & salt transport (T) or not (F) ln_diaptr = F Overturning heat & salt transport ln_ptrcomp = T T & S zonal mean and meridional stream function ln_diaznl = F Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = F Frequency of computation nn_fptr = 1 Frequency of outputs nn_fwri = 15 dia_hsb_init : check the heat and salt budgets ~~~~~~~~ dia_hsb_init : check the heat and salt budgets ~~~~~~~~~~~~ Namelist namhsb : set hsb parameters Switch for hsb diagnostic (T) or not (F) ln_diahsb = F Euler time step switch is 0 AAAAAAAA rot_rep : geographic <--> stretched ~~~~~ coordinate transformation sbc_ice_lim : update ocean surface boudary condition ~~~~~~~~~~~ via Louvain la Neuve Ice Model (LIM-3) time stepping file : namelist_ice_ref open ok unit = 23 status = OLD form = FORMATTED access = SEQUENTIAL file : namelist_ice_cfg open ok unit = 24 status = OLD form = FORMATTED access = SEQUENTIAL file : output.namelist.ice open ok unit = 25 status = UNKNOWN form = FORMATTED access = SEQUENTIAL ice_run : ice share parameters for dynamics/advection/thermo of sea-ice ~~~~~~ switch for ice dynamics (1) or not (0) ln_limdyn = T maximum ice concentration = 0.999000000000000 atmospheric drag over sea ice = 1.400000000000000E-003 atmospheric drag over ocean = 1.000000000000000E-003 Several ice points in the ice or not in ocean.output = F Diagnose heat/salt budget or not ln_limdiahsb = F Output heat/salt budget or not ln_limdiaout = F lim_thd : Ice Thermodynamics ~~~~~~~ Namelist of ice parameters for ice thermodynamic computation maximum melting at the bottom hmelt = -0.150000000000000 ice thick. for lateral accretion hiccrit = 0.100000000000000 Frazil ice thickness as a function of wind or not fraz_swi = 0 Maximum proportion of frazil ice collecting at bottom maxfrazb = 0.000000000000000E+000 Thresold relative drift speed for collection of frazil vfrazb = 0.416666700000000 Squeezing coefficient for collection of frazil Cfrazb = 5.00000000000000 minimum ice thickness hiclim = 0.100000000000000 numerical carac. of the scheme for diffusion in ice thickness of the surf. layer in temp. computation hnzst = 0.100000000000000 switch for snow sublimation (=1) or not (=0) parsub = 1.00000000000000 coefficient for ice-lead partition of snowfall betas = 0.600000000000000 extinction radiation parameter in sea ice (1.0) kappa_i = 1.00000000000000 maximal n. of iter. for heat diffusion computation nconv_i_thd = 50.0000000000000 maximal err. on T for heat diffusion computation maxer_i_thd = 1.000000000000000E-004 switch for comp. of thermal conductivity in the ice thcon_i_swi = 1 check heat conservation in the ice/snow con_i = F lim_thd_sal_init : Ice parameters for salinity ~~~~~~~~~~~~~~~~ switch for salinity num_sal : 2 bulk salinity value if num_sal = 1 : 4.00000000000000 restoring salinity for GD : 5.00000000000000 restoring time for GD : 1728000.00000000 restoring salinity for flushing : 2.00000000000000 restoring time for flushing : 864000.000000000 Maximum tolerated ice salinity : 20.0000000000000 Minimum tolerated ice salinity : 0.100000000000000 1st salinity for salinity profile : 3.50000000000000 2nd salinity for salinity profile : 4.50000000000000 lim_msh : LIM-3 sea-ice model, mesh initialization ~~~~~~~ the model domain is entirely in the southern hemisphere: njeq = 147 lim_itd_ini : Initialization of ice thickness distribution ~~~~~~~~~~~~ Number of ice types jpm = 1 Number of ice categories jpl = 5 Ice type 1 Number of thickness categories 5 Thickness category boundaries 1 5 Ice type vector 1 1 1 1 1 Thickness category boundaries independently of ice type hi_max 0.000000000000000E+000 0.629671477879617 1.32762233171754 2.24678625463993 3.84485861256699 7.22256121057878 Thickness category boundaries inside ice types Type number 1 hi_max_typ : 0.000000000000000E+000 0.629671477879617 1.32762233171754 2.24678625463993 3.84485861256699 7.22256121057878 lim_istate : Ice initialization ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ lim_istate_init : ice parameters inititialisation ~~~~~~~~~~~~~~~ initialization with ice (T) or not (F) ln_limini = F threshold water temp. for initial sea-ice thres_sst = 0.000000000000000E+000 initial snow thickness in the north hts_ini_n = 0.300000000000000 initial snow thickness in the south hts_ini_s = 0.300000000000000 initial ice thickness in the north hti_ini_n = 1.00000000000000 initial ice thickness in the south hti_ini_s = 1.00000000000000 initial ice concentr. in the north ati_ini_n = 0.900000000000000 initial ice concentr. in the north ati_ini_s = 0.900000000000000 initial ice salinity in the north smi_ini_n = 6.30100000000000 initial ice salinity in the south smi_ini_s = 6.30100000000000 initial ice/snw temp in the north tmi_ini_n = 270.000000000000 initial ice/snw temp in the south tmi_ini_s = 270.000000000000 lim_sbc_init : LIM-3 sea-ice - surface boundary condition ~~~~~~~~~~~~~ albedo : set albedo parameters ~~~~~~~ Namelist namsbc_alb : albedo correction for snow and ice albedo rn_cloud = 6.000000000000000E-002 albedo of melting ice in the arctic and antarctic rn_albice = 0.530000000000000 coefficients for linear rn_alphd = 0.800000000000000 interpolation used to compute albedo rn_alphdi = 0.720000000000000 between two extremes values (Pyane, 1972) rn_alphc = 0.650000000000000 lim_dyn_init : ice parameters for ice dynamics ~~~~~~~~~~~~ tolerance parameter epsd = 9.999999999999999E-021 relaxation constant om = 0.500000000000000 drag coefficient for oceanic stress cw = 5.000000000000000E-003 turning angle for oceanic stress angvg = 0.000000000000000E+000 first bulk-rheology parameter pstar = 20000.0000000000 second bulk-rhelogy parameter c_rhg = 20.0000000000000 creep limit creepl = 1.000000000000000E-012 eccentricity of the elliptical yield curve ecc = 2.00000000000000 horizontal diffusivity coeff. for sea-ice ahi0 = 350.000000000000 number of iterations for subcycling nevp = 300 timescale for elastic waves telast = 9600.00000000000 coefficient for the solution of int. stresses alphaevp = 1.00000000000000 min ice thickness for rheology calculations hminrhg = 1.000000000000000E-003 lim_itd_me_init : ice parameters for mechanical ice redistribution ~~~~~~~~~~~~~~~ Switch choosing the ice redistribution scheme ridge_scheme_swi 0 Fraction of shear energy contributing to ridging Cs 0.500000000000000 Ratio of ridging work to PotEner change in ridging Cf 17.0000000000000 Fraction of snow volume conserved during ridging fsnowrdg 0.500000000000000 Fraction of snow volume conserved during ridging fsnowrft 0.500000000000000 Fraction of total ice coverage contributing to ridging Gstar 0.150000000000000 Equivalent to G* for an exponential part function astar 5.000000000000000E-002 Quantity playing a role in max ridged ice thickness Hstar 100.000000000000 Rafting of ice sheets or not raft_swi 1 Parmeter thickness (threshold between ridge-raft) hparmeter 0.750000000000000 Rafting hyperbolic tangent coefficient Craft 5.00000000000000 Initial porosity of ridges ridge_por 0.300000000000000 Switch for part. function (0) linear (1) exponential partfun_swi 1 Switch for including brine volume in ice strength comp. brinstren_swi 0 lim_trp : Ice transport ~~~~~~~~~~~~