CNRS - NERC - Met OFFICE - MERCATOR-ocean - INGV - CMCC NEMO team Ocean General Circulation Model version 3.6 (2015) mynode : mpi initialisation ~~~~~~ Namelist nammpp mpi send type cn_mpi_send = I size in bytes of exported buffer nn_buffer = 0 processor grid extent in i jpni = 22 processor grid extent in j jpnj = 22 number of local domains jpnij = 360 avoid use of mpi_allgather at the north fold ln_nnogather = T Immediate non-blocking send (isend) AAAAAAAA nemo_ctl: Control prints & Benchmark ~~~~~~~ Namelist namctl run control (for debugging) ln_ctl = F level of print nn_print = 0 Start i indice for SUM control nn_ictls = 0 End i indice for SUM control nn_ictle = 0 Start j indice for SUM control nn_jctls = 0 End j indice for SUM control nn_jctle = 0 number of proc. following i nn_isplt = 1 number of proc. following j nn_jsplt = 1 benchmark parameter (0/1) nn_bench = 0 timing activated (0/1) nn_timing = 0 namcfg : configuration initialization through namelist read ~~~~~~~ Namelist namcfg configuration name cp_cfg = orca configuration zoom name cp_cfz = no zoom configuration resolution jp_cfg = 1 1st lateral dimension ( >= jpi ) jpidta = 362 2nd " " ( >= jpj ) jpjdta = 332 3nd " " jpkdta = 75 1st dimension of global domain in i jpiglo = 362 2nd - - in j jpjglo = 332 left bottom i index of the zoom (in data domain) jpizoom = 1 left bottom j index of the zoom (in data domain) jpizoom = 1 lateral cond. type (between 0 and 6) jperio = 6 use file attribute if exists as i/p j-start ln_use_jattr = F mpp_init : Message Passing MPI ~~~~~~~~ iom_nf90_open ~~~ open existing file: bathy_meter.nc in REA D mode ---> bathy_meter.nc OK ===>>> : W A R N I N G =============== iom_nf90_getatt: no attribute open_ocean_jstart found read Bathymetry (rec: 1) in bathy_meter.nc ok iom_close ~~~ close file: bathy_meter.nc ok mpp_init2: defines mpp subdomains ~~~~~~ ---------------------- iresti= 8 irestj= 22 jpni= 22 jpnj= 22 sum ilci(i,1)= 362.000000000000 jpiglo= 362 sum ilcj(1,j)= 332.000000000000 jpjglo= 332 **************************************************************** * * * * * 22 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * 345 * 346 * * * * * * **************************************************************** * * * * * 21 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * 333 * * * * * * **************************************************************** * * * * * 20 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 19 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * 306 * * * * * * **************************************************************** * * * * * 18 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * 289 * * * * * * **************************************************************** * * * * * 17 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * 271 * 272 * * * * * * **************************************************************** * * * * * 16 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * 251 * 252 * * * * * * **************************************************************** * * * * * 15 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 231 * 232 * 233 * 234 * * * * * * **************************************************************** * * * * * 14 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 210 * 211 * 212 * 213 * * * * * * **************************************************************** * * * * * 13 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 189 * 190 * 191 * 192 * * * * * * **************************************************************** * * * * * 12 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 168 * 169 * 170 * 171 * * * * * * **************************************************************** * * * * * 11 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 146 * 147 * 148 * 149 * * * * * * **************************************************************** * * * * * 10 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 124 * 125 * 126 * 127 * * * * * * **************************************************************** * * * * * 9 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 102 * 103 * 104 * 105 * * * * * * **************************************************************** * * * * * 8 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 80 * 81 * 82 * 83 * * * * * * **************************************************************** * * * * * 7 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 58 * 59 * 60 * 61 * * * * * * **************************************************************** * * * * * 6 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 36 * 37 * 38 * 39 * * * * * * **************************************************************** * * * * * 5 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 17 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 4 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 3 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 2 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 1 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** 1 2 3 4 **************************************************************** * * * * * 22 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 347 * 348 * 349 * 350 * * * * * * **************************************************************** * * * * * 21 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 334 * 335 * 336 * 337 * * * * * * **************************************************************** * * * * * 20 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 321 * 322 * 323 * 324 * * * * * * **************************************************************** * * * * * 19 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 307 * 308 * 309 * 310 * * * * * * **************************************************************** * * * * * 18 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 290 * 291 * 292 * 293 * * * * * * **************************************************************** * * * * * 17 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 273 * 274 * 275 * 276 * * * * * * **************************************************************** * * * * * 16 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 253 * 254 * 255 * 256 * * * * * * **************************************************************** * * * * * 15 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 235 * 236 * 237 * 238 * * * * * * **************************************************************** * * * * * 14 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 214 * 215 * 216 * 217 * * * * * * **************************************************************** * * * * * 13 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 193 * 194 * 195 * 196 * * * * * * **************************************************************** * * * * * 12 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 172 * 173 * 174 * 175 * * * * * * **************************************************************** * * * * * 11 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 150 * 151 * 152 * 153 * * * * * * **************************************************************** * * * * * 10 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 128 * 129 * 130 * 131 * * * * * * **************************************************************** * * * * * 9 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 106 * 107 * 108 * 109 * * * * * * **************************************************************** * * * * * 8 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 84 * 85 * 86 * 87 * * * * * * **************************************************************** * * * * * 7 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 62 * 63 * 64 * 65 * * * * * * **************************************************************** * * * * * 6 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 40 * 41 * 42 * 43 * * * * * * **************************************************************** * * * * * 5 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * 18 * 19 * 20 * 21 * * * * * * **************************************************************** * * * * * 4 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * 5 * 6 * 7 * * * * * * **************************************************************** * * * * * 3 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * 0 * 1 * 2 * * * * * * **************************************************************** * * * * * 2 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 1 * 19 x 17 * 19 x 17 * 19 x 17 * 19 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** 5 6 7 8 **************************************************************** * * * * * 22 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 351 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 21 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 338 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 20 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * 325 * * * * * * **************************************************************** * * * * * 19 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 311 * 312 * -1 * 313 * * * * * * **************************************************************** * * * * * 18 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 294 * 295 * -1 * 296 * * * * * * **************************************************************** * * * * * 17 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 277 * 278 * -1 * 279 * * * * * * **************************************************************** * * * * * 16 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 257 * 258 * 259 * 260 * * * * * * **************************************************************** * * * * * 15 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 239 * 240 * 241 * 242 * * * * * * **************************************************************** * * * * * 14 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 218 * 219 * 220 * 221 * * * * * * **************************************************************** * * * * * 13 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 197 * 198 * 199 * 200 * * * * * * **************************************************************** * * * * * 12 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 176 * 177 * 178 * 179 * * * * * * **************************************************************** * * * * * 11 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 154 * 155 * 156 * 157 * * * * * * **************************************************************** * * * * * 10 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 132 * 133 * 134 * 135 * * * * * * **************************************************************** * * * * * 9 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 110 * 111 * 112 * 113 * * * * * * **************************************************************** * * * * * 8 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 88 * 89 * 90 * 91 * * * * * * **************************************************************** * * * * * 7 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 66 * 67 * 68 * 69 * * * * * * **************************************************************** * * * * * 6 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 44 * 45 * 46 * 47 * * * * * * **************************************************************** * * * * * 5 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 22 * 23 * 24 * 25 * * * * * * **************************************************************** * * * * * 4 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 8 * 9 * 10 * 11 * * * * * * **************************************************************** * * * * * 3 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 2 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 1 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** 9 10 11 12 **************************************************************** * * * * * 22 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 352 * 353 * 354 * 355 * * * * * * **************************************************************** * * * * * 21 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * 339 * 340 * 341 * * * * * * **************************************************************** * * * * * 20 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 326 * 327 * 328 * 329 * * * * * * **************************************************************** * * * * * 19 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 314 * 315 * 316 * 317 * * * * * * **************************************************************** * * * * * 18 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 297 * 298 * 299 * 300 * * * * * * **************************************************************** * * * * * 17 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 280 * 281 * 282 * 283 * * * * * * **************************************************************** * * * * * 16 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 261 * 262 * 263 * 264 * * * * * * **************************************************************** * * * * * 15 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 243 * 244 * 245 * 246 * * * * * * **************************************************************** * * * * * 14 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 222 * 223 * 224 * 225 * * * * * * **************************************************************** * * * * * 13 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 201 * -1 * 202 * 203 * * * * * * **************************************************************** * * * * * 12 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 180 * -1 * 181 * 182 * * * * * * **************************************************************** * * * * * 11 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 158 * 159 * 160 * 161 * * * * * * **************************************************************** * * * * * 10 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 136 * 137 * 138 * 139 * * * * * * **************************************************************** * * * * * 9 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 114 * 115 * 116 * 117 * * * * * * **************************************************************** * * * * * 8 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 92 * 93 * 94 * 95 * * * * * * **************************************************************** * * * * * 7 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 70 * 71 * 72 * 73 * * * * * * **************************************************************** * * * * * 6 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 48 * 49 * 50 * 51 * * * * * * **************************************************************** * * * * * 5 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 26 * 27 * 28 * 29 * * * * * * **************************************************************** * * * * * 4 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 12 * 13 * 14 * 15 * * * * * * **************************************************************** * * * * * 3 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * 3 * 4 * * * * * * **************************************************************** * * * * * 2 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 1 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** 13 14 15 16 **************************************************************** * * * * * 22 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 356 * 357 * 358 * 359 * * * * * * **************************************************************** * * * * * 21 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 342 * 343 * 344 * -1 * * * * * * **************************************************************** * * * * * 20 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 330 * 331 * 332 * -1 * * * * * * **************************************************************** * * * * * 19 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 318 * 319 * 320 * -1 * * * * * * **************************************************************** * * * * * 18 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 301 * 302 * 303 * 304 * * * * * * **************************************************************** * * * * * 17 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 284 * 285 * 286 * 287 * * * * * * **************************************************************** * * * * * 16 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 265 * 266 * 267 * 268 * * * * * * **************************************************************** * * * * * 15 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 247 * -1 * -1 * 248 * * * * * * **************************************************************** * * * * * 14 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 226 * 227 * 228 * -1 * * * * * * **************************************************************** * * * * * 13 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 204 * 205 * 206 * 207 * * * * * * **************************************************************** * * * * * 12 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 183 * 184 * 185 * 186 * * * * * * **************************************************************** * * * * * 11 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 162 * 163 * 164 * 165 * * * * * * **************************************************************** * * * * * 10 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 140 * 141 * 142 * 143 * * * * * * **************************************************************** * * * * * 9 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 118 * 119 * 120 * 121 * * * * * * **************************************************************** * * * * * 8 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 96 * 97 * 98 * 99 * * * * * * **************************************************************** * * * * * 7 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 74 * 75 * 76 * 77 * * * * * * **************************************************************** * * * * * 6 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 52 * 53 * 54 * 55 * * * * * * **************************************************************** * * * * * 5 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 30 * 31 * 32 * 33 * * * * * * **************************************************************** * * * * * 4 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * 16 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 3 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 2 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** * * * * * 1 * 18 x 17 * 18 x 17 * 18 x 17 * 18 x 17 * * -1 * -1 * -1 * -1 * * * * * * **************************************************************** 17 18 19 20 ******************************** * * * 22 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 21 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 20 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 19 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 18 * 18 x 17 * 18 x 17 * * 305 * -1 * * * * ******************************** * * * 17 * 18 x 17 * 18 x 17 * * 288 * -1 * * * * ******************************** * * * 16 * 18 x 17 * 18 x 17 * * 269 * 270 * * * * ******************************** * * * 15 * 18 x 17 * 18 x 17 * * 249 * 250 * * * * ******************************** * * * 14 * 18 x 17 * 18 x 17 * * 229 * 230 * * * * ******************************** * * * 13 * 18 x 17 * 18 x 17 * * 208 * 209 * * * * ******************************** * * * 12 * 18 x 17 * 18 x 17 * * 187 * 188 * * * * ******************************** * * * 11 * 18 x 17 * 18 x 17 * * 166 * 167 * * * * ******************************** * * * 10 * 18 x 17 * 18 x 17 * * 144 * 145 * * * * ******************************** * * * 9 * 18 x 17 * 18 x 17 * * 122 * 123 * * * * ******************************** * * * 8 * 18 x 17 * 18 x 17 * * 100 * 101 * * * * ******************************** * * * 7 * 18 x 17 * 18 x 17 * * 78 * 79 * * * * ******************************** * * * 6 * 18 x 17 * 18 x 17 * * 56 * 57 * * * * ******************************** * * * 5 * 18 x 17 * 18 x 17 * * 34 * 35 * * * * ******************************** * * * 4 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 3 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 2 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** * * * 1 * 18 x 17 * 18 x 17 * * -1 * -1 * * * * ******************************** 21 22 nproc = 0 nowe = -1 noea = 1 nono = 5 noso = -1 nbondi = -1 nbondj = -1 npolj = 0 nperio = 0 nlci = 19 nlcj = 17 nimpp = 86 njmpp = 31 nreci = 2 npse = -1 nrecj = 2 npsw = -1 jpreci = 1 npne = 6 jprecj = 1 npnw = -1 mpp_init2 : North fold boundary prepared for jpni >1 mpp_init_ioipsl : iloc = 19 17 ~~~~~~~~~~~~~~~ iabsf = 86 31 ihals = 1 1 ihale = 1 1 iom_nf90_open ~~~ open existing file: coordinates_xios.nc i n READ mode ---> coordinates_xios.nc OK read nav_lon_grid_T (rec: 1) in coordinates_xios.nc ok read nav_lat_grid_T (rec: 1) in coordinates_xios.nc ok read area_grid_T (rec: 1) in coordinates_xios.nc ok read nav_lon_grid_U (rec: 1) in coordinates_xios.nc ok read nav_lat_grid_U (rec: 1) in coordinates_xios.nc ok read area_grid_U (rec: 1) in coordinates_xios.nc ok read nav_lon_grid_V (rec: 1) in coordinates_xios.nc ok read nav_lat_grid_V (rec: 1) in coordinates_xios.nc ok read area_grid_V (rec: 1) in coordinates_xios.nc ok read nav_lon_grid_W (rec: 1) in coordinates_xios.nc ok read nav_lat_grid_W (rec: 1) in coordinates_xios.nc ok read area_grid_W (rec: 1) in coordinates_xios.nc ok read bounds_lon_grid_T (rec: 1) in coordinates_xios.nc ok read bounds_lat_grid_T (rec: 1) in coordinates_xios.nc ok read bounds_lon_grid_U (rec: 1) in coordinates_xios.nc ok read bounds_lat_grid_U (rec: 1) in coordinates_xios.nc ok read bounds_lon_grid_V (rec: 1) in coordinates_xios.nc ok read bounds_lat_grid_V (rec: 1) in coordinates_xios.nc ok read bounds_lon_grid_W (rec: 1) in coordinates_xios.nc ok read bounds_lat_grid_W (rec: 1) in coordinates_xios.nc ok iom_close ~~~ close file: coordinates_xios.nc ok phy_cst : initialization of ocean parameters and constants ~~~~~~~ Domain info dimension of model Local domain Global domain Data domain jpi : 19 jpiglo : 362 jpidta : 362 jpj : 17 jpjglo : 332 jpjdta : 332 jpk : 75 jpk : 75 jpkdta : 75 jpij : 323 mpp local domain info (mpp) jpni : 22 jpreci : 1 jpnj : 22 jprecj : 1 jpnij : 360 lateral domain boundary condition type : jperio = 6 Constants mathematical constant rpi = 3.14159265358979 day rday = 86400.0000000000 s sideral year rsiyea = 31558149.0101107 s sideral day rsiday = 86164.0996559118 s omega omega = 7.292115083046062E-005 s^-1 nb of months per year raamo = 12.0000000000000 months nb of hours per day rjjhh = 24.0000000000000 hours nb of minutes per hour rhhmm = 60.0000000000000 mn nb of seconds per minute rmmss = 60.0000000000000 s earth radius ra = 6371229.00000000 m gravity grav = 9.80665000000000 m/s^2 triple point of temperature rtt = 273.160000000000 K freezing point of water rt0 = 273.150000000000 K melting point of snow rt0_snow = 273.150000000000 K melting point of ice rt0_ice = 273.150000000000 K reference density and heat capacity now defined in eosbn2.f90 thermal conductivity of pure ice = 2.03439600000000 J/s/m/K fresh ice specific heat = 2067.00000000000 J/kg/K latent heat of fusion of fresh ice / snow = 334000.000000000 J/kg latent heat of subl. of fresh ice / snow = 2834000.00000000 J/kg volumetric latent heat fusion of snow = 110220000.000000 J/m^3 density of sea ice = 917.000000000000 kg/m^3 density of snow = 330.000000000000 kg/m^3 emissivity of snow or ice = 0.970000000000000 salinity of ice = 6.00000000000000 psu salinity of sea = 34.7000000000000 psu latent heat of evaporation (water) = 2500000.00000000 J/m^3 correction factor for solar radiation = 0.900000000000000 von Karman constant = 0.400000000000000 Stefan-Boltzmann constant = 5.670000000000000E-008 J/s/m^2/K^4 conversion: degre ==> radian rad = 1.745329251994330E-002 smallest real computer value rsmall = 1.110223024625157E-016 eos_init : equation of state ~~~~~~~~ Namelist nameos : set eos parameters flag for eq. of state and N^2 nn_eos = -1 model uses Conservative Temperature Important: model must be initialized with CT and SA fields use of TEOS-10 equation of state (cons. temp. and abs. salinity) volumic mass of reference rau0 = 1026.00000000000 kg/m^3 1. / rau0 r1_rau0 = 9.746588693957114E-004 m^3/kg ocean specific heat rcp = 3991.86795711963 J/Kelvin rau0 * rcp rau0_rcp = 4095656.52400474 1. / ( rau0 * rcp ) r1_rau0_rcp = 2.441610994816036E-007 dom_cfg : set the ocean configuration ~~~~~~~ ocean model configuration used : cp_cfg = orca jp_cfg = 1 global domain lateral boundaries jperio= 6, cyclic east-west and north fold with F-point pivot dom_glo : domain: data / local ~~~~~~~ data input domain : jpidta = 362 jpjdta = 332 jpkdta = 75 global or zoom domain: jpiglo = 362 jpjglo = 332 jpk = 75 local domain : jpi = 19 jpj = 17 jpk = 75 south-west indices jpizoom = 1 jpjzoom = 1 conversion local ==> data i-index domain 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 conversion data ==> local i-index domain starting index 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ending index 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 conversion local ==> data j-index domain 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 conversion data ==> local j-index domain starting index 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 ending index 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 17 zoom flags : lzoom = F (T = zoom, F = global ) lzoom_e = F (T = forced closed east boundary) lzoom_w = F (T = forced closed west boundary) lzoom_s = F (T = forced closed South boundary) lzoom_n = F (T = forced closed North boundary) dom_init : domain initialization ~~~~~~~~ dom_nam : domain initialization through namelist read ~~~~~~~ Namelist namrun job number nn_no = 0 experiment name for output cn_exp = cm6111-ndg-std file prefix restart input cn_ocerst_in= restartopa restart input directory cn_ocerst_indir= . file prefix restart output cn_ocerst_out= restart restart output directory cn_ocerst_outdir= . restart logical ln_rstart = T restart but zero for current for T ln_rstart_ts = F start with forward time step nn_euler = 1 control of time step nn_rstctl = 2 number of the first time step nn_it000 = 9855713 number of the last time step nn_itend = 9856672 initial calendar date aammjj nn_date0 = 20000401 leap year calendar (0/1) nn_leapy = 1 initial state output nn_istate = 0 frequency of restart file nn_stock = 9856672 frequency of output file nn_write = 5475 multi file dimgout ln_dimgnnn = F mask land points ln_mskland = T additional CF standard metadata ln_cfmeta = T overwrite an existing file ln_clobber = F NetCDF chunksize (bytes) nn_chunksz = 0 The IOIPSL calendar is "gregorian", i.e. leap year Namelist namdom : space & time domain flag read/compute bathymetry nn_bathy = 1 Depth (if =0 bathy=jpkm1) rn_bathy = 0.000000000000000E+000 min depth of the ocean (>0) or rn_hmin = 20.0000000000000 min number of ocean level (<0) minimum thickness of partial rn_e3zps_min = 20.0000000000000 (m) step level rn_e3zps_rat = 0.100000000000000 create mesh/mask file(s) nn_msh = 0 = 0 no file created = 1 mesh_mask = 2 mesh and mask = 3 mesh_hgr, msh_zgr and mask ocean time step rn_rdt = 2700.00000000000 asselin time filter parameter rn_atfp = 0.100000000000000 acceleration of converge nn_acc = 0 nn_acc=1: surface tracer rdt rn_rdtmin = 28800.0000000000 bottom tracer rdt rdtmax = 28800.0000000000 depth of transition rn_rdth = 800.000000000000 suppression of closed seas (=0) nn_closea = 1 online coarsening of dynamical fields ln_crs = F type of horizontal mesh jphgr_msh = 0 longitude of first raw and column T-point ppglam0 = 999999.000000000 latitude of first raw and column T-point ppgphi0 = 999999.000000000 zonal grid-spacing (degrees) ppe1_deg = 999999.000000000 meridional grid-spacing (degrees) ppe2_deg = 999999.000000000 zonal grid-spacing (degrees) ppe1_m = 999999.000000000 meridional grid-spacing (degrees) ppe2_m = 999999.000000000 ORCA r4, r2 and r05 coefficients ppsur = -3958.95137127683 ppa0 = 103.953009600000 ppa1 = 2.41595126900000 ppkth = 15.3510137000000 ppacr = 7.00000000000000 Minimum vertical spacing ppdzmin = 999999.000000000 Maximum depth pphmax = 999999.000000000 Use double tanf function for vertical coordinates ldbletanh = T Double tanh function parameters ppa2 = 100.760928500000 ppkth2 = 48.0298937200000 ppacr2 = 13.0000000000000 Namelist namcla cross land advection nn_cla = 0 dom_clo : closed seas ~~~~~~~ dom_hgr : define the horizontal mesh from ithe following par_oce parameters ~~~~~~~ type of horizontal mesh jphgr_msh = 0 position of the first row and ppglam0 = 999999.000000000 column grid-point (degrees) ppgphi0 = 999999.000000000 zonal grid-spacing (degrees) ppe1_deg = 999999.000000000 meridional grid-spacing (degrees) ppe2_deg = 999999.000000000 zonal grid-spacing (meters) ppe1_m = 999999.000000000 meridional grid-spacing (meters) ppe2_m = 999999.000000000 curvilinear coordinate on the sphere read in "coordinate" file hgr_read : read the horizontal coordinates ~~~~~~~~ jpiglo = 362 jpjglo = 332 jpk = 75 iom_nf90_open ~~~ open existing file: coordinates.nc in REA D mode ---> coordinates.nc OK read glamt (rec: 1) in coordinates.nc ok read glamu (rec: 1) in coordinates.nc ok read glamv (rec: 1) in coordinates.nc ok read glamf (rec: 1) in coordinates.nc ok read gphit (rec: 1) in coordinates.nc ok read gphiu (rec: 1) in coordinates.nc ok read gphiv (rec: 1) in coordinates.nc ok read gphif (rec: 1) in coordinates.nc ok read e1t (rec: 1) in coordinates.nc ok read e1u (rec: 1) in coordinates.nc ok read e1v (rec: 1) in coordinates.nc ok read e1f (rec: 1) in coordinates.nc ok read e2t (rec: 1) in coordinates.nc ok read e2u (rec: 1) in coordinates.nc ok read e2v (rec: 1) in coordinates.nc ok read e2f (rec: 1) in coordinates.nc ok iom_close ~~~ close file: coordinates.nc ok orca_r1: Gibraltar : e2u reduced to 20 km orca_r1: Bhosporus : e2u reduced to 10 km orca_r1: Lombok : e1v reduced to 10 km orca_r1: Sumba : e1v reduced to 8 km orca_r1: Ombai : e1v reduced to 13 km orca_r1: Timor Passage : e1v reduced to 20 km orca_r1: W Halmahera : e1v reduced to 30 km orca_r1: E Halmahera : e1v reduced to 50 km dom_zgr : vertical coordinate ~~~~~~~ Namelist namzgr : set vertical coordinate z-coordinate - full steps ln_zco = F z-coordinate - partial steps ln_zps = T s- or hybrid z-s-coordinate ln_sco = F ice shelf cavities ln_isfcav = F zgr_z : Reference vertical z-coordinates ~~~~~~~ Value of coefficients for vertical mesh: zsur = -3958.95137127683 za0 = 103.953009600000 za1 = 2.41595126900000 zkth = 15.3510137000000 zacr = 7.00000000000000 (Double tanh za2 = 100.760928500000 parameters) zkth2= 48.0298937200000 zacr2= 13.0000000000000 Reference z-coordinate depth and scale factors: level gdept_1d gdepw_1d e3t_1d e3w_1d 1 0.51 0.00 1.02 1.00 2 1.56 1.02 1.08 1.05 3 2.67 2.10 1.15 1.11 4 3.86 3.25 1.23 1.19 5 5.14 4.49 1.34 1.28 6 6.54 5.83 1.47 1.40 7 8.09 7.30 1.63 1.55 8 9.82 8.93 1.83 1.73 9 11.77 10.77 2.08 1.95 10 13.99 12.85 2.37 2.22 11 16.53 15.22 2.71 2.53 12 19.43 17.93 3.11 2.90 13 22.76 21.04 3.56 3.33 14 26.56 24.60 4.05 3.80 15 30.87 28.65 4.59 4.31 16 35.74 33.24 5.15 4.86 17 41.18 38.39 5.73 5.44 18 47.21 44.12 6.33 6.03 19 53.85 50.45 6.95 6.64 20 61.11 57.40 7.58 7.26 21 69.02 64.98 8.24 7.91 22 77.61 73.23 8.94 8.59 23 86.93 82.18 9.70 9.32 24 97.04 91.88 10.53 10.11 25 108.03 102.42 11.46 10.98 26 120.00 113.89 12.50 11.96 27 133.08 126.39 13.68 13.07 28 147.41 140.07 15.01 14.32 29 163.16 155.10 16.54 15.75 30 180.55 171.64 18.27 17.38 31 199.79 189.92 20.25 19.23 32 221.14 210.18 22.50 21.34 33 244.89 232.70 25.05 23.74 34 271.36 257.76 27.94 26.45 35 300.89 285.72 31.19 29.52 36 333.86 316.92 34.83 32.96 37 370.69 351.77 38.89 36.81 38 411.79 390.68 43.39 41.09 39 457.63 434.09 48.35 45.81 40 508.64 482.46 53.76 51.00 41 565.29 536.23 59.62 56.63 42 628.03 595.87 65.92 62.72 43 697.26 661.81 72.61 69.22 44 773.37 734.43 79.66 76.10 45 856.68 814.11 87.00 83.30 46 947.45 901.12 94.56 90.76 47 1045.85 995.69 102.26 98.40 48 1151.99 1097.95 110.01 106.14 49 1265.86 1207.96 117.71 113.87 50 1387.38 1325.67 125.29 121.52 51 1516.36 1450.95 132.64 129.00 52 1652.57 1583.58 139.71 136.22 53 1795.67 1723.28 146.43 143.12 54 1945.30 1869.69 152.75 149.64 55 2101.03 2022.42 158.64 155.75 56 2262.42 2181.04 164.08 161.41 57 2429.03 2345.10 169.06 166.62 58 2600.38 2514.14 173.58 171.37 59 2776.04 2687.70 177.67 175.68 60 2955.57 2865.35 181.33 179.55 61 3138.56 3046.66 184.60 183.01 62 3324.64 3231.24 187.50 186.09 63 3513.45 3418.72 190.06 188.82 64 3704.66 3608.77 192.31 191.22 65 3897.98 3801.07 194.29 193.34 66 4093.16 3995.35 196.02 195.19 67 4289.95 4191.37 197.53 196.80 68 4488.15 4388.89 198.84 198.21 69 4687.58 4587.73 199.98 199.43 70 4888.07 4787.70 200.97 200.49 71 5089.48 4988.67 201.83 201.41 72 5291.68 5190.49 202.57 202.21 73 5494.58 5393.05 203.20 202.90 74 5698.06 5596.25 203.75 203.49 75 5902.06 5800.00 204.23 204.00 zgr_bat : defines level and meter bathymetry ~~~~~~~ iom_nf90_open ~~~ open existing file: bathy_meter.nc in REA D mode ---> bathy_meter.nc OK read Bathymetry (rec: 1) in bathy_meter.nc ok iom_close ~~~ close file: bathy_meter.nc ok Minimum ocean depth: 24.5959881615022 minimum number of ocean levels : 13 zgr_zps : z-coordinate with partial steps ~~~~~~~ mbathy is recomputed : bathy_level file is NOT used zgr_bat_ctl : check the bathymetry ~~~~~~~~~~~ suppress isolated ocean grid points ----------------------------------- the number of ocean level at grid-point (i,j) = 8 16 is changed from 40 to 39 the number of ocean level at grid-point (i,j) = 11 16 is changed from 44 to 43 1104 ocean grid points suppressed mbathy set to 0 along east and west boundary: nperio = 0 maximum number of ocean level = 45 < jpk-1 you can decrease jpk to 46 zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ~~~~~~~~~~~~~ zgr_top_level : ocean top k-index of T-, U-, V- and W-levels ~~~~~~~~~~~~~ dommsk : ocean mask ~~~~~~ Namelist namlbc lateral momentum boundary cond. rn_shlat = 0.000000000000000E+000 consistency with analytical form ln_vorlat = F ocean lateral free-slip orca_r1: increase friction near the following straits : Gibraltar Bhosporus Makassar (Top) Lombok Ombai Timor Passage West Halmahera East Halmahera dom_vvl_init : Variable volume activated ~~~~~~~~~~~~ dom_vvl_ctl : choice/control of the variable vertical coordinate ~~~~~~~~~~~ Namelist nam_vvl : chose a vertical coordinate zstar ln_vvl_zstar = T ztilde ln_vvl_ztilde = F layer ln_vvl_layer = F ztilde as zstar ln_vvl_ztilde_as_zstar = F ztilde near the equator ln_vvl_zstar_at_eqtor = F Namelist nam_vvl : thickness diffusion coefficient rn_ahe3 = 0.000000000000000E+000 Namelist nam_vvl : maximum e3t deformation fractional change rn_zdef_max = 0.900000000000000 Namelist nam_vvl : z-tilde to zstar restoration timescale (days) rn_rst_e3t = 30.0000000000000 Namelist nam_vvl : z-tilde cutoff frequency of low-pass filter (days ) rn_lf_cutoff = 5.00000000000000 Namelist nam_vvl : debug prints ln_vvl_dbg = T zstar vertical coordinate is used rst_read : read oce NetCDF restart file ~~~~~~~~ iom_nf90_open ~~~ open existing file: ./restartopa_0000.nc in READ mode ---> ./restartopa_0000.nc OK read sshn (rec: 1) in ./restartopa_0000.nc ok read fse3t_b (rec: 1) in ./restartopa_0000.nc ok read fse3t_n (rec: 1) in ./restartopa_0000.nc ok dom_vvl_rst : fse3t_b and fse3t_n found in restart files dom_stp : time stepping setting ~~~~~~~ synchronous time stepping dynamics and tracer time step = 0.750000000000000 hours istate_ini : Initialization of the dynamics and tracers ~~~~~~~~~~ dta_tsd_init : Temperature & Salinity data ~~~~~~~~~~~~ Namelist namtsd Initialisation of ocean T & S with T &S input data ln_tsd_init = T damping of ocean T & S toward T &S input data ln_tsd_tradmp = F ===>>> : W A R N I N G =============== dta_tsd_init: ocean restart and T & S data intialisation, we keep the restart T & S values and set ln_tsd_init to FALSE read ub (rec: 1) in ./restartopa_0000.nc ok read vb (rec: 1) in ./restartopa_0000.nc ok read tb (rec: 1) in ./restartopa_0000.nc ok read sb (rec: 1) in ./restartopa_0000.nc ok read rotb (rec: 1) in ./restartopa_0000.nc ok read hdivb (rec: 1) in ./restartopa_0000.nc ok read sshb (rec: 1) in ./restartopa_0000.nc ok read un (rec: 1) in ./restartopa_0000.nc ok read vn (rec: 1) in ./restartopa_0000.nc ok read tn (rec: 1) in ./restartopa_0000.nc ok read sn (rec: 1) in ./restartopa_0000.nc ok read sshn (rec: 1) in ./restartopa_0000.nc ok read rotn (rec: 1) in ./restartopa_0000.nc ok read hdivn (rec: 1) in ./restartopa_0000.nc ok read rhop (rec: 1) in ./restartopa_0000.nc ok *** Info read in restart : previous time-step : 9855712 *** restart option nrstdt = 2 : calendar parameters read in restart *** Info used values : date ndastp : 20000331 number of elapsed days since the begining of run : 307991.000000000 =======>> 1/2 time step before the start of the run DATE Y/M/D = 2000/ 3/31 nsec_day: 85050 nsec_week: 430650 ======>> time-step = 9855713 New day, DATE Y/M/D = 2000/04/01 nday_year = 092 nsec_year = 7863750 nsec_month = 1350 nsec_day = 1350 nsec_week = 433350 sbc_init : surface boundary condition setting ~~~~~~~~ Namelist namsbc (partly overwritten with CPP key setting) frequency update of sbc (and ice) nn_fsbc = 2 Type of sbc : analytical formulation ln_ana = F flux formulation ln_flx = F CLIO bulk formulation ln_blk_clio = F CORE bulk formulation ln_blk_core = F MFS bulk formulation ln_blk_mfs = F ocean-atmosphere coupled formulation ln_cpl = T forced-coupled mixed formulation ln_mixcpl = F OASIS coupling (with atm or sas) lk_oasis = T components of your executable nn_components = 0 Multicategory heat flux formulation (LIM3) nn_limflx = 2 Misc. options of sbc : Patm gradient added in ocean & ice Eqs. ln_apr_dyn = F ice management in the sbc (=0/1/2/3) nn_ice = 3 ice-ocean embedded/levitating (=0/1/2) nn_ice_embd = 1 daily mean to diurnal cycle qsr ln_dm2dc = F runoff / runoff mouths ln_rnf = F iceshelf formulation nn_isf = 3 Sea Surface Restoring on SST and/or SSS ln_ssr = F FreshWater Budget control (=0/1/2) nn_fwb = 0 closed sea (=0/1) (set in namdom) nn_closea = 1 n. of iterations if land-sea-mask applied nn_lsm = 0 Redistribute a single flux over categories (nn_limflx = 2) pure coupled formulation sbc_cpl_init : namsbc_cpl namelist ~~~~~~~~~~~~ received fields (mutiple ice categogies) 10m wind module = coupled (no) stress module = none (no) surface stress = mixed oce-ice (no) - referential = cartesian - orientation = eastward-northward - mesh = U,V non-solar heat flux sensitivity = coupled (no) solar heat flux = conservative (no) non-solar heat flux = conservative (no) freshwater budget = conservative (no) runoffs = coupled (no) calving = coupled (no) iceberg = coupled (no) ice shelf = coupled (no) sea ice heat fluxes = ( ) atm co2 = none (no) sent fields (multiple ice categories) surface temperature = weighted oce and ice (no) albedo = weighted ice (no) ice/snow thickness = none (no) surface current = mixed oce-ice (no) - referential = cartesian - orientation = eastward-northward - mesh = T oce co2 flux = none (no) nn_cplmodel = 1 ln_usecplmask = F runoffs received from oasis -> force ln_rnf = T iceshelf received from oasis cpl_define : initialization in coupled ocean/atmosphere case ~~~~~~~~~~~~~~~~~ sbc_ssm : sea surface mean fields ~~~~~~~ read ssu_m (rec: 1) in ./restartopa_0000.nc ok read ssv_m (rec: 1) in ./restartopa_0000.nc ok read sst_m (rec: 1) in ./restartopa_0000.nc ok read sss_m (rec: 1) in ./restartopa_0000.nc ok read ssh_m (rec: 1) in ./restartopa_0000.nc ok read e3t_m (rec: 1) in ./restartopa_0000.nc ok read frq_m (rec: 1) in ./restartopa_0000.nc ok ~~~~~~~ mean fields read in the ocean restart file sbc_isf: heat flux of the ice shelf ~~~~~~~~~ sbcisf : ~~~~~~~~ nn_isf = 3 nn_isfblk = 1 rn_hisf_tbl = 30.0000000000000 ln_divisf = T nn_gammablk = 1 rn_tfri2 = 0.000000000000000E+000 iom_nf90_open ~~~ open existing file: runoff-icb_DaiTrenber th_Depoorter_eORCA1_JD.nc in READ mode ---> runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc OK read sodepmax_isf (rec: 1) in runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc ok iom_close ~~~ close file: runoff-icb_DaiTrenberth_Depoorter _eORCA1_JD.nc ok iom_nf90_open ~~~ open existing file: runoff-icb_DaiTrenber th_Depoorter_eORCA1_JD.nc in READ mode ---> runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc OK read sodepmin_isf (rec: 1) in runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc ok iom_close ~~~ close file: runoff-icb_DaiTrenberth_Depoorter _eORCA1_JD.nc ok sbc_rnf : runoff ~~~~~~~ Namelist namsbc_rnf specific river mouths treatment ln_rnf_mouth = F river mouth additional Kz rn_avt_rnf = 1.000000000000000E-003 depth of river mouth additional mixing rn_hrnf = 15.0000000000000 multiplicative factor for runoff rn_rfact = 1.00000000000000 runoffs depth read in a file iom_nf90_open ~~~ open existing file: ./runoffs_eORCA1.0_de pths.nc in READ mode ---> ./runoffs_eORCA1.0_depths.nc OK read rodepth (rec: 1) in ./runoffs_eORCA1.0_depths.nc ok iom_close ~~~ close file: ./runoffs_eORCA1.0_depths.nc ok No specific treatment at river mouths sbc_ice_lim : update ocean surface boudary condition ~~~~~~~~~~~ via Louvain la Neuve Ice Model (LIM-3) time stepping file : namelist_ice_ref open ok unit = 20 status = OLD form = FORMATTED access = SEQUENTIAL file : namelist_ice_cfg open ok unit = 21 status = OLD form = FORMATTED access = SEQUENTIAL file : output.namelist.ice open ok unit = 22 status = UNKNOWN form = FORMATTED access = SEQUENTIAL ice_run : ice share parameters for dynamics/advection/thermo of sea-ice ~~~~~~ number of ice categories = 5 number of ice layers = 2 number of snow layers = 1 switch for ice dynamics (1) or not (0) ln_limdyn = T maximum ice concentration for NH = 0.997000000000000 maximum ice concentration for SH = 0.950000000000000 Diagnose heat/salt budget or not ln_limdiahsb = F Output heat/salt budget or not ln_limdiaout = T control prints in ocean.out for (i,j)=(iiceprt,jiceprt) = F i-index for control prints (ln_icectl=true) = 10 j-index for control prints (ln_icectl=true) = 10 activate sea ice nudging (T) or not (F) ln_limndg = T ice_itd : ice cat distribution ~~~~~~ shape of ice categories distribution nn_catbnd = 2 mean ice thickness in the domain (only active if nn_catbnd=2) rn_himean = 2.00000000000000 max ice thickness in the last category jpl rn_himax_bot = 99.0000000000000 lim_itd_init : Initialization of ice cat distribution ~~~~~~~~~~~~ Thickness category boundaries hi_max 0.000000000000000E+000 0.454016742415552 1.12931382308363 2.14145898894797 3.67059640991554 99.0000000000000 lim_hdf_init : Ice horizontal diffusion ~~~~~~~~~~~ horizontal diffusivity calculation nn_ahi0 = -1 horizontal diffusivity coeff. (orca2 grid) rn_ahi0_ref = 350.000000000000 convergence check frequency of the Crant-Nicholson scheme nn_convfrq = 5 No sea-ice diffusion applied lim_thd : Ice Thermodynamics ~~~~~~~ Namelist of ice parameters for ice thermodynamic computation ice thick. for lateral accretion rn_hnewice = 0.100000000000000 Frazil ice thickness as a function of wind or not ln_frazil = F Maximum proportion of frazil ice collecting at bottom rn_maxfrazb = 1.00000000000000 Thresold relative drift speed for collection of frazil rn_vfrazb = 0.417000000000000 Squeezing coefficient for collection of frazil rn_Cfrazb = 5.00000000000000 minimum ice thickness rn_himin = 0.100000000000000 numerical carac. of the scheme for diffusion in ice coefficient for ice-lead partition of snowfall rn_betas = 0.660000000000000 extinction radiation parameter in sea ice rn_kappa_i = 1.00000000000000 maximal n. of iter. for heat diffusion computation nn_conv_dif = 50.0000000000000 maximal err. on T for heat diffusion computation rn_terr_dif = 1.000000000000000E-004 switch for comp. of thermal conductivity in the ice nn_ice_thcon = 1 thermal conductivity of the snow rn_cdsn = 0.500000000000000 check heat conservation in the ice/snow con_i = F virtual ITD mono-category parameterizations (1) or not nn_monocat = 0 iterate the surface non-solar flux (T) or not (F) ln_it_qnsice = T lim_thd_sal_init : Ice parameters for salinity ~~~~~~~~~~~~~~~~ switch for salinity nn_icesal = 2 bulk salinity value if nn_icesal = 1 = 4.00000000000000 restoring salinity for GD = 5.00000000000000 restoring time for GD = 1730000.00000000 restoring salinity for flushing = 2.00000000000000 restoring time for flushing = 864000.000000000 Maximum tolerated ice salinity = 20.0000000000000 Minimum tolerated ice salinity = 0.100000000000000 lim_msh : LIM-3 sea-ice model, mesh initialization ~~~~~~~ the model domain is entirely in the southern hemisphere: njeq = 17 lim_itd_me_init : ice parameters for mechanical ice redistribution ~~~~~~~~~~~~~~~ Fraction of shear energy contributing to ridging rn_cs = 0.500000000000000 Fraction of snow volume conserved during ridging rn_fsnowrdg = 0.500000000000000 Fraction of snow volume conserved during ridging rn_fsnowrft = 0.500000000000000 Fraction of total ice coverage contributing to ridging rn_gstar = 0.150000000000000 Equivalent to G* for an exponential part function rn_astar = 3.000000000000000E-002 Quantity playing a role in max ridged ice thickness rn_hstar = 25.0000000000000 Rafting of ice sheets or not ln_rafting = T Parmeter thickness (threshold between ridge-raft) rn_hraft = 0.750000000000000 Rafting hyperbolic tangent coefficient rn_craft = 5.00000000000000 Initial porosity of ridges rn_por_rdg = 0.300000000000000 Switch for part. function (0) linear (1) exponential nn_partfun = 1 lim_rst_read : read ice NetCDF restart file ~~~~~~~~~~~~~ iom_nf90_open ~~~ open existing file: ./restart_ice_in_0000 .nc in READ mode ---> ./restart_ice_in_0000.nc OK read ice restart file at time step : 9855712.00000000 in any case we force it to nit000 - 1 : 9855712 read v_i_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read v_s_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read smv_i_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read oa_i_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read a_i_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read t_su_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_sl1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il2_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read v_i_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read v_s_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read smv_i_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read oa_i_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read a_i_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read t_su_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_sl1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il2_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read v_i_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read v_s_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read smv_i_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read oa_i_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read a_i_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read t_su_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_sl1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il2_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read v_i_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read v_s_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read smv_i_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read oa_i_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read a_i_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read t_su_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_sl1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il2_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read v_i_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read v_s_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read smv_i_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read oa_i_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read a_i_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read t_su_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_sl1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read tempt_il2_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read u_ice (rec: 1) in ./restart_ice_in_0000.nc ok read v_ice (rec: 1) in ./restart_ice_in_0000.nc ok read stress1_i (rec: 1) in ./restart_ice_in_0000.nc ok read stress2_i (rec: 1) in ./restart_ice_in_0000.nc ok read stress12_i (rec: 1) in ./restart_ice_in_0000.nc ok read snwice_mass (rec: 1) in ./restart_ice_in_0000.nc ok read snwice_mass_b (rec: 1) in ./restart_ice_in_0000.nc ok read sxice_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syice_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxice_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syyice_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyice_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsn_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sysn_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsn_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syysn_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysn_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxa_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sya_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxa_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syya_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxya_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxc0_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syc0_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxc0_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syyc0_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyc0_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsal_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sysal_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsal_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syysal_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysal_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxage_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syage_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxage_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syyage_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyage_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxice_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syice_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxice_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syyice_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyice_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsn_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sysn_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsn_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syysn_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysn_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxa_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sya_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxa_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syya_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxya_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxc0_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syc0_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxc0_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syyc0_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyc0_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsal_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sysal_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsal_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syysal_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysal_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxage_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syage_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxage_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syyage_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyage_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxice_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syice_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxice_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syyice_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyice_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsn_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sysn_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsn_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syysn_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysn_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxa_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sya_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxa_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syya_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxya_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxc0_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syc0_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxc0_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syyc0_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyc0_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsal_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sysal_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsal_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syysal_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysal_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxage_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syage_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxage_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syyage_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyage_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxice_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syice_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxice_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syyice_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyice_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsn_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sysn_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsn_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syysn_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysn_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxa_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sya_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxa_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syya_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxya_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxc0_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syc0_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxc0_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syyc0_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyc0_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsal_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sysal_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsal_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syysal_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysal_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxage_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syage_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxage_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syyage_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyage_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxice_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syice_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxice_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syyice_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyice_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsn_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sysn_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsn_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syysn_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysn_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxa_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sya_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxa_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syya_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxya_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxc0_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syc0_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxc0_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syyc0_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyc0_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxsal_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sysal_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxsal_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syysal_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxysal_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxage_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syage_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxage_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syyage_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxyage_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxopw (rec: 1) in ./restart_ice_in_0000.nc ok read syopw (rec: 1) in ./restart_ice_in_0000.nc ok read sxxopw (rec: 1) in ./restart_ice_in_0000.nc ok read syyopw (rec: 1) in ./restart_ice_in_0000.nc ok read sxyopw (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il1_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il2_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il2_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il2_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il2_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il2_htc1 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il1_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il2_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il2_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il2_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il2_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il2_htc2 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il1_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il2_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il2_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il2_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il2_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il2_htc3 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il1_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il2_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il2_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il2_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il2_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il2_htc4 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il1_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxe_il2_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sye_il2_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxxe_il2_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read syye_il2_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok read sxye_il2_htc5 (rec: 1) in ./restart_ice_in_0000.nc ok lim_sbc_init : LIM-3 sea-ice - surface boundary condition ~~~~~~~~~~~~~ lim_diahsb_init : check the heat and salt budgets ~~~~~~~~~~~~ ~~~~~~~ lim_diahsb_rst read at it= 0 date= 20000401 ~~~~~~~ read vol_loc_ini (rec: 1) in ./restart_ice_in_0000.nc ok read tem_loc_ini (rec: 1) in ./restart_ice_in_0000.nc ok read sal_loc_ini (rec: 1) in ./restart_ice_in_0000.nc ok ===>>> : W A R N I N G =============== W A R N I N G: end of record or file while reading namelist namicendg in refer ence namelist iostat = -1 lim_ndg_init : initialization of sea ice nudging ~~~~~~~ Namelist namicendg : SI relaxation term (Yes=1, based on SIC) nn_icendg = 1 (Yes=2, based on SIC and SIT) surface temperature restoration under ice = 1 dQ/dSIT (restoring coeff for SI grwth) rn_dqdt_mlt = 3500.00000000000 W/m2 dQ/dSIT (restoring coeff for SI melt) rn_dqdt_frz = 3500.00000000000 W/m2 dQ/dSST (restoring coeff for SST freezing) rn_dqdt_sst = -100.000000000000 W/m2 SI nudging activation threshold rn_ndg_thld = 1.000000000000000E-002 mask of the sea ice nudging (Yes=1) nn_msk = 0 lim_ndg_init : SI restoring target SIC data ~~~~~~~~~~~~ namicendg Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./sic_restore_data variable name: siconc frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: dyn_nept_init : Simplified Neptune module ~~~~~~~~~~~~~ --> Reading namelist namdyn_nept parameters: ln_neptsimp = F zdf_init: vertical physics ~~~~~~~~ Namelist namzdf : set vertical mixing mixing parameters vertical eddy viscosity rn_avm0 = 1.200000000000000E-004 vertical eddy diffusivity rn_avt0 = 1.200000000000000E-005 constant background or profile nn_avb = 0 horizontal variation for avtb nn_havtb = 0 time splitting / backward scheme ln_zdfexp = F number of time step nn_zdfexp = 3 enhanced vertical diffusion ln_zdfevd = T applied on momentum (=1/0) nn_evdm = 0 vertical coefficient for evd rn_avevd = 100.000000000000 non-penetrative convection (npc) ln_zdfnpc = F npc call frequency nn_npc = 1 npc print frequency nn_npcp = 365 vertical mixing option : TKE dependent eddy coefficients convection : use enhanced vertical dif. scheme use the 1.5 turbulent closure zdf_bfr_init : momentum bottom friction ~~~~~~~~~~~~~ Namelist nam_bfr : set bottom friction parameters quadratic bottom friction friction coef. rn_bfri2 = 1.000000000000000E-003 Max. coef. (log case) rn_bfri2_max = 0.100000000000000 background tke rn_bfeb2 = 2.500000000000000E-003 log formulation ln_bfr2d = F bottom roughness rn_bfrz0 [m] = 3.000000000000000E-003 implicit bottom friction switch ln_bfrimp = T zdf_tke_init : New tidal mixing case: force rn_emin = 1.e-10 and rmxl_min = 1. e-3 ~~~~~~~~~~~~ zdf_tke_init : tke turbulent closure scheme - initialisation ~~~~~~~~~~~~ Namelist namzdf_tke : set tke mixing parameters coef. to compute avt rn_ediff = 0.100000000000000 Kolmogoroff dissipation coef. rn_ediss = 0.700000000000000 tke surface input coef. rn_ebb = 67.8300000000000 minimum value of tke rn_emin = 1.000000000000000E-010 surface minimum value of tke rn_emin0 = 1.000000000000000E-004 background shear (>0) rn_bshear = 9.999999999999999E-021 mixing length type nn_mxl = 2 prandl number flag nn_pdl = 1 surface mixing length = F(stress) or not ln_mxl0 = T type of scaling under sea-ice nn_mxl0 = 2 surface mixing length minimum value rn_mxl0 = 4.000000000000000E-002 flag to take into acc. Langmuir circ. ln_lc = T coef to compute verticla velocity of LC rn_lc = 0.200000000000000 test param. to add tke induced by wind nn_etau = 0 flag for computation of exp. tke profile nn_htau = 1 fraction of en which pene. the thermocline rn_efr = 5.000000000000000E-002 critical Richardson nb with your parameters ri_cri = 0.222222222222222 use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min read en (rec: 1) in ./restartopa_0000.nc ok read avt (rec: 1) in ./restartopa_0000.nc ok read avm (rec: 1) in ./restartopa_0000.nc ok read avmu (rec: 1) in ./restartopa_0000.nc ok read avmv (rec: 1) in ./restartopa_0000.nc ok read dissl (rec: 1) in ./restartopa_0000.nc ok zdf_tmx_init : internal wave-driven mixing ~~~~~~~~~~~~ Namelist namzdf_tmx_new : set wave-driven mixing parameters Pycnocline-intensified diss. scales as N (=1) or N^2 (=2) = 2 Variable (T) or constant (F) mixing efficiency = T Differential internal wave-driven mixing (T) or not (F) = T Force the background value applied to avm & avt in TKE to be everywhere the viscous molecular value & a very small diffusive value, resp. iom_nf90_open ~~~ open existing file: mixing_power_bot.nc i n READ mode ---> mixing_power_bot.nc OK read field (rec: 1) in mixing_power_bot.nc ok iom_close ~~~ close file: mixing_power_bot.nc ok iom_nf90_open ~~~ open existing file: mixing_power_pyc.nc i n READ mode ---> mixing_power_pyc.nc OK read field (rec: 1) in mixing_power_pyc.nc ok iom_close ~~~ close file: mixing_power_pyc.nc ok iom_nf90_open ~~~ open existing file: mixing_power_cri.nc i n READ mode ---> mixing_power_cri.nc OK read field (rec: 1) in mixing_power_cri.nc ok iom_close ~~~ close file: mixing_power_cri.nc ok iom_nf90_open ~~~ open existing file: decay_scale_bot.nc in READ mode ---> decay_scale_bot.nc OK read field (rec: 1) in decay_scale_bot.nc ok iom_close ~~~ close file: decay_scale_bot.nc ok iom_nf90_open ~~~ open existing file: decay_scale_cri.nc in READ mode ---> decay_scale_cri.nc OK read field (rec: 1) in decay_scale_cri.nc ok iom_close ~~~ close file: decay_scale_cri.nc ok High-mode wave-breaking energy: 0.523104138124428 TW Pycnocline-intensifed wave-breaking energy: 0.560652161246472 TW Critical slope wave-breaking energy: 0.118846065884477 TW zdf_ddm : double diffusive mixing ~~~~~~~ Namelist namzdf_ddm : set dd mixing parameter maximum avs for dd mixing rn_avts = 1.000000000000000E-004 heat/salt buoyancy flux ratio rn_hsbfr = 1.60000000000000 ldf_tra_init : lateral tracer physics ~~~~~~~~~~~~ Namelist namtra_ldf : lateral mixing parameters (type, direction, coefficien ts) laplacian operator ln_traldf_lap = T bilaplacian operator ln_traldf_bilap = F iso-level ln_traldf_level = F horizontal (geopotential) ln_traldf_hor = F iso-neutral ln_traldf_iso = T iso-neutral (Griffies) ln_traldf_grif = F Griffies strmfn diagnostics ln_traldf_gdia = F lateral eddy diffusivity rn_aht_0 = 1000.00000000000 background hor. diffusivity rn_ahtb_0 = 0.000000000000000E+000 eddy induced velocity coef. rn_aeiv_0 = 1000.00000000000 maximum isoppycnal slope rn_slpmax = 1.000000000000000E-002 pure lateral mixing in ML ln_triad_iso = F lateral mixing on bottom ln_botmix_grif = F tracer mixing coef. = F( latitude, longitude) harmonic tracer diffusion (default) ldf_tra_c2d : 2D eddy diffusivity and eddy ~~~~~~~~~~~ -- induced velocity coefficients ldf_dyn : lateral momentum physics ~~~~~~~ Namelist namdyn_ldf : set lateral mixing parameters laplacian operator ln_dynldf_lap = T bilaplacian operator ln_dynldf_bilap = F iso-level ln_dynldf_level = F horizontal (geopotential) ln_dynldf_hor = T iso-neutral ln_dynldf_iso = F horizontal laplacian eddy viscosity rn_ahm_0_lap = 20000.0000000000 background viscosity rn_ahmb_0 = 0.000000000000000E+000 horizontal bilaplacian eddy viscosity rn_ahm_0_blp = 0.000000000000000E+000 upper limit for laplacian eddy visc rn_ahm_m_lap = 40000.0000000000 upper limit for bilap eddy viscosity rn_ahm_m_blp = -1000000000000.00 momentum mixing coef. = F( latitude, longitude, depth) harmonic momentum diff. (default) ldf_dyn_c3d : 3D lateral eddy viscosity coefficient ~~~~~~~~~~~ laplacian operator: ahm proportional to e1 maximum grid-spacing = 111198.923448546 maximum value for ahm = 20000.0000000000 ahm profile : jk ahm depth t-level 1 1.0000 0.5000 2 0.9999 1.5496 3 0.9999 2.6608 4 0.9998 3.8486 5 0.9997 5.1317 6 0.9997 6.5332 7 0.9996 8.0813 8 0.9995 9.8099 9 0.9994 11.7589 10 0.9993 13.9741 11 0.9991 16.5062 12 0.9989 19.4084 13 0.9987 22.7342 14 0.9985 26.5331 15 0.9983 30.8479 16 0.9980 35.7125 17 0.9976 41.1516 18 0.9972 47.1829 19 0.9968 53.8210 20 0.9963 61.0822 21 0.9958 68.9897 22 0.9952 77.5772 23 0.9945 86.8927 24 0.9937 97.0012 25 0.9928 107.9858 26 0.9918 119.9503 27 0.9907 133.0200 28 0.9894 147.3432 29 0.9878 163.0931 30 0.9860 180.4691 31 0.9839 199.6984 32 0.9813 221.0377 33 0.9783 244.7739 34 0.9745 271.2251 35 0.9699 300.7404 36 0.9641 333.6988 37 0.9568 370.5065 38 0.9475 411.5933 39 0.9355 457.4060 40 0.9200 508.4012 41 0.8997 565.0351 42 0.8732 627.7513 43 0.8386 696.9682 44 0.7942 773.0642 45 0.7387 856.3641 46 0.6721 947.1256 47 0.5972 1045.5281 48 0.5195 1151.6649 49 0.4462 1265.5388 50 0.3839 1387.0616 51 0.3360 1516.0589 52 0.3024 1652.2771 53 0.2804 1795.3952 54 0.2670 1945.0372 55 0.2592 2100.7868 56 0.2549 2262.2009 57 0.2525 2428.8236 58 0.2513 2600.1976 59 0.2506 2775.8744 60 0.2503 2955.4225 61 0.2501 3138.4329 62 0.2501 3324.5236 63 0.2500 3513.3418 64 0.2500 3704.5650 65 0.2500 3897.9012 66 0.2500 4093.0878 67 0.2500 4289.8902 68 0.2500 4488.1000 69 0.2500 4687.5332 70 0.2500 4888.0278 71 0.2500 5089.4418 72 0.2500 5291.6509 73 0.2500 5494.5469 74 0.2500 5698.0358 75 0.2500 5902.0358 ahm profile : jk ahm depth t-level 1 1.0000 0.5000 2 0.9999 1.5496 3 0.9999 2.6608 4 0.9998 3.8486 5 0.9997 5.1317 6 0.9997 6.5332 7 0.9996 8.0813 8 0.9995 9.8099 9 0.9994 11.7589 10 0.9993 13.9741 11 0.9991 16.5062 12 0.9989 19.4084 13 0.9987 22.7342 14 0.9985 26.5331 15 0.9983 30.8479 16 0.9980 35.7125 17 0.9976 41.1516 18 0.9972 47.1829 19 0.9968 53.8210 20 0.9963 61.0822 21 0.9958 68.9897 22 0.9952 77.5772 23 0.9945 86.8927 24 0.9937 97.0012 25 0.9928 107.9858 26 0.9918 119.9503 27 0.9907 133.0200 28 0.9894 147.3432 29 0.9878 163.0931 30 0.9860 180.4691 31 0.9839 199.6984 32 0.9813 221.0377 33 0.9783 244.7739 34 0.9745 271.2251 35 0.9699 300.7404 36 0.9641 333.6988 37 0.9568 370.5065 38 0.9475 411.5933 39 0.9355 457.4060 40 0.9200 508.4012 41 0.8997 565.0351 42 0.8732 627.7513 43 0.8386 696.9682 44 0.7942 773.0642 45 0.7387 856.3641 46 0.6721 947.1256 47 0.5972 1045.5281 48 0.5195 1151.6649 49 0.4462 1265.5388 50 0.3839 1387.0616 51 0.3360 1516.0589 52 0.3024 1652.2771 53 0.2804 1795.3952 54 0.2670 1945.0372 55 0.2592 2100.7868 56 0.2549 2262.2009 57 0.2525 2428.8236 58 0.2513 2600.1976 59 0.2506 2775.8744 60 0.2503 2955.4225 61 0.2501 3138.4329 62 0.2501 3324.5236 63 0.2500 3513.3418 64 0.2500 3704.5650 65 0.2500 3897.9012 66 0.2500 4093.0878 67 0.2500 4289.8902 68 0.2500 4488.1000 69 0.2500 4687.5332 70 0.2500 4888.0278 71 0.2500 5089.4418 72 0.2500 5291.6509 73 0.2500 5494.5469 74 0.2500 5698.0358 75 0.2500 5902.0358 ORCA R1, R2 or R4: overwrite the previous definition of ahm ================= ldfdyn_c3d_orca : 3D eddy viscosity coefficient ~~~~~~~~~~~~~~~ orca R1, R2 or R4 configuration: reduced in the surface Eq. strip iom_nf90_open ~~~ open existing file: ahmcoef.nc in READ mo de ---> ahmcoef.nc OK read icof (rec: 1) in ahmcoef.nc ok iom_close ~~~ close file: ahmcoef.nc ok 1D zcoef array ~~~~~~~~~~~~~~ jk zcoef 1 1.000 2 1.000 3 1.000 4 1.000 5 1.000 6 1.000 7 1.000 8 1.000 9 1.000 10 1.000 11 1.000 12 1.000 13 1.000 14 1.000 15 1.000 16 1.000 17 1.000 18 1.000 19 1.000 20 1.000 21 1.000 22 1.000 23 1.000 24 1.000 25 1.000 26 1.000 27 1.000 28 1.000 29 1.000 30 1.000 31 1.000 32 1.000 33 1.000 34 1.000 35 1.000 36 1.000 37 1.000 38 1.000 39 1.000 40 1.000 41 1.000 42 1.000 43 1.000 44 1.000 45 1.000 46 2.000 47 2.000 48 2.000 49 3.000 50 3.000 51 4.000 52 4.000 53 5.000 54 6.000 55 6.000 56 7.000 57 8.000 58 8.000 59 9.000 60 10.000 61 10.000 62 10.000 63 10.000 64 10.000 65 10.000 66 10.000 67 10.000 68 10.000 69 10.000 70 10.000 71 10.000 72 10.000 73 10.000 74 10.000 75 10.000 3D ahm1 array (k=1) 1 17 0.00 3D ahm2 array (k=1) 1 17 20.00 3D ahm2 array (k=jpk) 1 17 20.00 ldf_slp_init : direction of lateral mixing ~~~~~~~~~~~~ tra_qsr_init : penetration of the surface solar radiation ~~~~~~~~~~~~ Namelist namtra_qsr : set the parameter of penetration Light penetration (T) or not (F) ln_traqsr = T RGB (Red-Green-Blue) light penetration ln_qsr_rgb = F 2 band light penetration ln_qsr_2bd = F bio-model light penetration ln_qsr_bio = T light penetration for ice-model LIM3 ln_qsr_ice = T RGB : Chl data (=1/2) or cst value (=0) nn_chldta = 1 RGB & 2 bands: fraction of light (rn_si1) rn_abs = 0.580000000000000 RGB & 2 bands: shortess depth of extinction rn_si0 = 0.350000000000000 2 bands: longest depth of extinction rn_si1 = 23.0000000000000 bio-model light penetration read fraqsr_1lev (rec: 1) in ./restartopa_0000.nc ok tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating ~~~~~~~ Namelist nambbc : set bbc parameters Apply a geothermal heating at ocean bottom ln_trabbc = T type of geothermal flux nn_geoflx = 2 Constant geothermal flux value rn_geoflx_cst = 8.640000000000000E-002 *** variable geothermal heat flux tra_bbc_init : bottom temperature boundary condition ~~~~~~~~~~~~ nambbc Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./Goutorbe_ghflux.nc variable name: gh_flux frequency: -12.0000000000000 time interp: F climatology: T weights : ./weights_Goutorbe1_2_eorca1_bilinear.nc pairing : data type: yearly land/sea mask: iom_nf90_open ~~~ open existing file: ./Goutorbe_ghflux.nc in READ mode ---> ./Goutorbe_ghflux.nc OK iom_nf90_open ~~~ open existing file: ./Goutorbe_ghflux.nc in READ mode ---> ./Goutorbe_ghflux.nc OK iom_close ~~~ close file: ./Goutorbe_ghflux.nc ok iom_nf90_open ~~~ open existing file: ./weights_Goutorbe1_2 _eorca1_bilinear.nc in READ mode ---> ./weights_Goutorbe1_2_eorca1_bilinear.nc OK ===>>> : W A R N I N G =============== iom_nf90_getatt: no attribute ew_wrap found read src01 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read src02 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read src03 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read src04 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read wgt01 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read wgt02 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read wgt03 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok read wgt04 (rec: 1) in ./weights_Goutorbe1_2_eorca1_bilinear.nc ok iom_close ~~~ close file: ./weights_Goutorbe1_2_eorca1_bili near.nc ok read gh_flux (rec: 1) in ./Goutorbe_ghflux.nc ok fld_read: var gh_flux kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, record: 0001 (days 0.0000 <-> 366.0000) tra_bbl_init : bottom boundary layer initialisation ~~~~~~~~~~~~ Namelist nambbl : set bbl parameters diffusive bbl (=1) or not (=0) nn_bbl_ldf = 1 advective bbl (=1/2) or not (=0) nn_bbl_adv = 0 diffusive bbl coefficient rn_ahtbbl = 1000.00000000000 m2/s advective bbl coefficient rn_gambbl = 10.0000000000000 s tra_dmp_init : T and S newtonian relaxation ~~~~~~~ Namelist namtra_dmp : set relaxation parameters Apply relaxation or not ln_tradmp = F mixed layer damping option nn_zdmp = 0 Damping file name cn_resto = resto.nc tra_adv_init : choice/control of the tracer advection scheme ~~~~~~~~~~~ Namelist namtra_adv : chose a advection scheme for tracers 2nd order advection scheme ln_traadv_cen2 = F TVD advection scheme ln_traadv_tvd = T MUSCL advection scheme ln_traadv_muscl = F MUSCL2 advection scheme ln_traadv_muscl2 = F UBS advection scheme ln_traadv_ubs = F QUICKEST advection scheme ln_traadv_qck = F upstream scheme within muscl ln_traadv_msc_ups = F TVD advection scheme with zts ln_traadv_tvd_zts = F TVD scheme is used tra_adv_mle_init : mixed layer eddy (MLE) advection acting on tracers ~~~~~~~~~~~~~~~~ Namelist namtra_adv_mle : mixed layer eddy advection on tracers use mixed layer eddy (MLE, i.e. Fox-Kemper param) (T/F) ln_mle = T MLE type: =0 standard Fox-Kemper ; =1 new formulation nn_mle = 1 magnitude of the MLE (typical value: 0.06 to 0.08) rn_ce = 6.000000000000000E-002 scale of ML front (ML radius of deformation) (rn_mle=0) rn_lf = 5000.00000000000 m maximum time scale of MLE (rn_mle=0) rn_time = 172800.000000000 s reference latitude (degrees) of MLE coef. (rn_mle=1) rn_lat = 20.0000000000000 deg space interp. of MLD at u-(v-)pts (0=min,1=averaged,2=max) nn_mld_uv = 0 =1 no MLE in case of convection ; =0 always MLE nn_conv = 0 Density difference used to define ML for FK rn_rho_c_mle = 1.000000000000000E-002 Mixed Layer Eddy induced transport added to tracer advection New formulation ML buoyancy criteria = 9.558138401559454E-005 m/s2 associated ML density criteria defined in zdfmxl = 1.000000000000000E-002 kg/m3 tra_ldf_init : lateral tracer diffusive operator ~~~~~~~~~~~ Namelist namtra_ldf already read in ldftra module see ldf_tra_init report for lateral mixing parameters Rotated laplacian operator tra:ldf_ano : lateral diffusion acting on the full fields ~~~~~~~~~~~ tra_zdf_init : vertical tracer physics scheme ~~~~~~~~~~~ Implicit (euler backward) scheme dyn_adv_init : choice/control of the momentum advection scheme ~~~~~~~~~~~ Namelist namdyn_adv : chose a advection formulation & scheme for momentu m Vector/flux form (T/F) ln_dynadv_vec = T = 0 standard scheme ; =1 Hollingsworth scheme nn_dynkeg = 1 2nd order centred advection scheme ln_dynadv_cen2 = F 3rd order UBS advection scheme ln_dynadv_ubs = F Sub timestepping of vertical advection ln_dynzad_zts = F vector form : keg + zad + vor is used with Hollingsworth keg scheme dyn_vor_init : vorticity term : read namelist and control the consistency ~~~~~~~~~~~~ Namelist namdyn_vor : choice of the vorticity term scheme energy conserving scheme ln_dynvor_ene = F enstrophy conserving scheme ln_dynvor_ens = F mixed enstrophy/energy conserving scheme ln_dynvor_mix = F enstrophy and energy conserving scheme ln_dynvor_een = T enstrophy and energy conserving scheme (old) ln_dynvor_een_old= F Vector form advection : vorticity = Coriolis + relative vorticity vorticity scheme : energy and enstrophy conserving scheme dyn_ldf_init : Choice of the lateral diffusive operator on dynamics ~~~~~~~~~~~ Namelist nam_dynldf : set lateral mixing parameters (type, direction, co efficients) laplacian operator ln_dynldf_lap = T bilaplacian operator ln_dynldf_bilap = F iso-level ln_dynldf_level = F horizontal (geopotential) ln_dynldf_hor = T iso-neutral ln_dynldf_iso = F laplacian operator dyn_hpg_init : hydrostatic pressure gradient initialisation ~~~~~~~~~~~~ Namelist namdyn_hpg : choice of hpg scheme z-coord. - full steps ln_hpg_zco = F z-coord. - partial steps (interpolation) ln_hpg_zps = F s-coord. (standard jacobian formulation) ln_hpg_sco = T s-coord. (standard jacobian formulation) for isf ln_hpg_isf = F s-coord. (Density Jacobian: Cubic polynomial) ln_hpg_djc = F s-coord. (Pressure Jacobian: Cubic polynomial) ln_hpg_prj = F time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = F dyn_zdf_init : vertical dynamics physics scheme ~~~~~~~~~~~ Implicit (euler backward) scheme dyn_spg_init : choice of the surface pressure gradient scheme ~~~~~~~~~~~ Explicit free surface lk_dynspg_exp = F Free surface with time splitting lk_dynspg_ts = T Filtered free surface cst volume lk_dynspg_flt = F dyn_spg_ts : split-explicit free surface ~~~~~~~~~~ ln_ts_nn_auto=.true. Automatically set nn_baro Max. courant number allowed: 0.800000000000000 ln_bt_av=.true. => Time averaging over nn_baro time steps is on ln_bt_fw =.false.=> Centred integration of barotropic variables Time filter choice, nn_bt_flt: 1 Boxcar: width = nn_baro nn_baro = 37 Barotropic time step [s] is : 72.9729729729730 Maximum Courant number is : 0.781652123408619 free surface with time splitting scheme icbini : Namelist namberg ln_icebergs = F , NO icebergs used ~~~~~~~~ sto_par_init : stochastic parameterization ~~~~~~~~~~~~ Namelist namsto : stochastic parameterization restart stochastic parameters ln_rststo = F read seed of RNG from restart file ln_rstseed = T suffix of sto restart name (input) cn_storst_in = restart_sto suffix of sto restart name (output) cn_storst_out = restart_sto stochastic equation of state ln_sto_eos = F number of degrees of freedom nn_sto_eos = 1 random walk horz. std (in grid points) rn_eos_stdxy = 1.40000000000000 random walk vert. std (in grid points) rn_eos_stdz = 0.700000000000000 random walk tcor (in timesteps) rn_eos_tcor = 1440.00000000000 order of autoregressive processes nn_eos_ord = 1 passes of Laplacian filter nn_eos_flt = 0 limitation factor rn_eos_lim = 2.00000000000000 stochastic parameterization : trc_init : initial set up of the passive tracers ~~~~~~~ ===>>> : W A R N I N G =============== Coupling with passive tracers and used of diurnal cycle. Computation of a dai ly mean shortwave for some biogeochemical models) trc_nam : read the passive tracer namelists ~~~~~~~ Namelist : namtrc_run time step freq. for passive tracer nn_dttrc = 1 restart for passive tracer ln_rsttr = T control of time step for passive tracer nn_rsttr = 2 first time step for pass. trac. nittrc000 = 9855713 frequency of outputs for passive tracers nn_writetrc = 5475 Use euler integration for TRC (y/n) ln_top_euler = T trc_nam : read the passive tracer namelists ~~~~~~~ trc_nam_trp: read namelist for tracer transport ~~~~~~~~~~~ trc_adv_ctl : choice/control of the tracer advection scheme ~~~~~~~~~~~ Namelist namtrc_adv : chose a advection scheme for tracers 2nd order advection scheme ln_trcadv_cen2 = F TVD advection scheme ln_trcadv_tvd = F MUSCL advection scheme ln_trcadv_muscl = T MUSCL2 advection scheme ln_trcadv_muscl2 = F UBS advection scheme ln_trcadv_ubs = F QUICKEST advection scheme ln_trcadv_qck = F trc:ldf_ctl : lateral tracer diffusive operator ~~~~~~~~~~~ Namelist namtrc_ldf : set lateral mixing parameters (type, direction, coeffi cients) laplacian operator ln_trcldf_lap = T bilaplacian operator ln_trcldf_bilap = F iso-level ln_trcldf_level = F horizontal (geopotential) ln_trcldf_hor = F iso-neutral ln_trcldf_iso = T diffusivity coefficient rn_ahtrc_0 = 2000.00000000000 background hor. diffusivity rn_ahtrb_0 = 0.000000000000000E+000 enhanced zonal diffusivity rn_fact_lap = 15.0000000000000 Namelist namtrc_zdf : set vertical diffusion parameters time splitting / backward scheme ln_trczdf_exp = F number of time step nn_trczdf_exp = 3 Namelist namtrc_rad : treatment of negative concentrations correct artificially negative concen. or not ln_trcrad = T trc_dmp : Passive tracers newtonian damping ~~~~~~~ Namelist namtrc_dmp : set damping parameter mixed layer damping option nn_zdmp_tr = 1 (zoom: forced to 0) Restoration coeff file cn_resto_tr = resto_tr.nc trc_nam_ice : Read the namelist for trc_ice ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Sea ice tracers option (nn_ice_tr) : -1 trc_nam_pisces : read PISCES namelist ~~~~~~~~~~~~~~ MY_TRC not used CFCs trc_nam: Read namdates, namelist for CFC chemical model ~~~~~~~ initial calendar date (aammjj) for CFC ndate_beg = 300101 restoring time constant (year) nyear_res = 1932 Atmospheric CFC concentrations file clnamecfc = CFCs_CDIAC.dat Compute dynamics for CFC-11 lp_cfc11 = T Compute dynamics for CFC-12 lp_cfc12 = T Compute dynamics for SF6 lp_sf6 = T initial year (aa) nyear_beg = 30 C14 not used trc_nam_age: Read namage, namelist for Age passive tracer ~~~~~~~ depth over which age tracer reset to zero rn_age _depth = 10.0000000000000 recip of relax. timescale (s) for age tracer shallower than age_depth rn_age _kill_rate = -1.388888888888889E-004 Namelist : namtrc Read inputs data from file (y/n) ln_trcdta = T Damping of passive tracer (y/n) ln_trcdmp = F Restoring of tracer on closed seas ln_trcdmp_clo = T tracer nb : 1 short name : DIC tracer nb : 2 short name : Alkalini tracer nb : 3 short name : O2 tracer nb : 4 short name : CaCO3 tracer nb : 5 short name : PO4 tracer nb : 6 short name : POC tracer nb : 7 short name : Si tracer nb : 8 short name : PHY tracer nb : 9 short name : ZOO tracer nb : 10 short name : DOC tracer nb : 11 short name : PHY2 tracer nb : 12 short name : ZOO2 tracer nb : 13 short name : DSi tracer nb : 14 short name : Fer tracer nb : 15 short name : BFe tracer nb : 16 short name : GOC tracer nb : 17 short name : SFe tracer nb : 18 short name : DFe tracer nb : 19 short name : GSi tracer nb : 20 short name : NFe tracer nb : 21 short name : NCHL tracer nb : 22 short name : DCHL tracer nb : 23 short name : NO3 tracer nb : 24 short name : NH4 tracer nb : 25 short name : CFC11 tracer nb : 26 short name : CFC12 tracer nb : 27 short name : SF6 tracer nb : 28 short name : Age Read a restart file for passive tracer : restart_trc Passive Tracer time step rdttrc = 2700.00000000000 trc_rst_cal : read the TOP restart file for calendar ~~~~~~~~~~~~ iom_nf90_open ~~~ open existing file: ./restart_trc_0000.nc in READ mode ---> ./restart_trc_0000.nc OK *** Info read in restart : previous time-step : 9855712 *** restart option nn_rsttr = 2 : calendar parameters read in restart p4z_ini : PISCES biochemical model initialisation ~~~~~~~~~~~~~~ Namelist : nampisbio frequence pour la biologie nrdttrc = 1 POC sinking speed wsbio = 2.00000000000000 half saturation constant for mortality xkmort = 2.000000000000000E-007 Fe/C in zooplankton ferat3 = 1.000000000000000E-005 Big particles sinking speed wsbio2 = 30.0000000000000 Maximum number of iterations for POC niter1max = 1 Maximum number of iterations for GOC niter2max = 1 Namelist : nampisdmp Relaxation of tracer to glodap mean value ln_pisdmp = T Frequency of Relaxation nn_pisdmp = 11680 Namelist parameter for mass conservation checking ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Flag to check mass conservation of NO3/Si/TALK ln_check_mass = F Level corresponding to 100m depth 25 namelist : nampisopt ~~~~~~~~~~~~~~~~~ PAR as a variable fraction of SW ln_varpar = T Default value for the PAR fraction parlux = 0.430000000000000 initialize variable par fraction ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ p4z_opt_init : Variable PAR fraction ~~~~~~~~~~~~ nampisopt Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./par_fraction_gewex_clim90s00s_eORCA_R1.nc variable name: fr_par frequency: 24.0000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: iom_nf90_open ~~~ open existing file: par_fraction_gewex_cl im90s00s_eORCA_R1.nc in READ mode ---> par_fraction_gewex_clim90s00s_eORCA_R1.nc OK trc_oce_rgb : Initialisation of the optical look-up table ~~~~~~~~~~~ RGB longest depth of extinction r_si2 = 61.8046971569839 jc = 1 Chl = 1.000000000000000E-002 irgb = 1 jc = 2 Chl = 1.100000000000000E-002 irgb = 2 jc = 3 Chl = 1.300000000000000E-002 irgb = 3 jc = 4 Chl = 1.400000000000000E-002 irgb = 4 jc = 5 Chl = 1.600000000000000E-002 irgb = 5 jc = 6 Chl = 1.800000000000000E-002 irgb = 6 jc = 7 Chl = 2.000000000000000E-002 irgb = 7 jc = 8 Chl = 2.200000000000000E-002 irgb = 8 jc = 9 Chl = 2.500000000000000E-002 irgb = 9 jc = 10 Chl = 2.800000000000000E-002 irgb = 10 jc = 11 Chl = 3.200000000000000E-002 irgb = 11 jc = 12 Chl = 3.500000000000000E-002 irgb = 12 jc = 13 Chl = 4.000000000000000E-002 irgb = 13 jc = 14 Chl = 4.500000000000000E-002 irgb = 14 jc = 15 Chl = 5.000000000000000E-002 irgb = 15 jc = 16 Chl = 5.600000000000000E-002 irgb = 16 jc = 17 Chl = 6.300000000000000E-002 irgb = 17 jc = 18 Chl = 7.099999999999999E-002 irgb = 18 jc = 19 Chl = 7.900000000000000E-002 irgb = 19 jc = 20 Chl = 8.900000000000000E-002 irgb = 20 jc = 21 Chl = 0.100000000000000 irgb = 21 jc = 22 Chl = 0.112000000000000 irgb = 22 jc = 23 Chl = 0.126000000000000 irgb = 23 jc = 24 Chl = 0.141000000000000 irgb = 24 jc = 25 Chl = 0.158000000000000 irgb = 25 jc = 26 Chl = 0.178000000000000 irgb = 26 jc = 27 Chl = 0.200000000000000 irgb = 27 jc = 28 Chl = 0.224000000000000 irgb = 28 jc = 29 Chl = 0.251000000000000 irgb = 29 jc = 30 Chl = 0.282000000000000 irgb = 30 jc = 31 Chl = 0.316000000000000 irgb = 31 jc = 32 Chl = 0.355000000000000 irgb = 32 jc = 33 Chl = 0.398000000000000 irgb = 33 jc = 34 Chl = 0.447000000000000 irgb = 34 jc = 35 Chl = 0.501000000000000 irgb = 35 jc = 36 Chl = 0.562000000000000 irgb = 36 jc = 37 Chl = 0.631000000000000 irgb = 37 jc = 38 Chl = 0.708000000000000 irgb = 38 jc = 39 Chl = 0.794000000000000 irgb = 39 jc = 40 Chl = 0.891000000000000 irgb = 40 jc = 41 Chl = 1.00000000000000 irgb = 41 jc = 42 Chl = 1.12200000000000 irgb = 42 jc = 43 Chl = 1.25900000000000 irgb = 43 jc = 44 Chl = 1.41300000000000 irgb = 44 jc = 45 Chl = 1.58500000000000 irgb = 45 jc = 46 Chl = 1.77800000000000 irgb = 46 jc = 47 Chl = 1.99500000000000 irgb = 47 jc = 48 Chl = 2.23900000000000 irgb = 48 jc = 49 Chl = 2.51200000000000 irgb = 49 jc = 50 Chl = 2.81800000000000 irgb = 50 jc = 51 Chl = 3.16200000000000 irgb = 51 jc = 52 Chl = 3.54800000000000 irgb = 52 jc = 53 Chl = 3.98100000000000 irgb = 53 jc = 54 Chl = 4.46700000000000 irgb = 54 jc = 55 Chl = 5.01200000000000 irgb = 55 jc = 56 Chl = 5.62300000000000 irgb = 56 jc = 57 Chl = 6.31000000000000 irgb = 57 jc = 58 Chl = 7.07900000000000 irgb = 58 jc = 59 Chl = 7.94300000000000 irgb = 59 jc = 60 Chl = 8.91200000000000 irgb = 60 jc = 61 Chl = 10.0000000000000 irgb = 61 level of light extinction = 46 ref depth = 995.688512190334 m Namelist parameters for nutrient limitations, nampislim ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ mean rainratio caco3r = 0.300000000000000 NO3 half saturation of nanophyto concnno3 = 1.000000000000000E-006 NO3 half saturation of diatoms concdno3 = 3.000000000000000E-006 NH4 half saturation for phyto concnnh4 = 1.000000000000000E-007 NH4 half saturation for diatoms concdnh4 = 3.000000000000000E-007 half saturation constant for Si uptake xksi1 = 2.000000000000000E-006 half saturation constant for Si/C xksi2 = 2.000000000000000E-005 half-sat. of DOC remineralization xkdoc = 4.170000000000000E-004 Iron half saturation for nanophyto concnfer = 1.000000000000000E-009 Iron half saturation for diatoms concdfer = 3.000000000000000E-009 size ratio for nanophytoplankton xsizern = 3.00000000000000 size ratio for diatoms xsizerd = 3.00000000000000 NO3 half saturation of bacteria concbno3 = 2.000000000000000E-007 NH4 half saturation for bacteria concbnh4 = 2.000000000000000E-008 Minimum size criteria for diatoms xsizedia = 1.000000000000000E-006 Minimum size criteria for nanophyto xsizephy = 1.000000000000000E-006 Fe half saturation for bacteria concbfe = 9.999999999999999E-012 halk saturation constant for anoxia oxymin = 1.000000000000000E-006 optimal Fe quota for nano. qnfelim = 7.000000000000000E-006 Optimal Fe quota for diatoms qdfelim = 7.000000000000000E-006 Namelist parameters for phytoplankton growth, nampisprod ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Enable new parame. of production (T/F) ln_newprod = T mean Si/C ratio grosip = 0.159000000000000 P-I slope pislope = 2.00000000000000 Acclimation factor to low light xadap = 0.000000000000000E+000 excretion ratio of nanophytoplankton excret = 5.000000000000000E-002 excretion ratio of diatoms excret2 = 5.000000000000000E-002 basal respiration in phytoplankton bresp = 3.300000000000000E-002 Maximum Chl/C in phytoplankton chlcmin = 4.000000000000000E-003 P-I slope for diatoms pislope2 = 2.00000000000000 Minimum Chl/C in nanophytoplankton chlcnm = 3.300000000000000E-002 Minimum Chl/C in diatoms chlcdm = 5.000000000000000E-002 Maximum Fe/C in nanophytoplankton fecnm = 4.000000000000000E-005 Minimum Fe/C in diatoms fecdm = 4.000000000000000E-005 namelist : nampissbc ~~~~~~~~~~~~~~~~~ dust input from the atmosphere ln_dust = T Variable solubility of iron input ln_solub = T river input of nutrients ln_river = T atmospheric deposition of n ln_ndepo = T Fe input from sediments ln_ironsed = T Fe input from seaice ln_ironice = T fe input from hydrothermal vents ln_hydrofe = F coastal release of iron sedfeinput = 1.000000000000000E-009 solubility of the dust dustsolub = 2.000000000000000E-002 Mineral Fe content of the dust mfrac = 3.500000000000000E-002 Iron concentration in sea ice icefeinput = 1.500000000000000E-008 sinking speed of the dust wdust = 2.00000000000000 nitrogen fixation rate nitrfix = 1.000000000000000E-007 nitrogen fixation sensitivty to light diazolight = 50.0000000000000 fe half-saturation cste for diazotrophs concfediaz = 1.000000000000000E-010 Fe to 3He ratio assumed for vent iron supply hratio = 10000000.0000000 initialize dust input from atmosphere ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ p4z_sed_init : Atmospheric dust deposition ~~~~~~~~~~~~ nampissed Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./dust.orca.nc variable name: dust frequency: -1.00000000000000 time interp: T climatology: T weights : ./weights_lmd144142_bilin.nc pairing : data type: yearly land/sea mask: iom_nf90_open ~~~ open existing file: dust.orca.nc in READ mode ---> dust.orca.nc OK initialize variable solubility of Fe ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ p4z_sed_init : Solubility of atm. iron ~~~~~~~~~~~~ nampissed Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./Solubility_T62_Mahowald_eORCA_R1.nc variable name: solubility2 frequency: -12.0000000000000 time interp: F climatology: T weights : pairing : data type: yearly land/sea mask: iom_nf90_open ~~~ open existing file: Solubility_T62_Mahowa ld_eORCA_R1.nc in READ mode ---> Solubility_T62_Mahowald_eORCA_R1.nc OK iom_close ~~~ close file: Solubility_T62_Mahowald_eORCA_R1. nc ok p4z_sed_init : Input from river ~~~~~~~~~~~~ nampissed Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./river_global_news_eORCA_R1.nc variable name: riverdic frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./river_global_news_eORCA_R1.nc variable name: riverdoc frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./river_global_news_eORCA_R1.nc variable name: riverdin frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./river_global_news_eORCA_R1.nc variable name: riverdon frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./river_global_news_eORCA_R1.nc variable name: riverdip frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./river_global_news_eORCA_R1.nc variable name: riverdop frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./river_global_news_eORCA_R1.nc variable name: riverdsi frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdic (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdic (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdoc (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdoc (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdin (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdin (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdon (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdon (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdip (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdip (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdop (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdop (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok iom_nf90_open ~~~ open existing file: river_global_news_eOR CA_R1.nc in READ mode ---> river_global_news_eORCA_R1.nc OK read riverdsi (rec: 1) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 2) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 3) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 4) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 5) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 6) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 7) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 8) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 9) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 10) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 11) in river_global_news_eORCA_R1.nc ok read riverdsi (rec: 12) in river_global_news_eORCA_R1.nc ok iom_close ~~~ close file: river_global_news_eORCA_R1.nc ok initialize the nutrient input by dust from ndeposition.orca.nc ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ p4z_sed_init : Nutrient atmospheric depositon ~~~~~~~~~~~~ nampissed Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./ndeposition.orca.nc variable name: ndep frequency: -1.00000000000000 time interp: T climatology: T weights : ./weights_2d_bilin.nc pairing : data type: yearly land/sea mask: iom_nf90_open ~~~ open existing file: ndeposition.orca.nc i n READ mode ---> ndeposition.orca.nc OK computation of an island mask to enhance coastal supply of iron ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ iom_nf90_open ~~~ open existing file: pmarge_etopo_eORCA_R1 .nc in READ mode ---> pmarge_etopo_eORCA_R1.nc OK read bathy (rec: 1) in pmarge_etopo_eORCA_R1.nc ok iom_close ~~~ close file: pmarge_etopo_eORCA_R1.nc ok Level corresponding to 50m depth 18 53.8506372116824 iom_nf90_open ~~~ open existing file: ./dust.orca.nc in REA D mode ---> ./dust.orca.nc OK iom_nf90_open ~~~ open existing file: ./dust.orca.nc in REA D mode ---> ./dust.orca.nc OK iom_close ~~~ close file: ./dust.orca.nc ok iom_nf90_open ~~~ open existing file: ./weights_lmd144142_b ilin.nc in READ mode ---> ./weights_lmd144142_bilin.nc OK read src01 (rec: 1) in ./weights_lmd144142_bilin.nc ok read src02 (rec: 1) in ./weights_lmd144142_bilin.nc ok read src03 (rec: 1) in ./weights_lmd144142_bilin.nc ok read src04 (rec: 1) in ./weights_lmd144142_bilin.nc ok read wgt01 (rec: 1) in ./weights_lmd144142_bilin.nc ok read wgt02 (rec: 1) in ./weights_lmd144142_bilin.nc ok read wgt03 (rec: 1) in ./weights_lmd144142_bilin.nc ok read wgt04 (rec: 1) in ./weights_lmd144142_bilin.nc ok iom_close ~~~ close file: ./weights_lmd144142_bilin.nc ok read dust (rec: 3) in ./dust.orca.nc ok fld_init : time-interpolation for dust read previous record = 3 at time = 75.50 days weight file: ./weights_Goutorbe1_2_eorca1_bilinear.nc ddims: 360 180 numwgt: 4 jpiwgt: 23 jpjwgt: 6 botleft: 334 166 topright: 356 171 cyclical allocated weight file: ./weights_lmd144142_bilin.nc ddims: 144 143 numwgt: 4 jpiwgt: 10 jpjwgt: 5 botleft: 134 132 topright: 143 136 cyclical allocated read dust (rec: 4) in ./dust.orca.nc ok fld_read: var dust kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 weight file: ./weights_Goutorbe1_2_eorca1_bilinear.nc ddims: 360 180 numwgt: 4 jpiwgt: 23 jpjwgt: 6 botleft: 334 166 topright: 356 171 cyclical allocated weight file: ./weights_lmd144142_bilin.nc ddims: 144 143 numwgt: 4 jpiwgt: 10 jpjwgt: 5 botleft: 134 132 topright: 143 136 cyclical allocated iom_nf90_open ~~~ open existing file: ./Solubility_T62_Maho wald_eORCA_R1.nc in READ mode ---> ./Solubility_T62_Mahowald_eORCA_R1.nc OK read solubility2 (rec: 1) in ./Solubility_T62_Mahowald_eORCA_R1.nc ok fld_read: var solubility2 kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, record: 0001 (days 0.0000 <-> 366.0000) iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdic (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdic read previous record = 3 at time = 75.50 days iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdoc (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdoc read previous record = 3 at time = 75.50 days iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdin (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdin read previous record = 3 at time = 75.50 days iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdon (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdon read previous record = 3 at time = 75.50 days iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdip (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdip read previous record = 3 at time = 75.50 days iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdop (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdop read previous record = 3 at time = 75.50 days iom_nf90_open ~~~ open existing file: ./river_global_news_e ORCA_R1.nc in READ mode ---> ./river_global_news_eORCA_R1.nc OK read riverdsi (rec: 3) in ./river_global_news_eORCA_R1.nc ok fld_init : time-interpolation for riverdsi read previous record = 3 at time = 75.50 days weight file: ./weights_Goutorbe1_2_eorca1_bilinear.nc ddims: 360 180 numwgt: 4 jpiwgt: 23 jpjwgt: 6 botleft: 334 166 topright: 356 171 cyclical allocated weight file: ./weights_lmd144142_bilin.nc ddims: 144 143 numwgt: 4 jpiwgt: 10 jpjwgt: 5 botleft: 134 132 topright: 143 136 cyclical allocated read riverdic (rec: 4) in ./river_global_news_eORCA_R1.nc ok read riverdoc (rec: 4) in ./river_global_news_eORCA_R1.nc ok read riverdin (rec: 4) in ./river_global_news_eORCA_R1.nc ok read riverdon (rec: 4) in ./river_global_news_eORCA_R1.nc ok read riverdip (rec: 4) in ./river_global_news_eORCA_R1.nc ok read riverdop (rec: 4) in ./river_global_news_eORCA_R1.nc ok read riverdsi (rec: 4) in ./river_global_news_eORCA_R1.nc ok fld_read: var riverdic kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 fld_read: var riverdoc kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 fld_read: var riverdin kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 fld_read: var riverdon kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 fld_read: var riverdip kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 fld_read: var riverdop kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 fld_read: var riverdsi kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 iom_nf90_open ~~~ open existing file: ./ndeposition.orca.nc in READ mode ---> ./ndeposition.orca.nc OK iom_nf90_open ~~~ open existing file: ./ndeposition.orca.nc in READ mode ---> ./ndeposition.orca.nc OK iom_close ~~~ close file: ./ndeposition.orca.nc ok iom_nf90_open ~~~ open existing file: ./weights_2d_bilin.nc in READ mode ---> ./weights_2d_bilin.nc OK read src01 (rec: 1) in ./weights_2d_bilin.nc ok read src02 (rec: 1) in ./weights_2d_bilin.nc ok read src03 (rec: 1) in ./weights_2d_bilin.nc ok read src04 (rec: 1) in ./weights_2d_bilin.nc ok read wgt01 (rec: 1) in ./weights_2d_bilin.nc ok read wgt02 (rec: 1) in ./weights_2d_bilin.nc ok read wgt03 (rec: 1) in ./weights_2d_bilin.nc ok read wgt04 (rec: 1) in ./weights_2d_bilin.nc ok iom_close ~~~ close file: ./weights_2d_bilin.nc ok read ndep (rec: 3) in ./ndeposition.orca.nc ok fld_init : time-interpolation for ndep read previous record = 3 at time = 75.50 days weight file: ./weights_Goutorbe1_2_eorca1_bilinear.nc ddims: 360 180 numwgt: 4 jpiwgt: 23 jpjwgt: 6 botleft: 334 166 topright: 356 171 cyclical allocated weight file: ./weights_lmd144142_bilin.nc ddims: 144 143 numwgt: 4 jpiwgt: 10 jpjwgt: 5 botleft: 134 132 topright: 143 136 cyclical allocated weight file: ./weights_2d_bilin.nc ddims: 144 96 numwgt: 4 jpiwgt: 10 jpjwgt: 4 botleft: 62 6 topright: 71 9 cyclical allocated read ndep (rec: 4) in ./ndeposition.orca.nc ok fld_read: var ndep kt = 9855713 ( 91.0156 days), Y/M/D = 2000/04/01, records b/a: 0003/ 0004 (days 75.5000/ 106.0000) it_offset is : 0 Total input of elements from river supply ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N Supply : 37.6451048750217 TgN/yr Si Supply : 129.386583338715 TgSi/yr P Supply : 3.73248717459187 TgP/yr Alk Supply : 36.6109800087921 Teq/yr DIC Supply : 611.140107191845 TgC/yr Namelist parameters for Iron chemistry, nampisfer ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ enable complex iron chemistry scheme ln_fechem = F variable concentration of ligand ln_ligvar = F scavenging rate of Iron xlam1 = 5.000000000000000E-003 scavenging rate of Iron by dust xlamdust = 150.000000000000 ligand concentration in the ocean ligand = 6.000000000000000E-010 Namelist parameters for remineralization, nampisrem ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ remineralisation rate of POC xremip = 2.500000000000000E-002 remineralization rate of DOC xremik = 0.300000000000000 remineralization rate of Si xsirem = 3.000000000000000E-003 fast remineralization rate of Si xsiremlab = 3.000000000000000E-002 fraction of labile biogenic silica xsilab = 0.500000000000000 NH4 nitrification rate nitrif = 5.000000000000000E-002 Namelist parameters for phytoplankton mortality, nampismort ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ quadratic mortality of phytoplankton wchl = 1.000000000000000E-002 maximum quadratic mortality of diatoms wchld = 1.000000000000000E-002 maximum quadratic mortality of diatoms wchldm = 3.000000000000000E-002 phytoplankton mortality rate mprat = 1.000000000000000E-002 Diatoms mortality rate mprat2 = 1.000000000000000E-002 Namelist parameters for microzooplankton, nampiszoo ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ part of calcite not dissolved in microzoo guts part = 0.500000000000000 microzoo preference for POC xpref2c = 0.100000000000000 microzoo preference for nano xpref2p = 1.00000000000000 microzoo preference for diatoms xpref2d = 0.500000000000000 diatoms feeding threshold for microzoo xthreshdia = 1.000000000000000E-008 nanophyto feeding threshold for microzoo xthreshphy = 1.000000000000000E-008 poc feeding threshold for microzoo xthreshpoc = 1.000000000000000E-008 feeding threshold for microzooplankton xthresh = 3.000000000000000E-007 exsudation rate of microzooplankton resrat = 3.000000000000000E-002 microzooplankton mortality rate mzrat = 4.000000000000000E-003 maximal microzoo grazing rate grazrat = 3.00000000000000 non assimilated fraction of P by microzoo unass = 0.300000000000000 Efficicency of microzoo growth epsher = 0.300000000000000 Fraction of microzoo excretion as DOM sigma1 = 0.600000000000000 half sturation constant for grazing 1 xkgraz = 2.000000000000000E-005 Namelist parameters for mesozooplankton, nampismes ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ part of calcite not dissolved in mesozoo guts part2 = 0.750000000000000 mesozoo preference for phyto xprefc = 1.00000000000000 mesozoo preference for POC xprefp = 0.300000000000000 mesozoo preference for zoo xprefz = 1.00000000000000 mesozoo preference for poc xprefpoc = 0.300000000000000 microzoo feeding threshold for mesozoo xthresh2zoo = 1.000000000000000E-008 diatoms feeding threshold for mesozoo xthresh2dia = 1.000000000000000E-008 nanophyto feeding threshold for mesozoo xthresh2phy = 1.000000000000000E-008 poc feeding threshold for mesozoo xthresh2poc = 1.000000000000000E-008 feeding threshold for mesozooplankton xthresh2 = 3.000000000000000E-007 exsudation rate of mesozooplankton resrat2 = 5.000000000000000E-003 mesozooplankton mortality rate mzrat2 = 3.000000000000000E-002 maximal mesozoo grazing rate grazrat2 = 0.750000000000000 mesozoo flux feeding rate grazflux = 2000.00000000000 non assimilated fraction of P by mesozoo unass2 = 0.300000000000000 Efficicency of Mesozoo growth epsher2 = 0.350000000000000 Fraction of mesozoo excretion as DOM sigma2 = 0.600000000000000 half sturation constant for grazing 2 xkgraz2 = 2.000000000000000E-005 Namelist parameters for CaCO3 dissolution, nampiscal ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ diss. rate constant calcite (per month) kdca = 6.00000000000000 order of reaction for calcite dissolution nca = 1.00000000000000 Namelist parameters for air-sea exchange, nampisext ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Choice for reading in the atm pCO2 file or constant value, ln_co2int = F Constant Atmospheric pCO2 value atcco2 = 369.120000000000 Namelist nampisatm : Atmospheric Pressure as external forcing constant atmopsheric pressure (F) or from a file (T) ln_presatm = F Initialization of PISCES tracers done trc_ini_cfc: initialisation of CFC chemical model ~~~~~~~~~~~ Read annual atmospheric concentrations from formatted file : CFCs_CDIAC.dat ---> 115 years read Year c11NH c11SH c12NH c12SH SF6NH SF6SH 30 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 31 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 32 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 33 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 34 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 35 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 36 0.0000 0.0000 0.1000 0.0000 0.0000 0.0000 37 0.0000 0.0000 0.1000 0.1000 0.0000 0.0000 38 0.0000 0.0000 0.2000 0.1000 0.0000 0.0000 39 0.0000 0.0000 0.3000 0.2000 0.0000 0.0000 40 0.0000 0.0000 0.4000 0.3000 0.0000 0.0000 41 0.0000 0.0000 0.5000 0.4000 0.0000 0.0000 42 0.0000 0.0000 0.7000 0.5000 0.0000 0.0000 43 0.0000 0.0000 0.9000 0.7000 0.0000 0.0000 44 0.0000 0.0000 1.2000 0.9000 0.0000 0.0000 45 0.1000 0.0000 1.7000 1.2000 0.0000 0.0000 46 0.1000 0.1000 2.3000 1.7000 0.0000 0.0000 47 0.1000 0.1000 3.4000 2.4000 0.0000 0.0000 48 0.2000 0.1000 4.8000 3.4000 0.0000 0.0000 49 0.4000 0.2000 6.1000 4.7000 0.0000 0.0000 50 0.7000 0.4000 7.6000 6.0000 0.0000 0.0000 51 1.0100 0.7000 9.2000 7.4000 0.0000 0.0000 52 1.5100 1.0100 11.0000 9.0000 0.0000 0.0000 53 2.2100 1.5100 12.8000 10.7000 0.0400 0.0400 54 3.0200 2.2100 15.0000 12.6000 0.0400 0.0400 55 4.1200 3.0200 17.4000 14.7000 0.0400 0.0400 56 5.3300 4.0200 20.2000 17.1000 0.0400 0.0400 57 6.8300 5.2300 23.4000 19.9000 0.0500 0.0400 58 8.1400 6.5300 26.8000 23.0000 0.0500 0.0400 59 9.4500 7.8400 30.5000 26.3000 0.0500 0.0500 60 11.0600 9.1500 35.0000 30.1000 0.0500 0.0500 61 13.2700 10.8500 40.0000 34.4000 0.0600 0.0600 62 16.1800 13.0700 45.8000 39.4000 0.0700 0.0600 63 19.6000 15.7800 52.5000 45.1000 0.0800 0.0700 64 23.7200 19.2000 60.4000 51.8000 0.0900 0.0800 65 28.4400 23.1200 69.3000 59.5000 0.1100 0.1000 66 33.6700 27.6400 79.2000 68.2000 0.1300 0.1200 67 39.4000 32.6600 90.3000 77.9000 0.1500 0.1400 68 46.0300 38.2900 102.8000 88.8000 0.1800 0.1700 69 53.7700 44.8200 116.8000 101.1000 0.2100 0.1900 70 62.4100 52.2600 132.0000 114.7000 0.2300 0.2200 71 72.0600 60.7000 148.4000 129.6000 0.2600 0.2400 72 82.7100 69.9500 166.1000 145.7000 0.3000 0.2800 73 94.8700 80.4000 185.8000 163.3000 0.3400 0.3100 74 108.3400 92.1600 207.1000 182.5000 0.3800 0.3500 75 121.4100 104.7200 228.2000 202.9000 0.4400 0.4000 76 133.9700 117.0900 248.1000 223.2000 0.5000 0.4600 77 145.9300 129.3500 266.9000 242.7000 0.5800 0.5300 78 156.5800 140.8000 284.3000 261.2000 0.6600 0.6100 79 168.3400 148.7400 306.1000 273.5000 0.7600 0.7000 80 176.6800 159.3000 323.2000 292.3000 0.8800 0.8100 81 184.3200 167.8400 339.6000 308.8000 1.0000 0.9300 82 191.4600 176.0800 353.4000 325.5000 1.1300 1.0400 83 199.3000 184.5200 369.0000 342.6000 1.2700 1.1700 84 208.0400 192.4600 385.7000 359.4000 1.4000 1.2900 85 217.9900 202.0100 403.4000 378.2000 1.5500 1.4300 86 229.3500 211.3600 424.3000 396.5000 1.7100 1.5800 87 241.6100 222.2100 444.0000 416.3000 1.8800 1.7300 88 252.8600 233.2700 465.4000 435.8000 2.0500 1.8900 89 259.3000 242.1100 483.6000 454.4000 2.2200 2.0500 90 265.8300 251.0600 497.7000 472.7000 2.4100 2.2200 91 268.2400 256.6800 506.0000 487.3000 2.6200 2.4200 92 268.1400 260.8000 516.3000 498.3000 2.8500 2.6300 93 269.5500 262.5100 523.2000 507.0000 3.0900 2.8400 94 269.6500 263.7200 528.5000 514.8000 3.3300 3.0700 95 268.3400 263.2200 533.4000 521.0000 3.5800 3.3000 96 266.9300 262.9100 537.3000 526.5000 3.8600 3.5500 97 265.7300 262.0100 540.1000 530.8000 4.0700 3.8100 98 264.5200 261.0100 542.9000 534.3000 4.2500 4.0100 99 263.1200 259.9000 544.4000 537.2000 4.4800 4.2500 100 261.7100 258.2900 545.9000 539.0000 4.6800 4.4600 101 260.0000 256.9800 546.5000 540.6000 4.9000 4.6500 102 258.1900 255.0800 546.7000 541.3000 5.1400 4.8800 103 256.1800 253.2700 546.7000 541.6000 5.3700 5.0700 104 253.9700 251.3600 545.7000 541.5000 5.5800 5.3000 105 251.9600 249.1500 544.9000 540.7000 5.8000 5.5100 106 249.5500 247.3400 543.1000 539.8000 6.0400 5.7500 107 247.5400 245.0300 541.1000 538.1000 6.3300 6.0300 108 245.6300 243.1200 538.6000 536.2000 6.6200 6.2900 109 243.6100 241.0500 536.1200 533.4700 6.8900 6.5800 110 241.3100 239.1500 533.2600 531.0300 7.2100 6.8700 111 239.4000 236.8400 530.6000 528.5300 7.4600 7.1700 112 236.8900 234.7700 527.2000 526.0000 7.7300 7.4300 113 235.5700 232.8800 525.2200 523.2500 8.0900 7.7300 114 234.9700 231.4900 522.9700 521.0000 8.4000 8.1100 115 234.3600 230.0900 520.7100 518.7400 8.7200 8.5000 Initialization of CFC tracers done trc_ini_age: passive tracer age ~~~~~~~~~~~~~~ trc_ice_ini : Initialize sea ice tracer boundary condition ~~~~~~~~~~~~~ file : tracer.stat open ok unit = 33 status = REPLACE form = FORMATTED access = SEQUENTIAL trc_dta_init : Passive tracers Initial Conditions ~~~~~~~~~~~~~~ number of passive tracers to be initialize by data : 8 ===>>> : W A R N I N G =============== trc_dta_init: passive tracer data initialisation Input name of data file : TDIC differs from that of tracer : DIC Read IC file for tracer number : 1, name : TDIC, Multiplicative Scaling factor : 0.103E-05 Read IC file for tracer number : 2, name : Alkalini, Multiplicative Scaling factor : 0.103E-05 Read IC file for tracer number : 3, name : O2, Multiplicative Scaling factor : 0.446E-04 Read IC file for tracer number : 5, name : PO4, Multiplicative Scaling factor : 0.122E-03 Read IC file for tracer number : 7, name : Si, Multiplicative Scaling factor : 0.100E-05 Read IC file for tracer number : 10, name : DOC, Multiplicative Scaling factor : 0.100E-05 Read IC file for tracer number : 14, name : Fer, Multiplicative Scaling factor : 0.100E-05 Read IC file for tracer number : 23, name : NO3, Multiplicative Scaling factor : 0.760E-05 trc_dta_init : Passive tracer data ~~~~~~~~~~~~ namtrc Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./DIC_GLODAPv2.1_annual_eORCA_R1.nc variable name: TDIC frequency: -12.0000000000000 time interp: F climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./Alkalini_GLODAPv2.1_annual_eORCA_R1.nc variable name: Alkalini frequency: -12.0000000000000 time interp: F climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./O2_WOA2009_monthly_eORCA_R1.nc variable name: O2 frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./PO4_WOA2009_monthly_eORCA_R1.nc variable name: PO4 frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./Si_WOA2009_monthly_eORCA_R1.nc variable name: Si frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./DOC_PISCES_monthly_eORCA_R1.nc variable name: DOC frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./Fer_PISCES_monthly_eORCA_R1.nc variable name: Fer frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: root filename: ./NO3_WOA2009_monthly_eORCA_R1.nc variable name: NO3 frequency: -1.00000000000000 time interp: T climatology: T weights : pairing : data type: yearly land/sea mask: trc_rst_read : read data in the TOP restart file ~~~~~~~~~~~~ read TRNDIC (rec: 1) in ./restart_trc_0000.nc ok read TRNAlkalini (rec: 1) in ./restart_trc_0000.nc ok read TRNO2 (rec: 1) in ./restart_trc_0000.nc ok read TRNCaCO3 (rec: 1) in ./restart_trc_0000.nc ok read TRNPO4 (rec: 1) in ./restart_trc_0000.nc ok read TRNPOC (rec: 1) in ./restart_trc_0000.nc ok read TRNSi (rec: 1) in ./restart_trc_0000.nc ok read TRNPHY (rec: 1) in ./restart_trc_0000.nc ok read TRNZOO (rec: 1) in ./restart_trc_0000.nc ok read TRNDOC (rec: 1) in ./restart_trc_0000.nc ok read TRNPHY2 (rec: 1) in ./restart_trc_0000.nc ok read TRNZOO2 (rec: 1) in ./restart_trc_0000.nc ok read TRNDSi (rec: 1) in ./restart_trc_0000.nc ok read TRNFer (rec: 1) in ./restart_trc_0000.nc ok read TRNBFe (rec: 1) in ./restart_trc_0000.nc ok read TRNGOC (rec: 1) in ./restart_trc_0000.nc ok read TRNSFe (rec: 1) in ./restart_trc_0000.nc ok read TRNDFe (rec: 1) in ./restart_trc_0000.nc ok read TRNGSi (rec: 1) in ./restart_trc_0000.nc ok read TRNNFe (rec: 1) in ./restart_trc_0000.nc ok read TRNNCHL (rec: 1) in ./restart_trc_0000.nc ok read TRNDCHL (rec: 1) in ./restart_trc_0000.nc ok read TRNNO3 (rec: 1) in ./restart_trc_0000.nc ok read TRNNH4 (rec: 1) in ./restart_trc_0000.nc ok read TRNCFC11 (rec: 1) in ./restart_trc_0000.nc ok read TRNCFC12 (rec: 1) in ./restart_trc_0000.nc ok read TRNSF6 (rec: 1) in ./restart_trc_0000.nc ok read TRNAge (rec: 1) in ./restart_trc_0000.nc ok read TRBDIC (rec: 1) in ./restart_trc_0000.nc ok read TRBAlkalini (rec: 1) in ./restart_trc_0000.nc ok read TRBO2 (rec: 1) in ./restart_trc_0000.nc ok read TRBCaCO3 (rec: 1) in ./restart_trc_0000.nc ok read TRBPO4 (rec: 1) in ./restart_trc_0000.nc ok read TRBPOC (rec: 1) in ./restart_trc_0000.nc ok read TRBSi (rec: 1) in ./restart_trc_0000.nc ok read TRBPHY (rec: 1) in ./restart_trc_0000.nc ok read TRBZOO (rec: 1) in ./restart_trc_0000.nc ok read TRBDOC (rec: 1) in ./restart_trc_0000.nc ok read TRBPHY2 (rec: 1) in ./restart_trc_0000.nc ok read TRBZOO2 (rec: 1) in ./restart_trc_0000.nc ok read TRBDSi (rec: 1) in ./restart_trc_0000.nc ok read TRBFer (rec: 1) in ./restart_trc_0000.nc ok read TRBBFe (rec: 1) in ./restart_trc_0000.nc ok read TRBGOC (rec: 1) in ./restart_trc_0000.nc ok read TRBSFe (rec: 1) in ./restart_trc_0000.nc ok read TRBDFe (rec: 1) in ./restart_trc_0000.nc ok read TRBGSi (rec: 1) in ./restart_trc_0000.nc ok read TRBNFe (rec: 1) in ./restart_trc_0000.nc ok read TRBNCHL (rec: 1) in ./restart_trc_0000.nc ok read TRBDCHL (rec: 1) in ./restart_trc_0000.nc ok read TRBNO3 (rec: 1) in ./restart_trc_0000.nc ok read TRBNH4 (rec: 1) in ./restart_trc_0000.nc ok read TRBCFC11 (rec: 1) in ./restart_trc_0000.nc ok read TRBCFC12 (rec: 1) in ./restart_trc_0000.nc ok read TRBSF6 (rec: 1) in ./restart_trc_0000.nc ok read TRBAge (rec: 1) in ./restart_trc_0000.nc ok *** Total number of passive tracer jptra = 28 *** Total volume of ocean = 1.340033880867921E+018 *** Total inital content of all tracers tracer nb : 1 name : DIC initial content : 0.3121040896E+16 tracer nb : 2 name : Alkalini initial content : 0.3250888442E+16 tracer nb : 3 name : O2 initial content : 0.1945706953E+15 tracer nb : 4 name : CaCO3 initial content : 0.7001900381E+11 tracer nb : 5 name : PO4 initial content : 0.3539484377E+15 tracer nb : 6 name : POC initial content : 0.9525649834E+11 tracer nb : 7 name : Si initial content : 0.1226201308E+15 tracer nb : 8 name : PHY initial content : 0.5398812011E+11 tracer nb : 9 name : ZOO initial content : 0.2154511580E+11 tracer nb : 10 name : DOC initial content : 0.2633905930E+13 tracer nb : 11 name : PHY2 initial content : 0.1339595464E+11 tracer nb : 12 name : ZOO2 initial content : 0.3100611635E+11 tracer nb : 13 name : DSi initial content : 0.5563241368E+10 tracer nb : 14 name : Fer initial content : 0.7587123686E+09 tracer nb : 15 name : BFe initial content : 0.5599704035E+07 tracer nb : 16 name : GOC initial content : 0.1225416940E+12 tracer nb : 17 name : SFe initial content : 0.6204824121E+07 tracer nb : 18 name : DFe initial content : 0.3995349002E+06 tracer nb : 19 name : GSi initial content : 0.4642951910E+11 tracer nb : 20 name : NFe initial content : 0.1645925756E+07 tracer nb : 21 name : NCHL initial content : 0.1303799845E+11 tracer nb : 22 name : DCHL initial content : 0.4228124358E+10 tracer nb : 23 name : NO3 initial content : 0.3156891412E+15 tracer nb : 24 name : NH4 initial content : 0.1041873883E+12 tracer nb : 25 name : CFC11 initial content : 0.1375183195E+10 tracer nb : 26 name : CFC12 initial content : 0.3339042623E+09 tracer nb : 27 name : SF6 initial content : 0.1242997789E+06 tracer nb : 28 name : Age initial content : 0.7768181734E+21 dia_ptr_init : poleward transport and msf initialization ~~~~~~~~~~~~ Namelist namptr : set ptr parameters Poleward heat & salt transport (T) or not (F) ln_diaptr = T Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = T iom_nf90_open ~~~ open existing file: subbasins.nc in READ mode ---> subbasins.nc OK read atlmsk (rec: 1) in subbasins.nc ok read pacmsk (rec: 1) in subbasins.nc ok read indmsk (rec: 1) in subbasins.nc ok read basins (rec: 1) in subbasins.nc ok iom_close ~~~ close file: subbasins.nc ok dia_hsb_init ~~~~~~~~ check the heat and salt budgets (T) or not (F) ln_diahsb = T ~~~~~~~ dia_hsb_rst at it= 9855713 date= 20000401 ~~~~~~~ read ssh_ini (rec: 1) in ./restartopa_0000.nc ok read e3t_ini (rec: 1) in ./restartopa_0000.nc ok read hc_loc_ini (rec: 1) in ./restartopa_0000.nc ok read sc_loc_ini (rec: 1) in ./restartopa_0000.nc ok trd_init : Momentum/Tracers trends ~~~~~~~~~~ Namelist namtrd : set trends parameters global domain averaged dyn & tra trends ln_glo_trd = F U & V trends: 3D output ln_dyn_trd = F U & V trends: Mixed Layer averaged ln_dyn_mxl = F T & S trends: 3D output ln_tra_trd = T T & S trends: Mixed Layer averaged ln_tra_mxl = F Kinetic Energy trends ln_KE_trd = F Potential Energy trends ln_PE_trd = F Barotropic vorticity trends ln_vor_trd = F frequency of trends diagnostics (glo) nn_trd = 365 Euler time step switch is 1 AAAAAAAA