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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
ocean.output on Ticket #1403 – Attachment – NEMO

Ticket #1403: ocean.output

File ocean.output, 76.7 KB (added by ufla, 9 years ago)

Log file ocean.output for a test run that is interrupted by floating point overflow.

Line 
1 
2    CNRS - NERC - Met OFFICE - MERCATOR-ocean - INGV - CMCC
3                        NEMO team
4             Ocean General Circulation Model
5                   version 3.4  (2011)
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  =            1
41 
42 namcfg  : configuration initialization through namelist read
43 ~~~~~~~
44    Namelist namcfg
45       configuration name              cp_cfg      = orca
46       configuration zoom name         cp_cfz      =
47       configuration resolution        jp_cfg      =            1
48       1st lateral dimension ( >= jpi ) jpidta     =          362
49       2nd    "         "    ( >= jpj ) jpjdta     =          292
50       3nd    "         "               jpkdta     =           46
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 =            6
56 
57 mpp_init : Message Passing MPI
58 ~~~~~~~~
59 
60            defines mpp subdomains
61            ----------------------
62            iresti=           3  irestj=           2
63            jpni  =          21  jpnj  =           2
64 
65     ****************************************************************
66     *              *               *               *               *   
67   2 *   20  x147   *    20  x147   *    20  x147   *    19  x147   *   
68     *              *               *               *               *   
69     ****************************************************************
70     *              *               *               *               *   
71   1 *   20  x147   *    20  x147   *    20  x147   *    19  x147   *   
72     *              *               *               *               *   
73     ****************************************************************
74             1               2               3               4             
75 
76     ****************************************************************
77     *              *               *               *               *   
78   2 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
79     *              *               *               *               *   
80     ****************************************************************
81     *              *               *               *               *   
82   1 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
83     *              *               *               *               *   
84     ****************************************************************
85             5               6               7               8             
86 
87     ****************************************************************
88     *              *               *               *               *   
89   2 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
90     *              *               *               *               *   
91     ****************************************************************
92     *              *               *               *               *   
93   1 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
94     *              *               *               *               *   
95     ****************************************************************
96             9              10              11              12             
97 
98     ****************************************************************
99     *              *               *               *               *   
100   2 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
101     *              *               *               *               *   
102     ****************************************************************
103     *              *               *               *               *   
104   1 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
105     *              *               *               *               *   
106     ****************************************************************
107            13              14              15              16             
108 
109     ****************************************************************
110     *              *               *               *               *   
111   2 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
112     *              *               *               *               *   
113     ****************************************************************
114     *              *               *               *               *   
115   1 *   19  x147   *    19  x147   *    19  x147   *    19  x147   *   
116     *              *               *               *               *   
117     ****************************************************************
118            17              18              19              20             
119 
120     ****************
121     *              *   
122   2 *   19  x147   *   
123     *              *   
124     ****************
125     *              *   
126   1 *   19  x147   *   
127     *              *   
128     ****************
129            21             
130 
131  sum ilcit(i,1) =    362.000000000000       jpiglo =          362
132  sum ilcit(1,j) =    292.000000000000       jpjglo =          292
133 
134  nproc  =            0
135  nowe   =           20  noea   =             1
136  nono   =           21  noso   =           -21
137  nbondi =            0
138  nbondj =           -1
139  npolj  =            0
140  nperio =            0
141  nlci   =           20
142  nlcj   =          147
143  nimpp  =            1
144  njmpp  =            1
145  nbse   =            0  npse   =          -20
146  nbsw   =            0  npsw   =           -1
147  nbne   =            1  npne   =           22
148  nbnw   =            1  npnw   =           41
149 
150 mpp_init_ioipsl :   iloc  =           20         147
151 ~~~~~~~~~~~~~~~     iabsf =            1           1
152                     ihals =            0           0
153                     ihale =            1           1
154      file   :
155 timing.output                                                                 
156   open ok
157      unit   =           20
158      status = REPLACE
159      form   = FORMATTED
160      access = SEQUENTIAL
161 
162 
163  phy_cst : initialization of ocean parameters and constants
164  ~~~~~~~
165        Domain info
166           dimension of model
167                  Local domain      Global domain       Data domain
168               jpi     :      20   jpiglo  :     362   jpidta  :     362
169               jpj     :     147   jpjglo  :     292   jpjdta  :     292
170               jpk     :      46   jpk     :      46   jpkdta  :      46
171               jpij    :         2940
172           mpp local domain info (mpp)
173              jpni    :           21    jpreci  :            1
174              jpnj    :            2    jprecj  :            1
175              jpnij   :           42
176           lateral domain boundary condition type : jperio  =            6
177 
178        Constants
179 
180           mathematical constant                 rpi =    3.14159265358979     
181 
182           day                                rday   =    86400.0000000000     
183  s
184           sideral year                       rsiyea =    31558149.0101107     
185  s
186           sideral day                        rsiday =    86164.0996559118     
187  s
188           omega                              omega  =   7.292115083046062E-005
189  s^-1
190 
191           nb of months per year               raamo =    12.0000000000000     
192  months
193           nb of hours per day                 rjjhh =    24.0000000000000     
194  hours
195           nb of minutes per hour              rhhmm =    60.0000000000000     
196  mn
197           nb of seconds per minute            rmmss =    60.0000000000000     
198  s
199 
200           earth radius                         ra   =    6371229.00000000     
201  m
202           gravity                              grav =    9.80665000000000     
203  m/s^2
204 
205           triple point of temperature      rtt      =    273.160000000000     
206  K
207           freezing point of water          rt0      =    273.150000000000     
208  K
209           melting point of snow            rt0_snow =    273.160000000000     
210  K
211           melting point of ice             rt0_ice  =    273.160000000000     
212  K
213 
214           volumic mass of pure water          rauw  =    1000.00000000000     
215  kg/m^3
216           volumic mass of reference           rau0  =    1035.00000000000     
217  kg/m^3
218           1. / rau0                        r1_rau0  =   9.661835748792270E-004
219  m^3/kg
220           ocean specific heat                 rcp   =    4000.00000000000     
221  J/Kelvin
222           1. / ( rau0 * rcp )           r1_rau0_rcp =   2.415458937198068E-007
223 
224           thermal conductivity of the snow          =   0.310000000000000     
225  J/s/m/K
226           thermal conductivity of the ice           =    2.03439600000000     
227  J/s/m/K
228           fresh ice specific heat                   =    2067.00000000000     
229  J/kg/K
230           latent heat of fusion of fresh ice / snow =    334000.000000000     
231  J/kg
232           latent heat of subl.  of fresh ice / snow =    2834000.00000000     
233  J/kg
234           volumetric latent heat fusion of snow     =    110220000.000000     
235  J/m^3
236           density of sea ice                        =    917.000000000000     
237  kg/m^3
238           density of snow                           =    330.000000000000     
239  kg/m^3
240           emissivity of snow or ice                 =   0.970000000000000     
241           salinity of ice                           =    6.00000000000000     
242  psu
243           salinity of sea                           =    34.7000000000000     
244  psu
245           latent heat of evaporation (water)        =    2500000.00000000     
246  J/m^3
247           correction factor for solar radiation     =   0.900000000000000     
248           von Karman constant                       =   0.400000000000000     
249           Stefan-Boltzmann constant                 =   5.670000000000000E-008
250  J/s/m^2/K^4
251 
252           conversion: degre ==> radian          rad =   1.745329251994330E-002
253 
254           smallest real computer value       rsmall =   1.110223024625157E-016
255 
256 eos_init : equation of state
257 ~~~~~~~~
258           Namelist nameos : set eos parameters
259              flag for eq. of state and N^2  nn_eos   =            0
260              thermal exp. coef. (linear)    rn_alpha =   2.000000000000000E-004
261              saline  exp. coef. (linear)    rn_beta  =   7.700000000000000E-004
262 
263           use of Jackett & McDougall (1994) equation of state and
264                  McDougall (1987) Brunt-Vaisala frequency
265 
266 dom_cfg : set the ocean configuration
267 ~~~~~~~
268    ocean model configuration used :   cp_cfg =
269 orca                                                                           
270                                                                               
271                                                                               
272                     jp_cfg =            1
273    global domain lateral boundaries
274       jperio= 6, cyclic east-west and north fold with F-point pivot
275 
276 dom_glo : domain: data / local
277 ~~~~~~~
278           data input domain    : jpidta =          362  jpjdta =          292
279  jpkdta =           46
280           global or zoom domain: jpiglo =          362  jpjglo =          292
281  jpk    =           46
282           local domain         : jpi    =           20  jpj    =          147
283  jpk    =           46
284 
285           south-west indices    jpizoom =            1  jpjzoom =            1
286 
287           conversion local  ==> data i-index domain
288             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
289            20
290 
291           conversion data   ==> local  i-index domain
292              starting index
293             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
294            20  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
295            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
296            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
297            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
298            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
299            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
300            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
301            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
302            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
303            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
304            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
305            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
306            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
307            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
308            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
309            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
310            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
311            21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21  21
312            21
313              ending index
314             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
315            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
316            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
317            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
318            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
319            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
320            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
321            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
322            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
323            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
324            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
325            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
326            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
327            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
328            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
329            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
330            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
331            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
332            20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20  20
333            20
334 
335           conversion local  ==> data j-index domain
336             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
337            20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38
338            39  40  41  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57
339            58  59  60  61  62  63  64  65  66  67  68  69  70  71  72  73  74  75  76
340            77  78  79  80  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95
341            96  97  98  99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114
342           115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133
343           134 135 136 137 138 139 140 141 142 143 144 145 146 147
344 
345           conversion data  ==> local j-index domain
346              starting index
347             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
348            20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38
349            39  40  41  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57
350            58  59  60  61  62  63  64  65  66  67  68  69  70  71  72  73  74  75  76
351            77  78  79  80  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95
352            96  97  98  99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114
353           115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133
354           134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 148 148 148 148
355           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
356           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
357           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
358           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
359           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
360           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
361           148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148 148
362           148 148 148 148 148 148 148
363              ending index
364             1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19
365            20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38
366            39  40  41  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57
367            58  59  60  61  62  63  64  65  66  67  68  69  70  71  72  73  74  75  76
368            77  78  79  80  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95
369            96  97  98  99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114
370           115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133
371           134 135 136 137 138 139 140 141 142 143 144 145 146 147 147 147 147 147 147
372           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
373           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
374           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
375           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
376           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
377           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
378           147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147 147
379           147 147 147 147 147 147 147
380 
381           zoom flags :
382              lzoom   =  F  (T = zoom, F = global )
383              lzoom_e =  F  (T = forced closed east  boundary)
384              lzoom_w =  F  (T = forced closed west  boundary)
385              lzoom_s =  F  (T = forced closed South boundary)
386              lzoom_n =  F  (T = forced closed North boundary)
387 
388 dom_init : domain initialization
389 ~~~~~~~~
390 
391 dom_nam  : domain initialization through namelist read
392 ~~~~~~~
393    Namelist namrun
394       job number                      nn_no      =            0
395       experiment name for output      cn_exp     =
396 C36B                                                                           
397                                                                               
398                                                                               
399                   
400       file prefix restart input       cn_ocerst_in=
401 restart_oce                                                                   
402                                                                               
403                                                                               
404                   
405       file prefix restart output      cn_ocerst_out=
406 restart_oce                                                                   
407                                                                               
408                                                                               
409                   
410       restart logical                 ln_rstart  =  F
411       start with forward time step    nn_euler   =            1
412       control of time step            nn_rstctl  =            1
413       number of the first time step   nn_it000   =            1
414       number of the last time step    nn_itend   =          744
415       initial calendar date aammjj    nn_date0   =     19900101
416       leap year calendar (0/1)        nn_leapy   =            1
417       initial state output            nn_istate  =            0
418       frequency of restart file       nn_stock   =            0
419       frequency of output file        nn_write   =           -1
420       multi file dimgout              ln_dimgnnn =  F
421       mask land points                ln_mskland =  F
422       overwrite an existing file      ln_clobber =  F
423       NetCDF chunksize (bytes)        nn_chunksz =            0
424
425 ===>>> : W A R N I N G
426         ===============
427
428  ln_rstart =.FALSE., nn_euler is forced to 0                                   
429                                                                               
430                                                                               
431                   
432
433 ===>>> : W A R N I N G
434         ===============
435
436  nstock =            0 it is forced to          744                           
437                                                                               
438                                                                               
439                   
440    The IOIPSL calendar is "gregorian", i.e. leap year
441 
442    Namelist namdom : space & time domain
443       flag read/compute bathymetry      nn_bathy     =            1
444       Depth (if =0 bathy=jpkm1)         rn_bathy     =   0.000000000000000E+000
445       min depth of the ocean    (>0) or    rn_hmin   =   -3.00000000000000     
446       min number of ocean level (<0)       
447       minimum thickness of partial      rn_e3zps_min =    20.0000000000000     
448  (m)
449          step level                     rn_e3zps_rat =   0.100000000000000     
450       create mesh/mask file(s)          nn_msh       =            0
451            = 0   no file created           
452            = 1   mesh_mask                 
453            = 2   mesh and mask             
454            = 3   mesh_hgr, msh_zgr and mask
455       ocean time step                       rn_rdt    =
456   3600.00000000000     
457       asselin time filter parameter         rn_atfp   =
458  0.100000000000000     
459       acceleration of converge              nn_acc    =            0
460         nn_acc=1: surface tracer rdt        rn_rdtmin =
461   28800.0000000000     
462                   bottom  tracer rdt        rdtmax    =
463   28800.0000000000     
464                   depth of transition       rn_rdth   =
465   800.000000000000     
466       suppression of closed seas (=0)       nn_closea =            0
467       online coarsening of dynamical fields ln_crs    =  F
468       type of horizontal mesh jphgr_msh           =            0
469       longitude of first raw and column T-point ppglam0 =
470   999999.000000000     
471       latitude  of first raw and column T-point ppgphi0 =
472   999999.000000000     
473       zonal      grid-spacing (degrees) ppe1_deg        =
474   999999.000000000     
475       meridional grid-spacing (degrees) ppe2_deg        =
476   999999.000000000     
477       zonal      grid-spacing (degrees) ppe1_m          =
478   999999.000000000     
479       meridional grid-spacing (degrees) ppe2_m          =
480   999999.000000000     
481       ORCA r4, r2 and r05 coefficients  ppsur           =
482  -4762.96143546300     
483                                         ppa0            =
484   255.580490704400     
485                                         ppa1            =
486   245.581322324900     
487                                         ppkth           =
488   21.4333619793800     
489                                         ppacr           =
490   3.00000000000000     
491       Minimum vertical spacing ppdzmin                  =
492   999999.000000000     
493       Maximum depth pphmax                              =
494   999999.000000000     
495       Use double tanf function for vertical coordinates ldbletanh =  F
496       Double tanh function parameters ppa2              =
497   999999.000000000     
498                                       ppkth2            =
499   999999.000000000     
500                                       ppacr2            =
501   999999.000000000     
502 
503    Namelist namcla
504       cross land advection                 nn_cla    =            0
505 
506 dom_clo : closed seas
507 ~~~~~~~
508 
509 dom_hgr : define the horizontal mesh from ithe following par_oce parameters
510 ~~~~~~~      type of horizontal mesh           jphgr_msh =            0
511              position of the first row and     ppglam0  =
512   999999.000000000     
513              column grid-point (degrees)       ppgphi0  =
514   999999.000000000     
515              zonal      grid-spacing (degrees) ppe1_deg =
516   999999.000000000     
517              meridional grid-spacing (degrees) ppe2_deg =
518   999999.000000000     
519              zonal      grid-spacing (meters)  ppe1_m   =
520   999999.000000000     
521              meridional grid-spacing (meters)  ppe2_m   =
522   999999.000000000     
523 
524           curvilinear coordinate on the sphere read in "coordinate" file
525 
526 hgr_read : read the horizontal coordinates
527 ~~~~~~~~      jpiglo =          362  jpjglo =          292  jpk =           46
528                     iom_nf90_open ~~~ open existing file: coordinates.nc in REA
529 D mode
530                    ---> coordinates.nc OK
531           read glamt (rec:      1) in coordinates.nc ok
532           read glamu (rec:      1) in coordinates.nc ok
533           read glamv (rec:      1) in coordinates.nc ok
534           read glamf (rec:      1) in coordinates.nc ok
535           read gphit (rec:      1) in coordinates.nc ok
536           read gphiu (rec:      1) in coordinates.nc ok
537           read gphiv (rec:      1) in coordinates.nc ok
538           read gphif (rec:      1) in coordinates.nc ok
539           read e1t (rec:      1) in coordinates.nc ok
540           read e1u (rec:      1) in coordinates.nc ok
541           read e1v (rec:      1) in coordinates.nc ok
542           read e1f (rec:      1) in coordinates.nc ok
543           read e2t (rec:      1) in coordinates.nc ok
544           read e2u (rec:      1) in coordinates.nc ok
545           read e2v (rec:      1) in coordinates.nc ok
546           read e2f (rec:      1) in coordinates.nc ok
547                     iom_close ~~~ close file: coordinates.nc ok
548 
549              orca_r1: Gibraltar : e2u reduced to 20 km
550 
551              orca_r1: Bhosporus : e2u reduced to 10 km
552 
553              orca_r1: Lombok : e1v reduced to 10 km
554 
555              orca_r1: Sumba : e1v reduced to 8 km
556 
557              orca_r1: Ombai : e1v reduced to 13 km
558 
559              orca_r1: Timor Passage : e1v reduced to 20 km
560 
561              orca_r1: W Halmahera : e1v reduced to 30 km
562 
563              orca_r1: E Halmahera : e1v reduced to 50 km
564 
565           longitude and e1 scale factors
566           ------------------------------
567    1   72.50    73.00    72.50    73.00    22371.9996350562    22371.9996350562    22564.0248736808    22564.0248736808
568   11   82.50    83.00    82.50    83.00    22371.9996350562    22371.9996350562    22564.0248736808    22564.0248736808
569 
570           latitude and e2 scale factors
571           -----------------------------
572    1  -78.39   -78.39   -78.29   -78.29    22371.9996350562    22371.9996350562    22564.0248736808    22564.0248736808
573   11  -76.20   -76.20   -76.08   -76.08    26524.8299354813    26524.8299354813    26750.5179637307    26750.5179637307
574   21  -73.60   -73.60   -73.46   -73.46    31393.5671869958    31393.5671869958    31657.3897510458    31657.3897510458
575   31  -70.53   -70.53   -70.36   -70.36    37065.0717616678    37065.0717616678    37371.1271782116    37371.1271782116
576   41  -66.91   -66.91   -66.71   -66.71    43612.0739372158    43612.0739372158    43963.3187131693    43963.3187131693
577   51  -62.66   -62.66   -62.43   -62.43    51073.7561269995    51073.7561269995    51470.7880537348    51470.7880537348
578   61  -57.70   -57.70   -57.43   -57.43    59426.2991489708    59426.2991489708    59865.5913742379    59865.5913742379
579   71  -51.95   -51.95   -51.64   -51.64    68542.4202550136    68542.4202550136    69014.0413289406    69014.0413289406
580   81  -45.35   -45.35   -45.00   -45.00    78143.3450408712    78143.3450408712    78628.5295414200    78628.5295414200
581   91  -37.89   -37.89   -37.50   -37.50    87755.5230030943    87755.5230030943    88225.0221273672    88225.0221273672
582  101  -29.59   -29.59   -29.15   -29.15    96696.8748118213    96696.8748118213    97111.6439805149    97111.6439805149
583  111  -20.55   -20.55   -20.08   -20.08   104126.3062508076   104126.3062508076   104442.1936993819   104442.1936993819
584  121  -11.39   -11.39   -10.98   -10.98    90176.3801300774    90176.3801300774    87910.2293468842    87910.2293468842
585  131   -5.54    -5.54    -5.34    -5.34    44969.8813203525    44969.8813203525    43942.4027802381    43942.4027802381
586  141   -2.00    -2.00    -1.83    -1.83    37154.1706497673    37154.1706497673    37119.2489256383    37119.2489256383
587 
588 dom_zgr : vertical coordinate
589 ~~~~~~~
590           Namelist namzgr : set vertical coordinate
591              z-coordinate - full steps      ln_zco =  F
592              z-coordinate - partial steps   ln_zps =  T
593              s- or hybrid z-s-coordinate    ln_sco =  F
594 
595     zgr_z   : Reference vertical z-coordinates
596     ~~~~~~~
597            Value of coefficients for vertical mesh:
598                  zsur =   -4762.96143546300     
599                  za0  =    255.580490704400     
600                  za1  =    245.581322324900     
601                  zkth =    21.4333619793800     
602                  zacr =    3.00000000000000     
603 
604               Reference z-coordinate depth and scale factors:
605          level  gdept_1d  gdepw_1d  e3t_1d   e3w_1d 
606             1     5.00     0.00    10.00    10.00
607             2    15.00    10.00    10.00    10.00
608             3    25.00    20.00    10.00    10.00
609             4    35.01    30.00    10.01    10.00
610             5    45.01    40.01    10.01    10.01
611             6    55.03    50.02    10.02    10.02
612             7    65.06    60.04    10.04    10.03
613             8    75.13    70.09    10.09    10.06
614             9    85.25    80.18    10.17    10.12
615            10    95.49    90.35    10.33    10.24
616            11   105.97   100.69    10.65    10.47
617            12   116.90   111.36    11.27    10.91
618            13   128.70   122.65    12.47    11.77
619            14   142.20   135.16    14.78    13.43
620            15   158.96   150.03    19.23    16.65
621            16   181.96   169.42    27.66    22.78
622            17   216.65   197.37    43.26    34.30
623            18   272.48   241.13    70.88    55.21
624            19   364.30   312.74   116.11    90.99
625            20   511.53   429.72   181.55   146.43
626            21   732.20   611.89   261.03   220.35
627            22  1033.22   872.87   339.39   301.42
628            23  1405.70  1211.59   402.26   373.31
629            24  1830.89  1612.98   444.87   426.00
630            25  2289.77  2057.13   470.55   459.47
631            26  2768.24  2527.22   484.95   478.83
632            27  3257.48  3011.90   492.70   489.44
633            28  3752.44  3504.46   496.78   495.07
634            29  4250.40  4001.16   498.90   498.02
635            30  4749.91  4500.02   500.00   499.54
636            31  5250.23  5000.00   500.56   500.33
637            32  5750.95  5500.55   500.86   500.73
638            33  6251.89  6001.40   501.00   500.94
639            34  6752.94  6502.40   501.08   501.05
640            35  7254.04  7003.48   501.12   501.10
641            36  7755.17  7504.60   501.14   501.13
642            37  8256.32  8005.74   501.15   501.15
643            38  8757.47  8506.89   501.16   501.15
644            39  9258.63  9008.05   501.16   501.16
645            40  9759.79  9509.21   501.16   501.16
646            41 10260.95 10010.37   501.16   501.16
647            42 10762.11 10511.53   501.16   501.16
648            43 11263.27 11012.69   501.16   501.16
649            44 11764.43 11513.85   501.16   501.16
650            45 12265.60 12015.01   501.16   501.16
651            46 12766.76 12516.18   501.16   501.16
652 
653     zgr_bat : defines level and meter bathymetry
654     ~~~~~~~
655                     iom_nf90_open ~~~ open existing file: bathy_meter.nc in REA
656 D mode
657                    ---> bathy_meter.nc OK
658           read Bathymetry (rec:      1) in bathy_meter.nc ok
659                     iom_close ~~~ close file: bathy_meter.nc ok
660 Minimum ocean depth:    30.0032158075635       minimum number of ocean levels :
661             3
662 
663     zgr_zps : z-coordinate with partial steps
664     ~~~~~~~
665               mbathy is recomputed : bathy_level file is NOT used
666 
667     zgr_bat_ctl : check the bathymetry
668     ~~~~~~~~~~~
669 
670                    suppress isolated ocean grid points
671                    -----------------------------------
672  the number of ocean level at grid-point (i,j) =            10          41
673  is changed from           21  to           20
674  the number of ocean level at grid-point (i,j) =            11          42
675  is changed from           21  to           20
676  the number of ocean level at grid-point (i,j) =            19          43
677  is changed from           21  to           20
678  the number of ocean level at grid-point (i,j) =             5          52
679  is changed from           29  to           28
680  the number of ocean level at grid-point (i,j) =             8          68
681  is changed from           28  to           27
682  the number of ocean level at grid-point (i,j) =            13          99
683  is changed from           29  to           28
684  the number of ocean level at grid-point (i,j) =            10         100
685  is changed from           29  to           28
686  the number of ocean level at grid-point (i,j) =            16         122
687  is changed from           31  to           30
688              672  ocean grid points suppressed
689  mbathy set to 0 along east and west boundary: nperio =            0
690  maximum number of ocean level =           31  < jpk-1
691  you can decrease jpk to           32
692 
693     zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels
694     ~~~~~~~~~~~~~
695 
696 dommsk : ocean mask
697 ~~~~~~
698    Namelist namlbc
699       lateral momentum boundary cond.    rn_shlat  =    2.00000000000000     
700       consistency with analytical form   ln_vorlat =  F
701    ocean lateral  no-slip
702 
703    orca_r1: increase friction near the following straits :
704       Gibraltar
705       Bhosporus
706       Makassar (Top)
707       Lombok
708       Ombai
709       Timor Passage
710       West Halmahera
711       East Halmahera
712 
713 dom_stp : time stepping setting
714 ~~~~~~~
715                synchronous time stepping
716                dynamics and tracer time step =    1.00000000000000       hours
717 
718 dom_ctl : extrema of the masked scale factors
719 ~~~~~~~
720              e1t maxi:  111198.92 at i =     1 j=   147
721              e1t mini:   25207.90 at i =    92 j=     8
722              e2t maxi:  105405.52 at i =     1 j=   114
723              e2t mini:   24988.38 at i =   208 j=   283
724 
725 istate_ini : Initialization of the dynamics and tracers
726 ~~~~~~~~~~
727 
728 dta_tsd_init : Temperature & Salinity data
729 ~~~~~~~~~~~~
730    Namelist namtsd
731       Initialisation of ocean T & S with T &S input data   ln_tsd_init   =  T
732       damping of ocean T & S toward T &S input data        ln_tsd_tradmp =  F
733 
734 
735 dta_tsd : Temperature & Salinity data
736 ~~~~~~~
737           namtsd Namelist
738           list of files and frequency (>0: in hours ; <0 in months)
739                root filename: ./temperature variable name: votemper
740                frequency:   -1.00000000000000       time interp:  F
741  climatology:  T  weights    :  pairing    :  data type: monthly
742  land/sea mask:
743                root filename: ./salinity variable name: vosaline
744                frequency:   -1.00000000000000       time interp:  F
745  climatology:  T  weights    :  pairing    :  data type: monthly
746  land/sea mask:
747  *** Info used values :
748    date ndastp                                      :     19900100
749    number of elapsed days since the begining of run :   0.000000000000000E+000
750 
751 ==============>> 1/2 time step before the start of the run DATE Y/M/D =   1989/12/31  nsec_day: 84600
752  nsec_week: -1800
753======>> time-step =       1      New day, DATE Y/M/D = 1990/01/01      nday_year = 001
754         nsec_year =     1800   nsec_month =    1800   nsec_day =  1800
755   nsec_week =     1800
756                     iom_nf90_open ~~~ open existing file: ./temperature_m01.nc
757 in READ mode
758                    ---> ./temperature_m01.nc OK
759           read votemper (rec:      1) in ./temperature_m01.nc ok
760                     iom_nf90_open ~~~ open existing file: ./salinity_m01.nc in
761 READ mode
762                    ---> ./salinity_m01.nc OK
763           read vosaline (rec:      1) in ./salinity_m01.nc ok
764fld_read: var votemper kt =        1 (   0.0208 days), Y/M/D = 1990/01/01, record:   0001 (days    0.0000 <->   31.0000)
765fld_read: var vosaline kt =        1 (   0.0208 days), Y/M/D = 1990/01/01, record:   0001 (days    0.0000 <->   31.0000)
766  temperature Levitus
767 
768   level = 1
769
770
771        1
772 147  28.48
773 127  28.16
774 107  25.88
775  87  14.67
776  67   2.58
777  47   0.83
778  27   0.00
779   7   0.00
780   level =           23
781
782
783        1
784 147   7.51
785 127   6.54
786 107   7.44
787  87   6.20
788  67   2.12
789  47   1.03
790  27   0.00
791   7   0.00
792   level =           45
793
794
795        1
796 147   0.00
797 127   0.00
798 107   0.00
799  87   0.00
800  67   0.00
801  47   0.00
802  27   0.00
803   7   0.00
804 
805  salinity Levitus
806 
807   level = 1
808
809
810        1
811 147  35.15
812 127  34.36
813 107  35.42
814  87  34.89
815  67  33.88
816  47  33.67
817  27   0.00
818   7   0.00
819   level =           23
820
821
822        1
823 147  34.97
824 127  34.78
825 107  34.54
826  87  34.45
827  67  34.67
828  47  34.71
829  27   0.00
830   7   0.00
831   level =           45
832
833
834        1
835 147   0.00
836 127   0.00
837 107   0.00
838  87   0.00
839  67   0.00
840  47   0.00
841  27   0.00
842   7   0.00
843 
844 dta_tsd: deallocte T & S arrays as they are only use to initialize the run
845 
846 dyn_nept_init : Simplified Neptune module
847 ~~~~~~~~~~~~~
848  -->   Reading namelist namdyn_nept parameters:
849        ln_neptsimp          =  F
850 
851 
852 sbc_init : surface boundary condition setting
853 ~~~~~~~~
854         Namelist namsbc (partly overwritten with CPP key setting)
855            frequency update of sbc (and ice)             nn_fsbc     =
856           1
857            Type of sbc :
858               analytical formulation                     ln_ana      =  F
859               flux       formulation                     ln_flx      =  F
860               CLIO bulk  formulation                     ln_blk_clio =  F
861               CORE bulk  formulation                     ln_blk_core =  F
862               MFS  bulk  formulation                     ln_blk_mfs  =  F
863               coupled    formulation (T if key_sbc_cpl)  ln_cpl      =  T
864               Flux handling over ice categories          cn_iceflx   = linear
865            Misc. options of sbc :
866               Patm gradient added in ocean & ice Eqs.    ln_apr_dyn  =  F
867               ice management in the sbc (=0/1/2/3)       nn_ice      =
868           3
869               ice-ocean embedded/levitating (=0/1/2)     nn_ice_embd =
870           1
871               daily mean to diurnal cycle qsr            ln_dm2dc    =  F
872               runoff / runoff mouths                     ln_rnf      =  F
873               Sea Surface Restoring on SST and/or SSS    ln_ssr      =  F
874               FreshWater Budget control  (=0/1/2)        nn_fwb      =
875           0
876               closed sea (=0/1) (set in namdom)          nn_closea   =
877           0
878               n. of iterations if land-sea-mask applied  nn_lsm      =
879           0
880               Fluxes averaged over all ice categories         ln_iceflx_ave   
881 =  F
882               Fluxes distributed linearly over ice categories ln_iceflx_linear
883 =  T
884 
885               coupled formulation
886 
887 sbc_ssm : sea surface mean fields, nn_fsbc=1 : instantaneous values
888 ~~~~~~~
889 
890 sbc_cpl_init : namsbc_cpl namelist
891 ~~~~~~~~~~~~
892   received fields (mutiple ice categogies)
893       10m wind module                 = none (no)
894       stress module                   = none (no)
895       surface stress                  = oce and ice (no)
896                      - referential    = spherical                       
897                      - orientation    = eastward-northward             
898                      - mesh           = T                               
899       non-solar heat flux sensitivity = coupled (no)
900       solar heat flux                 = oce and ice (yes)
901       non-solar heat flux             = oce and ice (yes)
902       freshwater budget               = conservative (no)
903       runoffs                         = coupled (no)
904       calving                         = coupled (no)
905       sea ice heat fluxes             = none (no)
906       atm co2                         = none (no)
907   sent fields (multiple ice categories)
908       surface temperature             = oce and ice (yes)
909       albedo                          = ice (no)
910       ice/snow thickness              = ice and snow (yes)
911       surface current                 = none (no)
912                       - referential   = spherical                       
913                       - orientation   = eastward-northward             
914                       - mesh          = T                               
915       oce co2 flux                    = none (no)
916 
917 cpl_prism_define : initialization in coupled ocean/atmosphere case
918 ~~~~~~~~~~~~~~~~~
919 
920Define  1  1 OIceFrc.C001 for 20
921Define  1  2 OIceFrc.C002 for 20
922Define  1  3 OIceFrc.C003 for 20
923Define  1  4 OIceFrc.C004 for 20
924Define  1  5 OIceFrc.C005 for 20
925Define  2  1 O_SSTSST for 20
926Define  3  1 OTepIce.C001 for 20
927Define  3  2 OTepIce.C002 for 20
928Define  3  3 OTepIce.C003 for 20
929Define  3  4 OTepIce.C004 for 20
930Define  3  5 OTepIce.C005 for 20
931Define  5  1 OAlbIce for 20
932Define  7  1 OIceTck.C001 for 20
933Define  7  2 OIceTck.C002 for 20
934Define  7  3 OIceTck.C003 for 20
935Define  7  4 OIceTck.C004 for 20
936Define  7  5 OIceTck.C005 for 20
937Define  8  1 OSnwTck.C001 for 20
938Define  8  2 OSnwTck.C002 for 20
939Define  8  3 OSnwTck.C003 for 20
940Define  8  4 OSnwTck.C004 for 20
941Define  8  5 OSnwTck.C005 for 20
942Define  1  1 O_OTaux1 for 21
943Define  2  1 O_OTauy1 for 21
944Define  7  1 O_ITaux1 for 21
945Define  8  1 O_ITauy1 for 21
946Define 13  1 O_QsrOce for 21
947Define 14  1 O_QsrIce.C001 for 21
948Define 14  2 O_QsrIce.C002 for 21
949Define 14  3 O_QsrIce.C003 for 21
950Define 14  4 O_QsrIce.C004 for 21
951Define 14  5 O_QsrIce.C005 for 21
952Define 16  1 O_QnsOce for 21
953Define 17  1 O_QnsIce.C001 for 21
954Define 17  2 O_QnsIce.C002 for 21
955Define 17  3 O_QnsIce.C003 for 21
956Define 17  4 O_QnsIce.C004 for 21
957Define 17  5 O_QnsIce.C005 for 21
958Define 19  1 OTotRain for 21
959Define 20  1 OTotSnow for 21
960Define 21  1 OTotEvap for 21
961Define 22  1 OIceEvap for 21
962Define 27  1 O_dQnsdT for 21
963Define 28  1 O_Runoff for 21
964Define 29  1 OCalving for 21
965 
966 zdf_init: vertical physics
967 ~~~~~~~~
968    Namelist namzdf : set vertical mixing mixing parameters
969       vertical eddy viscosity             rn_avm0   =   1.200000000000000E-004
970       vertical eddy diffusivity           rn_avt0   =   1.200000000000000E-005
971       constant background or profile      nn_avb    =            0
972       horizontal variation for avtb       nn_havtb  =            0
973       time splitting / backward scheme    ln_zdfexp =  F
974       number of time step                 nn_zdfexp =            3
975       enhanced vertical diffusion         ln_zdfevd =  T
976          applied on momentum (=1/0)       nn_evdm   =            0
977       vertical coefficient for evd        rn_avevd  =    100.000000000000     
978       non-penetrative convection (npc)    ln_zdfnpc =  F
979       npc call  frequency                 nn_npc    =            1
980       npc print frequency                 nn_npcp   =          365
981 
982    vertical mixing option :
983       TKE dependent eddy coefficients
984 
985    convection :
986       use enhanced vertical dif. scheme
987       use the 1.5 turbulent closure
988 
989 zdf_bfr_init : momentum bottom friction
990 ~~~~~~~~~~~~~
991    Namelist nam_bfr : set bottom friction parameters
992       linear botton friction
993       friction coef.   rn_bfri1  =   4.000000000000000E-004
994       implicit bottom friction switch                ln_bfrimp  =  T
995 
996 zdf_tke_init : tke turbulent closure scheme - initialisation
997 ~~~~~~~~~~~~
998    Namelist namzdf_tke : set tke mixing parameters
999       coef. to compute avt                        rn_ediff  =
1000  0.100000000000000     
1001       Kolmogoroff dissipation coef.               rn_ediss  =
1002  0.700000000000000     
1003       tke surface input coef.                     rn_ebb    =
1004   67.8300000000000     
1005       minimum value of tke                        rn_emin   =
1006  1.000000000000000E-006
1007       surface minimum value of tke                rn_emin0  =
1008  1.000000000000000E-004
1009       background shear (>0)                       rn_bshear =
1010  9.999999999999999E-021
1011       mixing length type                          nn_mxl    =            2
1012       prandl number flag                          nn_pdl    =            1
1013       surface mixing length = F(stress) or not    ln_mxl0   =  T
1014       surface  mixing length minimum value        rn_mxl0   =
1015  4.000000000000000E-002
1016       flag to take into acc.  Langmuir circ.      ln_lc     =  T
1017       coef to compute verticla velocity of LC     rn_lc     =
1018  0.150000000000000     
1019       test param. to add tke induced by wind      nn_etau   =            1
1020       flag for computation of exp. tke profile    nn_htau   =            1
1021       fraction of en which pene. the thermocline  rn_efr    =
1022  5.000000000000000E-002
1023 
1024       critical Richardson nb with your parameters  ri_cri =
1025  0.222222222222222     
1026    use a surface mixing length = F(stress) :   set rn_mxl0 = rmxl_min
1027 
1028 zdf_ddm : double diffusive mixing
1029 ~~~~~~~
1030    Namelist namzdf_ddm : set dd mixing parameter
1031       maximum avs for dd mixing      rn_avts   =   1.000000000000000E-004
1032       heat/salt buoyancy flux ratio  rn_hsbfr  =    1.60000000000000     
1033 
1034 ldf_tra_init : lateral tracer physics
1035 ~~~~~~~~~~~~
1036    Namelist namtra_ldf : lateral mixing parameters (type, direction, coefficien
1037 ts)
1038       laplacian operator            ln_traldf_lap   =  T
1039       bilaplacian operator          ln_traldf_bilap =  F
1040       iso-level                     ln_traldf_level =  F
1041       horizontal (geopotential)     ln_traldf_hor   =  F
1042       iso-neutral                   ln_traldf_iso   =  T
1043       iso-neutral (Griffies)        ln_traldf_grif  =  F
1044       Griffies strmfn diagnostics   ln_traldf_gdia  =  F
1045       lateral eddy diffusivity      rn_aht_0        =    1000.00000000000     
1046       background hor. diffusivity   rn_ahtb_0       =   0.000000000000000E+000
1047       eddy induced velocity coef.   rn_aeiv_0       =    1000.00000000000     
1048       maximum isoppycnal slope      rn_slpmax       =   1.000000000000000E-002
1049       pure lateral mixing in ML     ln_triad_iso    =  F
1050       lateral mixing on bottom      ln_botmix_grif  =  F
1051 
1052           tracer mixing coef. = F( latitude, longitude)
1053           harmonic tracer diffusion (default)
1054 
1055  ldf_tra_c2d : 2D eddy diffusivity and eddy
1056  ~~~~~~~~~~~   --  induced velocity coefficients
1057 
1058 ldf_dyn : lateral momentum physics
1059 ~~~~~~~
1060    Namelist namdyn_ldf : set lateral mixing parameters
1061       laplacian operator                      ln_dynldf_lap   =  T
1062       bilaplacian operator                    ln_dynldf_bilap =  F
1063       iso-level                               ln_dynldf_level =  F
1064       horizontal (geopotential)               ln_dynldf_hor   =  T
1065       iso-neutral                             ln_dynldf_iso   =  F
1066       horizontal laplacian eddy viscosity     rn_ahm_0_lap    =
1067   20000.0000000000     
1068       background viscosity                    rn_ahmb_0       =
1069  0.000000000000000E+000
1070       horizontal bilaplacian eddy viscosity   rn_ahm_0_blp    =
1071  0.000000000000000E+000
1072       upper limit for laplacian eddy visc     rn_ahm_m_lap    =
1073   20000.0000000000     
1074       upper limit for bilap eddy viscosity    rn_ahm_m_blp    =
1075  -1000000000000.00     
1076    momentum mixing coef. = F( latitude, longitude, depth)
1077    harmonic momentum diff. (default)
1078 
1079 ldf_dyn_c3d : 3D lateral eddy viscosity coefficient
1080 ~~~~~~~~~~~
1081               laplacian operator: ahm proportional to e1
1082               maximum grid-spacing =    111198.923448546     
1083  maximum value for ahm =    20000.0000000000     
1084 
1085          ahm profile :
1086 
1087  jk      ahm         depth t-level
1088     1      1.0000      4.9999
1089     2      0.9994     15.0003
1090     3      0.9989     25.0018
1091     4      0.9982     35.0054
1092     5      0.9976     45.0133
1093     6      0.9970     55.0295
1094     7      0.9963     65.0618
1095     8      0.9956     75.1255
1096     9      0.9948     85.2504
1097    10      0.9940     95.4943
1098    11      0.9932    105.9699
1099    12      0.9923    116.8962
1100    13      0.9913    128.6979
1101    14      0.9901    142.1952
1102    15      0.9885    158.9606
1103    16      0.9861    181.9628
1104    17      0.9821    216.6479
1105    18      0.9745    272.4767
1106    19      0.9583    364.3030
1107    20      0.9192    511.5348
1108    21      0.8191    732.2009
1109    22      0.6067   1033.2173
1110    23      0.3759   1405.6975
1111    24      0.2766   1830.8850
1112    25      0.2544   2289.7679
1113    26      0.2506   2768.2423
1114    27      0.2501   3257.4789
1115    28      0.2500   3752.4422
1116    29      0.2500   4250.4012
1117    30      0.2500   4749.9133
1118    31      0.2500   5250.2266
1119    32      0.2500   5750.9524
1120    33      0.2500   6251.8903
1121    34      0.2500   6752.9371
1122    35      0.2500   7254.0399
1123    36      0.2500   7755.1714
1124    37      0.2500   8256.3176
1125    38      0.2500   8757.4714
1126    39      0.2500   9258.6291
1127    40      0.2500   9759.7888
1128    41      0.2500  10260.9496
1129    42      0.2500  10762.1108
1130    43      0.2500  11263.2724
1131    44      0.2500  11764.4340
1132    45      0.2500  12265.5958
1133    46      0.2500  12766.7575
1134 
1135          ahm profile :
1136 
1137  jk      ahm         depth t-level
1138     1      1.0000      4.9999
1139     2      0.9994     15.0003
1140     3      0.9989     25.0018
1141     4      0.9982     35.0054
1142     5      0.9976     45.0133
1143     6      0.9970     55.0295
1144     7      0.9963     65.0618
1145     8      0.9956     75.1255
1146     9      0.9948     85.2504
1147    10      0.9940     95.4943
1148    11      0.9932    105.9699
1149    12      0.9923    116.8962
1150    13      0.9913    128.6979
1151    14      0.9901    142.1952
1152    15      0.9885    158.9606
1153    16      0.9861    181.9628
1154    17      0.9821    216.6479
1155    18      0.9745    272.4767
1156    19      0.9583    364.3030
1157    20      0.9192    511.5348
1158    21      0.8191    732.2009
1159    22      0.6067   1033.2173
1160    23      0.3759   1405.6975
1161    24      0.2766   1830.8850
1162    25      0.2544   2289.7679
1163    26      0.2506   2768.2423
1164    27      0.2501   3257.4789
1165    28      0.2500   3752.4422
1166    29      0.2500   4250.4012
1167    30      0.2500   4749.9133
1168    31      0.2500   5250.2266
1169    32      0.2500   5750.9524
1170    33      0.2500   6251.8903
1171    34      0.2500   6752.9371
1172    35      0.2500   7254.0399
1173    36      0.2500   7755.1714
1174    37      0.2500   8256.3176
1175    38      0.2500   8757.4714
1176    39      0.2500   9258.6291
1177    40      0.2500   9759.7888
1178    41      0.2500  10260.9496
1179    42      0.2500  10762.1108
1180    43      0.2500  11263.2724
1181    44      0.2500  11764.4340
1182    45      0.2500  12265.5958
1183    46      0.2500  12766.7575
1184 
1185               ORCA R1, R2 or R4: overwrite the previous definition of ahm
1186               =================
1187 
1188 ldfdyn_c3d_orca : 3D eddy viscosity coefficient
1189 ~~~~~~~~~~~~~~~
1190         orca R1, R2 or R4 configuration: reduced in the surface Eq. strip
1191      file   :
1192 ahmcoef                                                                       
1193   open ok
1194      unit   =           21
1195      status = OLD
1196      form   = FORMATTED
1197      access = SEQUENTIAL
1198 
1199 
1200          1D zcoef array
1201          ~~~~~~~~~~~~~~
1202 
1203     jk        zcoef
1204      1        1.000
1205      2        1.000
1206      3        1.000
1207      4        1.000
1208      5        1.000
1209      6        1.000
1210      7        1.000
1211      8        1.000
1212      9        1.000
1213     10        1.000
1214     11        1.000
1215     12        1.000
1216     13        1.000
1217     14        1.000
1218     15        1.000
1219     16        1.000
1220     17        1.000
1221     18        1.000
1222     19        1.000
1223     20        1.000
1224     21        1.000
1225     22        2.000
1226     23        3.000
1227     24        5.000
1228     25        7.000
1229     26        9.000
1230     27       10.000
1231     28       10.000
1232     29       10.000
1233     30       10.000
1234     31       10.000
1235     32       10.000
1236     33       10.000
1237     34       10.000
1238     35       10.000
1239     36       10.000
1240     37       10.000
1241     38       10.000
1242     39       10.000
1243     40       10.000
1244     41       10.000
1245     42       10.000
1246     43       10.000
1247     44       10.000
1248     45       10.000
1249     46       10.000
1250 
1251          3D ahm1 array (k=1)
1252
1253
1254        1
1255 147   1.00
1256 127   4.41
1257 107  20.00
1258  87  20.00
1259  67  20.00
1260  47  20.00
1261  27  20.00
1262   7  20.00
1263 
1264          3D ahm2 array (k=1)
1265
1266
1267        1
1268 147   1.00
1269 127   4.08
1270 107  20.00
1271  87  20.00
1272  67  20.00
1273  47  20.00
1274  27  20.00
1275   7  20.00
1276 
1277          3D ahm2 array (k=jpk)
1278
1279
1280        1
1281 147  10.00
1282 127  20.00
1283 107  20.00
1284  87  20.00
1285  67  20.00
1286  47  20.00
1287  27  20.00
1288   7  20.00
1289 
1290 ldf_slp_init : direction of lateral mixing
1291 ~~~~~~~~~~~~
1292           Horizontal mixing in s-coordinate: slope = slope of s-surfaces
1293 
1294 tra_qsr_init : penetration of the surface solar radiation
1295 ~~~~~~~~~~~~
1296    Namelist namtra_qsr : set the parameter of penetration
1297       Light penetration (T) or not (F)         ln_traqsr  =  T
1298       RGB (Red-Green-Blue) light penetration   ln_qsr_rgb =  T
1299       2 band               light penetration   ln_qsr_2bd =  F
1300       bio-model            light penetration   ln_qsr_bio =  F
1301       light penetration for ice-model LIM3     ln_qsr_ice =  T
1302       RGB : Chl data (=1) or cst value (=0)    nn_chldta  =            0
1303       RGB & 2 bands: fraction of light (rn_si1)    rn_abs =
1304  0.580000000000000     
1305       RGB & 2 bands: shortess depth of extinction  rn_si0 =
1306  0.350000000000000     
1307       2 bands: longest depth of extinction         rn_si1 =
1308   23.0000000000000     
1309       light penetration for ice-model LIM3     ln_qsr_ice =  T
1310 
1311          R-G-B   light penetration - Constant Chlorophyll
1312 
1313 trc_oce_rgb : Initialisation of the optical look-up table
1314 ~~~~~~~~~~~
1315       RGB longest depth of extinction    r_si2 =    61.8046971569839     
1316     jc =           1   Chl =   1.000000000000000E-002   irgb =            1
1317     jc =           2   Chl =   1.100000000000000E-002   irgb =            2
1318     jc =           3   Chl =   1.300000000000000E-002   irgb =            3
1319     jc =           4   Chl =   1.400000000000000E-002   irgb =            4
1320     jc =           5   Chl =   1.600000000000000E-002   irgb =            5
1321     jc =           6   Chl =   1.800000000000000E-002   irgb =            6
1322     jc =           7   Chl =   2.000000000000000E-002   irgb =            7
1323     jc =           8   Chl =   2.200000000000000E-002   irgb =            8
1324     jc =           9   Chl =   2.500000000000000E-002   irgb =            9
1325     jc =          10   Chl =   2.800000000000000E-002   irgb =           10
1326     jc =          11   Chl =   3.200000000000000E-002   irgb =           11
1327     jc =          12   Chl =   3.500000000000000E-002   irgb =           12
1328     jc =          13   Chl =   4.000000000000000E-002   irgb =           13
1329     jc =          14   Chl =   4.500000000000000E-002   irgb =           14
1330     jc =          15   Chl =   5.000000000000000E-002   irgb =           15
1331     jc =          16   Chl =   5.600000000000000E-002   irgb =           16
1332     jc =          17   Chl =   6.300000000000000E-002   irgb =           17
1333     jc =          18   Chl =   7.099999999999999E-002   irgb =           18
1334     jc =          19   Chl =   7.900000000000000E-002   irgb =           19
1335     jc =          20   Chl =   8.900000000000000E-002   irgb =           20
1336     jc =          21   Chl =   0.100000000000000        irgb =           21
1337     jc =          22   Chl =   0.112000000000000        irgb =           22
1338     jc =          23   Chl =   0.126000000000000        irgb =           23
1339     jc =          24   Chl =   0.141000000000000        irgb =           24
1340     jc =          25   Chl =   0.158000000000000        irgb =           25
1341     jc =          26   Chl =   0.178000000000000        irgb =           26
1342     jc =          27   Chl =   0.200000000000000        irgb =           27
1343     jc =          28   Chl =   0.224000000000000        irgb =           28
1344     jc =          29   Chl =   0.251000000000000        irgb =           29
1345     jc =          30   Chl =   0.282000000000000        irgb =           30
1346     jc =          31   Chl =   0.316000000000000        irgb =           31
1347     jc =          32   Chl =   0.355000000000000        irgb =           32
1348     jc =          33   Chl =   0.398000000000000        irgb =           33
1349     jc =          34   Chl =   0.447000000000000        irgb =           34
1350     jc =          35   Chl =   0.501000000000000        irgb =           35
1351     jc =          36   Chl =   0.562000000000000        irgb =           36
1352     jc =          37   Chl =   0.631000000000000        irgb =           37
1353     jc =          38   Chl =   0.708000000000000        irgb =           38
1354     jc =          39   Chl =   0.794000000000000        irgb =           39
1355     jc =          40   Chl =   0.891000000000000        irgb =           40
1356     jc =          41   Chl =    1.00000000000000        irgb =           41
1357     jc =          42   Chl =    1.12200000000000        irgb =           42
1358     jc =          43   Chl =    1.25900000000000        irgb =           43
1359     jc =          44   Chl =    1.41300000000000        irgb =           44
1360     jc =          45   Chl =    1.58500000000000        irgb =           45
1361     jc =          46   Chl =    1.77800000000000        irgb =           46
1362     jc =          47   Chl =    1.99500000000000        irgb =           47
1363     jc =          48   Chl =    2.23900000000000        irgb =           48
1364     jc =          49   Chl =    2.51200000000000        irgb =           49
1365     jc =          50   Chl =    2.81800000000000        irgb =           50
1366     jc =          51   Chl =    3.16200000000000        irgb =           51
1367     jc =          52   Chl =    3.54800000000000        irgb =           52
1368     jc =          53   Chl =    3.98100000000000        irgb =           53
1369     jc =          54   Chl =    4.46700000000000        irgb =           54
1370     jc =          55   Chl =    5.01200000000000        irgb =           55
1371     jc =          56   Chl =    5.62300000000000        irgb =           56
1372     jc =          57   Chl =    6.31000000000000        irgb =           57
1373     jc =          58   Chl =    7.07900000000000        irgb =           58
1374     jc =          59   Chl =    7.94300000000000        irgb =           59
1375     jc =          60   Chl =    8.91200000000000        irgb =           60
1376     jc =          61   Chl =    10.0000000000000        irgb =           61
1377         level of light extinction =           25  ref depth =
1378   2527.21691545123       m
1379 
1380         Constant Chlorophyll concentration = 0.05
1381         fixed volume: light distribution computed one for all
1382 
1383 tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating
1384 ~~~~~~~   
1385    Namelist nambbc : set bbc parameters
1386       Apply a geothermal heating at ocean bottom   ln_trabbc     =  T
1387       type of geothermal flux                      nn_geoflx     =            1
1388       Constant geothermal flux value               rn_geoflx_cst =
1389  8.640000000000000E-002
1390 
1391       *** constant heat flux  =     8.640000000000000E-002
1392 
1393 tra_bbl_init : bottom boundary layer initialisation
1394 ~~~~~~~~~~~~
1395        Namelist nambbl : set bbl parameters
1396           diffusive bbl (=1)   or not (=0)    nn_bbl_ldf =            1
1397           advective bbl (=1/2) or not (=0)    nn_bbl_adv =            0
1398           diffusive bbl coefficient           rn_ahtbbl  =
1399   1000.00000000000       m2/s
1400           advective bbl coefficient           rn_gambbl  =
1401   10.0000000000000       s
1402 
1403 tra_dmp_init : T and S newtonian damping
1404 ~~~~~~~
1405    Namelist namtra_dmp : set damping parameter
1406       add a damping term or not       ln_tradmp =  F
1407       T and S damping option          nn_hdmp   =           -1
1408       mixed layer damping option      nn_zdmp   =            0
1409 (non-C1D zoom: forced to 0)
1410       surface time scale (days)       rn_surf   =    50.0000000000000     
1411       bottom time scale (days)        rn_bot    =    360.000000000000     
1412       depth of transition (meters)    rn_dep    =    800.000000000000     
1413       create a damping.coeff file     nn_file   =            0
1414 
1415 
1416 tra_adv_init : choice/control of the tracer advection scheme
1417 ~~~~~~~~~~~
1418    Namelist namtra_adv : chose a advection scheme for tracers
1419       2nd order advection scheme     ln_traadv_cen2    =  F
1420       TVD advection scheme           ln_traadv_tvd     =  T
1421       MUSCL  advection scheme        ln_traadv_muscl   =  F
1422       MUSCL2 advection scheme        ln_traadv_muscl2  =  F
1423       UBS    advection scheme        ln_traadv_ubs     =  F
1424       QUICKEST advection scheme      ln_traadv_qck     =  F
1425       upstream scheme within muscl   ln_traadv_msc_ups =  F
1426 
1427          TVD       scheme is used
1428 
1429 tra_adv_mle_init : mixed layer eddy (MLE) advection acting on tracers
1430 ~~~~~~~~~~~~~~~~
1431    Namelist namtra_adv_mle : mixed layer eddy advection on tracers
1432       use mixed layer eddy (MLE, i.e. Fox-Kemper param) (T/F)      ln_mle    =
1433 F
1434       MLE type: =0 standard Fox-Kemper ; =1 new formulation        nn_mle    =
1435           1
1436       magnitude of the MLE (typical value: 0.06 to 0.08)           rn_ce     =
1437  6.000000000000000E-002
1438       scale of ML front (ML radius of deformation) (rn_mle=0)      rn_lf     =
1439   5000.00000000000      m
1440       maximum time scale of MLE                    (rn_mle=0)      rn_time   =
1441   172800.000000000      s
1442       reference latitude (degrees) of MLE coef.    (rn_mle=1)      rn_lat    =
1443   20.0000000000000      deg
1444       space interp. of MLD at u-(v-)pts (0=min,1=averaged,2=max)   nn_mld_uv =
1445           0
1446       =1 no MLE in case of convection ; =0 always MLE              nn_conv   =
1447           0
1448       Density difference used to define ML for FK              rn_rho_c_mle  =
1449  1.000000000000000E-002
1450 
1451    Mixed Layer Eddy parametrisation NOT used
1452 
1453 tra_ldf_init : lateral tracer diffusive operator
1454 ~~~~~~~~~~~
1455    Namelist namtra_ldf already read in ldftra module
1456    see ldf_tra_init report for lateral mixing parameters
1457 
1458 
1459           Rotated laplacian operator
1460 
1461 tra:ldf_ano : lateral diffusion acting on the full fields
1462 ~~~~~~~~~~~
1463 
1464 tra_zdf_init : vertical tracer physics scheme
1465 ~~~~~~~~~~~
1466               Implicit (euler backward) scheme
1467 
1468 dyn_adv_init : choice/control of the momentum advection scheme
1469 ~~~~~~~~~~~
1470        Namelist namdyn_adv : chose a advection formulation & scheme for momentu
1471 m
1472           Vector/flux form (T/F)             ln_dynadv_vec  =  T
1473           2nd order centred advection scheme ln_dynadv_cen2 =  F
1474           3rd order UBS advection scheme     ln_dynadv_ubs  =  F
1475 
1476          vector form : keg + zad + vor is used
1477 
1478 dyn_vor_init : vorticity term : read namelist and control the consistency
1479 ~~~~~~~~~~~~
1480         Namelist namdyn_vor : choice of the vorticity term scheme
1481            energy    conserving scheme                ln_dynvor_ene =  F
1482            enstrophy conserving scheme                ln_dynvor_ens =  F
1483            mixed enstrophy/energy conserving scheme   ln_dynvor_mix =  F
1484            enstrophy and energy conserving scheme     ln_dynvor_een =  T
1485 
1486          Vector form advection : vorticity = Coriolis + relative vorticity
1487 
1488          vorticity scheme : energy and enstrophy conserving scheme
1489 
1490 dyn_ldf_init : Choice of the lateral diffusive operator on dynamics
1491 ~~~~~~~~~~~
1492        Namelist nam_dynldf : set lateral mixing parameters (type, direction, co
1493 efficients)
1494           laplacian operator          ln_dynldf_lap   =  T
1495           bilaplacian operator        ln_dynldf_bilap =  F
1496           iso-level                   ln_dynldf_level =  F
1497           horizontal (geopotential)   ln_dynldf_hor   =  T
1498           iso-neutral                 ln_dynldf_iso   =  F
1499 
1500               laplacian operator
1501 
1502 dyn_hpg_init : hydrostatic pressure gradient initialisation
1503 ~~~~~~~~~~~~
1504    Namelist namdyn_hpg : choice of hpg scheme
1505       z-coord. - full steps                             ln_hpg_zco    =  F
1506       z-coord. - partial steps (interpolation)          ln_hpg_zps    =  T
1507       s-coord. (standard jacobian formulation)          ln_hpg_sco    =  F
1508       s-coord. (Density Jacobian: Cubic polynomial)     ln_hpg_djc    =  F
1509       s-coord. (Pressure Jacobian: Cubic polynomial)    ln_hpg_prj    =  F
1510       time stepping: centered (F) or semi-implicit (T)  ln_dynhpg_imp =  F
1511 
1512 dyn_zdf_init : vertical dynamics physics scheme
1513 ~~~~~~~~~~~
1514               Implicit (euler backward) scheme
1515 
1516 dyn_spg_init : choice of the surface pressure gradient scheme
1517 ~~~~~~~~~~~
1518      Explicit free surface                  lk_dynspg_exp =  F
1519      Free surface with time splitting       lk_dynspg_ts  =  F
1520      Filtered free surface cst volume       lk_dynspg_flt =  T
1521 
1522      filtered free surface
1523      file   :
1524 solver.stat                                                                   
1525   open ok
1526      unit   =           22
1527      status = REPLACE
1528      form   = FORMATTED
1529      access = SEQUENTIAL
1530 
1531 
1532 solver_init : solver to compute the surface pressure gradient
1533 ~~~~~~~~~~~
1534    Namelist namsol : set solver parameters
1535       type of elliptic solver            nn_solv    =            1
1536       absolute/relative (0/1) precision  nn_sol_arp =            0
1537       minimum iterations for solver      nn_nmin    =          300
1538       maximum iterations for solver      nn_nmax    =          800
1539       frequency for test                 nn_nmod    =           10
1540       absolute precision of solver       rn_eps     =   1.000000000000000E-006
1541       absolute precision for SOR solver  rn_resmax  =   1.000000000000000E-010
1542       optimal coefficient of sor         rn_sor     =    1.92000000000000     
1543 
1544    a preconditioned conjugate gradient solver is used
1545 
1546 icbini :   Namelist namberg ln_icebergs = F , NO icebergs used
1547 ~~~~~~~~
1548 
1549 dia_ptr_init : poleward transport and msf initialization
1550 ~~~~~~~~~~~~
1551    Namelist namptr : set ptr parameters
1552       Poleward heat & salt transport (T) or not (F)      ln_diaptr  =  F
1553       Overturning heat & salt transport                  ln_ptrcomp =  T
1554       T & S zonal mean and meridional stream function    ln_diaznl  =  F
1555       Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins      ln_subbas  =  F
1556       Frequency of computation                           nn_fptr    =
1557           1
1558       Frequency of outputs                               nn_fwri    =
1559          15
1560 
1561 dia_hsb_init : check the heat and salt budgets
1562 ~~~~~~~~
1563 
1564 dia_hsb_init : check the heat and salt budgets
1565 ~~~~~~~~~~~~
1566    Namelist namhsb : set hsb parameters
1567       Switch for hsb diagnostic (T) or not (F)  ln_diahsb  =  F
1568 
1569 Euler time step switch is            0
1570
1571AAAAAAAA
1572
1573 
1574  rot_rep : geographic <--> stretched
1575  ~~~~~    coordinate transformation
1576 
1577 sbc_ice_lim : update ocean surface boudary condition
1578 ~~~~~~~~~~~   via Louvain la Neuve Ice Model (LIM-3) time stepping
1579      file   :
1580 namelist_ice_ref                                                               
1581   open ok
1582      unit   =           23
1583      status = OLD
1584      form   = FORMATTED
1585      access = SEQUENTIAL
1586 
1587      file   :
1588 namelist_ice_cfg                                                               
1589   open ok
1590      unit   =           24
1591      status = OLD
1592      form   = FORMATTED
1593      access = SEQUENTIAL
1594 
1595      file   :
1596 output.namelist.ice                                                           
1597   open ok
1598      unit   =           25
1599      status = UNKNOWN
1600      form   = FORMATTED
1601      access = SEQUENTIAL
1602 
1603 
1604 ice_run : ice share parameters for dynamics/advection/thermo of sea-ice
1605  ~~~~~~
1606    switch for ice dynamics (1) or not (0)      ln_limdyn   =  T
1607    maximum ice concentration                               =
1608  0.999000000000000     
1609    atmospheric drag over sea ice                           =
1610  1.400000000000000E-003
1611    atmospheric drag over ocean                             =
1612  1.000000000000000E-003
1613    Several ice points in the ice or not in ocean.output    =  F
1614    Diagnose heat/salt budget or not          ln_limdiahsb  =  F
1615    Output   heat/salt budget or not          ln_limdiaout  =  F
1616 
1617 lim_thd : Ice Thermodynamics
1618 ~~~~~~~
1619 
1620    Namelist of ice parameters for ice thermodynamic computation
1621       maximum melting at the bottom                           hmelt        =
1622 -0.150000000000000     
1623       ice thick. for lateral accretion                        hiccrit      =
1624  0.100000000000000     
1625       Frazil ice thickness as a function of wind or not       fraz_swi     =
1626           0
1627       Maximum proportion of frazil ice collecting at bottom   maxfrazb     =
1628  0.000000000000000E+000
1629       Thresold relative drift speed for collection of frazil  vfrazb       =
1630  0.416666700000000     
1631       Squeezing coefficient for collection of frazil          Cfrazb       =
1632   5.00000000000000     
1633       minimum ice thickness                                   hiclim       =
1634  0.100000000000000     
1635       numerical carac. of the scheme for diffusion in ice
1636       thickness of the surf. layer in temp. computation       hnzst        =
1637  0.100000000000000     
1638       switch for snow sublimation  (=1) or not (=0)           parsub       =
1639   1.00000000000000     
1640       coefficient for ice-lead partition of snowfall          betas        =
1641  0.600000000000000     
1642       extinction radiation parameter in sea ice (1.0)         kappa_i      =
1643   1.00000000000000     
1644       maximal n. of iter. for heat diffusion computation      nconv_i_thd  =
1645   50.0000000000000     
1646       maximal err. on T for heat diffusion computation        maxer_i_thd  =
1647  1.000000000000000E-004
1648       switch for comp. of thermal conductivity in the ice     thcon_i_swi  =
1649           1
1650       check heat conservation in the ice/snow                 con_i        =  F
1651 
1652 lim_thd_sal_init : Ice parameters for salinity
1653 ~~~~~~~~~~~~~~~~
1654  switch for salinity num_sal        :            2
1655  bulk salinity value if num_sal = 1 :    4.00000000000000     
1656  restoring salinity for GD          :    5.00000000000000     
1657  restoring time for GD              :    1728000.00000000     
1658  restoring salinity for flushing    :    2.00000000000000     
1659  restoring time for flushing        :    864000.000000000     
1660  Maximum tolerated ice salinity     :    20.0000000000000     
1661  Minimum tolerated ice salinity     :   0.100000000000000     
1662  1st salinity for salinity profile  :    3.50000000000000     
1663  2nd salinity for salinity profile  :    4.50000000000000     
1664 
1665 lim_msh : LIM-3 sea-ice model, mesh initialization
1666 ~~~~~~~
1667           the model domain is entirely in the southern hemisphere: njeq =
1668         147
1669 
1670 lim_itd_ini : Initialization of ice thickness distribution
1671 ~~~~~~~~~~~~
1672  Number of ice types jpm =                 1
1673  Number of ice categories jpl =            5
1674  Ice type            1
1675  Number of thickness categories            5
1676  Thickness category boundaries             1           5
1677 Ice type vector           1           1           1           1           1
1678 
1679  Thickness category boundaries independently of ice type
1680  hi_max   0.000000000000000E+000  0.629671477879617     
1681   1.32762233171754        2.24678625463993        3.84485861256699     
1682   7.22256121057878     
1683  Thickness category boundaries inside ice types
1684  Type number            1
1685  hi_max_typ :   0.000000000000000E+000  0.629671477879617     
1686   1.32762233171754        2.24678625463993        3.84485861256699     
1687   7.22256121057878     
1688 
1689 lim_istate : Ice initialization
1690 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1691 
1692 lim_istate_init : ice parameters inititialisation
1693 ~~~~~~~~~~~~~~~
1694    initialization with ice (T) or not (F)       ln_limini   =  F
1695    threshold water temp. for initial sea-ice    thres_sst  =
1696  0.000000000000000E+000
1697    initial snow thickness in the north          hts_ini_n  =
1698  0.300000000000000     
1699    initial snow thickness in the south          hts_ini_s  =
1700  0.300000000000000     
1701    initial ice thickness  in the north          hti_ini_n  =
1702   1.00000000000000     
1703    initial ice thickness  in the south          hti_ini_s  =
1704   1.00000000000000     
1705    initial ice concentr.  in the north          ati_ini_n  =
1706  0.900000000000000     
1707    initial ice concentr.  in the north          ati_ini_s  =
1708  0.900000000000000     
1709    initial  ice salinity  in the north          smi_ini_n  =
1710   6.30100000000000     
1711    initial  ice salinity  in the south          smi_ini_s  =
1712   6.30100000000000     
1713    initial  ice/snw temp  in the north          tmi_ini_n  =
1714   270.000000000000     
1715    initial  ice/snw temp  in the south          tmi_ini_s  =
1716   270.000000000000     
1717 
1718 lim_sbc_init : LIM-3 sea-ice - surface boundary condition
1719 ~~~~~~~~~~~~~   
1720 
1721 albedo : set albedo parameters
1722 ~~~~~~~
1723    Namelist namsbc_alb : albedo
1724       correction for snow and ice albedo                  rn_cloud  =
1725  6.000000000000000E-002
1726       albedo of melting ice in the arctic and antarctic   rn_albice =
1727  0.530000000000000     
1728       coefficients for linear                             rn_alphd  =
1729  0.800000000000000     
1730       interpolation used to compute albedo                rn_alphdi =
1731  0.720000000000000     
1732       between two extremes values (Pyane, 1972)           rn_alphc  =
1733  0.650000000000000     
1734 
1735 lim_dyn_init : ice parameters for ice dynamics
1736 ~~~~~~~~~~~~
1737    tolerance parameter                              epsd   =
1738  9.999999999999999E-021
1739    relaxation constant                              om     =
1740  0.500000000000000     
1741    drag coefficient for oceanic stress              cw     =
1742  5.000000000000000E-003
1743    turning angle for oceanic stress                 angvg  =
1744  0.000000000000000E+000
1745    first bulk-rheology parameter                    pstar  =
1746   20000.0000000000     
1747    second bulk-rhelogy parameter                    c_rhg  =
1748   20.0000000000000     
1749    creep limit                                      creepl =
1750  1.000000000000000E-012
1751    eccentricity of the elliptical yield curve       ecc    =
1752   2.00000000000000     
1753    horizontal diffusivity coeff. for sea-ice        ahi0   =
1754   350.000000000000     
1755    number of iterations for subcycling              nevp   =          300
1756    timescale for elastic waves                      telast =
1757   9600.00000000000     
1758    coefficient for the solution of int. stresses  alphaevp =
1759   1.00000000000000     
1760    min ice thickness for rheology calculations     hminrhg =
1761  1.000000000000000E-003
1762 
1763  lim_itd_me_init : ice parameters for mechanical ice redistribution
1764  ~~~~~~~~~~~~~~~
1765    Switch choosing the ice redistribution scheme           ridge_scheme_swi
1766           0
1767    Fraction of shear energy contributing to ridging        Cs             
1768  0.500000000000000     
1769    Ratio of ridging work to PotEner change in ridging      Cf             
1770   17.0000000000000     
1771    Fraction of snow volume conserved during ridging        fsnowrdg       
1772  0.500000000000000     
1773    Fraction of snow volume conserved during ridging        fsnowrft       
1774  0.500000000000000     
1775    Fraction of total ice coverage contributing to ridging  Gstar           
1776  0.150000000000000     
1777    Equivalent to G* for an exponential part function       astar           
1778  5.000000000000000E-002
1779    Quantity playing a role in max ridged ice thickness     Hstar           
1780   100.000000000000     
1781    Rafting of ice sheets or not                            raft_swi       
1782           1
1783    Parmeter thickness (threshold between ridge-raft)       hparmeter       
1784  0.750000000000000     
1785    Rafting hyperbolic tangent coefficient                  Craft           
1786   5.00000000000000     
1787    Initial porosity of ridges                              ridge_por       
1788  0.300000000000000     
1789    Switch for part. function (0) linear (1) exponential    partfun_swi     
1790           1
1791    Switch for including brine volume in ice strength comp. brinstren_swi   
1792           0
1793 
1794 lim_trp : Ice transport
1795 ~~~~~~~~~~~~