Topic #32: ocean.output

File ocean.output, 83.8 KB (added by mbareford, 2 years ago)

NEMO v3.6 r9301 ORCA2_LIM

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 =  F         
21 
22             Immediate non-blocking send (isend)                               
23 
24
25AAAAAAAA
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      =            2
48       1st lateral dimension ( >= jpi ) jpidta     =          182
49       2nd    "         "    ( >= jpj ) jpjdta     =          149
50       3nd    "         "               jpkdta     =           31
51       1st dimension of global domain in i jpiglo  =          182
52       2nd    -                  -    in j jpjglo  =          149
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=           6  irestj=           3
65 
66 jpni=           6  jpnj=           4
67 
68  sum ilcit(i,1)=   182.000000000000       jpiglo=         182
69  sum ilcit(1,j)=   149.000000000000       jpjglo=         149
70 
71 
72     ****************************************************************
73     *              *               *               *               *   
74   4 *   32  x 38   *    32  x 38   *    32  x 38   *    32  x 38   *   
75     *         18   *          19   *          20   *          21   *         
76     *              *               *               *               *   
77     ****************************************************************
78     *              *               *               *               *   
79   3 *   32  x 39   *    32  x 39   *    32  x 39   *    32  x 39   *   
80     *         12   *          13   *          14   *          15   *         
81     *              *               *               *               *   
82     ****************************************************************
83     *              *               *               *               *   
84   2 *   32  x 39   *    32  x 39   *    32  x 39   *    32  x 39   *   
85     *          6   *           7   *           8   *           9   *         
86     *              *               *               *               *   
87     ****************************************************************
88     *              *               *               *               *   
89   1 *   32  x 39   *    32  x 39   *    32  x 39   *    32  x 39   *   
90     *          0   *           1   *           2   *           3   *         
91     *              *               *               *               *   
92     ****************************************************************
93             1               2               3               4             
94 
95     ********************************
96     *              *               *   
97   4 *   32  x 38   *    32  x 38   *   
98     *         22   *          23   *         
99     *              *               *   
100     ********************************
101     *              *               *   
102   3 *   32  x 39   *    32  x 39   *   
103     *         16   *          17   *         
104     *              *               *   
105     ********************************
106     *              *               *   
107   2 *   32  x 39   *    32  x 39   *   
108     *         10   *          11   *         
109     *              *               *   
110     ********************************
111     *              *               *   
112   1 *   32  x 39   *    32  x 39   *   
113     *          4   *           5   *         
114     *              *               *   
115     ********************************
116             5               6             
117  nproc  =            0
118  nowe   =            5  noea   =             1
119  nono   =            6  noso   =            -6
120  nbondi =            0
121  nbondj =           -1
122  npolj  =            0
123  nperio =            0
124  nlci   =           32
125  nlcj   =           39
126  nimpp  =            1
127  njmpp  =            1
128  nreci  =            2  npse   =           -5
129  nrecj  =            2  npsw   =           -1
130  jpreci =            1  npne   =            7
131  jprecj =            1  npnw   =           11
132 
133  mpp_init : North fold boundary prepared for jpni >1
134 
135 mpp_init_ioipsl :   iloc  =           32          39
136 ~~~~~~~~~~~~~~~     iabsf =            1           1
137                     ihals =            0           0
138                     ihale =            1           1
139 
140  phy_cst : initialization of ocean parameters and constants
141  ~~~~~~~
142        Domain info
143           dimension of model
144                  Local domain      Global domain       Data domain
145               jpi     :      32   jpiglo  :     182   jpidta  :     182
146               jpj     :      39   jpjglo  :     149   jpjdta  :     149
147               jpk     :      31   jpk     :      31   jpkdta  :      31
148               jpij    :         1248
149           mpp local domain info (mpp)
150              jpni    :            6    jpreci  :            1
151              jpnj    :            4    jprecj  :            1
152              jpnij   :           24
153           lateral domain boundary condition type : jperio  =            4
154 
155        Constants
156 
157           mathematical constant                 rpi =    3.14159265358979     
158 
159           day                                rday   =    86400.0000000000     
160  s
161           sideral year                       rsiyea =    31558149.0101107     
162  s
163           sideral day                        rsiday =    86164.0996559118     
164  s
165           omega                              omega  =   7.292115083046062E-005
166  s^-1
167 
168           nb of months per year               raamo =    12.0000000000000     
169  months
170           nb of hours per day                 rjjhh =    24.0000000000000     
171  hours
172           nb of minutes per hour              rhhmm =    60.0000000000000     
173  mn
174           nb of seconds per minute            rmmss =    60.0000000000000     
175  s
176 
177           earth radius                         ra   =    6371229.00000000     
178  m
179           gravity                              grav =    9.80665000000000     
180  m/s^2
181 
182           triple point of temperature      rtt      =    273.160000000000     
183  K
184           freezing point of water          rt0      =    273.150000000000     
185  K
186           melting point of snow            rt0_snow =    273.150000000000     
187  K
188           melting point of ice             rt0_ice  =    273.050000000000     
189  K
190           reference density and heat capacity now defined in eosbn2.f90
191 
192           thermal conductivity of pure ice          =    2.03439600000000     
193  J/s/m/K
194           fresh ice specific heat                   =    2093.00000000000     
195  J/kg/K
196           latent heat of fusion of fresh ice / snow =    333700.000000000     
197  J/kg
198           density times specific heat for snow      =    690690.000000000     
199  J/m^3/K
200           density times specific heat for ice       =    1883700.00000000     
201  J/m^3/K
202           volumetric latent heat fusion of sea ice  =    300330000.000000     
203  J/m
204           latent heat of sublimation of snow        =    2800000.00000000     
205  J/kg
206           volumetric latent heat fusion of snow     =    110121000.000000     
207  J/m^3
208           density of sea ice                        =    900.000000000000     
209  kg/m^3
210           density of snow                           =    330.000000000000     
211  kg/m^3
212           emissivity of snow or ice                 =   0.970000000000000     
213           salinity of ice                           =    6.00000000000000     
214  psu
215           salinity of sea                           =    34.7000000000000     
216  psu
217           latent heat of evaporation (water)        =    2500000.00000000     
218  J/m^3
219           correction factor for solar radiation     =   0.900000000000000     
220           von Karman constant                       =   0.400000000000000     
221           Stefan-Boltzmann constant                 =   5.670000000000000E-008
222  J/s/m^2/K^4
223 
224           conversion: degre ==> radian          rad =   1.745329251994330E-002
225 
226           smallest real computer value       rsmall =   1.110223024625157E-016
227 
228 eos_init : equation of state
229 ~~~~~~~~
230           Namelist nameos : set eos parameters
231              flag for eq. of state and N^2  nn_eos   =           -1
232              model uses Conservative Temperature
233              Important: model must be initialized with CT and SA fields
234 
235           use of TEOS-10 equation of state (cons. temp. and abs. salinity)
236 
237           volumic mass of reference           rau0  =    1026.00000000000     
238  kg/m^3
239           1. / rau0                        r1_rau0  =   9.746588693957114E-004
240  m^3/kg
241           ocean specific heat                 rcp   =    3991.86795711963     
242  J/Kelvin
243           rau0 * rcp                       rau0_rcp =    4095656.52400474     
244           1. / ( rau0 * rcp )           r1_rau0_rcp =   2.441610994816036E-007
245 
246 dom_cfg : set the ocean configuration
247 ~~~~~~~
248    ocean model configuration used :   cp_cfg =
249 orca                                                                           
250                                                                               
251                                                                               
252                     jp_cfg =            2
253    global domain lateral boundaries
254       jperio= 4, cyclic east-west and north fold with T-point pivot
255 
256 dom_glo : domain: data / local
257 ~~~~~~~
258           data input domain    : jpidta =          182  jpjdta =          149
259  jpkdta =           31
260           global or zoom domain: jpiglo =          182  jpjglo =          149
261  jpk    =           31
262           local domain         : jpi    =           32  jpj    =           39
263  jpk    =           31
264 
265           south-west indices    jpizoom =            1  jpjzoom =            1
266 
267           conversion local  ==> data i-index domain
268             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
269            20  21  22  23  24  25  26  27  28  29  30  31  32
270 
271           conversion data   ==> local  i-index domain
272              starting index
273             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
274            20  21  22  23  24  25  26  27  28  29  30  31  32  33  33  33  33  33  33
275            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
276            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
277            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
278            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
279            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
280            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
281            33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33  33
282            33  33  33  33  33  33  33  33  33  33  33
283              ending index
284             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
285            20  21  22  23  24  25  26  27  28  29  30  31  32  32  32  32  32  32  32
286            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
287            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
288            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
289            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
290            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
291            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
292            32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32  32
293            32  32  32  32  32  32  32  32  32  32  32
294 
295           conversion local  ==> data j-index domain
296             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
297            20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38
298            39
299 
300           conversion data  ==> local j-index domain
301              starting index
302             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
303            20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38
304            39  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40
305            40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40
306            40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40
307            40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40
308            40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40
309            40  40  40  40  40  40  40  40  40  40  40  40  40  40  40  40
310              ending index
311             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
312            20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38
313            39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39
314            39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39
315            39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39
316            39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39
317            39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39
318            39  39  39  39  39  39  39  39  39  39  39  39  39  39  39  39
319 
320           zoom flags :
321              lzoom   =  F  (T = zoom, F = global )
322              lzoom_e =  F  (T = forced closed east  boundary)
323              lzoom_w =  F  (T = forced closed west  boundary)
324              lzoom_s =  F  (T = forced closed South boundary)
325              lzoom_n =  F  (T = forced closed North boundary)
326 
327 dom_init : domain initialization
328 ~~~~~~~~
329 
330 dom_nam  : domain initialization through namelist read
331 ~~~~~~~
332    Namelist namrun
333       job number                      nn_no      =            0
334       experiment name for output      cn_exp     =
335 ORCA2                                                                         
336                                                                               
337                                                                               
338                   
339       file prefix restart input       cn_ocerst_in=
340 restart                                                                       
341                                                                               
342                                                                               
343                   
344       restart input directory         cn_ocerst_indir=
345 .                                                                             
346                                                                               
347                                                                               
348                   
349       file prefix restart output      cn_ocerst_out=
350 restart                                                                       
351                                                                               
352                                                                               
353                   
354       restart output directory        cn_ocerst_outdir=
355 .                                                                             
356                                                                               
357                                                                               
358                   
359       restart logical                 ln_rstart  =  F
360       start with forward time step    nn_euler   =            1
361       control of time step            nn_rstctl  =            0
362       number of the first time step   nn_it000   =            1
363       number of the last time step    nn_itend   =         5475
364       initial calendar date aammjj    nn_date0   =        10101
365       leap year calendar (0/1)        nn_leapy   =            0
366       initial state output            nn_istate  =            0
367       frequency of restart file       nn_stock   =         5475
368       frequency of output file        nn_write   =         5475
369       multi file dimgout              ln_dimgnnn =  F
370       mask land points                ln_mskland =  F
371       additional CF standard metadata ln_cfmeta  =  F
372       overwrite an existing file      ln_clobber =  F
373       NetCDF chunksize (bytes)        nn_chunksz =            0
374
375 ===>>> : W A R N I N G
376         ===============
377
378  ln_rstart =.FALSE., nn_euler is forced to 0                                   
379                                                                               
380                                                                               
381                   
382    The IOIPSL calendar is "noleap", i.e. no leap year
383 
384    Namelist namdom : space & time domain
385       flag read/compute bathymetry      nn_bathy     =            1
386       Depth (if =0 bathy=jpkm1)         rn_bathy     =   0.000000000000000E+000
387       min depth of the ocean    (>0) or    rn_hmin   =   -3.00000000000000     
388       min number of ocean level (<0)       
389       minimum thickness of partial      rn_e3zps_min =    20.0000000000000     
390  (m)
391          step level                     rn_e3zps_rat =   0.100000000000000     
392       create mesh/mask file(s)          nn_msh       =            1
393            = 0   no file created           
394            = 1   mesh_mask                 
395            = 2   mesh and mask             
396            = 3   mesh_hgr, msh_zgr and mask
397       ocean time step                       rn_rdt    =
398   5760.00000000000     
399       asselin time filter parameter         rn_atfp   =
400  0.100000000000000     
401       acceleration of converge              nn_acc    =            0
402         nn_acc=1: surface tracer rdt        rn_rdtmin =
403   28800.0000000000     
404                   bottom  tracer rdt        rdtmax    =
405   28800.0000000000     
406                   depth of transition       rn_rdth   =
407   800.000000000000     
408       suppression of closed seas (=0)       nn_closea =            0
409       online coarsening of dynamical fields ln_crs    =  F
410       type of horizontal mesh jphgr_msh           =            0
411       longitude of first raw and column T-point ppglam0 =
412   999999.000000000     
413       latitude  of first raw and column T-point ppgphi0 =
414   999999.000000000     
415       zonal      grid-spacing (degrees) ppe1_deg        =
416   999999.000000000     
417       meridional grid-spacing (degrees) ppe2_deg        =
418   999999.000000000     
419       zonal      grid-spacing (degrees) ppe1_m          =
420   999999.000000000     
421       meridional grid-spacing (degrees) ppe2_m          =
422   999999.000000000     
423       ORCA r4, r2 and r05 coefficients  ppsur           =
424  -4762.96143546300     
425                                         ppa0            =
426   255.580490704400     
427                                         ppa1            =
428   245.581322324900     
429                                         ppkth           =
430   21.4333619793800     
431                                         ppacr           =
432   3.00000000000000     
433       Minimum vertical spacing ppdzmin                  =
434   999999.000000000     
435       Maximum depth pphmax                              =
436   999999.000000000     
437       Use double tanf function for vertical coordinates ldbletanh =  F
438       Double tanh function parameters ppa2              =
439   999999.000000000     
440                                       ppkth2            =
441   999999.000000000     
442                                       ppacr2            =
443   999999.000000000     
444 
445    Namelist namcla
446       cross land advection                 nn_cla    =            0
447 
448 dom_clo : closed seas
449 ~~~~~~~
450 
451 dom_hgr : define the horizontal mesh from ithe following par_oce parameters
452 ~~~~~~~      type of horizontal mesh           jphgr_msh =            0
453              position of the first row and     ppglam0  =
454   999999.000000000     
455              column grid-point (degrees)       ppgphi0  =
456   999999.000000000     
457              zonal      grid-spacing (degrees) ppe1_deg =
458   999999.000000000     
459              meridional grid-spacing (degrees) ppe2_deg =
460   999999.000000000     
461              zonal      grid-spacing (meters)  ppe1_m   =
462   999999.000000000     
463              meridional grid-spacing (meters)  ppe2_m   =
464   999999.000000000     
465 
466           curvilinear coordinate on the sphere read in "coordinate" file
467 
468 hgr_read : read the horizontal coordinates
469 ~~~~~~~~      jpiglo =          182  jpjglo =          149  jpk =           31
470                     iom_nf90_open ~~~ open existing file: coordinates.nc in REA
471 D mode
472                    ---> coordinates.nc OK
473           read glamt (rec:      1) in coordinates.nc ok
474           read glamu (rec:      1) in coordinates.nc ok
475           read glamv (rec:      1) in coordinates.nc ok
476           read glamf (rec:      1) in coordinates.nc ok
477           read gphit (rec:      1) in coordinates.nc ok
478           read gphiu (rec:      1) in coordinates.nc ok
479           read gphiv (rec:      1) in coordinates.nc ok
480           read gphif (rec:      1) in coordinates.nc ok
481           read e1t (rec:      1) in coordinates.nc ok
482           read e1u (rec:      1) in coordinates.nc ok
483           read e1v (rec:      1) in coordinates.nc ok
484           read e1f (rec:      1) in coordinates.nc ok
485           read e2t (rec:      1) in coordinates.nc ok
486           read e2u (rec:      1) in coordinates.nc ok
487           read e2v (rec:      1) in coordinates.nc ok
488           read e2f (rec:      1) in coordinates.nc ok
489                     iom_close ~~~ close file: coordinates.nc ok
490 
491              orca_r2: Gibraltar    : e2u reduced to 20 km
492 
493              orca_r2: Bab el Mandeb: e2u reduced to 30 km
494                                      e1v reduced to 18 km
495 
496              orca_r2: Danish Straits : e2u reduced to 10 km
497 
498           longitude and e1 scale factors
499           ------------------------------
500    1   78.00    79.00    78.00    79.00    45515.2539158887    45515.2539158887    46299.2119389419    46299.2119389419
501   11   98.00    99.00    98.00    99.00    45515.2539158887    45515.2539158887    46299.2119389419    46299.2119389419
502   21  118.00   119.00   118.00   119.00    45515.2539158887    45515.2539158887    46299.2119389419    46299.2119389419
503   31  138.00   139.00   138.00   139.00    45515.2539158887    45515.2539158887    46299.2119389419    46299.2119389419
504 
505           latitude and e2 scale factors
506           -----------------------------
507    1  -78.19   -78.19   -77.98   -77.98    45515.2539158899    45515.2539158899    46299.2119389409    46299.2119389409
508   11  -73.32   -73.32   -73.03   -73.03    63846.4647224122    63846.4647224122    64922.0377934163    64922.0377934163
509   21  -66.51   -66.51   -66.11   -66.11    88633.7296404815    88633.7296404815    90061.7329991986    90061.7329991986
510   31  -57.16   -57.16   -56.61   -56.61   120613.5995973325   120613.5995973325   122389.7264144751   122389.7264144751
511 
512 dom_zgr : vertical coordinate
513 ~~~~~~~
514           Namelist namzgr : set vertical coordinate
515              z-coordinate - full steps      ln_zco    =  F
516              z-coordinate - partial steps   ln_zps    =  T
517              s- or hybrid z-s-coordinate    ln_sco    =  F
518              ice shelf cavities             ln_isfcav =  F
519 
520     zgr_z   : Reference vertical z-coordinates
521     ~~~~~~~
522            Value of coefficients for vertical mesh:
523                  zsur =   -4762.96143546300     
524                  za0  =    255.580490704400     
525                  za1  =    245.581322324900     
526                  zkth =    21.4333619793800     
527                  zacr =    3.00000000000000     
528 
529               Reference z-coordinate depth and scale factors:
530          level  gdept_1d  gdepw_1d  e3t_1d   e3w_1d 
531             1     5.00     0.00    10.00    10.00
532             2    15.00    10.00    10.00    10.00
533             3    25.00    20.00    10.00    10.00
534             4    35.01    30.00    10.01    10.00
535             5    45.01    40.01    10.01    10.01
536             6    55.03    50.02    10.02    10.02
537             7    65.06    60.04    10.04    10.03
538             8    75.13    70.09    10.09    10.06
539             9    85.25    80.18    10.17    10.12
540            10    95.49    90.35    10.33    10.24
541            11   105.97   100.69    10.65    10.47
542            12   116.90   111.36    11.27    10.91
543            13   128.70   122.65    12.47    11.77
544            14   142.20   135.16    14.78    13.43
545            15   158.96   150.03    19.23    16.65
546            16   181.96   169.42    27.66    22.78
547            17   216.65   197.37    43.26    34.30
548            18   272.48   241.13    70.88    55.21
549            19   364.30   312.74   116.11    90.99
550            20   511.53   429.72   181.55   146.43
551            21   732.20   611.89   261.03   220.35
552            22  1033.22   872.87   339.39   301.42
553            23  1405.70  1211.59   402.26   373.31
554            24  1830.89  1612.98   444.87   426.00
555            25  2289.77  2057.13   470.55   459.47
556            26  2768.24  2527.22   484.95   478.83
557            27  3257.48  3011.90   492.70   489.44
558            28  3752.44  3504.46   496.78   495.07
559            29  4250.40  4001.16   498.90   498.02
560            30  4749.91  4500.02   500.00   499.54
561            31  5250.23  5000.00   500.56   500.33
562 
563     zgr_bat : defines level and meter bathymetry
564     ~~~~~~~
565                     iom_nf90_open ~~~ open existing file: bathy_meter.nc in REA
566 D mode
567                    ---> bathy_meter.nc OK
568           read Bathymetry (rec:      1) in bathy_meter.nc ok
569                     iom_close ~~~ close file: bathy_meter.nc ok
570 
571       orca_r2: Gibraltar strait open at i=         140  j=         102
572 
573              orca_r2: Bab el Mandeb strait open at i=         160  j=
574          88
575 Minimum ocean depth:    30.0032158075637       minimum number of ocean levels :
576             3
577 
578     zgr_zps : z-coordinate with partial steps
579     ~~~~~~~
580               mbathy is recomputed : bathy_level file is NOT used
581 
582     zgr_bat_ctl : check the bathymetry
583     ~~~~~~~~~~~
584 
585                    suppress isolated ocean grid points
586                    -----------------------------------
587  the number of ocean level at grid-point (i,j) =            12          28
588  is changed from           30  to           29
589  the number of ocean level at grid-point (i,j) =            14          29
590  is changed from           30  to           29
591              137  ocean grid points suppressed
592  mbathy set to 0 along east and west boundary: nperio =            0
593 
594     zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels
595     ~~~~~~~~~~~~~
596 
597     zgr_top_level : ocean top k-index of T-, U-, V- and W-levels
598     ~~~~~~~~~~~~~
599 
600 dommsk : ocean mask
601 ~~~~~~
602    Namelist namlbc
603       lateral momentum boundary cond.    rn_shlat  =    2.00000000000000     
604       consistency with analytical form   ln_vorlat =  F
605    ocean lateral  no-slip
606 
607 dom_stp : time stepping setting
608 ~~~~~~~
609                synchronous time stepping
610                dynamics and tracer time step =    1.60000000000000       hours
611 
612 dom_wri : create NetCDF mesh and mask information file(s)
613 ~~~~~~~
614                     iom_nf90_open ~~~ create new file: mesh_mask_0000.nc in WRI
615 TE mode
616                    ---> mesh_mask_0000.nc OK
617           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask define dimensio
618 n variables done
619           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask defined ok
620           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask write dimension
621  variables done
622           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmask written ok
623           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umask defined ok
624           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umask written ok
625           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmask defined ok
626           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmask written ok
627           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmask defined ok
628           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmask written ok
629           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmaskutil defined ok
630           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: tmaskutil written ok
631           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umaskutil defined ok
632           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: umaskutil written ok
633           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmaskutil defined ok
634           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: vmaskutil written ok
635           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmaskutil defined ok
636           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: fmaskutil written ok
637           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamt defined ok
638           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamt written ok
639           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamu defined ok
640           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamu written ok
641           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamv defined ok
642           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamv written ok
643           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamf defined ok
644           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: glamf written ok
645           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphit defined ok
646           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphit written ok
647           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiu defined ok
648           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiu written ok
649           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiv defined ok
650           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphiv written ok
651           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphif defined ok
652           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gphif written ok
653           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1t defined ok
654           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1t written ok
655           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1u defined ok
656           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1u written ok
657           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1v defined ok
658           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1v written ok
659           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1f defined ok
660           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e1f written ok
661           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2t defined ok
662           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2t written ok
663           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2u defined ok
664           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2u written ok
665           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2v defined ok
666           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2v written ok
667           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2f defined ok
668           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e2f written ok
669           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: ff defined ok
670           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: ff written ok
671           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: mbathy defined ok
672           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: mbathy written ok
673           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: misf defined ok
674           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: misf written ok
675           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: isfdraft defined ok
676           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: isfdraft written ok
677           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_0 defined ok
678           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_0 written ok
679           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3u_0 defined ok
680           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3u_0 written ok
681           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3v_0 defined ok
682           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3v_0 written ok
683           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_0 defined ok
684           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_0 written ok
685           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_0 defined ok
686           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_0 written ok
687           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepu defined ok
688           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepu written ok
689           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepv defined ok
690           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepv written ok
691           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_0 defined ok
692           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_0 written ok
693           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_1d defined ok
694           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdept_1d written ok
695           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_1d defined ok
696           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: gdepw_1d written ok
697           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_1d defined ok
698           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3t_1d written ok
699           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_1d defined ok
700           iom_nf90_rp0123d, file: mesh_mask_0000.nc, var: e3w_1d written ok
701                     iom_close ~~~ close file: mesh_mask_0000.nc ok
702 
703 dom_ctl : extrema of the masked scale factors
704 ~~~~~~~
705              e1t maxi:  222397.85 at i =     2 j=    74
706              e1t mini:   47096.05 at i =    51 j=     2
707              e2t maxi:  209504.00 at i =     5 j=    92
708              e2t mini:   47096.05 at i =    51 j=     2
709 
710 istate_ini : Initialization of the dynamics and tracers
711 ~~~~~~~~~~
712 
713 dta_tsd_init : Temperature & Salinity data
714 ~~~~~~~~~~~~
715    Namelist namtsd
716       Initialisation of ocean T & S with T &S input data   ln_tsd_init   =  T
717       damping of ocean T & S toward T &S input data        ln_tsd_tradmp =  T
718 
719 
720 dta_tsd : Temperature & Salinity data
721 ~~~~~~~
722           namtsd Namelist
723           list of files and frequency (>0: in hours ; <0 in months)
724                root filename: ./data_1m_potential_temperature_nomask
725  variable name: votemper
726                frequency:   -1.00000000000000       time interp:  T
727  climatology:  T  weights    :  pairing    :  data type: yearly 
728  land/sea mask:
729                root filename: ./data_1m_salinity_nomask variable name: vosaline
730                frequency:   -1.00000000000000       time interp:  T
731  climatology:  T  weights    :  pairing    :  data type: yearly 
732  land/sea mask:
733  *** Info used values :
734    date ndastp                                      :        10100
735    number of elapsed days since the begining of run :   0.000000000000000E+000
736 
737 =======>> 1/2 time step before the start of the run DATE Y/M/D =      0/12/31  nsec_day:   83520  nsec_week: -175680
738======>> time-step =       1      New day, DATE Y/M/D = 0001/01/01      nday_year = 001
739         nsec_year =     2880   nsec_month =    2880   nsec_day =  2880
740   nsec_week =  -169920
741                     iom_nf90_open ~~~ open existing file: ./data_1m_potential_t
742 emperature_nomask.nc in READ mode
743                    ---> ./data_1m_potential_temperature_nomask.nc OK
744           read votemper (rec:     12) in ./data_1m_potential_temperature_nomask.nc ok
745fld_init : time-interpolation for votemper read previous record =     12 at time =  -15.50 days
746                     iom_nf90_open ~~~ open existing file: ./data_1m_salinity_no
747 mask.nc in READ mode
748                    ---> ./data_1m_salinity_nomask.nc OK
749           read vosaline (rec:     12) in ./data_1m_salinity_nomask.nc ok
750fld_init : time-interpolation for vosaline read previous record =     12 at time =  -15.50 days
751           read votemper (rec:      1) in ./data_1m_potential_temperature_nomask.nc ok
752           read vosaline (rec:      1) in ./data_1m_salinity_nomask.nc ok
753fld_read: var votemper kt =        1 (   0.0333 days), Y/M/D = 0001/01/01, records b/a:   0012/  0001 (days  -15.5000/  15.5000)
754 it_offset is :            0
755fld_read: var vosaline kt =        1 (   0.0333 days), Y/M/D = 0001/01/01, records b/a:   0012/  0001 (days  -15.5000/  15.5000)
756 it_offset is :            0
757  temperature Levitus
758 
759   level = 1
760
761
762        1    21
763  39   7.01  9.01
764  19  -0.08  0.00
765   level =           15
766
767
768        1    21
769  39   5.35  8.16
770  19  -1.92  0.00
771   level =           30
772
773
774        1    21
775  39   0.00  0.00
776  19   0.00  0.00
777 
778  salinity Levitus
779 
780   level = 1
781
782
783        1    21
784  39  33.90 34.35
785  19  33.68  0.00
786   level =           15
787
788
789        1    21
790  39  34.10 34.46
791  19  34.40  0.00
792   level =           30
793
794
795        1    21
796  39   0.00  0.00
797  19   0.00  0.00
798 
799 
800 sbc_init : surface boundary condition setting
801 ~~~~~~~~
802         Namelist namsbc (partly overwritten with CPP key setting)
803            frequency update of sbc (and ice)             nn_fsbc     =
804           5
805            Type of sbc :
806               analytical formulation                     ln_ana      =  F
807               flux       formulation                     ln_flx      =  F
808               CLIO bulk  formulation                     ln_blk_clio =  F
809               CORE bulk  formulation                     ln_blk_core =  T
810               MFS  bulk  formulation                     ln_blk_mfs  =  F
811               ocean-atmosphere coupled formulation       ln_cpl      =  F
812               forced-coupled mixed formulation           ln_mixcpl   =  F
813               OASIS coupling (with atm or sas)           lk_oasis    =  F
814               components of your executable              nn_components =
815           0
816               Multicategory heat flux formulation (LIM3) nn_limflx   =
817          -1
818            Misc. options of sbc :
819               Patm gradient added in ocean & ice Eqs.    ln_apr_dyn  =  F
820               ice management in the sbc (=0/1/2/3)       nn_ice      =
821           2
822               ice-ocean embedded/levitating (=0/1/2)     nn_ice_embd =
823           1
824               daily mean to diurnal cycle qsr            ln_dm2dc    =  F
825               runoff / runoff mouths                     ln_rnf      =  T
826               iceshelf formulation                       nn_isf      =
827           0
828               Sea Surface Restoring on SST and/or SSS    ln_ssr      =  T
829               FreshWater Budget control  (=0/1/2)        nn_fwb      =
830           2
831               closed sea (=0/1) (set in namdom)          nn_closea   =
832           0
833               n. of iterations if land-sea-mask applied  nn_lsm      =
834           0
835               Use of per-category fluxes (nn_limflx = -1)
836 
837               CORE bulk formulation
838 
839 sbc_ssm : sea surface mean fields
840 ~~~~~~~
841           default initialisation of ss?_m arrays
842 
843 sbc_ssr : SST and/or SSS damping term
844 ~~~~~~~
845    Namelist namsbc_ssr :
846       SST restoring term (Yes=1)             nn_sstr     =            0
847       SSS damping term (Yes=1, salt flux)    nn_sssr     =            2
848                        (Yes=2, volume flux)
849       dQ/dT (restoring magnitude on SST)     rn_dqdt     =
850  -40.0000000000000       W/m2/K
851       dE/dS (restoring magnitude on SST)     rn_deds     =
852  -166.670000000000       mm/day
853       flag to bound erp term                 ln_sssr_bnd =  T
854       ABS(Max./Min.) erp threshold           rn_sssr_bnd =
855   4.00000000000000       mm/day
856 
857 sbc_ssr : SSS restoring term toward SSS data
858 ~~~~~~~
859           namsbc_ssr Namelist
860           list of files and frequency (>0: in hours ; <0 in months)
861                root filename: ./sss_data variable name: sss
862                frequency:   -1.00000000000000       time interp:  T
863  climatology:  T  weights    :  pairing    :  data type: yearly 
864  land/sea mask:
865 
866 sbc_rnf : runoff
867 ~~~~~~~
868    Namelist namsbc_rnf
869       specific river mouths treatment            ln_rnf_mouth =  T
870       river mouth additional Kz                  rn_avt_rnf   =
871  1.000000000000000E-003
872       depth of river mouth additional mixing     rn_hrnf      =
873   15.0000000000000     
874       multiplicative factor for runoff           rn_rfact     =
875   1.00000000000000     
876 
877           runoffs inflow read in a file
878 
879 sbc_rnf_init : read runoffs data
880 ~~~~~~~~~~~~
881           namsbc_rnf Namelist
882           list of files and frequency (>0: in hours ; <0 in months)
883                root filename: ./runoff_core_monthly variable name: sorunoff
884                frequency:   -1.00000000000000       time interp:  T
885  climatology:  T  weights    :  pairing    :  data type: yearly 
886  land/sea mask:
887 
888           Specific treatment used in vicinity of river mouths :
889              - Increase Kz in surface layers (if rn_hrnf > 0 )
890                by   1.000000000000000E-003  m2/s  over            2  w-levels
891              - set to zero SSS damping       (if ln_ssr=T)
892              - mixed upstream-centered       (if ln_traadv_cen2=T)
893 
894 rnf_mouth : river mouth mask
895 ~~~~~~~~~
896                     iom_nf90_open ~~~ open existing file: ./runoff_core_monthly
897 .nc in READ mode
898                    ---> ./runoff_core_monthly.nc OK
899           read socoefr0 (rec:      1) in ./runoff_core_monthly.nc ok
900                     iom_close ~~~ close file: ./runoff_core_monthly.nc ok
901 
902 dyn_nept_init : Simplified Neptune module
903 ~~~~~~~~~~~~~
904  -->   Reading namelist namdyn_nept parameters:
905        ln_neptsimp          =  F
906 
907 
908 zdf_init: vertical physics
909 ~~~~~~~~
910    Namelist namzdf : set vertical mixing mixing parameters
911       vertical eddy viscosity             rn_avm0   =   1.200000000000000E-004
912       vertical eddy diffusivity           rn_avt0   =   1.200000000000000E-005
913       constant background or profile      nn_avb    =            0
914       horizontal variation for avtb       nn_havtb  =            0
915       time splitting / backward scheme    ln_zdfexp =  F
916       number of time step                 nn_zdfexp =            3
917       enhanced vertical diffusion         ln_zdfevd =  T
918          applied on momentum (=1/0)       nn_evdm   =            0
919       vertical coefficient for evd        rn_avevd  =    100.000000000000     
920       non-penetrative convection (npc)    ln_zdfnpc =  F
921       npc call  frequency                 nn_npc    =            1
922       npc print frequency                 nn_npcp   =          365
923 
924    vertical mixing option :
925       TKE dependent eddy coefficients
926 
927    convection :
928       use enhanced vertical dif. scheme
929       use the 1.5 turbulent closure
930 
931 zdf_bfr_init : momentum bottom friction
932 ~~~~~~~~~~~~~
933    Namelist nam_bfr : set bottom friction parameters
934       linear botton friction
935       bottom friction coef.   rn_bfri1  =   4.000000000000000E-004
936       implicit bottom friction switch                ln_bfrimp  =  T
937 
938 zdf_tke_init : tke turbulent closure scheme - initialisation
939 ~~~~~~~~~~~~
940    Namelist namzdf_tke : set tke mixing parameters
941       coef. to compute avt                        rn_ediff  =
942  0.100000000000000     
943       Kolmogoroff dissipation coef.               rn_ediss  =
944  0.700000000000000     
945       tke surface input coef.                     rn_ebb    =
946   67.8300000000000     
947       minimum value of tke                        rn_emin   =
948  1.000000000000000E-006
949       surface minimum value of tke                rn_emin0  =
950  1.000000000000000E-004
951       background shear (>0)                       rn_bshear =
952  9.999999999999999E-021
953       mixing length type                          nn_mxl    =            2
954       prandl number flag                          nn_pdl    =            1
955       surface mixing length = F(stress) or not    ln_mxl0   =  T
956       surface  mixing length minimum value        rn_mxl0   =
957  4.000000000000000E-002
958       flag to take into acc.  Langmuir circ.      ln_lc     =  T
959       coef to compute verticla velocity of LC     rn_lc     =
960  0.150000000000000     
961       test param. to add tke induced by wind      nn_etau   =            1
962       flag for computation of exp. tke profile    nn_htau   =            1
963       fraction of en which pene. the thermocline  rn_efr    =
964  5.000000000000000E-002
965 
966       critical Richardson nb with your parameters  ri_cri =
967  0.222222222222222     
968    use a surface mixing length = F(stress) :   set rn_mxl0 = rmxl_min
969 
970 zdf_tmx_init : tidal mixing
971 ~~~~~~~~~~~~
972    Namelist namzdf_tmx : set tidal mixing parameters
973       Vertical decay scale for turbulence   =    500.000000000000     
974       Brunt-Vaisala frequency threshold     =   1.000000000000000E-008
975       Tidal dissipation efficiency          =   0.333000000000000     
976       Mixing efficiency                     =   0.200000000000000     
977       ITF specific parameterisation         =  T
978       ITF tidal dissipation efficiency      =    1.00000000000000     
979                     iom_nf90_open ~~~ open existing file: mask_itf.nc in READ m
980 ode
981                    ---> mask_itf.nc OK
982           read tmaskitf (rec:      1) in mask_itf.nc ok
983                     iom_close ~~~ close file: mask_itf.nc ok
984                     iom_nf90_open ~~~ open existing file: M2rowdrg.nc in READ m
985 ode
986                    ---> M2rowdrg.nc OK
987           read field (rec:      1) in M2rowdrg.nc ok
988                     iom_close ~~~ close file: M2rowdrg.nc ok
989                     iom_nf90_open ~~~ open existing file: K1rowdrg.nc in READ m
990 ode
991                    ---> K1rowdrg.nc OK
992           read field (rec:      1) in K1rowdrg.nc ok
993                     iom_close ~~~ close file: K1rowdrg.nc ok
994 
995 zdf_ddm : double diffusive mixing
996 ~~~~~~~
997    Namelist namzdf_ddm : set dd mixing parameter
998       maximum avs for dd mixing      rn_avts   =   1.000000000000000E-004
999       heat/salt buoyancy flux ratio  rn_hsbfr  =    1.60000000000000     
1000 
1001 ldf_tra_init : lateral tracer physics
1002 ~~~~~~~~~~~~
1003    Namelist namtra_ldf : lateral mixing parameters (type, direction, coefficien
1004 ts)
1005       laplacian operator            ln_traldf_lap   =  T
1006       bilaplacian operator          ln_traldf_bilap =  F
1007       iso-level                     ln_traldf_level =  F
1008       horizontal (geopotential)     ln_traldf_hor   =  F
1009       iso-neutral                   ln_traldf_iso   =  T
1010       iso-neutral (Griffies)        ln_traldf_grif  =  F
1011       Griffies strmfn diagnostics   ln_traldf_gdia  =  F
1012       lateral eddy diffusivity      rn_aht_0        =    2000.00000000000     
1013       background hor. diffusivity   rn_ahtb_0       =   0.000000000000000E+000
1014       eddy induced velocity coef.   rn_aeiv_0       =    2000.00000000000     
1015       maximum isoppycnal slope      rn_slpmax       =   1.000000000000000E-002
1016       pure lateral mixing in ML     ln_triad_iso    =  F
1017       lateral mixing on bottom      ln_botmix_grif  =  F
1018 
1019           tracer mixing coef. = F( latitude, longitude)
1020           harmonic tracer diffusion (default)
1021 
1022  ldf_tra_c2d : 2D eddy diffusivity and eddy
1023  ~~~~~~~~~~~   --  induced velocity coefficients
1024                ORCA R2 or R4 case
1025                Constant values used for eddy diffusivity coefficients
1026                Variation lat/lon only for eddy induced velocity coefficients
1027 
1028 ldf_dyn : lateral momentum physics
1029 ~~~~~~~
1030    Namelist namdyn_ldf : set lateral mixing parameters
1031       laplacian operator                      ln_dynldf_lap   =  T
1032       bilaplacian operator                    ln_dynldf_bilap =  F
1033       iso-level                               ln_dynldf_level =  F
1034       horizontal (geopotential)               ln_dynldf_hor   =  T
1035       iso-neutral                             ln_dynldf_iso   =  F
1036       horizontal laplacian eddy viscosity     rn_ahm_0_lap    =
1037   40000.0000000000     
1038       background viscosity                    rn_ahmb_0       =
1039  0.000000000000000E+000
1040       horizontal bilaplacian eddy viscosity   rn_ahm_0_blp    =
1041  0.000000000000000E+000
1042       upper limit for laplacian eddy visc     rn_ahm_m_lap    =
1043   40000.0000000000     
1044       upper limit for bilap eddy viscosity    rn_ahm_m_blp    =
1045  -1000000000000.00     
1046    momentum mixing coef. = F( latitude, longitude, depth)
1047    harmonic momentum diff. (default)
1048 
1049 ldf_dyn_c3d : 3D lateral eddy viscosity coefficient
1050 ~~~~~~~~~~~
1051               laplacian operator: ahm proportional to e1
1052               maximum grid-spacing =    222397.846897091     
1053  maximum value for ahm =    40000.0000000000     
1054 
1055          ahm profile :
1056 
1057  jk      ahm         depth t-level
1058     1      1.0000      4.9999
1059     2      0.9994     15.0003
1060     3      0.9989     25.0018
1061     4      0.9982     35.0054
1062     5      0.9976     45.0133
1063     6      0.9970     55.0295
1064     7      0.9963     65.0618
1065     8      0.9956     75.1255
1066     9      0.9948     85.2504
1067    10      0.9940     95.4943
1068    11      0.9932    105.9699
1069    12      0.9923    116.8962
1070    13      0.9913    128.6979
1071    14      0.9901    142.1952
1072    15      0.9885    158.9606
1073    16      0.9861    181.9628
1074    17      0.9821    216.6479
1075    18      0.9745    272.4767
1076    19      0.9583    364.3030
1077    20      0.9192    511.5348
1078    21      0.8191    732.2009
1079    22      0.6067   1033.2173
1080    23      0.3759   1405.6975
1081    24      0.2766   1830.8850
1082    25      0.2544   2289.7679
1083    26      0.2506   2768.2423
1084    27      0.2501   3257.4789
1085    28      0.2500   3752.4422
1086    29      0.2500   4250.4012
1087    30      0.2500   4749.9133
1088    31      0.2500   5250.2266
1089 
1090          ahm profile :
1091 
1092  jk      ahm         depth t-level
1093     1      1.0000      4.9999
1094     2      0.9994     15.0003
1095     3      0.9989     25.0018
1096     4      0.9982     35.0054
1097     5      0.9976     45.0133
1098     6      0.9970     55.0295
1099     7      0.9963     65.0618
1100     8      0.9956     75.1255
1101     9      0.9948     85.2504
1102    10      0.9940     95.4943
1103    11      0.9932    105.9699
1104    12      0.9923    116.8962
1105    13      0.9913    128.6979
1106    14      0.9901    142.1952
1107    15      0.9885    158.9606
1108    16      0.9861    181.9628
1109    17      0.9821    216.6479
1110    18      0.9745    272.4767
1111    19      0.9583    364.3030
1112    20      0.9192    511.5348
1113    21      0.8191    732.2009
1114    22      0.6067   1033.2173
1115    23      0.3759   1405.6975
1116    24      0.2766   1830.8850
1117    25      0.2544   2289.7679
1118    26      0.2506   2768.2423
1119    27      0.2501   3257.4789
1120    28      0.2500   3752.4422
1121    29      0.2500   4250.4012
1122    30      0.2500   4749.9133
1123    31      0.2500   5250.2266
1124 
1125               ORCA R1, R2 or R4: overwrite the previous definition of ahm
1126               =================
1127 
1128 ldfdyn_c3d_orca : 3D eddy viscosity coefficient
1129 ~~~~~~~~~~~~~~~
1130         orca R1, R2 or R4 configuration: reduced in the surface Eq. strip
1131                     iom_nf90_open ~~~ open existing file: ahmcoef.nc in READ mo
1132 de
1133                    ---> ahmcoef.nc OK
1134           read icof (rec:      1) in ahmcoef.nc ok
1135                     iom_close ~~~ close file: ahmcoef.nc ok
1136 
1137          1D zcoef array
1138          ~~~~~~~~~~~~~~
1139 
1140     jk        zcoef
1141      1        1.000
1142      2        1.000
1143      3        1.000
1144      4        1.000
1145      5        1.000
1146      6        1.000
1147      7        1.000
1148      8        1.000
1149      9        1.000
1150     10        1.000
1151     11        1.000
1152     12        1.000
1153     13        1.000
1154     14        1.000
1155     15        1.000
1156     16        1.000
1157     17        1.000
1158     18        1.000
1159     19        1.000
1160     20        1.000
1161     21        1.000
1162     22        2.000
1163     23        3.000
1164     24        5.000
1165     25        7.000
1166     26        9.000
1167     27       10.000
1168     28       10.000
1169     29       10.000
1170     30       10.000
1171     31       10.000
1172 
1173          3D ahm1 array (k=1)
1174
1175
1176        1    21
1177  39  40.00 40.00
1178  19  40.00 40.00
1179 
1180          3D ahm2 array (k=1)
1181
1182
1183        1    21
1184  39  40.00 40.00
1185  19  40.00 40.00
1186 
1187          3D ahm2 array (k=jpk)
1188
1189
1190        1    21
1191  39  40.00 40.00
1192  19  40.00 40.00
1193 
1194 ldf_slp_init : direction of lateral mixing
1195 ~~~~~~~~~~~~
1196 
1197 tra_qsr_init : penetration of the surface solar radiation
1198 ~~~~~~~~~~~~
1199    Namelist namtra_qsr : set the parameter of penetration
1200       Light penetration (T) or not (F)         ln_traqsr  =  T
1201       RGB (Red-Green-Blue) light penetration   ln_qsr_rgb =  T
1202       2 band               light penetration   ln_qsr_2bd =  F
1203       bio-model            light penetration   ln_qsr_bio =  F
1204       light penetration for ice-model LIM3     ln_qsr_ice =  T
1205       RGB : Chl data (=1/2) or cst value (=0)  nn_chldta  =            1
1206       RGB & 2 bands: fraction of light (rn_si1)    rn_abs =
1207  0.580000000000000     
1208       RGB & 2 bands: shortess depth of extinction  rn_si0 =
1209  0.350000000000000     
1210       2 bands: longest depth of extinction         rn_si1 =
1211   23.0000000000000     
1212 
1213          R-G-B   light penetration - 2D Chl data
1214 
1215 trc_oce_rgb : Initialisation of the optical look-up table
1216 ~~~~~~~~~~~
1217       RGB longest depth of extinction    r_si2 =    61.8046971569839     
1218     jc =           1   Chl =   1.000000000000000E-002   irgb =            1
1219     jc =           2   Chl =   1.100000000000000E-002   irgb =            2
1220     jc =           3   Chl =   1.300000000000000E-002   irgb =            3
1221     jc =           4   Chl =   1.400000000000000E-002   irgb =            4
1222     jc =           5   Chl =   1.600000000000000E-002   irgb =            5
1223     jc =           6   Chl =   1.800000000000000E-002   irgb =            6
1224     jc =           7   Chl =   2.000000000000000E-002   irgb =            7
1225     jc =           8   Chl =   2.200000000000000E-002   irgb =            8
1226     jc =           9   Chl =   2.500000000000000E-002   irgb =            9
1227     jc =          10   Chl =   2.800000000000000E-002   irgb =           10
1228     jc =          11   Chl =   3.200000000000000E-002   irgb =           11
1229     jc =          12   Chl =   3.500000000000000E-002   irgb =           12
1230     jc =          13   Chl =   4.000000000000000E-002   irgb =           13
1231     jc =          14   Chl =   4.500000000000000E-002   irgb =           14
1232     jc =          15   Chl =   5.000000000000000E-002   irgb =           15
1233     jc =          16   Chl =   5.600000000000000E-002   irgb =           16
1234     jc =          17   Chl =   6.300000000000000E-002   irgb =           17
1235     jc =          18   Chl =   7.099999999999999E-002   irgb =           18
1236     jc =          19   Chl =   7.900000000000000E-002   irgb =           19
1237     jc =          20   Chl =   8.900000000000000E-002   irgb =           20
1238     jc =          21   Chl =   0.100000000000000        irgb =           21
1239     jc =          22   Chl =   0.112000000000000        irgb =           22
1240     jc =          23   Chl =   0.126000000000000        irgb =           23
1241     jc =          24   Chl =   0.141000000000000        irgb =           24
1242     jc =          25   Chl =   0.158000000000000        irgb =           25
1243     jc =          26   Chl =   0.178000000000000        irgb =           26
1244     jc =          27   Chl =   0.200000000000000        irgb =           27
1245     jc =          28   Chl =   0.224000000000000        irgb =           28
1246     jc =          29   Chl =   0.251000000000000        irgb =           29
1247     jc =          30   Chl =   0.282000000000000        irgb =           30
1248     jc =          31   Chl =   0.316000000000000        irgb =           31
1249     jc =          32   Chl =   0.355000000000000        irgb =           32
1250     jc =          33   Chl =   0.398000000000000        irgb =           33
1251     jc =          34   Chl =   0.447000000000000        irgb =           34
1252     jc =          35   Chl =   0.501000000000000        irgb =           35
1253     jc =          36   Chl =   0.562000000000000        irgb =           36
1254     jc =          37   Chl =   0.631000000000000        irgb =           37
1255     jc =          38   Chl =   0.708000000000000        irgb =           38
1256     jc =          39   Chl =   0.794000000000000        irgb =           39
1257     jc =          40   Chl =   0.891000000000000        irgb =           40
1258     jc =          41   Chl =    1.00000000000000        irgb =           41
1259     jc =          42   Chl =    1.12200000000000        irgb =           42
1260     jc =          43   Chl =    1.25900000000000        irgb =           43
1261     jc =          44   Chl =    1.41300000000000        irgb =           44
1262     jc =          45   Chl =    1.58500000000000        irgb =           45
1263     jc =          46   Chl =    1.77800000000000        irgb =           46
1264     jc =          47   Chl =    1.99500000000000        irgb =           47
1265     jc =          48   Chl =    2.23900000000000        irgb =           48
1266     jc =          49   Chl =    2.51200000000000        irgb =           49
1267     jc =          50   Chl =    2.81800000000000        irgb =           50
1268     jc =          51   Chl =    3.16200000000000        irgb =           51
1269     jc =          52   Chl =    3.54800000000000        irgb =           52
1270     jc =          53   Chl =    3.98100000000000        irgb =           53
1271     jc =          54   Chl =    4.46700000000000        irgb =           54
1272     jc =          55   Chl =    5.01200000000000        irgb =           55
1273     jc =          56   Chl =    5.62300000000000        irgb =           56
1274     jc =          57   Chl =    6.31000000000000        irgb =           57
1275     jc =          58   Chl =    7.07900000000000        irgb =           58
1276     jc =          59   Chl =    7.94300000000000        irgb =           59
1277     jc =          60   Chl =    8.91200000000000        irgb =           60
1278     jc =          61   Chl =    10.0000000000000        irgb =           61
1279         level of light extinction =           25  ref depth =
1280   2527.21691545123       m
1281 
1282         Chlorophyll read in a file
1283 
1284 tra_qsr_init : Solar penetration function of read chlorophyll
1285 ~~~~~~~~~~~~
1286           namtra_qsr Namelist
1287           list of files and frequency (>0: in hours ; <0 in months)
1288                root filename: ./chlorophyll variable name: CHLA
1289                frequency:   -1.00000000000000       time interp:  T
1290  climatology:  T  weights    :  pairing    :  data type: yearly 
1291  land/sea mask:
1292 
1293 tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating
1294 ~~~~~~~   
1295    Namelist nambbc : set bbc parameters
1296       Apply a geothermal heating at ocean bottom   ln_trabbc     =  T
1297       type of geothermal flux                      nn_geoflx     =            2
1298       Constant geothermal flux value               rn_geoflx_cst =
1299  8.640000000000000E-002
1300 
1301       *** variable geothermal heat flux
1302 
1303 tra_bbc_init : bottom temperature boundary condition
1304 ~~~~~~~~~~~~
1305           nambbc Namelist
1306           list of files and frequency (>0: in hours ; <0 in months)
1307                root filename: ./geothermal_heating.nc variable name: heatflow
1308                frequency:   -12.0000000000000       time interp:  F
1309  climatology:  T  weights    :  pairing    :  data type: yearly 
1310  land/sea mask:
1311                     iom_nf90_open ~~~ open existing file: ./geothermal_heating.
1312 nc in READ mode
1313                    ---> ./geothermal_heating.nc OK
1314           read heatflow (rec:      1) in ./geothermal_heating.nc ok
1315fld_read: var heatflow kt =        1 (   0.0333 days), Y/M/D = 0001/01/01, record:   0001 (days    0.0000 <->  365.0000)
1316 
1317 tra_bbl_init : bottom boundary layer initialisation
1318 ~~~~~~~~~~~~
1319        Namelist nambbl : set bbl parameters
1320           diffusive bbl (=1)   or not (=0)    nn_bbl_ldf =            1
1321           advective bbl (=1/2) or not (=0)    nn_bbl_adv =            0
1322           diffusive bbl coefficient           rn_ahtbbl  =
1323   1000.00000000000       m2/s
1324           advective bbl coefficient           rn_gambbl  =
1325   10.0000000000000       s
1326 
1327 tra_dmp_init : T and S newtonian relaxation
1328 ~~~~~~~
1329    Namelist namtra_dmp : set relaxation parameters
1330       Apply relaxation   or not       ln_tradmp =  T
1331       mixed layer damping option      nn_zdmp   =            0
1332       Damping file name               cn_resto  =
1333 resto.nc                                                                       
1334                                                                               
1335                                           
1336 
1337    tracer damping as specified by mask
1338                     iom_nf90_open ~~~ open existing file: resto.nc in READ mode
1339                    ---> resto.nc OK
1340           read resto (rec:      1) in resto.nc ok
1341                     iom_close ~~~ close file: resto.nc ok
1342 
1343 tra_adv_init : choice/control of the tracer advection scheme
1344 ~~~~~~~~~~~
1345    Namelist namtra_adv : chose a advection scheme for tracers
1346       2nd order advection scheme     ln_traadv_cen2    =  F
1347       TVD advection scheme           ln_traadv_tvd     =  T
1348       MUSCL  advection scheme        ln_traadv_muscl   =  F
1349       MUSCL2 advection scheme        ln_traadv_muscl2  =  F
1350       UBS    advection scheme        ln_traadv_ubs     =  F
1351       QUICKEST advection scheme      ln_traadv_qck     =  F
1352       upstream scheme within muscl   ln_traadv_msc_ups =  F
1353       TVD advection scheme with zts  ln_traadv_tvd_zts =  F
1354 
1355          TVD       scheme is used
1356 
1357 tra_adv_mle_init : mixed layer eddy (MLE) advection acting on tracers
1358 ~~~~~~~~~~~~~~~~
1359    Namelist namtra_adv_mle : mixed layer eddy advection on tracers
1360       use mixed layer eddy (MLE, i.e. Fox-Kemper param) (T/F)      ln_mle    =
1361 T
1362       MLE type: =0 standard Fox-Kemper ; =1 new formulation        nn_mle    =
1363           1
1364       magnitude of the MLE (typical value: 0.06 to 0.08)           rn_ce     =
1365  6.000000000000000E-002
1366       scale of ML front (ML radius of deformation) (rn_mle=0)      rn_lf     =
1367   5000.00000000000      m
1368       maximum time scale of MLE                    (rn_mle=0)      rn_time   =
1369   172800.000000000      s
1370       reference latitude (degrees) of MLE coef.    (rn_mle=1)      rn_lat    =
1371   20.0000000000000      deg
1372       space interp. of MLD at u-(v-)pts (0=min,1=averaged,2=max)   nn_mld_uv =
1373           0
1374       =1 no MLE in case of convection ; =0 always MLE              nn_conv   =
1375           0
1376       Density difference used to define ML for FK              rn_rho_c_mle  =
1377  1.000000000000000E-002
1378 
1379    Mixed Layer Eddy induced transport added to tracer advection
1380    New formulation
1381 
1382       ML buoyancy criteria =   9.558138401559454E-005  m/s2
1383       associated ML density criteria defined in zdfmxl =
1384  1.000000000000000E-002 kg/m3
1385 
1386 tra_ldf_init : lateral tracer diffusive operator
1387 ~~~~~~~~~~~
1388    Namelist namtra_ldf already read in ldftra module
1389    see ldf_tra_init report for lateral mixing parameters
1390 
1391 
1392           Rotated laplacian operator
1393 
1394 tra:ldf_ano : lateral diffusion acting on the full fields
1395 ~~~~~~~~~~~
1396 
1397 tra_zdf_init : vertical tracer physics scheme
1398 ~~~~~~~~~~~
1399               Implicit (euler backward) scheme
1400 
1401 dyn_adv_init : choice/control of the momentum advection scheme
1402 ~~~~~~~~~~~
1403        Namelist namdyn_adv : chose a advection formulation & scheme for momentu
1404 m
1405           Vector/flux form (T/F)                           ln_dynadv_vec  =  T
1406           = 0 standard scheme  ; =1 Hollingsworth scheme   nn_dynkeg      =
1407           0
1408           2nd order centred advection scheme               ln_dynadv_cen2 =  F
1409           3rd order UBS advection scheme                   ln_dynadv_ubs  =  F
1410           Sub timestepping of vertical advection           ln_dynzad_zts  =  F
1411 
1412          vector form : keg + zad + vor is used
1413 with Centered standard keg scheme
1414 
1415 dyn_vor_init : vorticity term : read namelist and control the consistency
1416 ~~~~~~~~~~~~
1417         Namelist namdyn_vor : choice of the vorticity term scheme
1418            energy    conserving scheme                ln_dynvor_ene =  F
1419            enstrophy conserving scheme                ln_dynvor_ens =  F
1420            mixed enstrophy/energy conserving scheme   ln_dynvor_mix =  F
1421            enstrophy and energy conserving scheme     ln_dynvor_een =  T
1422            enstrophy and energy conserving scheme (old) ln_dynvor_een_old=  F
1423 
1424          Vector form advection : vorticity = Coriolis + relative vorticity
1425 
1426          vorticity scheme : energy and enstrophy conserving scheme
1427 
1428 dyn_ldf_init : Choice of the lateral diffusive operator on dynamics
1429 ~~~~~~~~~~~
1430        Namelist nam_dynldf : set lateral mixing parameters (type, direction, co
1431 efficients)
1432           laplacian operator          ln_dynldf_lap   =  T
1433           bilaplacian operator        ln_dynldf_bilap =  F
1434           iso-level                   ln_dynldf_level =  F
1435           horizontal (geopotential)   ln_dynldf_hor   =  T
1436           iso-neutral                 ln_dynldf_iso   =  F
1437 
1438               laplacian operator
1439 
1440 dyn_hpg_init : hydrostatic pressure gradient initialisation
1441 ~~~~~~~~~~~~
1442    Namelist namdyn_hpg : choice of hpg scheme
1443       z-coord. - full steps                             ln_hpg_zco    =  F
1444       z-coord. - partial steps (interpolation)          ln_hpg_zps    =  T
1445       s-coord. (standard jacobian formulation)          ln_hpg_sco    =  F
1446       s-coord. (standard jacobian formulation) for isf  ln_hpg_isf    =  F
1447       s-coord. (Density Jacobian: Cubic polynomial)     ln_hpg_djc    =  F
1448       s-coord. (Pressure Jacobian: Cubic polynomial)    ln_hpg_prj    =  F
1449       time stepping: centered (F) or semi-implicit (T)  ln_dynhpg_imp =  F
1450 
1451 dyn_zdf_init : vertical dynamics physics scheme
1452 ~~~~~~~~~~~
1453               Implicit (euler backward) scheme
1454 
1455 dyn_spg_init : choice of the surface pressure gradient scheme
1456 ~~~~~~~~~~~
1457      Explicit free surface                  lk_dynspg_exp =  F
1458      Free surface with time splitting       lk_dynspg_ts  =  F
1459      Filtered free surface cst volume       lk_dynspg_flt =  T
1460 
1461      filtered free surface
1462      file   :
1463 solver.stat                                                                   
1464   open ok
1465      unit   =           11
1466      status = REPLACE
1467      form   = FORMATTED
1468      access = SEQUENTIAL
1469 
1470 
1471 solver_init : solver to compute the surface pressure gradient
1472 ~~~~~~~~~~~
1473    Namelist namsol : set solver parameters
1474       type of elliptic solver            nn_solv    =            1
1475       absolute/relative (0/1) precision  nn_sol_arp =            0
1476       minimum iterations for solver      nn_nmin    =          300
1477       maximum iterations for solver      nn_nmax    =          800
1478       frequency for test                 nn_nmod    =           10
1479       absolute precision of solver       rn_eps     =   1.000000000000000E-006
1480       absolute precision for SOR solver  rn_resmax  =   1.000000000000000E-010
1481       optimal coefficient of sor         rn_sor     =    1.92000000000000     
1482 
1483    a preconditioned conjugate gradient solver is used
1484 
1485 icbini :   AGRIF is not compatible with namelist namberg : 
1486          definition of rn_initial_mass(nclasses) with nclasses as PARAMETER
1487  namelist namberg not read
1488 
1489 icbini :   Namelist namberg ln_icebergs = F , NO icebergs used
1490 ~~~~~~~~
1491
1492 ===>>> : W A R N I N G
1493         ===============
1494
1495 W A R N I N G:  end of record or file while reading namelist namsto in configur
1496 ation namelist iostat =    -1
1497 
1498 sto_par_init : stochastic parameterization
1499 ~~~~~~~~~~~~
1500    Namelist namsto : stochastic parameterization
1501       restart stochastic parameters           ln_rststo     =  F
1502       read seed of RNG from restart file      ln_rstseed    =  T
1503       suffix of sto restart name (input)      cn_storst_in  =
1504 restart_sto                     
1505       suffix of sto restart name (output)     cn_storst_out =
1506 restart_sto                     
1507       stochastic equation of state            ln_sto_eos    =  F
1508       number of degrees of freedom            nn_sto_eos    =            1
1509       random walk horz. std (in grid points)  rn_eos_stdxy  =
1510   1.40000000000000     
1511       random walk vert. std (in grid points)  rn_eos_stdz   =
1512  0.700000000000000     
1513       random walk tcor (in timesteps)         rn_eos_tcor   =
1514   1440.00000000000     
1515       order of autoregressive  processes      nn_eos_ord    =            1
1516       passes of Laplacian filter              nn_eos_flt    =            0
1517       limitation factor                       rn_eos_lim    =
1518   2.00000000000000     
1519 
1520    stochastic parameterization :
1521 
1522 dia_ptr_init : poleward transport and msf initialization
1523 ~~~~~~~~~~~~
1524    Namelist namptr : set ptr parameters
1525       Poleward heat & salt transport (T) or not (F)      ln_diaptr  =  F
1526       Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins      ln_subbas  =  F
1527 
1528 dia_hsb_init
1529 ~~~~~~~~
1530   check the heat and salt budgets (T) or not (F)       ln_diahsb =  F
1531
1532 ===>>> : W A R N I N G
1533         ===============
1534
1535 W A R N I N G:  end of record or file while reading namelist namtrd in configur
1536 ation namelist iostat =    -1
1537 
1538  trd_init : Momentum/Tracers trends
1539  ~~~~~~~~~~
1540    Namelist namtrd : set trends parameters
1541       global domain averaged dyn & tra trends   ln_glo_trd  =  F
1542       U & V trends: 3D output                   ln_dyn_trd  =  F
1543       U & V trends: Mixed Layer averaged        ln_dyn_mxl  =  F
1544       T & S trends: 3D output                   ln_tra_trd  =  F
1545       T & S trends: Mixed Layer averaged        ln_tra_mxl  =  F
1546       Kinetic   Energy trends                   ln_KE_trd   =  F
1547       Potential Energy trends                   ln_PE_trd   =  F
1548       Barotropic vorticity trends               ln_vor_trd  =  F
1549       frequency of trends diagnostics (glo)     nn_trd      =          365
1550 
1551 dia_obs_init : Observation diagnostic initialization
1552 ~~~~~~~~~~~~
1553           Namelist namobs : set observation diagnostic parameters
1554              Logical switch for T profile observations          ln_t3d =  F
1555              Logical switch for S profile observations          ln_s3d =  F
1556              Logical switch for ENACT insitu data set           ln_ena =  F
1557              Logical switch for Coriolis insitu data set        ln_cor =  F
1558              Logical switch for feedback insitu data set      ln_profb =  F
1559              Logical switch for SLA observations                ln_sla =  F
1560              Logical switch for AVISO SLA data                ln_sladt =  F
1561              Logical switch for feedback SLA data             ln_slafb =  F
1562              Logical switch for SSH observations                ln_ssh =  F
1563              Logical switch for SST observations                ln_sst =  F
1564              Logical switch for Reynolds observations        ln_reysst =  F
1565              Logical switch for GHRSST observations          ln_ghrsst =  F
1566              Logical switch for feedback SST data             ln_sstfb =  F
1567              Logical switch for night-time SST obs         ln_sstnight =  F
1568              Logical switch for SSS observations                ln_sss =  F
1569              Logical switch for Sea Ice observations         ln_seaice =  F
1570              Logical switch for velocity observations         ln_vel3d =  F
1571              Logical switch for velocity daily av. cur.    ln_velavcur =  F
1572              Logical switch for velocity high freq. cur.   ln_velhrcur =  F
1573              Logical switch for velocity daily av. ADCP   ln_velavadcp =  F
1574              Logical switch for velocity high freq. ADCP  ln_velhradcp =  F
1575              Logical switch for feedback velocity data        ln_velfb =  F
1576              Global distribtion of observations         ln_grid_global =  F
1577              Logical switch for obs grid search w/lookup table  ln_grid_search_
1578 lookup =  F
1579              Initial date in window YYYYMMDD.HHMMSS        dobsini =
1580   20000101.0000000     
1581              Final date in window YYYYMMDD.HHMMSS          dobsend =
1582   20010101.0000000     
1583              Type of vertical interpolation method          n1dint =
1584           0
1585              Type of horizontal interpolation method        n2dint =
1586           0
1587              Rejection of observations near land swithch    ln_nea =  F
1588              MSSH correction scheme                         nmsshc =
1589           0
1590              MDT  correction                               mdtcorr =
1591   1.61000000000000     
1592              MDT cutoff for computed correction          mdtcutoff =
1593   65.0000000000000     
1594              Logical switch for alt bias                ln_altbias =  F
1595              Logical switch for ignoring missing files   ln_ignmis =  T
1596              ENACT daily average types                             =
1597         820          -1          -1          -1          -1          -1
1598          -1          -1          -1          -1          -1          -1
1599          -1          -1          -1          -1          -1          -1
1600          -1          -1
1601
1602 ===>>> : W A R N I N G
1603         ===============
1604
1605  key_diaobs is activated but logical flags
1606  ln_t3d, ln_s3d, ln_sla, ln_ssh, ln_sst, ln_sss, ln_seaice, ln_vel3d are all se
1607 t to .FALSE.
1608 
1609 dia_obs : Call the observation operators           0
1610 ~~~~~~~
1611 
1612 asm_inc_init : Assimilation increment initialization :
1613 ~~~~~~~~~~~~
1614    Namelist namasm : set assimilation increment parameters
1615       Logical switch for writing out background state          ln_bkgwri =  F
1616       Logical switch for applying tracer increments            ln_trainc =  F
1617       Logical switch for applying velocity increments          ln_dyninc =  F
1618       Logical switch for applying SSH increments               ln_sshinc =  F
1619       Logical switch for Direct Initialization (DI)            ln_asmdin =  F
1620       Logical switch for applying sea ice increments        ln_seaiceinc =  F
1621       Logical switch for Incremental Analysis Updating (IAU)   ln_asmiau =  F
1622       Timestep of background in [0,nitend-nit000-1]            nitbkg    =
1623           0
1624       Timestep of background for DI in [0,nitend-nit000-1]     nitdin    =
1625           0
1626       Timestep of start of IAU interval in [0,nitend-nit000-1] nitiaustr =
1627           1
1628       Timestep of end of IAU interval in [0,nitend-nit000-1]   nitiaufin =
1629          15
1630       Type of IAU weighting function                           niaufn    =
1631           0
1632       Logical switch for ensuring that the sa > salfixmin      ln_salfix =  F
1633       Minimum salinity after applying the increments           salfixmin =
1634  -9999.00000000000     
1635 
1636    Time steps referenced to current cycle:
1637        iitrst      =            0
1638        nit000      =            1
1639        nitend      =         5475
1640        nitbkg_r    =            0
1641        nitdin_r    =            0
1642        nitiaustr_r =            1
1643        nitiaufin_r =           15
1644 
1645    Dates referenced to current cycle:
1646        ndastp         =        10101
1647        ndate0         =        10101
1648        iitend_date    =        11231
1649        iitbkg_date    =        10101
1650        iitdin_date    =        10101
1651        iitiaustr_date =        10101
1652        iitiaufin_date =        10101
1653
1654 ===>>> : W A R N I N G
1655         ===============
1656
1657  ln_trainc, ln_dyninc, ln_sshinc and ln_seaiceinc are set to .false. :
1658  The assimilation increments are not applied
1659 Euler time step switch is            0
1660
1661AAAAAAAA
1662
1663 
1664 ~~~~~~~   mean fields initialised to instantaneous values
1665
1666 ===>>> : W A R N I N G
1667         ===============
1668
1669 sbcmod time step rdt * nn_fsbc is NOT a submultiple of atmospheric forcing freq
1670 uency
1671
1672 ===>>> : W A R N I N G
1673         ===============
1674
1675 sbcmod time step rdt * nn_fsbc is NOT a submultiple of atmospheric forcing freq
1676 uency
1677
1678 ===>>> : W A R N I N G
1679         ===============
1680
1681 sbcmod time step rdt * nn_fsbc is NOT a submultiple of atmospheric forcing freq
1682 uency
1683
1684 ===>>> : W A R N I N G
1685         ===============
1686
1687 sbcmod time step rdt * nn_fsbc is NOT a submultiple of atmospheric forcing freq
1688 uency
1689 
1690 sbc_blk_core : flux formulation for ocean surface boundary condition
1691 ~~~~~~~~~~~~
1692           namsbc_core Namelist
1693           list of files and frequency (>0: in hours ; <0 in months)
1694                root filename: ./u_10.15JUNE2009_fill variable name: U_10_MOD
1695                frequency:    6.00000000000000       time interp:  F
1696  climatology:  T  weights    : ./weights_core_orca2_bicubic_noc.nc
1697  pairing    : Uwnd data type: yearly   land/sea mask:
1698                root filename: ./v_10.15JUNE2009_fill variable name: V_10_MOD
1699                frequency:    6.00000000000000       time interp:  F
1700  climatology:  T  weights    : ./weights_core_orca2_bicubic_noc.nc
1701  pairing    : Vwnd data type: yearly   land/sea mask:
1702                root filename: ./q_10.15JUNE2009_fill variable name: Q_10_MOD
1703                frequency:    6.00000000000000       time interp:  F
1704  climatology:  T  weights    : ./weights_core_orca2_bilinear_noc.nc
1705  pairing    :  data type: yearly   land/sea mask:
1706                root filename: ./ncar_rad.15JUNE2009_fill variable name:
1707 SWDN_MOD
1708                frequency:    24.0000000000000       time interp:  F
1709  climatology:  T  weights    : ./weights_core_orca2_bilinear_noc.nc
1710  pairing    :  data type: yearly   land/sea mask:
1711                root filename: ./ncar_rad.15JUNE2009_fill variable name:
1712 LWDN_MOD
1713                frequency:    24.0000000000000       time interp:  F
1714  climatology:  T  weights    : ./weights_core_orca2_bilinear_noc.nc
1715  pairing    :  data type: yearly   land/sea mask:
1716                root filename: ./t_10.15JUNE2009_fill variable name: T_10_MOD
1717                frequency:    6.00000000000000       time interp:  F
1718  climatology:  T  weights    : ./weights_core_orca2_bilinear_noc.nc
1719  pairing    :  data type: yearly   land/sea mask:
1720                root filename: ./ncar_precip.15JUNE2009_fill variable name:
1721 PRC_MOD1
1722                frequency:   -1.00000000000000       time interp:  F
1723  climatology:  T  weights    : ./weights_core_orca2_bilinear_noc.nc
1724  pairing    :  data type: yearly   land/sea mask:
1725                root filename: ./ncar_precip.15JUNE2009_fill variable name: SNOW
1726                frequency:   -1.00000000000000       time interp:  F
1727  climatology:  T  weights    : ./weights_core_orca2_bilinear_noc.nc
1728  pairing    :  data type: yearly   land/sea mask:
1729                     iom_nf90_open ~~~ open existing file: ./u_10.15JUNE2009_fil
1730 l.nc in READ mode
1731                    ---> ./u_10.15JUNE2009_fill.nc OK
1732                     iom_nf90_open ~~~ open existing file: ./u_10.15JUNE2009_fil
1733 l.nc in READ mode
1734                    ---> ./u_10.15JUNE2009_fill.nc OK
1735                     iom_close ~~~ close file: ./u_10.15JUNE2009_fill.nc ok
1736                     iom_nf90_open ~~~ open existing file: ./weights_core_orca2_
1737 bicubic_noc.nc in READ mode
1738                    ---> ./weights_core_orca2_bicubic_noc.nc OK
1739           read src01 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1740           read src02 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1741           read src03 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1742           read src04 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1743           read wgt01 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1744           read wgt02 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1745           read wgt03 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1746           read wgt04 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1747           read wgt05 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1748           read wgt06 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1749           read wgt07 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1750           read wgt08 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1751           read wgt09 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1752           read wgt10 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1753           read wgt11 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1754           read wgt12 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1755           read wgt13 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1756           read wgt14 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1757           read wgt15 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1758           read wgt16 (rec:      1) in ./weights_core_orca2_bicubic_noc.nc ok
1759                     iom_close ~~~ close file: ./weights_core_orca2_bicubic_noc.
1760 nc ok
1761           read U_10_MOD (rec:      1) in ./u_10.15JUNE2009_fill.nc ok
1762                     iom_nf90_open ~~~ open existing file: ./v_10.15JUNE2009_fil
1763 l.nc in READ mode
1764                    ---> ./v_10.15JUNE2009_fill.nc OK
1765           read V_10_MOD (rec:      1) in ./v_10.15JUNE2009_fill.nc ok
1766                     iom_nf90_open ~~~ open existing file: ./q_10.15JUNE2009_fil
1767 l.nc in READ mode
1768                    ---> ./q_10.15JUNE2009_fill.nc OK
1769                     iom_nf90_open ~~~ open existing file: ./q_10.15JUNE2009_fil
1770 l.nc in READ mode
1771                    ---> ./q_10.15JUNE2009_fill.nc OK
1772                     iom_close ~~~ close file: ./q_10.15JUNE2009_fill.nc ok
1773                     iom_nf90_open ~~~ open existing file: ./weights_core_orca2_
1774 bilinear_noc.nc in READ mode
1775                    ---> ./weights_core_orca2_bilinear_noc.nc OK
1776           read src01 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1777           read src02 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1778           read src03 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1779           read src04 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1780           read wgt01 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1781           read wgt02 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1782           read wgt03 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1783           read wgt04 (rec:      1) in ./weights_core_orca2_bilinear_noc.nc ok
1784                     iom_close ~~~ close file: ./weights_core_orca2_bilinear_noc
1785 .nc ok
1786           read Q_10_MOD (rec:      1) in ./q_10.15JUNE2009_fill.nc ok
1787                     iom_nf90_open ~~~ open existing file: ./ncar_rad.15JUNE2009
1788 _fill.nc in READ mode
1789                    ---> ./ncar_rad.15JUNE2009_fill.nc OK
1790           read SWDN_MOD (rec:      1) in ./ncar_rad.15JUNE2009_fill.nc ok
1791                     iom_nf90_open ~~~ open existing file: ./ncar_rad.15JUNE2009
1792 _fill.nc in READ mode
1793                    ---> ./ncar_rad.15JUNE2009_fill.nc OK
1794           read LWDN_MOD (rec:      1) in ./ncar_rad.15JUNE2009_fill.nc ok
1795                     iom_nf90_open ~~~ open existing file: ./t_10.15JUNE2009_fil
1796 l.nc in READ mode
1797                    ---> ./t_10.15JUNE2009_fill.nc OK
1798           read T_10_MOD (rec:      1) in ./t_10.15JUNE2009_fill.nc ok
1799                     iom_nf90_open ~~~ open existing file: ./ncar_precip.15JUNE2
1800 009_fill.nc in READ mode
1801                    ---> ./ncar_precip.15JUNE2009_fill.nc OK
1802           read PRC_MOD1 (rec:      1) in ./ncar_precip.15JUNE2009_fill.nc ok
1803                     iom_nf90_open ~~~ open existing file: ./ncar_precip.15JUNE2
1804 009_fill.nc in READ mode
1805                    ---> ./ncar_precip.15JUNE2009_fill.nc OK
1806           read SNOW (rec:      1) in ./ncar_precip.15JUNE2009_fill.nc ok
1807 
1808  rot_rep : geographic <--> stretched
1809  ~~~~~    coordinate transformation