Context Navigation

Back to Topic #27

Topic #27: ocean.output_v3.6

File ocean.output_v3.6, 115.7 KB (added by dandanouc, 7 years ago)
ocean.output

Line
1
2	CNRS - NERC - Met OFFICE - MERCATOR-ocean - INGV - CMCC
3	NEMO team
4	Ocean General Circulation Model
5	version 3.6 (2015)
6
7
8	mynode : mpi initialisation
9
10	~~~~~~
11
12	Namelist nammpp
13
14	mpi send type cn_mpi_send = I
15
16	size in bytes of exported buffer nn_buffer = 0
17
18	jpni, jpnj and jpnij will be calculated automatically
19
20	avoid use of mpi_allgather at the north fold ln_nnogather = T
21
22	Immediate non-blocking send (isend)
23
24
25	AAAAAAAA
26
27
28	nemo_ctl: Control prints & Benchmark
29	~~~~~~~
30	Namelist namctl
31	run control (for debugging) ln_ctl = F
32	level of print nn_print = 0
33	Start i indice for SUM control nn_ictls = 0
34	End i indice for SUM control nn_ictle = 0
35	Start j indice for SUM control nn_jctls = 0
36	End j indice for SUM control nn_jctle = 0
37	number of proc. following i nn_isplt = 1
38	number of proc. following j nn_jsplt = 1
39	benchmark parameter (0/1) nn_bench = 0
40	timing activated (0/1) nn_timing = 0
41
42	namcfg : configuration initialization through namelist read
43	~~~~~~~
44	Namelist namcfg
45	configuration name cp_cfg = orca
46	configuration zoom name cp_cfz = no zoom
47	configuration resolution jp_cfg = 1
48	1st lateral dimension ( >= jpi ) jpidta = 362
49	2nd " " ( >= jpj ) jpjdta = 292
50	3nd " " jpkdta = 75
51	1st dimension of global domain in i jpiglo = 362
52	2nd - - in j jpjglo = 292
53	left bottom i index of the zoom (in data domain) jpizoom = 1
54	left bottom j index of the zoom (in data domain) jpizoom = 1
55	lateral cond. type (between 0 and 6) jperio = 4
56	use file attribute if exists as i/p j-start ln_use_jattr = F
57
58	mpp_init : Message Passing MPI
59	~~~~~~~~
60
61	mpp_init: defines mpp subdomains
62	~~~~~~ ----------------------
63
64	iresti= 12 irestj= 2
65
66	jpni= 12 jpnj= 8
67
68	sum ilcit(i,1)= 362.000000000000 jpiglo= 362
69	sum ilcit(1,j)= 292.000000000000 jpjglo= 292
70
71
72	****************************************************************
73	* * * * *
74	8 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
75	* 84 * 85 * 86 * 87 *
76	* * * * *
77	****************************************************************
78	* * * * *
79	7 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
80	* 72 * 73 * 74 * 75 *
81	* * * * *
82	****************************************************************
83	* * * * *
84	6 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
85	* 60 * 61 * 62 * 63 *
86	* * * * *
87	****************************************************************
88	* * * * *
89	5 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
90	* 48 * 49 * 50 * 51 *
91	* * * * *
92	****************************************************************
93	* * * * *
94	4 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
95	* 36 * 37 * 38 * 39 *
96	* * * * *
97	****************************************************************
98	* * * * *
99	3 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
100	* 24 * 25 * 26 * 27 *
101	* * * * *
102	****************************************************************
103	* * * * *
104	2 * 32 x 39 * 32 x 39 * 32 x 39 * 32 x 39 *
105	* 12 * 13 * 14 * 15 *
106	* * * * *
107	****************************************************************
108	* * * * *
109	1 * 32 x 39 * 32 x 39 * 32 x 39 * 32 x 39 *
110	* 0 * 1 * 2 * 3 *
111	* * * * *
112	****************************************************************
113	1 2 3 4
114
115	****************************************************************
116	* * * * *
117	8 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
118	* 88 * 89 * 90 * 91 *
119	* * * * *
120	****************************************************************
121	* * * * *
122	7 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
123	* 76 * 77 * 78 * 79 *
124	* * * * *
125	****************************************************************
126	* * * * *
127	6 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
128	* 64 * 65 * 66 * 67 *
129	* * * * *
130	****************************************************************
131	* * * * *
132	5 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
133	* 52 * 53 * 54 * 55 *
134	* * * * *
135	****************************************************************
136	* * * * *
137	4 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
138	* 40 * 41 * 42 * 43 *
139	* * * * *
140	****************************************************************
141	* * * * *
142	3 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
143	* 28 * 29 * 30 * 31 *
144	* * * * *
145	****************************************************************
146	* * * * *
147	2 * 32 x 39 * 32 x 39 * 32 x 39 * 32 x 39 *
148	* 16 * 17 * 18 * 19 *
149	* * * * *
150	****************************************************************
151	* * * * *
152	1 * 32 x 39 * 32 x 39 * 32 x 39 * 32 x 39 *
153	* 4 * 5 * 6 * 7 *
154	* * * * *
155	****************************************************************
156	5 6 7 8
157
158	****************************************************************
159	* * * * *
160	8 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
161	* 92 * 93 * 94 * 95 *
162	* * * * *
163	****************************************************************
164	* * * * *
165	7 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
166	* 80 * 81 * 82 * 83 *
167	* * * * *
168	****************************************************************
169	* * * * *
170	6 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
171	* 68 * 69 * 70 * 71 *
172	* * * * *
173	****************************************************************
174	* * * * *
175	5 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
176	* 56 * 57 * 58 * 59 *
177	* * * * *
178	****************************************************************
179	* * * * *
180	4 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
181	* 44 * 45 * 46 * 47 *
182	* * * * *
183	****************************************************************
184	* * * * *
185	3 * 32 x 38 * 32 x 38 * 32 x 38 * 32 x 38 *
186	* 32 * 33 * 34 * 35 *
187	* * * * *
188	****************************************************************
189	* * * * *
190	2 * 32 x 39 * 32 x 39 * 32 x 39 * 32 x 39 *
191	* 20 * 21 * 22 * 23 *
192	* * * * *
193	****************************************************************
194	* * * * *
195	1 * 32 x 39 * 32 x 39 * 32 x 39 * 32 x 39 *
196	* 8 * 9 * 10 * 11 *
197	* * * * *
198	****************************************************************
199	9 10 11 12
200	nproc = 0
201	nowe = 11 noea = 1
202	nono = 12 noso = -12
203	nbondi = 0
204	nbondj = -1
205	npolj = 0
206	nperio = 0
207	nlci = 32
208	nlcj = 39
209	nimpp = 1
210	njmpp = 1
211	nreci = 2 npse = -11
212	nrecj = 2 npsw = -1
213	jpreci = 1 npne = 13
214	jprecj = 1 npnw = 23
215
216	mpp_init : North fold boundary prepared for jpni >1
217
218	mpp_init_ioipsl : iloc = 32 39
219	~~~~~~~~~~~~~~~ iabsf = 1 1
220	ihals = 0 0
221	ihale = 1 1
222
223	phy_cst : initialization of ocean parameters and constants
224	~~~~~~~
225	Domain info
226	dimension of model
227	Local domain Global domain Data domain
228	jpi : 32 jpiglo : 362 jpidta : 362
229	jpj : 39 jpjglo : 292 jpjdta : 292
230	jpk : 75 jpk : 75 jpkdta : 75
231	jpij : 1248
232	mpp local domain info (mpp)
233	jpni : 12 jpreci : 1
234	jpnj : 8 jprecj : 1
235	jpnij : 96
236	lateral domain boundary condition type : jperio = 4
237
238	Constants
239
240	mathematical constant rpi = 3.14159265358979
241
242	day rday = 86400.0000000000
243	s
244	sideral year rsiyea = 31558149.0101107
245	s
246	sideral day rsiday = 86164.0996559118
247	s
248	omega omega = 7.292115083046062E-005
249	s^-1
250
251	nb of months per year raamo = 12.0000000000000
252	months
253	nb of hours per day rjjhh = 24.0000000000000
254	hours
255	nb of minutes per hour rhhmm = 60.0000000000000
256	mn
257	nb of seconds per minute rmmss = 60.0000000000000
258	s
259
260	earth radius ra = 6371229.00000000
261	m
262	gravity grav = 9.80665000000000
263	m/s^2
264
265	triple point of temperature rtt = 273.160000000000
266	K
267	freezing point of water rt0 = 273.150000000000
268	K
269	melting point of snow rt0_snow = 273.150000000000
270	K
271	melting point of ice rt0_ice = 273.050000000000
272	K
273	reference density and heat capacity now defined in eosbn2.f90
274
275	thermal conductivity of pure ice = 2.03439600000000
276	J/s/m/K
277	fresh ice specific heat = 2093.00000000000
278	J/kg/K
279	latent heat of fusion of fresh ice / snow = 333700.000000000
280	J/kg
281	density times specific heat for snow = 690690.000000000
282	J/m^3/K
283	density times specific heat for ice = 1883700.00000000
284	J/m^3/K
285	volumetric latent heat fusion of sea ice = 300330000.000000
286	J/m
287	latent heat of sublimation of snow = 2800000.00000000
288	J/kg
289	volumetric latent heat fusion of snow = 110121000.000000
290	J/m^3
291	density of sea ice = 900.000000000000
292	kg/m^3
293	density of snow = 330.000000000000
294	kg/m^3
295	emissivity of snow or ice = 0.970000000000000
296	salinity of ice = 6.00000000000000
297	psu
298	salinity of sea = 34.7000000000000
299	psu
300	latent heat of evaporation (water) = 2500000.00000000
301	J/m^3
302	correction factor for solar radiation = 0.900000000000000
303	von Karman constant = 0.400000000000000
304	Stefan-Boltzmann constant = 5.670000000000000E-008
305	J/s/m^2/K^4
306
307	conversion: degre ==> radian rad = 1.745329251994330E-002
308
309	smallest real computer value rsmall = 1.110223024625157E-016
310
311	eos_init : equation of state
312	~~~~~~~~
313	Namelist nameos : set eos parameters
314	flag for eq. of state and N^2 nn_eos = -1
315	model uses Conservative Temperature
316	Important: model must be initialized with CT and SA fields
317
318	use of TEOS-10 equation of state (cons. temp. and abs. salinity)
319
320	volumic mass of reference rau0 = 1026.00000000000
321	kg/m^3
322	1. / rau0 r1_rau0 = 9.746588693957114E-004
323	m^3/kg
324	ocean specific heat rcp = 3991.86795711963
325	J/Kelvin
326	rau0 * rcp rau0_rcp = 4095656.52400474
327	1. / ( rau0 * rcp ) r1_rau0_rcp = 2.441610994816036E-007
328
329	dom_cfg : set the ocean configuration
330	~~~~~~~
331	ocean model configuration used : cp_cfg =
332	orca
333
334
335	jp_cfg = 1
336	global domain lateral boundaries
337	jperio= 4, cyclic east-west and north fold with T-point pivot
338
339	dom_glo : domain: data / local
340	~~~~~~~
341	data input domain : jpidta = 362 jpjdta = 292
342	jpkdta = 75
343	global or zoom domain: jpiglo = 362 jpjglo = 292
344	jpk = 75
345	local domain : jpi = 32 jpj = 39
346	jpk = 75
347
348	south-west indices jpizoom = 1 jpjzoom = 1
349
350	conversion local ==> data i-index domain
351	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
352	20 21 22 23 24 25 26 27 28 29 30 31 32
353
354	conversion data ==> local i-index domain
355	starting index
356	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
357	20 21 22 23 24 25 26 27 28 29 30 31 32 33 33 33 33 33 33
358	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
359	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
360	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
361	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
362	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
363	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
364	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
365	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
366	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
367	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
368	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
369	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
370	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
371	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
372	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
373	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
374	33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
375	33
376	ending index
377	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
378	20 21 22 23 24 25 26 27 28 29 30 31 32 32 32 32 32 32 32
379	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
380	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
381	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
382	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
383	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
384	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
385	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
386	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
387	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
388	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
389	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
390	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
391	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
392	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
393	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
394	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
395	32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
396	32
397
398	conversion local ==> data j-index domain
399	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
400	20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
401	39
402
403	conversion data ==> local j-index domain
404	starting index
405	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
406	20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
407	39 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
408	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
409	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
410	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
411	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
412	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
413	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
414	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
415	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
416	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
417	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
418	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
419	40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40
420	40 40 40 40 40 40 40
421	ending index
422	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
423	20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
424	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
425	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
426	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
427	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
428	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
429	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
430	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
431	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
432	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
433	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
434	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
435	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
436	39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39 39
437	39 39 39 39 39 39 39
438
439	zoom flags :
440	lzoom = F (T = zoom, F = global )
441	lzoom_e = F (T = forced closed east boundary)
442	lzoom_w = F (T = forced closed west boundary)
443	lzoom_s = F (T = forced closed South boundary)
444	lzoom_n = F (T = forced closed North boundary)
445
446	dom_init : domain initialization
447	~~~~~~~~
448
449	dom_nam : domain initialization through namelist read
450	~~~~~~~
451	Namelist namrun
452	job number nn_no = 0
453	experiment name for output cn_exp =
454	ORCA1
455
456
457
458	file prefix restart input cn_ocerst_in=
459	restart
460
461
462
463	restart input directory cn_ocerst_indir=
464	.
465
466
467
468	file prefix restart output cn_ocerst_out=
469	restart
470
471
472
473	restart output directory cn_ocerst_outdir=
474	.
475
476
477
478	restart logical ln_rstart = F
479	start with forward time step nn_euler = 1
480	control of time step nn_rstctl = 0
481	number of the first time step nn_it000 = 1
482	number of the last time step nn_itend = 8760
483	initial calendar date aammjj nn_date0 = 10101
484	leap year calendar (0/1) nn_leapy = 0
485	initial state output nn_istate = 0
486	frequency of restart file nn_stock = 8760
487	frequency of output file nn_write = 8760
488	multi file dimgout ln_dimgnnn = F
489	mask land points ln_mskland = F
490	additional CF standard metadata ln_cfmeta = F
491	overwrite an existing file ln_clobber = F
492	NetCDF chunksize (bytes) nn_chunksz = 0
493
494	===>>> : W A R N I N G
495	===============
496
497	ln_rstart =.FALSE., nn_euler is forced to 0
498
499
500
501	The IOIPSL calendar is "noleap", i.e. no leap year
502
503	Namelist namdom : space & time domain
504	flag read/compute bathymetry nn_bathy = 1
505	Depth (if =0 bathy=jpkm1) rn_bathy = 0.000000000000000E+000
506	min depth of the ocean (>0) or rn_hmin = 20.0000000000000
507	min number of ocean level (<0)
508	minimum thickness of partial rn_e3zps_min = 20.0000000000000
509	(m)
510	step level rn_e3zps_rat = 0.100000000000000
511	create mesh/mask file(s) nn_msh = 1
512	= 0 no file created
513	= 1 mesh_mask
514	= 2 mesh and mask
515	= 3 mesh_hgr, msh_zgr and mask
516	ocean time step rn_rdt =
517	3600.00000000000
518	asselin time filter parameter rn_atfp =
519	0.100000000000000
520	acceleration of converge nn_acc = 0
521	nn_acc=1: surface tracer rdt rn_rdtmin =
522	28800.0000000000
523	bottom tracer rdt rdtmax =
524	28800.0000000000
525	depth of transition rn_rdth =
526	800.000000000000
527	suppression of closed seas (=0) nn_closea = 0
528	online coarsening of dynamical fields ln_crs = F
529	type of horizontal mesh jphgr_msh = 0
530	longitude of first raw and column T-point ppglam0 =
531	999999.000000000
532	latitude of first raw and column T-point ppgphi0 =
533	999999.000000000
534	zonal grid-spacing (degrees) ppe1_deg =
535	999999.000000000
536	meridional grid-spacing (degrees) ppe2_deg =
537	999999.000000000
538	zonal grid-spacing (degrees) ppe1_m =
539	999999.000000000
540	meridional grid-spacing (degrees) ppe2_m =
541	999999.000000000
542	ORCA r4, r2 and r05 coefficients ppsur =
543	-3958.95137127683
544	ppa0 =
545	103.953009600000
546	ppa1 =
547	2.41595126900000
548	ppkth =
549	15.3510137000000
550	ppacr =
551	7.00000000000000
552	Minimum vertical spacing ppdzmin =
553	999999.000000000
554	Maximum depth pphmax =
555	999999.000000000
556	Use double tanf function for vertical coordinates ldbletanh = T
557	Double tanh function parameters ppa2 =
558	100.760928500000
559	ppkth2 =
560	48.0298937200000
561	ppacr2 =
562	13.0000000000000
563
564	Namelist namcla
565	cross land advection nn_cla = 0
566
567	dom_clo : closed seas
568	~~~~~~~
569
570	dom_hgr : define the horizontal mesh from ithe following par_oce parameters
571	~~~~~~~ type of horizontal mesh jphgr_msh = 0
572	position of the first row and ppglam0 =
573	999999.000000000
574	column grid-point (degrees) ppgphi0 =
575	999999.000000000
576	zonal grid-spacing (degrees) ppe1_deg =
577	999999.000000000
578	meridional grid-spacing (degrees) ppe2_deg =
579	999999.000000000
580	zonal grid-spacing (meters) ppe1_m =
581	999999.000000000
582	meridional grid-spacing (meters) ppe2_m =
583	999999.000000000
584
585	curvilinear coordinate on the sphere read in "coordinate" file
586
587	hgr_read : read the horizontal coordinates
588	~~~~~~~~ jpiglo = 362 jpjglo = 292 jpk = 75
589	iom_nf90_open ~~~ open existing file: coordinates.nc in REA
590	D mode
591	---> coordinates.nc OK
592	read glamt (rec: 1) in coordinates.nc ok
593	read glamu (rec: 1) in coordinates.nc ok
594	read glamv (rec: 1) in coordinates.nc ok
595	read glamf (rec: 1) in coordinates.nc ok
596	read gphit (rec: 1) in coordinates.nc ok
597	read gphiu (rec: 1) in coordinates.nc ok
598	read gphiv (rec: 1) in coordinates.nc ok
599	read gphif (rec: 1) in coordinates.nc ok
600	read e1t (rec: 1) in coordinates.nc ok
601	read e1u (rec: 1) in coordinates.nc ok
602	read e1v (rec: 1) in coordinates.nc ok
603	read e1f (rec: 1) in coordinates.nc ok
604	read e2t (rec: 1) in coordinates.nc ok
605	read e2u (rec: 1) in coordinates.nc ok
606	read e2v (rec: 1) in coordinates.nc ok
607	read e2f (rec: 1) in coordinates.nc ok
608	iom_close ~~~ close file: coordinates.nc ok
609
610	orca_r1: Gibraltar : e2u reduced to 20 km
611
612	orca_r1: Bhosporus : e2u reduced to 10 km
613
614	orca_r1: Lombok : e1v reduced to 10 km
615
616	orca_r1: Sumba : e1v reduced to 8 km
617
618	orca_r1: Ombai : e1v reduced to 13 km
619
620	orca_r1: Timor Passage : e1v reduced to 20 km
621
622	orca_r1: W Halmahera : e1v reduced to 30 km
623
624	orca_r1: E Halmahera : e1v reduced to 50 km
625
626	longitude and e1 scale factors
627	------------------------------
628	1 72.50 73.00 72.50 73.00 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808
629	11 82.50 83.00 82.50 83.00 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808
630	21 92.50 93.00 92.50 93.00 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808
631	31 102.50 103.00 102.50 103.00 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808
632
633	latitude and e2 scale factors
634	-----------------------------
635	1 -78.39 -78.39 -78.29 -78.29 22371.9996350562 22371.9996350562 22564.0248736808 22564.0248736808
636	11 -76.20 -76.20 -76.08 -76.08 26524.8299354813 26524.8299354813 26750.5179637307 26750.5179637307
637	21 -73.60 -73.60 -73.46 -73.46 31393.5671869958 31393.5671869958 31657.3897510458 31657.3897510458
638	31 -70.53 -70.53 -70.36 -70.36 37065.0717616678 37065.0717616678 37371.1271782116 37371.1271782116
639
640	dom_zgr : vertical coordinate
641	~~~~~~~
642	Namelist namzgr : set vertical coordinate
643	z-coordinate - full steps ln_zco = F
644	z-coordinate - partial steps ln_zps = T
645	s- or hybrid z-s-coordinate ln_sco = F
646	ice shelf cavities ln_isfcav = F
647
648	zgr_z : Reference vertical z-coordinates
649	~~~~~~~
650	Value of coefficients for vertical mesh:
651	zsur = -3958.95137127683
652	za0 = 103.953009600000
653	za1 = 2.41595126900000
654	zkth = 15.3510137000000
655	zacr = 7.00000000000000
656	(Double tanh za2 = 100.760928500000
657	parameters) zkth2= 48.0298937200000
658	zacr2= 13.0000000000000
659
660	Reference z-coordinate depth and scale factors:
661	level gdept_1d gdepw_1d e3t_1d e3w_1d
662	1 0.51 0.00 1.02 1.00
663	2 1.56 1.02 1.08 1.05
664	3 2.67 2.10 1.15 1.11
665	4 3.86 3.25 1.23 1.19
666	5 5.14 4.49 1.34 1.28
667	6 6.54 5.83 1.47 1.40
668	7 8.09 7.30 1.63 1.55
669	8 9.82 8.93 1.83 1.73
670	9 11.77 10.77 2.08 1.95
671	10 13.99 12.85 2.37 2.22
672	11 16.53 15.22 2.71 2.53
673	12 19.43 17.93 3.11 2.90
674	13 22.76 21.04 3.56 3.33
675	14 26.56 24.60 4.05 3.80
676	15 30.87 28.65 4.59 4.31
677	16 35.74 33.24 5.15 4.86
678	17 41.18 38.39 5.73 5.44
679	18 47.21 44.12 6.33 6.03
680	19 53.85 50.45 6.95 6.64
681	20 61.11 57.40 7.58 7.26
682	21 69.02 64.98 8.24 7.91
683	22 77.61 73.23 8.94 8.59
684	23 86.93 82.18 9.70 9.32
685	24 97.04 91.88 10.53 10.11
686	25 108.03 102.42 11.46 10.98
687	26 120.00 113.89 12.50 11.96
688	27 133.08 126.39 13.68 13.07
689	28 147.41 140.07 15.01 14.32
690	29 163.16 155.10 16.54 15.75
691	30 180.55 171.64 18.27 17.38
692	31 199.79 189.92 20.25 19.23
693	32 221.14 210.18 22.50 21.34
694	33 244.89 232.70 25.05 23.74
695	34 271.36 257.76 27.94 26.45
696	35 300.89 285.72 31.19 29.52
697	36 333.86 316.92 34.83 32.96
698	37 370.69 351.77 38.89 36.81
699	38 411.79 390.68 43.39 41.09
700	39 457.63 434.09 48.35 45.81
701	40 508.64 482.46 53.76 51.00
702	41 565.29 536.23 59.62 56.63
703	42 628.03 595.87 65.92 62.72
704	43 697.26 661.81 72.61 69.22
705	44 773.37 734.43 79.66 76.10
706	45 856.68 814.11 87.00 83.30
707	46 947.45 901.12 94.56 90.76
708	47 1045.85 995.69 102.26 98.40
709	48 1151.99 1097.95 110.01 106.14
710	49 1265.86 1207.96 117.71 113.87
711	50 1387.38 1325.67 125.29 121.52
712	51 1516.36 1450.95 132.64 129.00
713	52 1652.57 1583.58 139.71 136.22
714	53 1795.67 1723.28 146.43 143.12
715	54 1945.30 1869.69 152.75 149.64
716	55 2101.03 2022.42 158.64 155.75
717	56 2262.42 2181.04 164.08 161.41
718	57 2429.03 2345.10 169.06 166.62
719	58 2600.38 2514.14 173.58 171.37
720	59 2776.04 2687.70 177.67 175.68
721	60 2955.57 2865.35 181.33 179.55
722	61 3138.56 3046.66 184.60 183.01
723	62 3324.64 3231.24 187.50 186.09
724	63 3513.45 3418.72 190.06 188.82
725	64 3704.66 3608.77 192.31 191.22
726	65 3897.98 3801.07 194.29 193.34
727	66 4093.16 3995.35 196.02 195.19
728	67 4289.95 4191.37 197.53 196.80
729	68 4488.15 4388.89 198.84 198.21
730	69 4687.58 4587.73 199.98 199.43
731	70 4888.07 4787.70 200.97 200.49
732	71 5089.48 4988.67 201.83 201.41
733	72 5291.68 5190.49 202.57 202.21
734	73 5494.58 5393.05 203.20 202.90
735	74 5698.06 5596.25 203.75 203.49
736	75 5902.06 5800.00 204.23 204.00
737
738	zgr_bat : defines level and meter bathymetry
739	~~~~~~~
740	iom_nf90_open ~~~ open existing file: bathy_meter.nc in REA
741	D mode
742	---> bathy_meter.nc OK
743	read Bathymetry (rec: 1) in bathy_meter.nc ok
744	iom_close ~~~ close file: bathy_meter.nc ok
745	Minimum ocean depth: 24.5959881615020 minimum number of ocean levels :
746	13
747
748	zgr_zps : z-coordinate with partial steps
749	~~~~~~~
750	mbathy is recomputed : bathy_level file is NOT used
751
752	zgr_bat_ctl : check the bathymetry
753	~~~~~~~~~~~
754
755	suppress isolated ocean grid points
756	-----------------------------------
757	1064 ocean grid points suppressed
758	mbathy set to 0 along east and west boundary: nperio = 0
759	maximum number of ocean level = 43 < jpk-1
760	you can decrease jpk to 44
761
762	zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels
763	~~~~~~~~~~~~~
764
765	zgr_top_level : ocean top k-index of T-, U-, V- and W-levels
766	~~~~~~~~~~~~~
767
768	dommsk : ocean mask
769	~~~~~~
770	Namelist namlbc
771	lateral momentum boundary cond. rn_shlat = 0.000000000000000E+000
772	consistency with analytical form ln_vorlat = F
773	ocean lateral free-slip
774
775	orca_r1: increase friction near the following straits :
776	Gibraltar
777	Bhosporus
778	Makassar (Top)
779	Lombok
780	Ombai
781	Timor Passage
782	West Halmahera
783	East Halmahera
784
785	dom_stp : time stepping setting
786	~~~~~~~
787	synchronous time stepping
788	dynamics and tracer time step = 1.00000000000000 hours
789
790	dom_wri : create NetCDF mesh and mask information file(s)
791	~~~~~~~
792	iom_nf90_open ~~~ create new file: mesh_mask_0000.nc in WRI
793	TE mode
794	---> mesh_mask_0000.nc OK
795	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask define dimensio
796	n variables done
797	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask defined ok
798	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask write dimension
799	variables done
800	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask written ok
801	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umask defined ok
802	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umask written ok
803	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmask defined ok
804	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmask written ok
805	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmask defined ok
806	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmask written ok
807	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmaskutil defined ok
808	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmaskutil written ok
809	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umaskutil defined ok
810	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umaskutil written ok
811	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmaskutil defined ok
812	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmaskutil written ok
813	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmaskutil defined ok
814	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmaskutil written ok
815	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamt defined ok
816	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamt written ok
817	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamu defined ok
818	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamu written ok
819	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamv defined ok
820	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamv written ok
821	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamf defined ok
822	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamf written ok
823	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphit defined ok
824	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphit written ok
825	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiu defined ok
826	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiu written ok
827	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiv defined ok
828	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiv written ok
829	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphif defined ok
830	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphif written ok
831	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1t defined ok
832	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1t written ok
833	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1u defined ok
834	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1u written ok
835	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1v defined ok
836	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1v written ok
837	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1f defined ok
838	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1f written ok
839	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2t defined ok
840	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2t written ok
841	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2u defined ok
842	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2u written ok
843	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2v defined ok
844	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2v written ok
845	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2f defined ok
846	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2f written ok
847	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: ff defined ok
848	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: ff written ok
849	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: mbathy defined ok
850	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: mbathy written ok
851	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: misf defined ok
852	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: misf written ok
853	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: isfdraft defined ok
854	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: isfdraft written ok
855	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_0 defined ok
856	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_0 written ok
857	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3u_0 defined ok
858	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3u_0 written ok
859	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3v_0 defined ok
860	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3v_0 written ok
861	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_0 defined ok
862	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_0 written ok
863	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_0 defined ok
864	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_0 written ok
865	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepu defined ok
866	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepu written ok
867	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepv defined ok
868	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepv written ok
869	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_0 defined ok
870	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_0 written ok
871	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_1d defined ok
872	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_1d written ok
873	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_1d defined ok
874	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_1d written ok
875	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_1d defined ok
876	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_1d written ok
877	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_1d defined ok
878	iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_1d written ok
879	iom_close ~~~ close file: mesh_mask_0000.nc ok
880
881	dom_ctl : extrema of the masked scale factors
882	~~~~~~~
883	e1t maxi: 111198.92 at i = 2 j= 147
884	e1t mini: 25207.90 at i = 92 j= 8
885	e2t maxi: 105405.52 at i = 2 j= 114
886	e2t mini: 24988.38 at i = 208 j= 283
887
888	istate_ini : Initialization of the dynamics and tracers
889	~~~~~~~~~~
890
891	dta_tsd_init : Temperature & Salinity data
892	~~~~~~~~~~~~
893	Namelist namtsd
894	Initialisation of ocean T & S with T &S input data ln_tsd_init = T
895	damping of ocean T & S toward T &S input data ln_tsd_tradmp = F
896
897
898	dta_tsd : Temperature & Salinity data
899	~~~~~~~
900	namtsd Namelist
901	list of files and frequency (>0: in hours ; <0 in months)
902	root filename: ./data_1m_potential_temperature_nomask
903	variable name: votemper
904	frequency: -1.00000000000000 time interp: F
905	climatology: T weights : pairing : data type: yearly
906	land/sea mask:
907	root filename: ./data_1m_salinity_nomask variable name: vosaline
908	frequency: -1.00000000000000 time interp: F
909	climatology: T weights : pairing : data type: yearly
910	land/sea mask:
911	*** Info used values :
912	date ndastp : 10100
913	number of elapsed days since the begining of run : 0.000000000000000E+000
914
915	=======>> 1/2 time step before the start of the run DATE Y/M/D = 0/12/31 nsec_day: 84600 nsec_week: -174600
916	======>> time-step = 1 New day, DATE Y/M/D = 0001/01/01 nday_year = 001
917	nsec_year = 1800 nsec_month = 1800 nsec_day = 1800
918	nsec_week = -171000
919	iom_nf90_open ~~~ open existing file: ./data_1m_potential_t
920	emperature_nomask.nc in READ mode
921	---> ./data_1m_potential_temperature_nomask.nc OK
922	read votemper (rec: 1) in ./data_1m_potential_temperature_nomask.nc ok
923	iom_nf90_open ~~~ open existing file: ./data_1m_salinity_no
924	mask.nc in READ mode
925	---> ./data_1m_salinity_nomask.nc OK
926	read vosaline (rec: 1) in ./data_1m_salinity_nomask.nc ok
927	fld_read: var votemper kt = 1 ( 0.0208 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 31.0000)
928	fld_read: var vosaline kt = 1 ( 0.0208 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 31.0000)
929	temperature Levitus
930
931	level = 1
932
933
934	1 21
935	39 0.79 0.00
936	19 0.00 0.00
937	level = 37
938
939
940	1 21
941	39 -1.15 0.00
942	19 0.00 0.00
943	level = 74
944
945
946	1 21
947	39 0.00 0.00
948	19 0.00 0.00
949
950	salinity Levitus
951
952	level = 1
953
954
955	1 21
956	39 33.69 0.00
957	19 0.00 0.00
958	level = 37
959
960
961	1 21
962	39 34.54 0.00
963	19 0.00 0.00
964	level = 74
965
966
967	1 21
968	39 0.00 0.00
969	19 0.00 0.00
970
971	dta_tsd: deallocte T & S arrays as they are only use to initialize the run
972
973	sbc_init : surface boundary condition setting
974	~~~~~~~~
975	Namelist namsbc (partly overwritten with CPP key setting)
976	frequency update of sbc (and ice) nn_fsbc =
977	3
978	Type of sbc :
979	analytical formulation ln_ana = F
980	flux formulation ln_flx = F
981	CLIO bulk formulation ln_blk_clio = F
982	CORE bulk formulation ln_blk_core = T
983	MFS bulk formulation ln_blk_mfs = F
984	ocean-atmosphere coupled formulation ln_cpl = F
985	forced-coupled mixed formulation ln_mixcpl = F
986	OASIS coupling (with atm or sas) lk_oasis = F
987	components of your executable nn_components =
988	0
989	Multicategory heat flux formulation (LIM3) nn_limflx =
990	-1
991	Misc. options of sbc :
992	Patm gradient added in ocean & ice Eqs. ln_apr_dyn = F
993	ice management in the sbc (=0/1/2/3) nn_ice =
994	2
995	ice-ocean embedded/levitating (=0/1/2) nn_ice_embd =
996	1
997	daily mean to diurnal cycle qsr ln_dm2dc = T
998	runoff / runoff mouths ln_rnf = T
999	iceshelf formulation nn_isf =
1000	0
1001	Sea Surface Restoring on SST and/or SSS ln_ssr = T
1002	FreshWater Budget control (=0/1/2) nn_fwb =
1003	1
1004	closed sea (=0/1) (set in namdom) nn_closea =
1005	0
1006	n. of iterations if land-sea-mask applied nn_lsm =
1007	0
1008	Use of per-category fluxes (nn_limflx = -1)
1009
1010	CORE bulk formulation
1011
1012	sbc_ssm : sea surface mean fields
1013	~~~~~~~
1014	default initialisation of ss?_m arrays
1015
1016	sbc_ssr : SST and/or SSS damping term
1017	~~~~~~~
1018	Namelist namsbc_ssr :
1019	SST restoring term (Yes=1) nn_sstr = 0
1020	SSS damping term (Yes=1, salt flux) nn_sssr = 2
1021	(Yes=2, volume flux)
1022	dQ/dT (restoring magnitude on SST) rn_dqdt =
1023	-40.0000000000000 W/m2/K
1024	dE/dS (restoring magnitude on SST) rn_deds =
1025	-166.670000000000 mm/day
1026	flag to bound erp term ln_sssr_bnd = F
1027	ABS(Max./Min.) erp threshold rn_sssr_bnd =
1028	4.00000000000000 mm/day
1029
1030	sbc_ssr : SSS restoring term toward SSS data
1031	~~~~~~~
1032	namsbc_ssr Namelist
1033	list of files and frequency (>0: in hours ; <0 in months)
1034	root filename: ./sss_data variable name: sss
1035	frequency: -1.00000000000000 time interp: T
1036	climatology: T weights : pairing : data type: yearly
1037	land/sea mask:
1038
1039	sbc_rnf : runoff
1040	~~~~~~~
1041	Namelist namsbc_rnf
1042	specific river mouths treatment ln_rnf_mouth = T
1043	river mouth additional Kz rn_avt_rnf =
1044	1.000000000000000E-003
1045	depth of river mouth additional mixing rn_hrnf =
1046	15.0000000000000
1047	multiplicative factor for runoff rn_rfact =
1048	1.00000000000000
1049
1050	runoffs inflow read in a file
1051
1052	sbc_rnf_init : read runoffs data
1053	~~~~~~~~~~~~
1054	namsbc_rnf Namelist
1055	list of files and frequency (>0: in hours ; <0 in months)
1056	root filename: ./runoff_1m_ORCA1 variable name: sorunoff
1057	frequency: -1.00000000000000 time interp: T
1058	climatology: T weights : pairing : data type: yearly
1059	land/sea mask:
1060
1061	Specific treatment used in vicinity of river mouths :
1062	- Increase Kz in surface layers (if rn_hrnf > 0 )
1063	by 1.000000000000000E-003 m2/s over 10 w-levels
1064	- set to zero SSS damping (if ln_ssr=T)
1065	- mixed upstream-centered (if ln_traadv_cen2=T)
1066
1067	rnf_mouth : river mouth mask
1068	~~~~~~~~~
1069	iom_nf90_open ~~~ open existing file: ./runoff_1m_ORCA1.nc
1070	in READ mode
1071	---> ./runoff_1m_ORCA1.nc OK
1072	read socoeff (rec: 1) in ./runoff_1m_ORCA1.nc ok
1073	iom_close ~~~ close file: ./runoff_1m_ORCA1.nc ok
1074
1075	dyn_nept_init : Simplified Neptune module
1076	~~~~~~~~~~~~~
1077	--> Reading namelist namdyn_nept parameters:
1078	ln_neptsimp = F
1079
1080
1081	zdf_init: vertical physics
1082	~~~~~~~~
1083	Namelist namzdf : set vertical mixing mixing parameters
1084	vertical eddy viscosity rn_avm0 = 1.200000000000000E-004
1085	vertical eddy diffusivity rn_avt0 = 1.200000000000000E-005
1086	constant background or profile nn_avb = 0
1087	horizontal variation for avtb nn_havtb = 0
1088	time splitting / backward scheme ln_zdfexp = F
1089	number of time step nn_zdfexp = 3
1090	enhanced vertical diffusion ln_zdfevd = T
1091	applied on momentum (=1/0) nn_evdm = 0
1092	vertical coefficient for evd rn_avevd = 100.000000000000
1093	non-penetrative convection (npc) ln_zdfnpc = F
1094	npc call frequency nn_npc = 1
1095	npc print frequency nn_npcp = 365
1096
1097	vertical mixing option :
1098	TKE dependent eddy coefficients
1099
1100	convection :
1101	use enhanced vertical dif. scheme
1102	use the 1.5 turbulent closure
1103
1104	zdf_bfr_init : momentum bottom friction
1105	~~~~~~~~~~~~~
1106	Namelist nam_bfr : set bottom friction parameters
1107	quadratic bottom friction
1108	friction coef. rn_bfri2 = 1.000000000000000E-003
1109	Max. coef. (log case) rn_bfri2_max = 0.100000000000000
1110	background tke rn_bfeb2 = 2.500000000000000E-003
1111	log formulation ln_bfr2d = F
1112	bottom roughness rn_bfrz0 [m] = 3.000000000000000E-003
1113	implicit bottom friction switch ln_bfrimp = T
1114
1115	zdf_tke_init : tke turbulent closure scheme - initialisation
1116	~~~~~~~~~~~~
1117	Namelist namzdf_tke : set tke mixing parameters
1118	coef. to compute avt rn_ediff =
1119	0.100000000000000
1120	Kolmogoroff dissipation coef. rn_ediss =
1121	0.700000000000000
1122	tke surface input coef. rn_ebb =
1123	67.8300000000000
1124	minimum value of tke rn_emin =
1125	1.000000000000000E-006
1126	surface minimum value of tke rn_emin0 =
1127	1.000000000000000E-004
1128	background shear (>0) rn_bshear =
1129	9.999999999999999E-021
1130	mixing length type nn_mxl = 2
1131	prandl number flag nn_pdl = 1
1132	surface mixing length = F(stress) or not ln_mxl0 = T
1133	surface mixing length minimum value rn_mxl0 =
1134	4.000000000000000E-002
1135	flag to take into acc. Langmuir circ. ln_lc = T
1136	coef to compute verticla velocity of LC rn_lc =
1137	0.150000000000000
1138	test param. to add tke induced by wind nn_etau = 1
1139	flag for computation of exp. tke profile nn_htau = 1
1140	fraction of en which pene. the thermocline rn_efr =
1141	5.000000000000000E-002
1142
1143	critical Richardson nb with your parameters ri_cri =
1144	0.222222222222222
1145	use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min
1146
1147	zdf_tmx_init : tidal mixing
1148	~~~~~~~~~~~~
1149	Namelist namzdf_tmx : set tidal mixing parameters
1150	Vertical decay scale for turbulence = 500.000000000000
1151	Brunt-Vaisala frequency threshold = 1.000000000000000E-008
1152	Tidal dissipation efficiency = 0.333000000000000
1153	Mixing efficiency = 0.200000000000000
1154	ITF specific parameterisation = T
1155	ITF tidal dissipation efficiency = 1.00000000000000
1156	iom_nf90_open ~~~ open existing file: mask_itf.nc in READ m
1157	ode
1158	---> mask_itf.nc OK
1159	read tmaskitf (rec: 1) in mask_itf.nc ok
1160	iom_close ~~~ close file: mask_itf.nc ok
1161	iom_nf90_open ~~~ open existing file: M2rowdrg.nc in READ m
1162	ode
1163	---> M2rowdrg.nc OK
1164	read field (rec: 1) in M2rowdrg.nc ok
1165	iom_close ~~~ close file: M2rowdrg.nc ok
1166	iom_nf90_open ~~~ open existing file: K1rowdrg.nc in READ m
1167	ode
1168	---> K1rowdrg.nc OK
1169	read field (rec: 1) in K1rowdrg.nc ok
1170	iom_close ~~~ close file: K1rowdrg.nc ok
1171
1172	zdf_ddm : double diffusive mixing
1173	~~~~~~~
1174	Namelist namzdf_ddm : set dd mixing parameter
1175	maximum avs for dd mixing rn_avts = 1.000000000000000E-004
1176	heat/salt buoyancy flux ratio rn_hsbfr = 1.60000000000000
1177
1178	ldf_tra_init : lateral tracer physics
1179	~~~~~~~~~~~~
1180	Namelist namtra_ldf : lateral mixing parameters (type, direction, coefficien
1181	ts)
1182	laplacian operator ln_traldf_lap = T
1183	bilaplacian operator ln_traldf_bilap = F
1184	iso-level ln_traldf_level = F
1185	horizontal (geopotential) ln_traldf_hor = F
1186	iso-neutral ln_traldf_iso = T
1187	iso-neutral (Griffies) ln_traldf_grif = F
1188	Griffies strmfn diagnostics ln_traldf_gdia = F
1189	lateral eddy diffusivity rn_aht_0 = 1000.00000000000
1190	background hor. diffusivity rn_ahtb_0 = 0.000000000000000E+000
1191	eddy induced velocity coef. rn_aeiv_0 = 1000.00000000000
1192	maximum isoppycnal slope rn_slpmax = 1.000000000000000E-002
1193	pure lateral mixing in ML ln_triad_iso = F
1194	lateral mixing on bottom ln_botmix_grif = F
1195
1196	tracer mixing coef. = F( latitude, longitude)
1197	harmonic tracer diffusion (default)
1198
1199	ldf_tra_c2d : 2D eddy diffusivity and eddy
1200	~~~~~~~~~~~ -- induced velocity coefficients
1201
1202	ldf_dyn : lateral momentum physics
1203	~~~~~~~
1204	Namelist namdyn_ldf : set lateral mixing parameters
1205	laplacian operator ln_dynldf_lap = T
1206	bilaplacian operator ln_dynldf_bilap = F
1207	iso-level ln_dynldf_level = F
1208	horizontal (geopotential) ln_dynldf_hor = T
1209	iso-neutral ln_dynldf_iso = F
1210	horizontal laplacian eddy viscosity rn_ahm_0_lap =
1211	20000.0000000000
1212	background viscosity rn_ahmb_0 =
1213	0.000000000000000E+000
1214	horizontal bilaplacian eddy viscosity rn_ahm_0_blp =
1215	0.000000000000000E+000
1216	upper limit for laplacian eddy visc rn_ahm_m_lap =
1217	40000.0000000000
1218	upper limit for bilap eddy viscosity rn_ahm_m_blp =
1219	-1000000000000.00
1220	momentum mixing coef. = F( latitude, longitude, depth)
1221	harmonic momentum diff. (default)
1222
1223	ldf_dyn_c3d : 3D lateral eddy viscosity coefficient
1224	~~~~~~~~~~~
1225	laplacian operator: ahm proportional to e1
1226	maximum grid-spacing = 111198.923448546
1227	maximum value for ahm = 20000.0000000000
1228
1229	ahm profile :
1230
1231	jk ahm depth t-level
1232	1 1.0000 0.5058
1233	2 0.9999 1.5559
1234	3 0.9999 2.6677
1235	4 0.9998 3.8563
1236	5 0.9997 5.1404
1237	6 0.9997 6.5430
1238	7 0.9996 8.0925
1239	8 0.9995 9.8228
1240	9 0.9994 11.7737
1241	10 0.9993 13.9910
1242	11 0.9991 16.5253
1243	12 0.9989 19.4298
1244	13 0.9987 22.7576
1245	14 0.9985 26.5583
1246	15 0.9983 30.8746
1247	16 0.9980 35.7402
1248	17 0.9976 41.1800
1249	18 0.9972 47.2119
1250	19 0.9968 53.8506
1251	20 0.9963 61.1128
1252	21 0.9958 69.0217
1253	22 0.9951 77.6112
1254	23 0.9945 86.9294
1255	24 0.9937 97.0413
1256	25 0.9928 108.0303
1257	26 0.9918 120.0000
1258	27 0.9907 133.0758
1259	28 0.9893 147.4062
1260	29 0.9878 163.1645
1261	30 0.9860 180.5499
1262	31 0.9839 199.7900
1263	32 0.9813 221.1412
1264	33 0.9782 244.8906
1265	34 0.9745 271.3564
1266	35 0.9698 300.8875
1267	36 0.9640 333.8628
1268	37 0.9567 370.6885
1269	38 0.9474 411.7938
1270	39 0.9355 457.6256
1271	40 0.9199 508.6399
1272	41 0.8996 565.2923
1273	42 0.8731 628.0260
1274	43 0.8385 697.2586
1275	44 0.7940 773.3683
1276	45 0.7384 856.6789
1277	46 0.6719 947.4479
1278	47 0.5970 1045.8543
1279	48 0.5192 1151.9912
1280	49 0.4460 1265.8614
1281	50 0.3838 1387.3770
1282	51 0.3359 1516.3636
1283	52 0.3023 1652.5684
1284	53 0.2804 1795.6708
1285	54 0.2670 1945.2955
1286	55 0.2592 2101.0265
1287	56 0.2549 2262.4216
1288	57 0.2525 2429.0252
1289	58 0.2513 2600.3805
1290	59 0.2506 2776.0393
1291	60 0.2503 2955.5704
1292	61 0.2501 3138.5649
1293	62 0.2501 3324.6408
1294	63 0.2500 3513.4456
1295	64 0.2500 3704.6567
1296	65 0.2500 3897.9819
1297	66 0.2500 4093.1587
1298	67 0.2500 4289.9524
1299	68 0.2500 4488.1546
1300	69 0.2500 4687.5811
1301	70 0.2500 4888.0698
1302	71 0.2500 5089.4786
1303	72 0.2500 5291.6832
1304	73 0.2500 5494.5753
1305	74 0.2500 5698.0608
1306	75 0.2500 5902.0578
1307
1308	ahm profile :
1309
1310	jk ahm depth t-level
1311	1 1.0000 0.5058
1312	2 0.9999 1.5559
1313	3 0.9999 2.6677
1314	4 0.9998 3.8563
1315	5 0.9997 5.1404
1316	6 0.9997 6.5430
1317	7 0.9996 8.0925
1318	8 0.9995 9.8228
1319	9 0.9994 11.7737
1320	10 0.9993 13.9910
1321	11 0.9991 16.5253
1322	12 0.9989 19.4298
1323	13 0.9987 22.7576
1324	14 0.9985 26.5583
1325	15 0.9983 30.8746
1326	16 0.9980 35.7402
1327	17 0.9976 41.1800
1328	18 0.9972 47.2119
1329	19 0.9968 53.8506
1330	20 0.9963 61.1128
1331	21 0.9958 69.0217
1332	22 0.9951 77.6112
1333	23 0.9945 86.9294
1334	24 0.9937 97.0413
1335	25 0.9928 108.0303
1336	26 0.9918 120.0000
1337	27 0.9907 133.0758
1338	28 0.9893 147.4062
1339	29 0.9878 163.1645
1340	30 0.9860 180.5499
1341	31 0.9839 199.7900
1342	32 0.9813 221.1412
1343	33 0.9782 244.8906
1344	34 0.9745 271.3564
1345	35 0.9698 300.8875
1346	36 0.9640 333.8628
1347	37 0.9567 370.6885
1348	38 0.9474 411.7938
1349	39 0.9355 457.6256
1350	40 0.9199 508.6399
1351	41 0.8996 565.2923
1352	42 0.8731 628.0260
1353	43 0.8385 697.2586
1354	44 0.7940 773.3683
1355	45 0.7384 856.6789
1356	46 0.6719 947.4479
1357	47 0.5970 1045.8543
1358	48 0.5192 1151.9912
1359	49 0.4460 1265.8614
1360	50 0.3838 1387.3770
1361	51 0.3359 1516.3636
1362	52 0.3023 1652.5684
1363	53 0.2804 1795.6708
1364	54 0.2670 1945.2955
1365	55 0.2592 2101.0265
1366	56 0.2549 2262.4216
1367	57 0.2525 2429.0252
1368	58 0.2513 2600.3805
1369	59 0.2506 2776.0393
1370	60 0.2503 2955.5704
1371	61 0.2501 3138.5649
1372	62 0.2501 3324.6408
1373	63 0.2500 3513.4456
1374	64 0.2500 3704.6567
1375	65 0.2500 3897.9819
1376	66 0.2500 4093.1587
1377	67 0.2500 4289.9524
1378	68 0.2500 4488.1546
1379	69 0.2500 4687.5811
1380	70 0.2500 4888.0698
1381	71 0.2500 5089.4786
1382	72 0.2500 5291.6832
1383	73 0.2500 5494.5753
1384	74 0.2500 5698.0608
1385	75 0.2500 5902.0578
1386
1387	ORCA R1, R2 or R4: overwrite the previous definition of ahm
1388	=================
1389
1390	ldfdyn_c3d_orca : 3D eddy viscosity coefficient
1391	~~~~~~~~~~~~~~~
1392	orca R1, R2 or R4 configuration: reduced in the surface Eq. strip
1393	iom_nf90_open ~~~ open existing file: ahmcoef.nc in READ mo
1394	de
1395	---> ahmcoef.nc OK
1396	read icof (rec: 1) in ahmcoef.nc ok
1397	iom_close ~~~ close file: ahmcoef.nc ok
1398
1399	1D zcoef array
1400	~~~~~~~~~~~~~~
1401
1402	jk zcoef
1403	1 1.000
1404	2 1.000
1405	3 1.000
1406	4 1.000
1407	5 1.000
1408	6 1.000
1409	7 1.000
1410	8 1.000
1411	9 1.000
1412	10 1.000
1413	11 1.000
1414	12 1.000
1415	13 1.000
1416	14 1.000
1417	15 1.000
1418	16 1.000
1419	17 1.000
1420	18 1.000
1421	19 1.000
1422	20 1.000
1423	21 1.000
1424	22 1.000
1425	23 1.000
1426	24 1.000
1427	25 1.000
1428	26 1.000
1429	27 1.000
1430	28 1.000
1431	29 1.000
1432	30 1.000
1433	31 1.000
1434	32 1.000
1435	33 1.000
1436	34 1.000
1437	35 1.000
1438	36 1.000
1439	37 1.000
1440	38 1.000
1441	39 1.000
1442	40 1.000
1443	41 1.000
1444	42 1.000
1445	43 1.000
1446	44 1.000
1447	45 1.000
1448	46 2.000
1449	47 2.000
1450	48 2.000
1451	49 3.000
1452	50 3.000
1453	51 4.000
1454	52 4.000
1455	53 5.000
1456	54 6.000
1457	55 6.000
1458	56 7.000
1459	57 8.000
1460	58 8.000
1461	59 9.000
1462	60 10.000
1463	61 10.000
1464	62 10.000
1465	63 10.000
1466	64 10.000
1467	65 10.000
1468	66 10.000
1469	67 10.000
1470	68 10.000
1471	69 10.000
1472	70 10.000
1473	71 10.000
1474	72 10.000
1475	73 10.000
1476	74 10.000
1477	75 10.000
1478
1479	3D ahm1 array (k=1)
1480
1481
1482	1 21
1483	39 20.00 20.00
1484	19 20.00 20.00
1485
1486	3D ahm2 array (k=1)
1487
1488
1489	1 21
1490	39 20.00 20.00
1491	19 20.00 20.00
1492
1493	3D ahm2 array (k=jpk)
1494
1495
1496	1 21
1497	39 20.00 20.00
1498	19 20.00 20.00
1499
1500	ldf_slp_init : direction of lateral mixing
1501	~~~~~~~~~~~~
1502
1503	tra_qsr_init : penetration of the surface solar radiation
1504	~~~~~~~~~~~~
1505	Namelist namtra_qsr : set the parameter of penetration
1506	Light penetration (T) or not (F) ln_traqsr = T
1507	RGB (Red-Green-Blue) light penetration ln_qsr_rgb = T
1508	2 band light penetration ln_qsr_2bd = F
1509	bio-model light penetration ln_qsr_bio = F
1510	light penetration for ice-model LIM3 ln_qsr_ice = T
1511	RGB : Chl data (=1/2) or cst value (=0) nn_chldta = 1
1512	RGB & 2 bands: fraction of light (rn_si1) rn_abs =
1513	0.580000000000000
1514	RGB & 2 bands: shortess depth of extinction rn_si0 =
1515	0.350000000000000
1516	2 bands: longest depth of extinction rn_si1 =
1517	23.0000000000000
1518
1519	R-G-B light penetration - 2D Chl data
1520
1521	trc_oce_rgb : Initialisation of the optical look-up table
1522	~~~~~~~~~~~
1523	RGB longest depth of extinction r_si2 = 61.8046971569839
1524	jc = 1 Chl = 1.000000000000000E-002 irgb = 1
1525	jc = 2 Chl = 1.100000000000000E-002 irgb = 2
1526	jc = 3 Chl = 1.300000000000000E-002 irgb = 3
1527	jc = 4 Chl = 1.400000000000000E-002 irgb = 4
1528	jc = 5 Chl = 1.600000000000000E-002 irgb = 5
1529	jc = 6 Chl = 1.800000000000000E-002 irgb = 6
1530	jc = 7 Chl = 2.000000000000000E-002 irgb = 7
1531	jc = 8 Chl = 2.200000000000000E-002 irgb = 8
1532	jc = 9 Chl = 2.500000000000000E-002 irgb = 9
1533	jc = 10 Chl = 2.800000000000000E-002 irgb = 10
1534	jc = 11 Chl = 3.200000000000000E-002 irgb = 11
1535	jc = 12 Chl = 3.500000000000000E-002 irgb = 12
1536	jc = 13 Chl = 4.000000000000000E-002 irgb = 13
1537	jc = 14 Chl = 4.500000000000000E-002 irgb = 14
1538	jc = 15 Chl = 5.000000000000000E-002 irgb = 15
1539	jc = 16 Chl = 5.600000000000000E-002 irgb = 16
1540	jc = 17 Chl = 6.300000000000000E-002 irgb = 17
1541	jc = 18 Chl = 7.099999999999999E-002 irgb = 18
1542	jc = 19 Chl = 7.900000000000000E-002 irgb = 19
1543	jc = 20 Chl = 8.900000000000000E-002 irgb = 20
1544	jc = 21 Chl = 0.100000000000000 irgb = 21
1545	jc = 22 Chl = 0.112000000000000 irgb = 22
1546	jc = 23 Chl = 0.126000000000000 irgb = 23
1547	jc = 24 Chl = 0.141000000000000 irgb = 24
1548	jc = 25 Chl = 0.158000000000000 irgb = 25
1549	jc = 26 Chl = 0.178000000000000 irgb = 26
1550	jc = 27 Chl = 0.200000000000000 irgb = 27
1551	jc = 28 Chl = 0.224000000000000 irgb = 28
1552	jc = 29 Chl = 0.251000000000000 irgb = 29
1553	jc = 30 Chl = 0.282000000000000 irgb = 30
1554	jc = 31 Chl = 0.316000000000000 irgb = 31
1555	jc = 32 Chl = 0.355000000000000 irgb = 32
1556	jc = 33 Chl = 0.398000000000000 irgb = 33
1557	jc = 34 Chl = 0.447000000000000 irgb = 34
1558	jc = 35 Chl = 0.501000000000000 irgb = 35
1559	jc = 36 Chl = 0.562000000000000 irgb = 36
1560	jc = 37 Chl = 0.631000000000000 irgb = 37
1561	jc = 38 Chl = 0.708000000000000 irgb = 38
1562	jc = 39 Chl = 0.794000000000000 irgb = 39
1563	jc = 40 Chl = 0.891000000000000 irgb = 40
1564	jc = 41 Chl = 1.00000000000000 irgb = 41
1565	jc = 42 Chl = 1.12200000000000 irgb = 42
1566	jc = 43 Chl = 1.25900000000000 irgb = 43
1567	jc = 44 Chl = 1.41300000000000 irgb = 44
1568	jc = 45 Chl = 1.58500000000000 irgb = 45
1569	jc = 46 Chl = 1.77800000000000 irgb = 46
1570	jc = 47 Chl = 1.99500000000000 irgb = 47
1571	jc = 48 Chl = 2.23900000000000 irgb = 48
1572	jc = 49 Chl = 2.51200000000000 irgb = 49
1573	jc = 50 Chl = 2.81800000000000 irgb = 50
1574	jc = 51 Chl = 3.16200000000000 irgb = 51
1575	jc = 52 Chl = 3.54800000000000 irgb = 52
1576	jc = 53 Chl = 3.98100000000000 irgb = 53
1577	jc = 54 Chl = 4.46700000000000 irgb = 54
1578	jc = 55 Chl = 5.01200000000000 irgb = 55
1579	jc = 56 Chl = 5.62300000000000 irgb = 56
1580	jc = 57 Chl = 6.31000000000000 irgb = 57
1581	jc = 58 Chl = 7.07900000000000 irgb = 58
1582	jc = 59 Chl = 7.94300000000000 irgb = 59
1583	jc = 60 Chl = 8.91200000000000 irgb = 60
1584	jc = 61 Chl = 10.0000000000000 irgb = 61
1585	level of light extinction = 44 ref depth =
1586	814.105680447617 m
1587
1588	Chlorophyll read in a file
1589
1590	tra_qsr_init : Solar penetration function of read chlorophyll
1591	~~~~~~~~~~~~
1592	namtra_qsr Namelist
1593	list of files and frequency (>0: in hours ; <0 in months)
1594	root filename: ./chlorophyll variable name: CHLA
1595	frequency: -1.00000000000000 time interp: T
1596	climatology: T weights : pairing : data type: yearly
1597	land/sea mask:
1598
1599	tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating
1600	~~~~~~~
1601	Namelist nambbc : set bbc parameters
1602	Apply a geothermal heating at ocean bottom ln_trabbc = F
1603	type of geothermal flux nn_geoflx = 2
1604	Constant geothermal flux value rn_geoflx_cst =
1605	8.640000000000000E-002
1606
1607	*** no geothermal heat flux
1608
1609	tra_bbl_init : bottom boundary layer initialisation
1610	~~~~~~~~~~~~
1611	Namelist nambbl : set bbl parameters
1612	diffusive bbl (=1) or not (=0) nn_bbl_ldf = 1
1613	advective bbl (=1/2) or not (=0) nn_bbl_adv = 0
1614	diffusive bbl coefficient rn_ahtbbl =
1615	1000.00000000000 m2/s
1616	advective bbl coefficient rn_gambbl =
1617	10.0000000000000 s
1618
1619	tra_dmp_init : T and S newtonian relaxation
1620	~~~~~~~
1621	Namelist namtra_dmp : set relaxation parameters
1622	Apply relaxation or not ln_tradmp = F
1623	mixed layer damping option nn_zdmp = 0
1624	Damping file name cn_resto =
1625	resto.nc
1626
1627
1628
1629
1630	tra_adv_init : choice/control of the tracer advection scheme
1631	~~~~~~~~~~~
1632	Namelist namtra_adv : chose a advection scheme for tracers
1633	2nd order advection scheme ln_traadv_cen2 = F
1634	TVD advection scheme ln_traadv_tvd = T
1635	MUSCL advection scheme ln_traadv_muscl = F
1636	MUSCL2 advection scheme ln_traadv_muscl2 = F
1637	UBS advection scheme ln_traadv_ubs = F
1638	QUICKEST advection scheme ln_traadv_qck = F
1639	upstream scheme within muscl ln_traadv_msc_ups = F
1640	TVD advection scheme with zts ln_traadv_tvd_zts = F
1641
1642	TVD scheme is used
1643
1644	tra_adv_mle_init : mixed layer eddy (MLE) advection acting on tracers
1645	~~~~~~~~~~~~~~~~
1646	Namelist namtra_adv_mle : mixed layer eddy advection on tracers
1647	use mixed layer eddy (MLE, i.e. Fox-Kemper param) (T/F) ln_mle =
1648	T
1649	MLE type: =0 standard Fox-Kemper ; =1 new formulation nn_mle =
1650	1
1651	magnitude of the MLE (typical value: 0.06 to 0.08) rn_ce =
1652	6.000000000000000E-002
1653	scale of ML front (ML radius of deformation) (rn_mle=0) rn_lf =
1654	5000.00000000000 m
1655	maximum time scale of MLE (rn_mle=0) rn_time =
1656	172800.000000000 s
1657	reference latitude (degrees) of MLE coef. (rn_mle=1) rn_lat =
1658	20.0000000000000 deg
1659	space interp. of MLD at u-(v-)pts (0=min,1=averaged,2=max) nn_mld_uv =
1660	0
1661	=1 no MLE in case of convection ; =0 always MLE nn_conv =
1662	0
1663	Density difference used to define ML for FK rn_rho_c_mle =
1664	1.000000000000000E-002
1665
1666	Mixed Layer Eddy induced transport added to tracer advection
1667	New formulation
1668
1669	ML buoyancy criteria = 9.558138401559454E-005 m/s2
1670	associated ML density criteria defined in zdfmxl =
1671	1.000000000000000E-002 kg/m3
1672
1673	tra_ldf_init : lateral tracer diffusive operator
1674	~~~~~~~~~~~
1675	Namelist namtra_ldf already read in ldftra module
1676	see ldf_tra_init report for lateral mixing parameters
1677
1678
1679	Rotated laplacian operator
1680
1681	tra:ldf_ano : lateral diffusion acting on the full fields
1682	~~~~~~~~~~~
1683
1684	tra_zdf_init : vertical tracer physics scheme
1685	~~~~~~~~~~~
1686	Implicit (euler backward) scheme
1687
1688	dyn_adv_init : choice/control of the momentum advection scheme
1689	~~~~~~~~~~~
1690	Namelist namdyn_adv : chose a advection formulation & scheme for momentu
1691	m
1692	Vector/flux form (T/F) ln_dynadv_vec = T
1693	= 0 standard scheme ; =1 Hollingsworth scheme nn_dynkeg =
1694	0
1695	2nd order centred advection scheme ln_dynadv_cen2 = F
1696	3rd order UBS advection scheme ln_dynadv_ubs = F
1697	Sub timestepping of vertical advection ln_dynzad_zts = F
1698
1699	vector form : keg + zad + vor is used
1700	with Centered standard keg scheme
1701
1702	dyn_vor_init : vorticity term : read namelist and control the consistency
1703	~~~~~~~~~~~~
1704	Namelist namdyn_vor : choice of the vorticity term scheme
1705	energy conserving scheme ln_dynvor_ene = F
1706	enstrophy conserving scheme ln_dynvor_ens = F
1707	mixed enstrophy/energy conserving scheme ln_dynvor_mix = F
1708	enstrophy and energy conserving scheme ln_dynvor_een = T
1709	enstrophy and energy conserving scheme (old) ln_dynvor_een_old= F
1710
1711	Vector form advection : vorticity = Coriolis + relative vorticity
1712
1713	vorticity scheme : energy and enstrophy conserving scheme
1714
1715	dyn_ldf_init : Choice of the lateral diffusive operator on dynamics
1716	~~~~~~~~~~~
1717	Namelist nam_dynldf : set lateral mixing parameters (type, direction, co
1718	efficients)
1719	laplacian operator ln_dynldf_lap = T
1720	bilaplacian operator ln_dynldf_bilap = F
1721	iso-level ln_dynldf_level = F
1722	horizontal (geopotential) ln_dynldf_hor = T
1723	iso-neutral ln_dynldf_iso = F
1724
1725	laplacian operator
1726
1727	dyn_hpg_init : hydrostatic pressure gradient initialisation
1728	~~~~~~~~~~~~
1729	Namelist namdyn_hpg : choice of hpg scheme
1730	z-coord. - full steps ln_hpg_zco = F
1731	z-coord. - partial steps (interpolation) ln_hpg_zps = F
1732	s-coord. (standard jacobian formulation) ln_hpg_sco = T
1733	s-coord. (standard jacobian formulation) for isf ln_hpg_isf = F
1734	s-coord. (Density Jacobian: Cubic polynomial) ln_hpg_djc = F
1735	s-coord. (Pressure Jacobian: Cubic polynomial) ln_hpg_prj = F
1736	time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = F
1737
1738	dyn_zdf_init : vertical dynamics physics scheme
1739	~~~~~~~~~~~
1740	Implicit (euler backward) scheme
1741
1742	dyn_spg_init : choice of the surface pressure gradient scheme
1743	~~~~~~~~~~~
1744	Explicit free surface lk_dynspg_exp = F
1745	Free surface with time splitting lk_dynspg_ts = F
1746	Filtered free surface cst volume lk_dynspg_flt = T
1747
1748	filtered free surface
1749	file :
1750	solver.stat
1751	open ok
1752	unit = 20
1753	status = REPLACE
1754	form = FORMATTED
1755	access = SEQUENTIAL
1756
1757
1758	solver_init : solver to compute the surface pressure gradient
1759	~~~~~~~~~~~
1760	Namelist namsol : set solver parameters
1761	type of elliptic solver nn_solv = 1
1762	absolute/relative (0/1) precision nn_sol_arp = 0
1763	minimum iterations for solver nn_nmin = 300
1764	maximum iterations for solver nn_nmax = 800
1765	frequency for test nn_nmod = 10
1766	absolute precision of solver rn_eps = 1.000000000000000E-006
1767	absolute precision for SOR solver rn_resmax = 1.000000000000000E-010
1768	optimal coefficient of sor rn_sor = 1.92000000000000
1769
1770	a preconditioned conjugate gradient solver is used
1771
1772	icbini : Namelist namberg ln_icebergs = F , NO icebergs used
1773	~~~~~~~~
1774
1775	===>>> : W A R N I N G
1776	===============
1777
1778	W A R N I N G: end of record or file while reading namelist namsto in configur
1779	ation namelist iostat = -1
1780
1781	sto_par_init : stochastic parameterization
1782	~~~~~~~~~~~~
1783	Namelist namsto : stochastic parameterization
1784	restart stochastic parameters ln_rststo = F
1785	read seed of RNG from restart file ln_rstseed = T
1786	suffix of sto restart name (input) cn_storst_in =
1787	restart_sto
1788	suffix of sto restart name (output) cn_storst_out =
1789	restart_sto
1790	stochastic equation of state ln_sto_eos = F
1791	number of degrees of freedom nn_sto_eos = 1
1792	random walk horz. std (in grid points) rn_eos_stdxy =
1793	1.40000000000000
1794	random walk vert. std (in grid points) rn_eos_stdz =
1795	0.700000000000000
1796	random walk tcor (in timesteps) rn_eos_tcor =
1797	1440.00000000000
1798	order of autoregressive processes nn_eos_ord = 1
1799	passes of Laplacian filter nn_eos_flt = 0
1800	limitation factor rn_eos_lim =
1801	2.00000000000000
1802
1803	stochastic parameterization :
1804
1805	dia_ptr_init : poleward transport and msf initialization
1806	~~~~~~~~~~~~
1807	Namelist namptr : set ptr parameters
1808	Poleward heat & salt transport (T) or not (F) ln_diaptr = T
1809	Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = T
1810
1811	dia_hsb_init
1812	~~~~~~~~
1813	check the heat and salt budgets (T) or not (F) ln_diahsb = F
1814
1815	===>>> : W A R N I N G
1816	===============
1817
1818	W A R N I N G: end of record or file while reading namelist namtrd in configur
1819	ation namelist iostat = -1
1820
1821	trd_init : Momentum/Tracers trends
1822	~~~~~~~~~~
1823	Namelist namtrd : set trends parameters
1824	global domain averaged dyn & tra trends ln_glo_trd = F
1825	U & V trends: 3D output ln_dyn_trd = F
1826	U & V trends: Mixed Layer averaged ln_dyn_mxl = F
1827	T & S trends: 3D output ln_tra_trd = F
1828	T & S trends: Mixed Layer averaged ln_tra_mxl = F
1829	Kinetic Energy trends ln_KE_trd = F
1830	Potential Energy trends ln_PE_trd = F
1831	Barotropic vorticity trends ln_vor_trd = F
1832	frequency of trends diagnostics (glo) nn_trd = 365
1833	Euler time step switch is 0
1834
1835	AAAAAAAA
1836
1837
1838	~~~~~~~ mean fields initialised to instantaneous values
1839
1840	===>>> : W A R N I N G
1841	===============
1842
1843	sbc_blk_core: ln_dm2dc is taking care of the temporal interpolation of daily qs
1844	r
1845	==> We force time interpolation = .false. for qsr
1846
1847	sbc_blk_core : flux formulation for ocean surface boundary condition
1848	~~~~~~~~~~~~
1849	namsbc_core Namelist
1850	list of files and frequency (>0: in hours ; <0 in months)
1851	root filename: ./U_10_MOD_15JUNE2009_unmasked variable name:
1852	U_10_MOD
1853	frequency: 6.00000000000000 time interp: F
1854	climatology: T weights : ./weights_core_bicubic_orca1.nc pairing :
1855	Uwnd data type: yearly land/sea mask:
1856	root filename: ./V_10_MOD_15JUNE2009_unmasked variable name:
1857	V_10_MOD
1858	frequency: 6.00000000000000 time interp: F
1859	climatology: T weights : ./weights_core_bicubic_orca1.nc pairing :
1860	Vwnd data type: yearly land/sea mask:
1861	root filename: ./Q_10_MOD_15JUNE2009_unmasked variable name:
1862	Q_10_MOD
1863	frequency: 6.00000000000000 time interp: F
1864	climatology: T weights : ./weights_core_bilinear_orca1.nc pairing :
1865	data type: yearly land/sea mask:
1866	root filename: ./SWDN_MOD_15JUNE2009_unmasked variable name:
1867	SWDN_MOD
1868	frequency: 24.0000000000000 time interp: F
1869	climatology: T weights : ./weights_core_bilinear_orca1.nc pairing :
1870	data type: yearly land/sea mask:
1871	root filename: ./LWDN_MOD_15JUNE2009_unmasked variable name:
1872	LWDN_MOD
1873	frequency: 24.0000000000000 time interp: T
1874	climatology: T weights : ./weights_core_bilinear_orca1.nc pairing :
1875	data type: yearly land/sea mask:
1876	root filename: ./T_10_MOD_15JUNE2009_unmasked variable name:
1877	T_10_MOD
1878	frequency: 6.00000000000000 time interp: F
1879	climatology: T weights : ./weights_core_bilinear_orca1.nc pairing :
1880	data type: yearly land/sea mask:
1881	root filename: ./TPRECIP_15JUNE2009_unmasked variable name:
1882	TPRECIP
1883	frequency: -1.00000000000000 time interp: T
1884	climatology: T weights : ./weights_core_bilinear_orca1.nc pairing :
1885	data type: yearly land/sea mask:
1886	root filename: ./SNOW_15JUNE2009_unmasked variable name: SNOW
1887	frequency: -1.00000000000000 time interp: T
1888	climatology: T weights : ./weights_core_bilinear_orca1.nc pairing :
1889	data type: yearly land/sea mask:
1890	iom_nf90_open ~~~ open existing file: ./LWDN_MOD_15JUNE2009
1891	_unmasked.nc in READ mode
1892	---> ./LWDN_MOD_15JUNE2009_unmasked.nc OK
1893	iom_nf90_open ~~~ open existing file: ./LWDN_MOD_15JUNE2009
1894	_unmasked.nc in READ mode
1895	---> ./LWDN_MOD_15JUNE2009_unmasked.nc OK
1896	iom_close ~~~ close file: ./LWDN_MOD_15JUNE2009_unmasked.nc
1897	ok
1898	iom_nf90_open ~~~ open existing file: ./weights_core_biline
1899	ar_orca1.nc in READ mode
1900	---> ./weights_core_bilinear_orca1.nc OK
1901
1902	===>>> : W A R N I N G
1903	===============
1904
1905	iom_nf90_getatt: no attribute ew_wrap found
1906	read src01 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1907	read src02 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1908	read src03 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1909	read src04 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1910	read wgt01 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1911	read wgt02 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1912	read wgt03 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1913	read wgt04 (rec: 1) in ./weights_core_bilinear_orca1.nc ok
1914	iom_close ~~~ close file: ./weights_core_bilinear_orca1.nc
1915	ok
1916	read LWDN_MOD (rec: 365) in ./LWDN_MOD_15JUNE2009_unmasked.nc ok
1917	fld_init : time-interpolation for LWDN_MOD read previous record = 365 at time = -0.50 days
1918	iom_nf90_open ~~~ open existing file: ./TPRECIP_15JUNE2009_
1919	unmasked.nc in READ mode
1920	---> ./TPRECIP_15JUNE2009_unmasked.nc OK
1921	read TPRECIP (rec: 12) in ./TPRECIP_15JUNE2009_unmasked.nc ok
1922	fld_init : time-interpolation for TPRECIP read previous record = 12 at time = -15.50 days
1923	iom_nf90_open ~~~ open existing file: ./SNOW_15JUNE2009_unm
1924	asked.nc in READ mode
1925	---> ./SNOW_15JUNE2009_unmasked.nc OK
1926	read SNOW (rec: 12) in ./SNOW_15JUNE2009_unmasked.nc ok
1927	fld_init : time-interpolation for SNOW read previous record = 12 at time = -15.50 days
1928	weight file: ./weights_core_bilinear_orca1.nc
1929	ddims: 192 94
1930	numwgt: 4
1931	jpiwgt: 19
1932	jpjwgt: 7
1933	botleft: 39 6
1934	topright: 57 12
1935	cyclical
1936	allocated
1937	iom_nf90_open ~~~ open existing file: ./U_10_MOD_15JUNE2009
1938	_unmasked.nc in READ mode
1939	---> ./U_10_MOD_15JUNE2009_unmasked.nc OK
1940	iom_nf90_open ~~~ open existing file: ./U_10_MOD_15JUNE2009
1941	_unmasked.nc in READ mode
1942	---> ./U_10_MOD_15JUNE2009_unmasked.nc OK
1943	iom_close ~~~ close file: ./U_10_MOD_15JUNE2009_unmasked.nc
1944	ok
1945	iom_nf90_open ~~~ open existing file: ./weights_core_bicubi
1946	c_orca1.nc in READ mode
1947	---> ./weights_core_bicubic_orca1.nc OK
1948
1949	===>>> : W A R N I N G
1950	===============
1951
1952	iom_nf90_getatt: no attribute ew_wrap found
1953	read src01 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1954	read src02 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1955	read src03 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1956	read src04 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1957	read wgt01 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1958	read wgt02 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1959	read wgt03 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1960	read wgt04 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1961	read wgt05 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1962	read wgt06 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1963	read wgt07 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1964	read wgt08 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1965	read wgt09 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1966	read wgt10 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1967	read wgt11 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1968	read wgt12 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1969	read wgt13 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1970	read wgt14 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1971	read wgt15 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1972	read wgt16 (rec: 1) in ./weights_core_bicubic_orca1.nc ok
1973	iom_close ~~~ close file: ./weights_core_bicubic_orca1.nc o
1974	k
1975	read U_10_MOD (rec: 1) in ./U_10_MOD_15JUNE2009_unmasked.nc ok
1976	iom_nf90_open ~~~ open existing file: ./V_10_MOD_15JUNE2009
1977	_unmasked.nc in READ mode
1978	---> ./V_10_MOD_15JUNE2009_unmasked.nc OK
1979	read V_10_MOD (rec: 1) in ./V_10_MOD_15JUNE2009_unmasked.nc ok
1980	iom_nf90_open ~~~ open existing file: ./Q_10_MOD_15JUNE2009
1981	_unmasked.nc in READ mode
1982	---> ./Q_10_MOD_15JUNE2009_unmasked.nc OK
1983	read Q_10_MOD (rec: 1) in ./Q_10_MOD_15JUNE2009_unmasked.nc ok
1984	iom_nf90_open ~~~ open existing file: ./SWDN_MOD_15JUNE2009
1985	_unmasked.nc in READ mode
1986	---> ./SWDN_MOD_15JUNE2009_unmasked.nc OK
1987	read SWDN_MOD (rec: 1) in ./SWDN_MOD_15JUNE2009_unmasked.nc ok
1988	read LWDN_MOD (rec: 1) in ./LWDN_MOD_15JUNE2009_unmasked.nc ok
1989	iom_nf90_open ~~~ open existing file: ./T_10_MOD_15JUNE2009
1990	_unmasked.nc in READ mode
1991	---> ./T_10_MOD_15JUNE2009_unmasked.nc OK
1992	read T_10_MOD (rec: 1) in ./T_10_MOD_15JUNE2009_unmasked.nc ok
1993	read TPRECIP (rec: 1) in ./TPRECIP_15JUNE2009_unmasked.nc ok
1994	read SNOW (rec: 1) in ./SNOW_15JUNE2009_unmasked.nc ok
1995
1996	rot_rep : geographic <--> stretched
1997	~~~~~ coordinate transformation
1998	fld_read: vector pair (U_10_MOD, V_10_MOD) rotated on to model grid
1999	fld_read: var U_10_MOD kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2000	fld_read: var V_10_MOD kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2001	fld_read: var Q_10_MOD kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2002	fld_read: var SWDN_MOD kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 1.0000)
2003	fld_read: var LWDN_MOD kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, records b/a: 0365/ 0001 (days -0.5000/ 0.5000)
2004	it_offset is : 0
2005	fld_read: var T_10_MOD kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2006	fld_read: var TPRECIP kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2007	it_offset is : 0
2008	fld_read: var SNOW kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2009	it_offset is : 0
2010
2011	sbc_dcy : introduce diurnal cycle from daily mean qsr
2012	~~~~~~~
2013
2014
2015	sbc_ice_lim_2 : update ocean surface boudary condition
2016	~~~~~~~~~~~~~ via Louvain la Neuve Ice Model (LIM) time stepping
2017
2018	ice_init_2 : LIM-2 sea-ice - initialization
2019	~~~~~~~~~~~
2020	file :
2021	namelist_ice_ref
2022	open ok
2023	unit = 21
2024	status = OLD
2025	form = FORMATTED
2026	access = SEQUENTIAL
2027
2028	file :
2029	namelist_ice_cfg
2030	open ok
2031	unit = 22
2032	status = OLD
2033	form = FORMATTED
2034	access = SEQUENTIAL
2035
2036	file :
2037	output.namelist.ice
2038	open ok
2039	unit = 23
2040	status = UNKNOWN
2041	form = FORMATTED
2042	access = SEQUENTIAL
2043
2044
2045	ice_run : ice share parameters for dynamics/advection/thermo of sea-ice
2046	~~~~~~
2047	switch for ice dynamics (1) or not (0) ln_limdyn = T
2048	Ice damping ln_limdmp = F
2049	minimum fraction for leads in the NH (SH) acrit(1/2) =
2050	1.000000000000000E-006 1.000000000000000E-006
2051	computation of temp. in snow (=0) or not (=9999) hsndif =
2052	0.000000000000000E+000
2053	computation of temp. in ice (=0) or not (=9999) hicdif =
2054	0.000000000000000E+000
2055
2056	lim_msh_2 : LIM 2.0 sea-ice model, mesh initialization
2057	~~~~~~~~~
2058	the model domain is entirely in the southern hemisphere: njeq =
2059	39
2060
2061	lim_istate_init_2 : ice parameters inititialisation
2062	~~~~~~~~~~~~~~~~~
2063	threshold water temp. for initial sea-ice ttest =
2064	2.00000000000000
2065	initial snow thickness in the north hninn =
2066	0.500000000000000
2067	initial ice thickness in the north hginn =
2068	3.00000000000000
2069	initial leads area in the north alinn =
2070	5.000000000000000E-002
2071	initial snow thickness in the south hnins =
2072	0.100000000000000
2073	initial ice thickness in the south hgins =
2074	1.00000000000000
2075	initial leads area in the south alins =
2076	0.100000000000000
2077	Ice state initialization using input file ln_limini = F
2078
2079	lim_sbc_init_2 : LIM-2 sea-ice - surface boundary condition
2080	~~~~~~~~~~~~~~~
2081	EVP rheology - C-grid case
2082
2083	albedo : set albedo parameters
2084	~~~~~~~
2085	Namelist namsbc_alb : albedo
2086	choose the albedo parameterization nn_ice_alb =
2087	1
2088	albedo of dry snow rn_alb_sdry =
2089	0.850000000000000
2090	albedo of melting snow rn_alb_smlt =
2091	0.750000000000000
2092	albedo of dry ice rn_alb_idry =
2093	0.600000000000000
2094	albedo of bare puddled ice rn_alb_imlt =
2095	0.500000000000000
2096
2097	lim_dyn_init_2: ice parameters for ice dynamics
2098	~~~~~~~~~~~~~~
2099	tolerance parameter epsd =
2100	9.999999999999999E-021
2101	coefficient for semi-implicit coriolis alpha =
2102	0.500000000000000
2103	diffusion constant for dynamics dm =
2104	600.000000000000
2105	number of sub-time steps for relaxation nbiter = 1
2106	maximum number of iterations for relaxation nbitdr = 100
2107	relaxation constant om =
2108	0.500000000000000
2109	maximum value for the residual of relaxation resl =
2110	5.000000000000000E-005
2111	drag coefficient for oceanic stress cw =
2112	5.000000000000000E-003
2113	turning angle for oceanic stress angvg =
2114	0.000000000000000E+000 degrees
2115	first bulk-rheology parameter pstar =
2116	10000.0000000000
2117	second bulk-rhelogy parameter c_rhg =
2118	20.0000000000000
2119	minimun value for viscosity etamn =
2120	0.000000000000000E+000
2121	creep limit rn_creepl =
2122	1.000000000000000E-008
2123	eccentricity of the elliptical yield curve rn_ecc =
2124	2.00000000000000
2125	horizontal diffusivity coeff. for sea-ice ahi0 =
2126	350.000000000000
2127	number of iterations for subcycling nn_nevp= 120
2128	timescale for elastic waves telast = 9600
2129	coefficient for the solution of int. stresses alphaevp =
2130	1.00000000000000
2131
2132	lim_trp_init_2 : Ice parameters for advection
2133	~~~~~~~~~~~~~~
2134	boundary conditions (0. no-slip, 1. free-slip) bound =
2135	0.000000000000000E+000
2136
2137	lim_thd_init_2: ice parameters for ice thermodynamic computation
2138	~~~~~~~~~~~~~~
2139	maximum melting at the bottom hmelt =
2140	-0.150000000000000
2141	ice thick. for lateral accretion in NH (SH) hiccrit(1/2) =
2142	0.300000000000000 0.300000000000000
2143	ice thick. corr. to max. energy stored in brine pocket hicmin =
2144	0.200000000000000
2145	minimum ice thickness hiclim =
2146	5.000000000000000E-002
2147	maximum lead fraction amax =
2148	0.999000000000000
2149	energy stored in brine pocket (=1) or not (=0) swiqst =
2150	1.00000000000000
2151	numerical carac. of the scheme for diffusion in ice
2152	Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) sbeta =
2153	1.00000000000000
2154	percentage of energy used for lateral ablation parlat =
2155	0.000000000000000E+000
2156	slope of distr. for Hakkinen-Mellor lateral melting hakspl =
2157	0.500000000000000
2158	slope of distribution for Hibler lateral melting hibspl =
2159	0.500000000000000
2160	exponent for leads-closure rate exld =
2161	2.00000000000000
2162	coefficient for diffusions of ice and snow hakdif =
2163	1.00000000000000
2164	threshold thick. for comp. of eq. thermal conductivity zhth =
2165	0.200000000000000
2166	thickness of the surf. layer in temp. computation hnzst =
2167	0.100000000000000
2168	switch for snow sublimation (=1) or not (=0) parsub =
2169	1.00000000000000
2170	coefficient for snow density when snow ice formation alphs =
2171	1.00000000000000
2172
2173	===>>> : W A R N I N G
2174	===============
2175
2176	W A R N I N G: end of record or file while reading namelist namiceout in confi
2177	guration namelist iostat = -1
2178
2179	lim_wri_init_2 : Ice parameters for outputs
2180	~~~~~~~~~~~~~~
2181	number of fields to be stored noumef = 19
2182	title name unit Saving (1/0)
2183	multiplicative constant additive constant
2184	Snow thickness isnowthi m 1
2185	1.00000000000000 0.000000000000000E+000
2186	Ice thickness iicethic m 1
2187	1.00000000000000 0.000000000000000E+000
2188	Ice produced iiceprod m/kt 1
2189	1.00000000000000 0.000000000000000E+000
2190	Ice concentration ileadfra % 1
2191	-1.00000000000000 1.00000000000000
2192	Ice temperature iicetemp C 1
2193	1.00000000000000 -273.150000000000
2194	Oceanic flux at the ice base ioceflxb w/m2 1
2195	1.00000000000000 0.000000000000000E+000
2196	Ice velocity u iicevelu m/s 1
2197	1.00000000000000 0.000000000000000E+000
2198	Ice velocity v iicevelv m/s 1
2199	1.00000000000000 0.000000000000000E+000
2200	Sea surface temperature isstempe C 1
2201	1.00000000000000 -273.150000000000
2202	Sea surface salinity isssalin PSU 1
2203	1.00000000000000 0.000000000000000E+000
2204	Total flux at ocean surface iocetflx w/m2 1
2205	1.00000000000000 0.000000000000000E+000
2206	Solar flux at ocean surface iocesflx w/m2 1
2207	1.00000000000000 0.000000000000000E+000
2208	Non-solar flux at ocean surface iocwnsfl w/m2 1
2209	1.00000000000000 0.000000000000000E+000
2210	Salt flux at ocean surface iocesafl kg/m2/kt 1
2211	1.00000000000000 0.000000000000000E+000
2212	Wind stress u iocestru Pa 1
2213	1.00000000000000 0.000000000000000E+000
2214	Wind stress v iocestrv Pa 1
2215	1.00000000000000 0.000000000000000E+000
2216	Solar flux at ice/ocean surface iicesflx w/m2 1
2217	1.00000000000000 0.000000000000000E+000
2218	Non-solar flux at ice/ocean surface iicenflx w/m2 1
2219	1.00000000000000 0.000000000000000E+000
2220	Snow precipitation isnowpre kg/day 1
2221	1.00000000000000 0.000000000000000E+000
2222	iom_nf90_open ~~~ open existing file: ./runoff_1m_ORCA1.nc
2223	in READ mode
2224	---> ./runoff_1m_ORCA1.nc OK
2225	read sorunoff (rec: 12) in ./runoff_1m_ORCA1.nc ok
2226	fld_init : time-interpolation for sorunoff read previous record = 12 at time = -15.50 days
2227	weight file: ./weights_core_bilinear_orca1.nc
2228	ddims: 192 94
2229	numwgt: 4
2230	jpiwgt: 19
2231	jpjwgt: 7
2232	botleft: 39 6
2233	topright: 57 12
2234	cyclical
2235	allocated
2236	weight file: ./weights_core_bicubic_orca1.nc
2237	ddims: 192 94
2238	numwgt: 16
2239	jpiwgt: 19
2240	jpjwgt: 7
2241	botleft: 39 6
2242	topright: 57 12
2243	cyclical
2244	allocated
2245	read sorunoff (rec: 1) in ./runoff_1m_ORCA1.nc ok
2246	fld_read: var sorunoff kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2247	it_offset is : 0
2248	nit000-1 runoff forcing fields set to nit000
2249	iom_nf90_open ~~~ open existing file: ./sss_data.nc in READ
2250	mode
2251	---> ./sss_data.nc OK
2252	read sss (rec: 12) in ./sss_data.nc ok
2253	fld_init : time-interpolation for sss read previous record = 12 at time = -15.50 days
2254	weight file: ./weights_core_bilinear_orca1.nc
2255	ddims: 192 94
2256	numwgt: 4
2257	jpiwgt: 19
2258	jpjwgt: 7
2259	botleft: 39 6
2260	topright: 57 12
2261	cyclical
2262	allocated
2263	weight file: ./weights_core_bicubic_orca1.nc
2264	ddims: 192 94
2265	numwgt: 16
2266	jpiwgt: 19
2267	jpjwgt: 7
2268	botleft: 39 6
2269	topright: 57 12
2270	cyclical
2271	allocated
2272	read sss (rec: 1) in ./sss_data.nc ok
2273	fld_read: var sss kt = 1 ( 0.0625 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2274	it_offset is : 0
2275
2276	sbc_fwb : FreshWater Budget correction
2277	~~~~~~~
2278	instantaneously set to zero
2279	nit000-1 surface forcing fields set to nit000
2280
2281	zdf_bfr : Set bottom friction coefficient (non-linear case)
2282	~~~~~~~~
2283
2284	zdf_evd : Enhanced Vertical Diffusion (evd)
2285	~~~~~~~
2286
2287
2288	N Total power consumption by av_tide : ztpc =
2289	3.027013904413500E-009 TW
2290	N Total power consumption by zavt_itf: ztpc =
2291	3.027013904413500E-009 TW
2292
2293	zdf_mxl : mixed layer depth
2294	~~~~~~~
2295
2296	ldf_eiv : eddy induced velocity coefficients
2297	~~~~~~~
2298
2299	ssh_nxt : after sea surface height
2300	~~~~~~~
2301
2302	div_cur : horizontal velocity divergence and
2303	~~~~~~~ relative vorticity
2304
2305	wzv : now vertical velocity
2306	~~~~~
2307
2308	Date 0 used : 1 YEAR 1 MONTH 1 DAY
2309	1 Julian day : 365.000000000000
2310	indexes of zoom = 1 32 1 39
2311	limit storage in depth = 75
2312	file :
2313	date.file
2314	open ok
2315	unit = 24
2316	status = REPLACE
2317	form = FORMATTED
2318	access = SEQUENTIAL
2319
2320	Name of NETCDF file ORCA1_1y_00010101_00011231_grid_T
2321	Name of NETCDF file ORCA1_1y_00010101_00011231_grid_U
2322	Name of NETCDF file ORCA1_1y_00010101_00011231_grid_V
2323	Name of NETCDF file ORCA1_1y_00010101_00011231_grid_W
2324
2325	End of NetCDF Initialization
2326
2327	tra_sbc : TRAcer Surface Boundary Condition
2328	~~~~~~~
2329
2330	tra_qsr : penetration of the surface solar radiation
2331	~~~~~~~
2332	iom_nf90_open ~~~ open existing file: ./chlorophyll.nc in R
2333	EAD mode
2334	---> ./chlorophyll.nc OK
2335	read CHLA (rec: 12) in ./chlorophyll.nc ok
2336	fld_init : time-interpolation for CHLA read previous record = 12 at time = -15.50 days
2337	weight file: ./weights_core_bilinear_orca1.nc
2338	ddims: 192 94
2339	numwgt: 4
2340	jpiwgt: 19
2341	jpjwgt: 7
2342	botleft: 39 6
2343	topright: 57 12
2344	cyclical
2345	allocated
2346	weight file: ./weights_core_bicubic_orca1.nc
2347	ddims: 192 94
2348	numwgt: 16
2349	jpiwgt: 19
2350	jpjwgt: 7
2351	botleft: 39 6
2352	topright: 57 12
2353	cyclical
2354	allocated
2355	read CHLA (rec: 1) in ./chlorophyll.nc ok
2356	fld_read: var CHLA kt = 1 ( 0.0208 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2357	it_offset is : 0
2358
2359	trabbl:bbl : Compute bbl velocities and diffusive coefficients in TRA
2360	~~~~~~~~~~
2361
2362	tra_adv_eiv : eddy induced advection on TRA :
2363	~~~~~~~~~~~ add to velocity fields the eiv component
2364
2365	tra_adv_tvd : TVD advection scheme on TRA
2366	~~~~~~~~~~~
2367
2368	tra_ldf_iso : rotated laplacian diffusion operator on TRA
2369	~~~~~~~~~~~
2370
2371	tra_zdf_imp : implicit vertical mixing on TRA
2372	~~~~~~~~~~~
2373
2374	tra_nxt : achieve the time stepping by Asselin filter and array swap
2375	~~~~~~~
2376
2377	dyn_keg : kinetic energy gradient trend, scheme number= 0
2378	~~~~~~~
2379
2380	dyn_zad : arakawa advection scheme
2381
2382	dyn:vor_een : vorticity term: energy and enstrophy conserving scheme
2383	~~~~~~~~~~~
2384
2385	dyn_ldf : iso-level harmonic (laplacian) operator
2386	~~~~~~~
2387
2388	dyn:hpg_sco : hydrostatic pressure gradient trend
2389	~~~~~~~~~~~ s-coordinate case, OPA original scheme used
2390
2391	dyn_zdf_imp : vertical momentum diffusion implicit operator
2392	~~~~~~~~~~~
2393
2394	dyn_spg_flt : surface pressure gradient trend
2395	~~~~~~~~~~~ (free surface constant volume case)
2396
2397	dyn_nxt : time stepping
2398	~~~~~~~
2399
2400	ssh_swp : Asselin time filter and swap of sea surface height
2401	~~~~~~~
2402
2403	stp_ctl : time-stepping control
2404	~~~~~~~
2405	file :
2406	time.step
2407	open ok
2408	unit = 24
2409	status = REPLACE
2410	form = FORMATTED
2411	access = SEQUENTIAL
2412
2413	==>> time-step= 1 abs(U) max: 0.398008369514494
2414	==>> time-step= 1 SSS min: 5.28036190306677
2415	fld_read: var CHLA kt = 2 ( 0.0625 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2416	it_offset is : 0
2417	fld_read: var CHLA kt = 3 ( 0.1042 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2418	it_offset is : 0
2419	fld_read: var U_10_MOD kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2420	fld_read: var V_10_MOD kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2421	fld_read: var Q_10_MOD kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2422	fld_read: var SWDN_MOD kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 1.0000)
2423	fld_read: var LWDN_MOD kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, records b/a: 0365/ 0001 (days -0.5000/ 0.5000)
2424	it_offset is : 0
2425	fld_read: var T_10_MOD kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 0.2500)
2426	fld_read: var TPRECIP kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2427	it_offset is : 0
2428	fld_read: var SNOW kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2429	it_offset is : 0
2430	fld_read: var sorunoff kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2431	it_offset is : 0
2432	fld_read: var sss kt = 4 ( 0.1875 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2433	it_offset is : 0
2434	fld_read: var CHLA kt = 4 ( 0.1458 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2435	it_offset is : 0
2436	fld_read: var CHLA kt = 5 ( 0.1875 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2437	it_offset is : 0
2438	fld_read: var CHLA kt = 6 ( 0.2292 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2439	it_offset is : 0
2440	read U_10_MOD (rec: 2) in ./U_10_MOD_15JUNE2009_unmasked.nc ok
2441	read V_10_MOD (rec: 2) in ./V_10_MOD_15JUNE2009_unmasked.nc ok
2442	read Q_10_MOD (rec: 2) in ./Q_10_MOD_15JUNE2009_unmasked.nc ok
2443	read T_10_MOD (rec: 2) in ./T_10_MOD_15JUNE2009_unmasked.nc ok
2444	fld_read: var U_10_MOD kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, record: 0002 (days 0.2500 <-> 0.5000)
2445	fld_read: var V_10_MOD kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, record: 0002 (days 0.2500 <-> 0.5000)
2446	fld_read: var Q_10_MOD kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, record: 0002 (days 0.2500 <-> 0.5000)
2447	fld_read: var SWDN_MOD kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, record: 0001 (days 0.0000 <-> 1.0000)
2448	fld_read: var LWDN_MOD kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, records b/a: 0365/ 0001 (days -0.5000/ 0.5000)
2449	it_offset is : 0
2450	fld_read: var T_10_MOD kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, record: 0002 (days 0.2500 <-> 0.5000)
2451	fld_read: var TPRECIP kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2452	it_offset is : 0
2453	fld_read: var SNOW kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2454	it_offset is : 0
2455	lim_trp_2 : violation of cfl criterion the 1 th day, cfl =
2456	0.556601947087279
2457	fld_read: var sorunoff kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2458	it_offset is : 0
2459	fld_read: var sss kt = 7 ( 0.3125 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2460	it_offset is : 0
2461	fld_read: var CHLA kt = 7 ( 0.2708 days), Y/M/D = 0001/01/01, records b/a: 0012/ 0001 (days -15.5000/ 15.5000)
2462	it_offset is : 0
2463
2464	===>>> : E R R O R
2465	===========
2466
2467	stpctl: the zonal velocity is larger than 20 m/s
2468	======
2469	kt= 7 max abs(U): 42.49 , i j k: 77 291 17
2470
2471	output of last fields in numwso
2472
2473	===>>> : E R R O R
2474	===========
2475
2476	step: indic < 0
2477
2478	dia_wri_state : single instantaneous ocean state
2479	~~~~~~~~~~~~~ and forcing fields file created
2480	and named :output.abort .nc
2481
2482	AAAAAAAA
2483
2484
2485	===>>> : E R R O R
2486	===========
2487
2488	1 error have been found
2489	iom_close ~~~ close file: ./data_1m_potential_temperature_n
2490	omask.nc ok
2491	iom_close ~~~ close file: ./data_1m_salinity_nomask.nc ok
2492	iom_close ~~~ close file: ./LWDN_MOD_15JUNE2009_unmasked.nc
2493	ok
2494	iom_close ~~~ close file: ./TPRECIP_15JUNE2009_unmasked.nc
2495	ok
2496	iom_close ~~~ close file: ./SNOW_15JUNE2009_unmasked.nc ok
2497	iom_close ~~~ close file: ./U_10_MOD_15JUNE2009_unmasked.nc
2498	ok
2499	iom_close ~~~ close file: ./V_10_MOD_15JUNE2009_unmasked.nc
2500	ok
2501	iom_close ~~~ close file: ./Q_10_MOD_15JUNE2009_unmasked.nc
2502	ok
2503	iom_close ~~~ close file: ./SWDN_MOD_15JUNE2009_unmasked.nc
2504	ok
2505	iom_close ~~~ close file: ./T_10_MOD_15JUNE2009_unmasked.nc
2506	ok
2507	iom_close ~~~ close file: ./runoff_1m_ORCA1.nc ok
2508	iom_close ~~~ close file: ./sss_data.nc ok
2509	iom_close ~~~ close file: ./chlorophyll.nc ok

Download in other formats:

Original Format

New URL for NEMO forge! http://forge.nemo-ocean.eu

Context Navigation

Topic #27: ocean.output_v3.6

Download in other formats: