Changeset 561 for CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM
- Timestamp:
- 02/18/09 12:30:42 (15 years ago)
- Location:
- CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM
- Files:
-
- 3 added
- 4 edited
-
Bands_96x95x19_3prc.dat (added)
-
gcm.def_LMD9695 (added)
-
namcouple_ORCA2xLMD144142 (modified) (24 diffs)
-
namcouple_ORCA2xLMD9671 (modified) (19 diffs)
-
namcouple_ORCA2xLMD9695 (added)
-
namelist_ORCA2 (modified) (11 diffs)
-
namelist_ice_ORCA2 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namcouple_ORCA2xLMD144142
r409 r561 29 29 # 30 30 $CHANNEL 31 31 <lib_mpi> 32 32 <nproc_atm> <nproc_atm> 33 33 <nproc_oce> <nproc_oce> … … 59 59 # 60 60 $RUNTIME 61 61 <total_time> 62 62 $END 63 63 ############################################################################### … … 67 67 # 68 68 $INIDATE 69 69 <initial_date> 70 70 $END 71 71 ############################################################################### … … 108 108 # Field 1 : Weighted sea surface temperature (o->a 1) 109 109 # 110 O_SSTSST SISUTESW 1 86400 5 sstoc.nc EXPORTED110 O_SSTSST SISUTESW 1 86400 5 sstoc.nc <output_mode> 111 111 torc tlmd LAG=5760 112 112 P 2 P 0 113 113 LOCTRANS CHECKIN MOZAIC CHECKOUT REVERSE 114 # LOCTRANS: AVERAGE to average value over coupling period 114 115 AVERAGE 115 116 # CHECKIN: indicate computation of global, land and sea field integrals. … … 126 127 # Field 2 : Sea ice extent (o->a 2) 127 128 # 128 OIceFrac SIICECOV 2 86400 5 sstoc.nc EXPORTED129 OIceFrac SIICECOV 2 86400 5 sstoc.nc <output_mode> 129 130 torc tlmd LAG=5760 130 131 P 2 P 0 … … 144 145 # Field 3 : Weighted Sea Ice Temperature (o->a 3) 145 146 # 146 O_TepIce SIICTEMW 34 86400 5 sstoc.nc EXPORTED147 O_TepIce SIICTEMW 34 86400 5 sstoc.nc <output_mode> 147 148 torc tlmd LAG=5760 148 149 P 2 P 0 … … 161 162 # Field 4 : Weighted Sea ice albedo (o->a 4) 162 163 # 163 O_AlbIce SIICEALW 17 86400 5 sstoc.nc EXPORTED164 O_AlbIce SIICEALW 17 86400 5 sstoc.nc <output_mode> 164 165 torc tlmd LAG=5760 165 166 P 2 P 0 … … 182 183 # Field 5 : wind stress along X axis (a->o tau 1) 183 184 # 184 COTAUXXU O_OTaux1 23 86400 6 flxat.nc EXPORTED185 COTAUXXU O_OTaux1 23 86400 6 flxat.nc <output_mode> 185 186 tlmd uorc LAG=1800 186 187 P 0 P 2 … … 202 203 # Field 6 : stress along Y axis (a->o tau 1) 203 204 # 204 COTAUYYU O_OTauy1 23 86400 6 flxat.nc EXPORTED205 COTAUYYU O_OTauy1 23 86400 6 flxat.nc <output_mode> 205 206 tlmd uorc LAG=1800 206 207 P 0 P 2 … … 222 223 # Field 7 :wind stress along Z axis (a->o tau 1) 223 224 # 224 COTAUZZU O_OTauz1 23 86400 6 flxat.nc EXPORTED225 COTAUZZU O_OTauz1 23 86400 6 flxat.nc <output_mode> 225 226 tlmd uorc LAG=1800 226 227 P 0 P 2 … … 242 243 # Field 8 : wind stress along X axis 2 (a->o tau 2) 243 244 # 244 COTAUXXV O_OTaux2 24 86400 6 flxat.nc EXPORTED245 COTAUXXV O_OTaux2 24 86400 6 flxat.nc <output_mode> 245 246 tlmd vorc LAG=1800 246 247 P 0 P 2 … … 262 263 # Field 9 : wind stress along Y axis 2 (a->o tau 2) 263 264 # 264 COTAUYYV O_OTauy2 24 86400 6 flxat.nc EXPORTED265 COTAUYYV O_OTauy2 24 86400 6 flxat.nc <output_mode> 265 266 tlmd vorc LAG=1800 266 267 P 0 P 2 … … 282 283 # Field 10 : wind stress along Z axis 2 (a->o tau 2) 283 284 # 284 COTAUZZV O_OTauz2 24 86400 6 flxat.nc EXPORTED285 COTAUZZV O_OTauz2 24 86400 6 flxat.nc <output_mode> 285 286 tlmd vorc LAG=1800 286 287 P 0 P 2 … … 299 300 INT=1 300 301 ## 301 ##########################################################################302 # Field 11 : 10-m wind speed (a->o tau 2)303 #304 COWINDSP O_Wind10 36 86400 6 flxat.nc EXPORTED305 tlmd torc LAG=1800306 P 0 P 2307 INVERT CHECKIN MASK EXTRAP INTERP CHECKOUT308 # Reverse309 NORSUD WSTEST310 # CHECKIN: indicate computation of global, land and sea field integrals.311 INT=1312 # Mask value313 888888.314 # Extrapolation method315 NINENN 2 1 1316 # Interpolation method317 BILINEAR Z SCALAR318 # CHECKOUT: indicate computation of global, land and sea field integrals.319 INT=1320 #321 302 ######################################################################### 322 # Field 1 2 : tprecip = ( rain_ice + snow_ice ) - evap_ice(a->o flx 7)323 # 324 CO PEFWAT OIPr-Sub 25 86400 5 flxat.nc EXPORTED325 tlmd torc LAG=1800 326 P 0 P 2 327 INVERT CHECKIN MOZAIC BLASNEWCHECKOUT303 # Field 11 : emp_ice = evap_ice - ( rain_ice + snow_ice ) (a->o flx 7) 304 # 305 COEMPSIC OSubMPre 25 86400 4 flxat.nc <output_mode> 306 tlmd torc LAG=1800 307 P 0 P 2 308 INVERT CHECKIN MOZAIC CHECKOUT 328 309 # Reverse 329 310 NORSUD WSTEST … … 332 313 # Interpolation method ou parametres mozaic 333 314 wa2o.flx 92 2 32 334 # change of sign 335 -1.0 0 336 # CHECKOUT: indicate computation of global, land and sea field integrals. 337 INT=1 338 # 339 ########################################################################## 340 # Field 13 : sprecip = snow_ice - evap_ice (o->a flx 6) 341 # 342 COPEFICE OISn-Sub 25 86400 5 flxat.nc EXPORTED 343 tlmd torc LAG=1800 344 P 0 P 2 345 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 315 # CHECKOUT: indicate computation of global, land and sea field integrals. 316 INT=1 317 # 318 ########################################################################## 319 # Field 12 : sprecip = evap_ice - snow_ice (o->a flx 6) 320 # 321 CONESOPR OISubMSn 25 86400 4 flxat.nc <output_mode> 322 tlmd torc LAG=1800 323 P 0 P 2 324 INVERT CHECKIN MOZAIC CHECKOUT 346 325 # Reverse 347 326 NORSUD WSTEST … … 350 329 # Interpolation method ou parametres mozaic 351 330 wa2o.flx 92 2 32 352 # change of sign 353 -1.0 0 354 # CHECKOUT: indicate computation of global, land and sea field integrals. 355 INT=1 356 # 357 ########################################################################## 358 # Field 14: emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 359 # 360 COTOSPSU OOEv-OPr 28 86400 4 flxat.nc EXPORTED 331 # CHECKOUT: indicate computation of global, land and sea field integrals. 332 INT=1 333 # 334 ########################################################################## 335 # Field 13 : emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 336 # 337 COEMPOCE OOEvaMPr 28 86400 4 flxat.nc <output_mode> 361 338 tlmd torc LAG=1800 362 339 P 0 P 2 … … 368 345 # Interpolation method ou parametres mozaic 369 346 wa2o.flx 92 2 32 370 # change of sign 371 ###-1.0 0 372 # CHECKOUT: indicate computation of global, land and sea field integrals. 373 INT=1 347 # CHECKOUT: indicate computation of global, land and sea field integrals. 348 INT=1 349 # 350 ########################################################################## 351 # Field 14 : Evaporation Ice (a->o flx 8) 352 # 353 COICEVAP OIceEvap 27 86400 4 flxat.nc <output_mode> 354 tlmd torc LAG=1800 355 P 0 P 2 356 INVERT CHECKIN MOZAIC CHECKOUT 357 # Reverse 358 NORSUD WSTEST 359 # CHECKIN: indicate computation of global, land and sea field integrals. 360 INT=1 361 # Interpolation method ou parametres mozaic 362 wa2o.flx 92 2 32 363 # CHECKOUT: indicate computation of global, land and sea field integrals. 364 INT=1 365 # 374 366 ########################################################################## 375 367 # Field 15 : solar heat flux on ocean (a->o flx 2) 376 368 # 377 COS WFLDO O_QsrOce 7 86400 4 flxat.nc EXPORTED369 COSHFOCE O_QsrOce 7 86400 4 flxat.nc <output_mode> 378 370 tlmd torc LAG=1800 379 371 P 0 P 2 … … 387 379 # CHECKOUT: indicate computation of global, land and sea field integrals. 388 380 INT=1 381 # 389 382 ######################################################################### 390 383 # Field 16 : Non solar heat flux on ocean (a->o flx 4) 391 384 # 392 CONSF LDO O_QnsOce 6 86400 4 flxat.nc EXPORTED385 CONSFOCE O_QnsOce 6 86400 4 flxat.nc <output_mode> 393 386 tlmd torc LAG=1800 394 387 P 0 P 2 … … 402 395 # CHECKOUT: indicate computation of global, land and sea field integrals. 403 396 INT=1 397 # 404 398 ########################################################################## 405 399 # Field 17 : solar heat flux on ice (a->o flx 1) 406 400 # 407 COSHF LIC O_QsrIce 7 86400 4 flxat.nc EXPORTED401 COSHFICE O_QsrIce 7 86400 4 flxat.nc <output_mode> 408 402 tlmd torc LAG=1800 409 403 P 0 P 2 … … 421 415 # Field 18 : Non solar heat flux on ice (a->o flx 3) 422 416 # 423 CONSF LIC O_QnsIce 6 86400 4 flxat.nc EXPORTED417 CONSFICE O_QnsIce 6 86400 4 flxat.nc <output_mode> 424 418 tlmd torc LAG=1800 425 419 P 0 P 2 … … 437 431 # Field 19 : Non solar heat flux derivative (o->a flx 5) 438 432 # 439 CODFLXDT O_dQnsdT 35 86400 4 flxat.nc EXPORTED433 CODFLXDT O_dQnsdT 35 86400 4 flxat.nc <output_mode> 440 434 tlmd torc LAG=1800 441 435 P 0 P 2 … … 449 443 # CHECKOUT: indicate computation of global, land and sea field integrals. 450 444 INT=1 445 # 451 446 ########################################################################### 452 447 # Field 20 : iceberg calving (a->o flx 12) 453 448 # 454 C RWOCEIS OCalving 32 86400 5 flxat.nc EXPORTED455 tlmd torc LAG=1800 456 P 0 P 2 457 INVERT CHECKIN MOZAIC BLASNEWCHECKOUT449 COCALVIN OCalving 32 86400 4 flxat.nc <output_mode> 450 tlmd torc LAG=1800 451 P 0 P 2 452 INVERT CHECKIN MOZAIC CHECKOUT 458 453 # Reverse 459 454 NORSUD WSTEST … … 462 457 # Interpolation method ou parametres mozaic 463 458 wa2o.cal 98 5 71 464 # change of sign 465 -1.0 0 466 # CHECKOUT: indicate computation of global, land and sea field integrals. 467 INT=1 459 # CHECKOUT: indicate computation of global, land and sea field integrals. 460 INT=1 461 # 468 462 ############################################################################ 469 463 # Field 21 : run-off (a->o flx 11) 470 464 # 471 C RWOCERD O_Runoff 32 86400 5 flxat.nc EXPORTED465 COLIQRUN O_Runoff 32 86400 5 flxat.nc <output_mode> 472 466 tlmd torc LAG=1800 473 467 P 0 P 2 … … 481 475 # change units from m/s to kg/s, change of sign 482 476 # weights convert from kg/s to kg/m^2/s 483 -1000.0 0484 # CHECKOUT: indicate computation of global, land and sea field integrals. 485 INT=1 486 # 487 ############################################################################ 488 # 489 $END 477 1000.0 0 478 # CHECKOUT: indicate computation of global, land and sea field integrals. 479 INT=1 480 # 481 ############################################################################ 482 # 483 $END -
CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namcouple_ORCA2xLMD9671
r396 r561 55 55 $END 56 56 ############################################################################### 57 # 57 # 58 58 # RUNTIME: total simulated time for the actual run in seconds (<I8) 59 59 # … … 64 64 # 65 65 # INIDATE (I8) 66 # initial date of the run (YYYYMMDD) 66 # initial date of the run (YYYYMMDD) 67 67 # 68 68 $INIDATE … … 105 105 # OCEAN --->>> ATMOS 106 106 # -------------------- 107 # Field 1 : Sea surface temperature (o->a 1) 108 O_SSTSST SISUTESW 1 86400 5 sstoc.nc EXPORTED 107 # Field 1 : Weighted Sea surface temperature (o->a 1) 108 # 109 O_SSTSST SISUTESW 1 86400 5 sstoc.nc <output_mode> 109 110 torc tlmd LAG=5760 110 111 P 2 P 0 111 112 LOCTRANS CHECKIN MOZAIC CHECKOUT REVERSE 113 # LOCTRANS: AVERAGE to average value over coupling period 112 114 AVERAGE 113 115 # CHECKIN: indicate computation of global, land and sea field integrals. … … 121 123 NORSUD WSTEST 122 124 # 123 #124 ############################################################################125 #126 125 ############################################################################ 127 126 # Field 2 : Sea ice extent (o->a 2) 128 127 # 129 OIceFrac SIICECOV 2 86400 5 sstoc.nc EXPORTED128 OIceFrac SIICECOV 2 86400 5 sstoc.nc <output_mode> 130 129 torc tlmd LAG=5760 131 130 P 2 P 0 … … 143 142 # 144 143 ############################################################################ 145 # Field 3 : Sea Ice Temperature (o->a 4)146 # 147 O_TepIce SIICTEMW 34 86400 5 sstoc.nc EXPORTED144 # Field 3 : Weighted Sea Ice Temperature (o->a 3) 145 # 146 O_TepIce SIICTEMW 34 86400 5 sstoc.nc <output_mode> 148 147 torc tlmd LAG=5760 149 148 P 2 P 0 … … 160 159 # 161 160 ############################################################################ 162 # Field 4 : Sea ice albedo (o->a 3)163 # 164 O_AlbIce SIICEALW 17 86400 5 sstoc.nc EXPORTED161 # Field 4 : Weighted Sea ice albedo (o->a 4) 162 # 163 O_AlbIce SIICEALW 17 86400 5 sstoc.nc <output_mode> 165 164 torc tlmd LAG=5760 166 165 P 2 P 0 … … 182 181 # ------------------------- 183 182 ############################################################################ 184 # Field 17 : wind stress along X axis (a->o tau 1) 185 # 186 COTAUXXU O_OTaux1 23 86400 6 flxat.nc EXPORTED 183 # 184 # Field 5 : wind stress along X axis (a->o tau 1) 185 # 186 COTAUXXU O_OTaux1 23 86400 6 flxat.nc <output_mode> 187 187 tlmd uorc LAG=1800 188 188 P 0 P 2 … … 203 203 ########################################################################### 204 204 # 205 # 206 # Field 18 : stress along Y axis (a->o tau 1) 207 # 208 COTAUYYU O_OTauy1 23 86400 6 flxat.nc EXPORTED 205 # Field 6 : stress along Y axis (a->o tau 1) 206 # 207 COTAUYYU O_OTauy1 23 86400 6 flxat.nc <output_mode> 209 208 tlmd uorc LAG=1800 210 209 P 0 P 2 … … 225 224 ########################################################################### 226 225 # 227 # 228 # Field 19 :wind stress along Z axis (a->o tau 1) 229 # 230 COTAUZZU O_OTauz1 23 86400 6 flxat.nc EXPORTED 226 # Field 7 :wind stress along Z axis (a->o tau 1) 227 # 228 COTAUZZU O_OTauz1 23 86400 6 flxat.nc <output_mode> 231 229 tlmd uorc LAG=1800 232 230 P 0 P 2 … … 247 245 ########################################################################## 248 246 # 249 # Field 20: wind stress along X axis 2 (a->o tau 2)250 # 251 COTAUXXV O_OTaux2 24 86400 6 flxat.nc EXPORTED247 # Field 8 : wind stress along X axis 2 (a->o tau 2) 248 # 249 COTAUXXV O_OTaux2 24 86400 6 flxat.nc <output_mode> 252 250 tlmd vorc LAG=1800 253 251 P 0 P 2 … … 268 266 ########################################################################## 269 267 # 270 # Field 21: wind stress along Y axis 2 (a->o tau 2)271 # 272 COTAUYYV O_OTauy2 24 86400 6 flxat.nc EXPORTED268 # Field 9 : wind stress along Y axis 2 (a->o tau 2) 269 # 270 COTAUYYV O_OTauy2 24 86400 6 flxat.nc <output_mode> 273 271 tlmd vorc LAG=1800 274 272 P 0 P 2 … … 289 287 ########################################################################## 290 288 # 291 # Field 22: wind stress along Z axis 2 (a->o tau 2)292 # 293 COTAUZZV O_OTauz2 24 86400 6 flxat.nc EXPORTED289 # Field 10 : wind stress along Z axis 2 (a->o tau 2) 290 # 291 COTAUZZV O_OTauz2 24 86400 6 flxat.nc <output_mode> 294 292 tlmd vorc LAG=1800 295 293 P 0 P 2 … … 307 305 # CHECKOUT: indicate computation of global, land and sea field integrals. 308 306 INT=1 309 ## 307 # 308 ######################################################################### 309 # 310 # Field 11 : emp_ice = evap_ice - ( rain_ice + snow_ice ) (a->o flx 7) 311 # 312 COEMPSIC OSubMPre 25 86400 4 flxat.nc <output_mode> 313 tlmd torc LAG=1800 314 P 0 P 2 315 INVERT CHECKIN MOZAIC CHECKOUT 316 # Reverse 317 NORSUD WSTEST 318 # CHECKIN: indicate computation of global, land and sea field integrals. 319 INT=1 320 # Interpolation method ou parametres mozaic 321 wa2o.flx 92 2 14 322 # CHECKOUT: indicate computation of global, land and sea field integrals. 323 INT=1 324 # 310 325 ########################################################################## 311 326 # 312 # Field 11 : 10-m wind speed (a->o tau 2) 313 # 314 COWINDSP O_Wind10 36 86400 6 flxat.nc EXPORTED 315 tlmd torc LAG=1800 316 P 0 P 2 317 INVERT CHECKIN MASK EXTRAP INTERP CHECKOUT 318 # Reverse 319 NORSUD WSTEST 320 # CHECKIN: indicate computation of global, land and sea field integrals. 321 INT=1 322 # Mask value 323 888888. 324 # Extrapolation method 325 NINENN 2 1 1 326 # Interpolation method 327 BILINEAR Z SCALAR 327 # Field 12 : sprecip = evap_ice - snow_ice (o->a flx 6) 328 # 329 CONESOPR OISubMSn 25 86400 4 flxat.nc <output_mode> 330 tlmd torc LAG=1800 331 P 0 P 2 332 INVERT CHECKIN MOZAIC CHECKOUT 333 # Reverse 334 NORSUD WSTEST 335 # CHECKIN: indicate computation of global, land and sea field integrals. 336 INT=1 337 # Interpolation method ou parametres mozaic 338 wa2o.flx 92 2 14 328 339 # CHECKOUT: indicate computation of global, land and sea field integrals. 329 340 INT=1 … … 331 342 ######################################################################### 332 343 # 333 # Field 10 : tprecip = ( rain_ice + snow_ice ) - evap_ice 334 ### (a->o flx 7) 335 # 336 COPEFWAT OIPr-Sub 25 86400 5 flxat.nc EXPORTED 344 # Field 13 : emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 345 # 346 COEMPOCE OOEvaMPr 29 86400 4 flxat.nc <output_mode> 337 347 tlmd torc LAG=1800 338 348 P 0 P 2 339 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 340 # Reverse 341 NORSUD WSTEST 342 # CHECKIN: indicate computation of global, land and sea field integrals. 343 INT=1 344 # Interpolation method ou parametres mozaic 345 wa2o.flx 92 2 14 346 # change of sign 347 -1.0 0 349 INVERT CHECKIN MOZAIC CHECKOUT 350 # Reverse 351 NORSUD WSTEST 352 # CHECKIN: indicate computation of global, land and sea field integrals. 353 INT=1 354 # Interpolation method ou parametres mozaic 355 wa2o.flx 92 2 14 348 356 # CHECKOUT: indicate computation of global, land and sea field integrals. 349 357 INT=1 … … 351 359 ########################################################################## 352 360 # 353 # Field 1 1 : sprecip = snow_ice - evap_ice (o->a flx 6)354 # 355 CO PEFICE OISn-Sub 25 86400 5 flxat.nc EXPORTED361 # Field 14 : Evaporation Ice (a->o flx 8) 362 # 363 COICEVAP OIceEvap 41 86400 4 flxat.nc <output_mode> 356 364 tlmd torc LAG=1800 357 365 P 0 P 2 358 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 359 # Reverse 360 NORSUD WSTEST 361 # CHECKIN: indicate computation of global, land and sea field integrals. 362 INT=1 363 # Interpolation method ou parametres mozaic 364 wa2o.flx 92 2 14 365 # change of sign 366 -1.0 0 367 # CHECKOUT: indicate computation of global, land and sea field integrals. 368 INT=1 366 INVERT CHECKIN MOZAIC CHECKOUT 367 # Reverse 368 NORSUD WSTEST 369 # CHECKIN: indicate computation of global, land and sea field integrals. 370 INT=1 371 # Interpolation method ou parametres mozaic 372 wa2o.flx 92 2 14 373 # CHECKOUT: indicate computation of global, land and sea field integrals. 374 INT=1 375 # 376 ########################################################################### 377 # 378 # Field 15 : solar heat flux on ocean (a->o flx 2) 379 # 380 COSHFOCE O_QsrOce 7 86400 4 flxat.nc <output_mode> 381 tlmd torc LAG=1800 382 P 0 P 2 383 INVERT CHECKIN MOZAIC CHECKOUT 384 # Reverse 385 NORSUD WSTEST 386 # CHECKIN: indicate computation of global, land and sea field integrals. 387 INT=1 388 # Interpolation method or mozaic parameters 389 wa2o.flx 92 2 14 390 # CHECKOUT: indicate computation of global, land and sea field integrals. 391 INT=1 392 # 369 393 ######################################################################### 370 394 # 371 # Field 12: emp = emp_oce = evap_oce - 372 ### ( rain_oce + snow_oce ) (a->o flx 9) 373 # 374 COTOSPSU OOEv-OPr 28 86400 4 flxat.nc EXPORTED 375 tlmd torc LAG=1800 376 P 0 P 2 377 INVERT CHECKIN MOZAIC CHECKOUT 378 # Reverse 379 NORSUD WSTEST 380 # CHECKIN: indicate computation of global, land and sea field integrals. 381 INT=1 382 # Interpolation method ou parametres mozaic 383 wa2o.flx 92 2 14 384 # change of sign 385 ###-1.0 0 386 # CHECKOUT: indicate computation of global, land and sea field integrals. 387 INT=1 388 ########################################################################## 389 # 390 # Field 13 : Evaporation Ice (a->o flx 8) 391 # 392 #####COICEVAP OIceEvap 27 86400 4 flxat.nc EXPORTED 393 #####tlmd torc LAG=1800 394 #####P 0 P 2 395 ###INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 396 #####INVERT CHECKIN MOZAIC CHECKOUT 397 # Reverse 398 #####NORSUD WSTEST 399 # CHECKIN: indicate computation of global, land and sea field integrals. 400 ##### INT=1 401 # Interpolation method ou parametres mozaic 402 #####wa2o.flx 92 2 14 403 # change of sign 404 ###-1.0 0 405 # CHECKOUT: indicate computation of global, land and sea field integrals. 406 ##### INT=1 395 # Field 16 : Non solar heat flux on ocean (a->o flx 4) 396 # 397 CONSFOCE O_QnsOce 6 86400 4 flxat.nc <output_mode> 398 tlmd torc LAG=1800 399 P 0 P 2 400 INVERT CHECKIN MOZAIC CHECKOUT 401 # Reverse 402 NORSUD WSTEST 403 # CHECKIN: indicate computation of global, land and sea field integrals. 404 INT=1 405 # Interpolation method ou parametres mozaic 406 wa2o.flx 92 2 14 407 # CHECKOUT: indicate computation of global, land and sea field integrals. 408 INT=1 407 409 # 408 410 ########################################################################### 409 # Field 6 : solar heat flux on ocean (a->o flx 2) 410 # 411 COSWFLDO O_QsrOce 7 86400 4 flxat.nc EXPORTED 412 tlmd torc LAG=1800 413 P 0 P 2 414 INVERT CHECKIN MOZAIC CHECKOUT 415 # Reverse 416 NORSUD WSTEST 417 # CHECKIN: indicate computation of global, land and sea field integrals. 418 INT=1 419 # Interpolation method or mozaic parameters 420 wa2o.flx 92 2 14 421 # CHECKOUT: indicate computation of global, land and sea field integrals. 422 INT=1 423 ######################################################################### 424 # 425 # Field 8 : Non solar heat flux on ocean (a->o flx 4) 426 # 427 CONSFLDO O_QnsOce 6 86400 4 flxat.nc EXPORTED 428 tlmd torc LAG=1800 429 P 0 P 2 430 INVERT CHECKIN MOZAIC CHECKOUT 431 # Reverse 432 NORSUD WSTEST 433 # CHECKIN: indicate computation of global, land and sea field integrals. 434 INT=1 435 # Interpolation method ou parametres mozaic 436 wa2o.flx 92 2 14 437 # CHECKOUT: indicate computation of global, land and sea field integrals. 438 INT=1 439 # 440 ########################################################################### 441 # Field 5 : solar heat flux on ice (a->o flx 1) 442 # 443 COSHFLIC O_QsrIce 7 86400 4 flxat.nc EXPORTED 411 # 412 # Field 17 : solar heat flux on ice (a->o flx 1) 413 # 414 COSHFICE O_QsrIce 7 86400 4 flxat.nc <output_mode> 444 415 tlmd torc LAG=1800 445 416 P 0 P 2 … … 454 425 INT=1 455 426 # 456 #457 427 ############################################################################ 458 428 # 459 # Field 7: Non solar heat flux on ice (a->o flx 3)460 # 461 CONSF LIC O_QnsIce 6 86400 4 flxat.nc EXPORTED429 # Field 18 : Non solar heat flux on ice (a->o flx 3) 430 # 431 CONSFICE O_QnsIce 6 86400 4 flxat.nc <output_mode> 462 432 tlmd torc LAG=1800 463 433 P 0 P 2 … … 474 444 ############################################################################ 475 445 # 476 # Field 9 : Non solar heat flux derivative (o->a flx 5)477 # 478 CODFLXDT O_dQnsdT 35 86400 4 flxat.nc EXPORTED446 # Field 19 : Non solar heat flux derivative (o->a flx 5) 447 # 448 CODFLXDT O_dQnsdT 35 86400 4 flxat.nc <output_mode> 479 449 tlmd torc LAG=1800 480 450 P 0 P 2 … … 491 461 ######################################################################### 492 462 # 493 # Field 16: iceberg calving (a->o flx 12)494 # 495 C RWOCEIS OCalving 32 86400 5 flxat.nc EXPORTED463 # Field 20 : iceberg calving (a->o flx 12) 464 # 465 COCALVIN OCalving 36 86400 4 flxat.nc <output_mode> 496 466 tlmd torc LAG=1800 497 467 P 0 P 2 468 INVERT CHECKIN MOZAIC CHECKOUT 469 # Reverse 470 NORSUD WSTEST 471 # CHECKIN: indicate computation of global, land and sea field integrals. 472 INT=1 473 # Interpolation method ou parametres mozaic 474 wa2o.cal 98 5 36 475 # CHECKOUT: indicate computation of global, land and sea field integrals. 476 INT=1 477 # 478 ######################################################################### 479 # 480 # Field 21 : liquid run-off (river + direct) (a->o flx 11) 481 # 482 COLIQRUN O_Runoff 32 86400 5 flxat.nc <output_mode> 483 tlmd torc LAG=1800 484 P 0 P 2 498 485 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 499 486 # Reverse … … 502 489 INT=1 503 490 # Interpolation method ou parametres mozaic 504 wa2o.cal 98 5 36505 # change of sign506 -1.0 0507 # CHECKOUT: indicate computation of global, land and sea field integrals.508 INT=1509 #########################################################################510 #511 # Field 15 : river run-off (a->o flx 11)512 #513 CRWOCERD O_Runoff 32 86400 5 flxat.nc EXPORTED514 tlmd torc LAG=1800515 P 0 P 2516 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT517 # Reverse518 NORSUD WSTEST519 # CHECKIN: indicate computation of global, land and sea field integrals.520 INT=1521 # Interpolation method ou parametres mozaic522 491 wa2o.run 94 4 158 523 # change units from m/s to kg/s , change of sign492 # change units from m/s to kg/s 524 493 # weights convert from kg/s to kg/m^2/s 525 -1000.0 0 526 # CHECKOUT: indicate computation of global, land and sea field integrals. 527 INT=1 528 # 529 ######################################################################### 530 # 531 # Field 14 : direct run-off (a->o flx 10) 532 # 533 ##CRWOCECD FRWOCECD 32 86400 5 flxat.nc EXPORTED 534 ##tlmd torc LAG=1800 535 ##P 0 P 2 536 ##INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 537 # Reverse 538 ##NORSUD WSTEST 539 # CHECKIN: indicate computation of global, land and sea field integrals. 540 ## INT=1 541 # Interpolation method ou parametres mozaic 542 ##wa2o.run 94 4 158 543 # change units from m/s to kg/s, change of sign 544 # weights convert from kg/s to kg/m^2/s 545 ##-1000.0 0 546 # CHECKOUT: indicate computation of global, land and sea field integrals. 547 ## INT=1 548 # 549 ########################################################################## 494 1000.0 0 495 # CHECKOUT: indicate computation of global, land and sea field integrals. 496 INT=1 497 # 550 498 ########################################################################### 551 499 # -
CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namelist_ORCA2
r396 r561 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun , nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl)2 !! NEMO/OPA : 1 - run manager (namrun) 3 3 !! namelists 2 - Domain (nam_zgr, nam_zgr_sco, namdom) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc )6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, namtide) 7 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 8 !! 6 - Tracer (nameos, nam_traadv, nam_traldf, namtdp) … … 10 10 !! 8 - Verical physics (namzdf, namnpc, namric, namtke, namkpp, namddm) 11 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol )12 !! 9 - miscellaneous (namsol, nam_mpp, nam_mpp_dyndist, namctl) 13 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical definition 15 ! ******* use .true. or .false. and NOT .TRUE. or .FALSE. 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. 16 15 17 16 !!====================================================================== … … 19 18 !!====================================================================== 20 19 !! namrun parameters of the run 21 !! nam_mpp Massively Parallel Processing 22 !! nam_mpp_dyndist 23 !! nam_ctl Control prints & Benchmark 24 !!====================================================================== 25 26 !----------------------------------------------------------------------- 27 ! namrun parameters of the run 28 !----------------------------------------------------------------------- 29 ! no job number 30 ! cexper experience name for vairmer format 31 ! ln_rstart boolean term for restart (true or false) 32 ! nrstdt restart control = 0 restart, do not control nit000 in the restart file. 33 ! ! = 1 restart, control nit000 in the restart file. Do not 34 ! ! use the date in the restart file (use ndate0 in namelist) 35 ! ! = 2 restart, control nit000 in the restart file, use the date 36 ! ! in the restart file. ndate0 in the namelist is ignored. 37 ! nit000 number of the first time step 38 ! nitend number of the last time step 39 ! ndate0 initial calendar date aammjj 40 ! nleapy Leap year calendar (0/1) 41 ! ninist initial state output flag (0/1) 42 ! nstock frequency of restart file 43 ! nwrite frequency of OUTPUT file 44 ! ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 45 &namrun 46 no = 0 47 cexper = "ORCA2" 48 ln_rstart = .false. 49 nrstdt = 1 50 nit000 = 76 51 nitend = 150 52 ndate0 =18600106 53 nleapy = 30 54 ninist = 0 55 nstock = 5400 56 nwrite = 1 57 ln_dimgnnn = .false. 58 / 59 !----------------------------------------------------------------------- 60 ! nam_mpp Massively Parallel Processing 61 !----------------------------------------------------------------------- 62 ! c_mpi_send mpi send/recieve type 63 ! = 'S' : standard blocking send 64 ! = 'B' : buffer blocking send 65 ! = 'I' : immediate non-blocking send 66 &nam_mpp 67 c_mpi_send = 'S' 68 / 69 !----------------------------------------------------------------------- 70 ! nam_mpp_dyndist Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 71 !----------------------------------------------------------------------- 72 ! jpni number of processors following i 73 ! jpnj number of processors following j 74 ! jpnij number of local domains 75 &nam_mpp_dyndist 76 jpni = 1 77 jpnj = 1 78 jpnij = 1 79 / 80 !----------------------------------------------------------------------- 81 ! nam_ctl Control prints & Benchmark 82 !----------------------------------------------------------------------- 83 ! ln_ctl trends control print (expensive!) 84 ! nprint level of print (0 no print) 85 ! nictls start i indice to make the control SUM (very usefull to compare mono- 86 ! nictle end i indice to make the control SUM (-versus multi processor runs) 87 ! njctls start j indice to make the control SUM (very usefull to compare mono- 88 ! njctle end j indice to make the control SUM (-versus multi processor runs) 89 ! nisplt number of processors following i 90 ! njsplt number of processors following j 91 ! nbench Bench parameter (0/1): CAUTION it must be zero except for bench 92 ! for which we don't care about physical meaning of the results 93 ! nbit_cmp bit comparison mode parameter (0/1): enables bit comparison between 94 ! single and multiple processor runs. 95 &namctl 96 ln_ctl = .false. 97 nprint = 0 98 nictls = 0 99 nictle = 0 100 njctls = 0 101 njctle = 0 102 isplt = 1 103 jsplt = 1 104 nbench = 0 105 nbit_cmp = 0 106 / 107 20 !!====================================================================== 21 22 !----------------------------------------------------------------------- 23 &namrun ! parameters of the run 24 !----------------------------------------------------------------------- 25 no = 0 ! job number 26 cexper = "ORCA2" ! experience name 27 cn_ocerst_in = "restartopa"! suffix of ocean restart name (input) 28 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 29 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 30 nrstdt = 1 ! restart control = 0 nit000 is not compared to the restart file value 31 ! = 1 use ndate0 in namelist (not the value in the restart file) 32 ! = 2 calendar parameters read in the restart file 33 nit000 = 76 ! first time step 34 nitend = 150 ! last time step 35 ndate0 =18600106 ! initial calendar date yymmdd (used if nrstdt=1) 36 nleapy = 30 ! Leap year calendar (1) or not (0) 37 ninist = 0 ! output the initial state (1) or not (0) 38 nstock = 5400 ! frequency of creation of a restart file 39 nwrite = 1 ! frequency of write in the output file 40 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 41 / 108 42 !!====================================================================== 109 43 !! *** Domain namelists *** … … 115 49 116 50 !----------------------------------------------------------------------- 117 ! nam_zgr vertical coordinate 118 !----------------------------------------------------------------------- 119 ! ln_zco z-coordinate - full steps (T/F) ("key_zco" may also be defined) 120 ! ln_zps z-coordinate - partial steps (T/F) 121 ! ln_sco s- or hybrid z-s-coordinate (T/F) 122 &nam_zgr 123 ln_zco = .false. 124 ln_zps = .true. 125 ln_sco = .false. 126 / 127 !----------------------------------------------------------------------- 128 ! nam_zgr_sco s-coordinate or hybrid z-s-coordinate 129 !----------------------------------------------------------------------- 130 ! sbot_min minimum depth of s-bottom surface (>0) (m) 131 ! sbot_max maximum depth of s-bottom surface (= ocean depth) (>0) (m) 132 ! theta surface control parameter (0<=theta<=20) 133 ! thetb bottom control parameter (0<=thetb<= 1) 134 ! r_max maximum cut-off r-value allowed (0<r_max<1) 135 &nam_zgr_sco 136 sbot_min = 300. 137 sbot_max = 5250. 138 theta = 6.0 139 thetb = 0.75 140 r_max = 0.15 141 / 142 !----------------------------------------------------------------------- 143 ! namdom space and time domain (bathymetry, mesh, timestep) 144 !----------------------------------------------------------------------- 145 ! ntopo = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 146 ! e3zps_min the thickness of the partial step is set larger than the 147 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 148 ! nmsh =1 create a mesh file (coordinates, scale factors, masks) 149 ! nacc acceleration of convergence method = 0, not used, rdt = rdttra 150 ! = 1, used, rdt < rdttra(k) 151 ! atfp asselin time filter parameter 152 ! rdt time step for the dynamics (and tracer if nacc=0) 153 ! rdtmin minimum time step on tracers 154 ! rdtmax maximum time step on tracers 155 ! rdth depth variation of tracer time step 156 ! rdtbt barotropic time step (for the time splitting algorithm) 157 ! nclosea =0 remove the closed sea from the global domain (orca configuration) 158 ! =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 159 &namdom 160 ntopo = 1 161 e3zps_min = 5. 162 e3zps_rat = 0.1 163 nmsh = 0 164 nacc = 0 165 atfp = 0.1 166 rdt = 5760. 167 rdtmin = 5760. 168 rdtmax = 5760. 169 rdth = 800. 170 rdtbt = 90. 171 nclosea = 0 172 / 173 51 &nam_zgr ! vertical coordinate 52 !----------------------------------------------------------------------- 53 ln_zco = .true. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) 54 ln_zps = .false. ! z-coordinate - partial steps (T/F) 55 ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) 56 / 57 !----------------------------------------------------------------------- 58 &nam_zgr_sco ! s-coordinate or hybrid z-s-coordinate 59 !----------------------------------------------------------------------- 60 sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) 61 sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 62 theta = 6.0 ! surface control parameter (0<=theta<=20) 63 thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) 64 r_max = 0.15 ! maximum cut-off r-value allowed (0<r_max<1) 65 / 66 !----------------------------------------------------------------------- 67 &namdom ! space and time domain (bathymetry, mesh, timestep) 68 !----------------------------------------------------------------------- 69 ntopo = 1 ! compute (=0) or read(=1) the bathymetry file 70 e3zps_min = 5. ! the thickness of the partial step is set larger than the minimum 71 e3zps_rat = 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 72 nmsh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 73 nacc = 0 ! =1 acceleration of convergence method used, rdt < rdttra(k) 74 ! =0, no acceleration, rdt = rdttra 75 atfp = 0.1 ! asselin time filter parameter 76 rdt = 5760. ! time step for the dynamics (and tracer if nacc=0) 77 rdtmin = 5760. ! minimum time step on tracers (used if nacc=1) 78 rdtmax = 5760. ! maximum time step on tracers (used if nacc=1) 79 rdth = 800. ! depth variation of tracer time step (used if nacc=1) 80 nn_baro = 64 ! number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 81 nclosea = 1 ! = 0 no closed sea in the model domain 82 ! = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 83 / 174 84 !!====================================================================== 175 85 !! *** Surface Boundary Condition namelists *** … … 188 98 189 99 !----------------------------------------------------------------------- 190 ! namsbc Surface Boundary Condition (surface module) 191 !----------------------------------------------------------------------- 192 ! nn_fsbc frequency of surface boundary condition computation 193 ! (= the frequency of sea-ice model call) 194 ! ln_ana analytical formulation (fill namsbc_ana ) 195 ! ln_flx flux formulation (fill namsbc_flx ) 196 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) 197 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) 198 ! ln_cpl Coupled formulation (fill namsbc_cpl ) 199 ! nn_ice =0 no ice boundary condition , 200 ! =1 observed ice-cover , 201 ! =2 ice-model used ("key_lim3" or "key_lim2) 202 ! nn_ico_cpl ice-ocean coupling : =0 LIM-3 old case 203 ! =1 stresses computed using now ocean velocity 204 ! =2 combination of 0 and 1 cases 205 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 206 ! ln_ssr Sea Surface Restoring on T and/or S (fill namsbc_ssr) 207 ! nn_fwb FreshWater Budget: =0 no control , 208 ! =1 annual global mean of e-p-r set to zero , 209 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 210 &namsbc 211 nn_fsbc = 5 212 ln_ana = .false. 213 ln_flx = .false. 214 ln_blk_clio = .false. 215 ln_blk_core = .false. 216 ln_cpl = .true. 217 nn_ice = 2 218 nn_ico_cpl = 0 219 ln_dm2dc = .false. 220 ln_rnf = .false. 221 ln_ssr = .false. 222 nn_fwb = 0 223 / 224 !----------------------------------------------------------------------- 225 ! namsbc_ana analytical surface boundary condition 226 !----------------------------------------------------------------------- 227 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 228 ! rn_utau0 uniform value for the i-stress 229 ! rn_vtau0 uniform value for the j-stress 230 ! rn_q0 uniform value for the total heat flux 231 ! rn_qsr0 uniform value for the solar radiation 232 ! rn_emp0 uniform value for the freswater budget (E-P) 233 &namsbc_ana 234 nn_tau000 = 0 235 rn_utau0 = 0.5 236 rn_vtau0 = 0.e0 237 rn_q0 = 0.e0 238 rn_qsr0 = 0.e0 239 rn_emp0 = 0.e0 240 / 241 !----------------------------------------------------------------------- 242 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 243 !----------------------------------------------------------------------- 244 ! cn_dir directory in which the model is executed 100 &namsbc ! Surface Boundary Condition (surface module) 101 !----------------------------------------------------------------------- 102 nn_fsbc = 5 ! frequency of surface boundary condition computation 103 ! (= the frequency of sea-ice model call) 104 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 105 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 106 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 107 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 108 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl ) 109 nn_ice = 2 ! =0 no ice boundary condition , 110 ! =1 use observed ice-cover , 111 ! =2 ice-model used ("key_lim3" or "key_lim2) 112 nn_ico_cpl = 0 ! ice-ocean coupling : =0 each nn_fsbc 113 ! =1 stresses recomputed each ocean time step ("key_lim3" only) 114 ! =2 combination of 0 and 1 cases ("key_lim3" only) 115 ln_dm2dc = .false. ! daily mean to diurnal cycle short wave (qsr) 116 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 117 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 118 nn_fwb = 0 ! FreshWater Budget: =0 unchecked , 119 ! =1 global mean of e-p-r set to zero at each nn_fsbc time step , 120 ! =2 annual global mean of e-p-r set to zero 121 / 122 !----------------------------------------------------------------------- 123 &namsbc_ana ! analytical surface boundary condition 124 !----------------------------------------------------------------------- 125 nn_tau000 = 0 ! gently increase the stress over the first ntau_rst time-steps 126 rn_utau0 = 0.5 ! uniform value for the i-stress 127 rn_vtau0 = 0.e0 ! uniform value for the j-stress 128 rn_q0 = 0.e0 ! uniform value for the total heat flux 129 rn_qsr0 = 0.e0 ! uniform value for the solar radiation 130 rn_emp0 = 0.e0 ! uniform value for the freswater budget (E-P) 131 / 132 !----------------------------------------------------------------------- 133 &namsbc_flx ! surface boundary condition : flux formulation 134 !----------------------------------------------------------------------- 135 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 136 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 137 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. , 'yearly' 138 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. , 'yearly' 139 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. , 'yearly' 140 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. , 'yearly' 141 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. , 'yearly' 245 142 ! 246 ! THE ORDER OF THE FILES MATTER: 247 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 248 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 249 &namsbc_flx 250 cn_dir = './' 251 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 252 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 253 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 254 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 255 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 143 cn_dir = './' ! root directory for the location of the flux files 256 144 / 257 145 !----------------------------------------------------------------------- 258 ! namsbc_clio CLIO bulk formulea 259 !----------------------------------------------------------------------- 260 ! cn_dir directory in which the model is executed 146 &namsbc_clio ! namsbc_clio CLIO bulk formulea 147 !----------------------------------------------------------------------- 148 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 149 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 150 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. , 'yearly' 151 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. , 'yearly' 152 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. , 'yearly' 153 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. , 'yearly' 154 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. , 'yearly' 155 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. , 'yearly' 156 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. , 'yearly' 261 157 ! 262 ! THE ORDER OF THE FILES MATTER: 263 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 264 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 265 &namsbc_clio 266 cn_dir = './' 267 sn_utau = 'taux_1m' , -12. , 'sozotaux', .false. , 1 , 0 268 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .false. , 1 , 0 269 sn_wndm = 'flx' , -12. , 'socliowi', .false. , 1 , 0 270 sn_tair = 'flx' , -12. , 'socliot1', .false. , 1 , 0 271 sn_humi = 'flx' , -12. , 'socliohu', .false. , 1 , 0 272 sn_ccov = 'flx' , -12. , 'socliocl', .false. , 1 , 0 273 sn_prec = 'flx' , -12. , 'socliopl', .false. , 1 , 0 274 / 275 !----------------------------------------------------------------------- 276 ! namsbc_core CORE bulk formulea 277 !----------------------------------------------------------------------- 278 ! cn_dir directory in which the model is executed 279 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 280 ! alpha_precip multiplication factor for precipitation (total & snow) 158 cn_dir = './' ! root directory for the location of the bulk files are 159 / 160 !----------------------------------------------------------------------- 161 &namsbc_core ! namsbc_core CORE bulk formulea 162 !----------------------------------------------------------------------- 163 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 164 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 165 sn_wndi = 'u10' , 24. , 'U_10_MOD' , .false. , .true. , 'yearly' 166 sn_wndj = 'v10' , 24. , 'V_10_MOD' , .false. , .true. , 'yearly' 167 sn_qsr = 'rad' , 24. , 'SWDN_MOD' , .false. , .true. , 'yearly' 168 sn_qlw = 'rad' , 24. , 'LWDN_MOD' , .false. , .true. , 'yearly' 169 sn_tair = 't10' , 24. , 'T_10_MOD' , .false. , .true. , 'yearly' 170 sn_humi = 'q10' , 24. , 'Q_10_MOD' , .false. , .true. , 'yearly' 171 sn_prec = 'precip' , -1. , 'RAIN' , .false. , .true. , 'yearly' 172 sn_snow = 'precip' , -1. , 'SNOW' , .false. , .true. , 'yearly' 281 173 ! 282 ! THE ORDER OF THE FILES MATTER: 283 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 284 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 285 &namsbc_core 286 cn_dir = './' 287 ln_2m = .false. 288 alpha_precip = 1. 289 sn_wndi = 'u10' , 24. , 'U_10_MOD', .false. , 1 , 0 290 sn_wndj = 'v10' , 24. , 'V_10_MOD', .false. , 1 , 0 291 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .false. , 1 , 0 292 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .false. , 1 , 0 293 sn_tair = 't10' , 24. , 'T_10_MOD', .false. , 1 , 0 294 sn_humi = 'q10' , 24. , 'Q_10_MOD', .false. , 1 , 0 295 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 296 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 297 / 298 !----------------------------------------------------------------------- 299 ! namsbc_cpl coupled ocean/atmosphere model ("key_coupled") 300 !----------------------------------------------------------------------- 301 &namsbc_cpl 174 cn_dir = './' ! root directory for the location of the bulk files 175 ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 176 alpha_precip= 1. ! multiplicative factor for precipitation (total & snow) 177 / 178 !----------------------------------------------------------------------- 179 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 180 !----------------------------------------------------------------------- 302 181 ! SEND 303 182 cn_snd_temperature= 'weighted oce and ice' ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 304 183 cn_snd_albedo = 'weighted ice' ! 'none' 'weighted ice' 'mixed oce-ice' 305 184 cn_snd_thickness = 'none' ! 'none' 'weighted ice and snow' 306 cn_snd_c urrent_1= 'none' ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice'307 cn_snd_c urrent_2= 'spherical' ! 'spherical' 'cartesian'308 cn_snd_c urrent_3= 'eastward-northward' ! 'eastward-northward' or 'local grid'309 cn_snd_c urrent_4 = 'T' ! 'T''U,V'185 cn_snd_crt_nature = 'none' ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice' 186 cn_snd_crt_refere = 'spherical' ! 'spherical' 'cartesian' 187 cn_snd_crt_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 188 cn_snd_crt_grid = 'T' ! 'T' 310 189 ! RECEIVE 311 cn_rcv_w10m = ' coupled'! 'none' 'coupled'312 cn_rcv_ stress_1 = 'oce only'! 'oce only' 'oce and ice' 'mixed oce-ice'313 cn_rcv_ stress_2= 'cartesian' ! 'spherical' 'cartesian'314 cn_rcv_ stress_3= 'eastward-northward' ! 'eastward-northward' or 'local grid'315 cn_rcv_ stress_4= 'U,V' ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V'190 cn_rcv_w10m = 'none' ! 'none' 'coupled' 191 cn_rcv_tau_nature = 'mixed oce-ice' ! 'oce only' 'oce and ice' 'mixed oce-ice' 192 cn_rcv_tau_refere = 'cartesian' ! 'spherical' 'cartesian' 193 cn_rcv_tau_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 194 cn_rcv_tau_grid = 'U,V' ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V' 316 195 cn_rcv_dqnsdt = 'coupled' ! 'none' 'coupled' 317 196 cn_rcv_qsr = 'oce and ice' ! 'conservative' 'oce and ice' 'mixed oce-ice' 318 197 cn_rcv_qns = 'oce and ice' ! 'conservative' 'oce and ice' 'mixed oce-ice' 319 cn_rcv_emp = 'oce and ice' ! 'conservative' 'oce and ice' 'mixed oce-ice'198 cn_rcv_emp = 'oce and ice' ! 'conservative' 'oce and ice' 'mixed oce-ice' 320 199 cn_rcv_rnf = 'coupled' ! 'coupled' 'climato' 'mixed' 321 200 cn_rcv_cal = 'coupled' ! 'none' 'coupled' 322 201 / 323 202 !----------------------------------------------------------------------- 324 ! namqsr penetrative solar radiation 325 !----------------------------------------------------------------------- 326 ! ln_traqsr : penetrative solar radiation (T) or not (F) (Default=T) 327 ! rabs fraction of qsr associated with xsi1 328 ! xsi1 first depth of extinction 329 ! xsi2 second depth of extinction 330 &namqsr 331 ln_traqsr = .true. 332 rabs = 0.58 333 xsi1 = 0.35 334 xsi2 = 23.0 335 / 336 !----------------------------------------------------------------------- 337 ! namsbc_rnf runoffs namelist surface boundary condition 338 !----------------------------------------------------------------------- 339 !cn_dir Root directory for location of ssr files 340 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 341 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 342 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 343 !sn_rnf informations about the runoff file to be read 344 !sn_cnf informations about the runoff mouth file to be read 345 ! THE ORDER OF THE FILES MATTER: 346 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 347 ! ! ! (if <0 months) ! name ! (logical) ! 348 &namsbc_rnf 349 cn_dir = './' 350 nn_runoff = 2 351 rn_hrnf = 0.e0 352 rn_avt_rnf = 1.e-3 353 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .true. , 1 , 0 354 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .false. , 1 , 0 355 / 356 !----------------------------------------------------------------------- 357 ! namsbc_ssr surface boundary condition : sea surface restoring 358 !----------------------------------------------------------------------- 359 ! cn_dir directory in which the model is executed 360 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 361 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 362 ! (=1, salt flux, concentration/dillution emps only) 363 ! (=2, volume flux, both emp and emps are updated) 364 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 365 ! deds magnitude of the damping on salinity [mm/day] 366 ! THE ORDER OF THE FILES MATTER: 367 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 368 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 369 &namsbc_ssr 370 cn_dir = './' 371 nn_sstr = 0 372 nn_sssr = 0 373 dqdt = -40. 374 deds = 27.7 375 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 376 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 203 &namqsr ! penetrative solar radiation 204 !----------------------------------------------------------------------- 205 ln_traqsr = .true. ! penetrative solar radiation (T) or not (F) 206 rabs = 0.58 ! fraction of qsr associated with xsi1 207 xsi1 = 0.35 ! first depth of extinction 208 xsi2 = 23.0 ! second depth of extinction 209 / 210 !----------------------------------------------------------------------- 211 &namsbc_rnf ! runoffs namelist surface boundary condition 212 !----------------------------------------------------------------------- 213 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 214 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 215 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. , 'yearly' 216 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. , 'yearly' 217 ! 218 cn_dir = './' ! root directory for the location of the runoff files 219 ln_rnf_emp = .true. ! runoffs included into precipitation field (T) or into a file (F) 220 ln_rnf_mouth = .true. ! specific treatment at rivers mouths 221 rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used 222 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 223 / 224 !----------------------------------------------------------------------- 225 &namsbc_ssr ! surface boundary condition : sea surface restoring 226 !----------------------------------------------------------------------- 227 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 228 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 229 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly' 230 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' 231 ! 232 cn_dir = './' ! root directory for the location of the runoff files 233 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 234 nn_sssr = 1 ! add a damping term in the surface freshwater flux (=1) or not (=0) 235 dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 236 deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu] 377 237 / 378 238 !----------------------------------------------------------------------- 379 ! namalb albedo parameters 380 !----------------------------------------------------------------------- 381 ! cgren correction of the snow or ice albedo to take into account 382 ! albice albedo of melting ice in the arctic and antarctic 383 ! alphd coefficients for linear interpolation used to compute albedo 384 ! between two extremes values (Pyane, 1972) 385 ! alphc " " 386 ! alphdi " " 387 &namalb 388 cgren = 0.06 389 albice = 0.5 390 alphd = 0.80 391 alphc = 0.65 392 alphdi = 0.72 393 / 394 239 &namalb ! albedo parameters 240 !----------------------------------------------------------------------- 241 cgren = 0.0 ! correction of the snow or ice albedo to take into account the 242 albice = 0.5 ! albedo of melting ice in the arctic and antarctic 243 alphd = 0.80 ! coefficients for linear interpolation used to 244 alphc = 0.65 ! compute albedo between two extremes values 245 alphdi = 0.72 ! (Pyane, 1972) 246 / 395 247 !!====================================================================== 396 248 !! *** Lateral boundary condition *** … … 399 251 !! namcla cross land advection 400 252 !! namobc open boundaries parameters ("key_obc") 401 !!====================================================================== 402 403 !----------------------------------------------------------------------- 404 ! namlbc lateral momentum boundary condition 405 !----------------------------------------------------------------------- 406 ! shlat shlat = 0 , free slip ; 0 < shlat < 2 , partial slip 407 ! shlat = 2 , no slip ; 2 < shlat , strong slip 408 &namlbc 409 shlat = 2. 410 / 411 !----------------------------------------------------------------------- 412 ! namcla cross land advection 413 !----------------------------------------------------------------------- 414 ! n_cla advection between 2 ocean pts separates by land 415 &namcla 416 n_cla = 1 417 / 418 !----------------------------------------------------------------------- 419 ! namobc open boundaries parameters (#ifdef key_obc) 420 !----------------------------------------------------------------------- 421 ! nobc_dta = 0 the obc data are equal to the initial state 422 ! = 1 the obc data are read in 'obc .dta' files 423 ! rdpeob time relaxation (days) for the east open boundary 424 ! rdpwob time relaxation (days) for the west open boundary 425 ! rdpnob time relaxation (days) for the north open boundary 426 ! rdpsob time relaxation (days) for the south open boundary 427 ! zbsic1 barotropic stream function on isolated coastline 1 428 ! zbsic2 barotropic stream function on isolated coastline 2 429 ! zbsic3 barotropic stream function on isolated coastline 3 430 ! ln_obc_clim climatological obc data files (default T) 431 ! ln_vol_cst total volume conserved 432 &namobc 433 nobc_dta = 0 434 rdpein = 1. 435 rdpwin = 1. 436 rdpnin = 30. 437 rdpsin = 1. 438 rdpeob = 1500. 439 rdpwob = 15. 440 rdpnob = 150. 441 rdpsob = 15. 442 zbsic1 = 140.e+6 443 zbsic2 = 1.e+6 444 zbsic3 = 0. 445 ln_obc_clim = .true. 446 ln_vol_cst = .false. 447 / 448 253 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 254 !! nambdy Unstructured open boundaries ("key_bdy") 255 !! namtide Tidal forcing at open boundaries ("key_bdy_tides") 256 !!====================================================================== 257 258 !----------------------------------------------------------------------- 259 &namlbc ! lateral momentum boundary condition 260 !----------------------------------------------------------------------- 261 shlat = 2. ! shlat = 0 : free slip 262 ! 0 < shlat < 2 : partial slip 263 ! shlat = 2 : no slip 264 ! 2 < shlat : strong slip 265 / 266 !----------------------------------------------------------------------- 267 &namcla ! cross land advection 268 !----------------------------------------------------------------------- 269 n_cla = 0 ! advection between 2 ocean pts separates by land 270 / 271 !----------------------------------------------------------------------- 272 &namobc ! open boundaries parameters ("key_obc") 273 !----------------------------------------------------------------------- 274 nobc_dta = 1 ! = 0 the obc data are equal to the initial state 275 ! = 1 the obc data are read in 'obc.dta' files 276 cffile = 'annual' ! set to annual if obc datafile hold 1 year of data 277 ! set to monthly if obc datafile hold 1 month of data 278 rdpein = 1. ! ??? 279 rdpwin = 1. ! ??? 280 rdpnin = 1. ! ??? 281 rdpsin = 1. ! ??? 282 rdpeob = 3000. ! time relaxation (days) for the east open boundary 283 rdpwob = 15. ! " " " west " 284 rdpnob = 3000. ! " " " north " 285 rdpsob = 15. ! " " " south " 286 zbsic1 = 140.e+6 ! barotropic stream function on first isolated coastline 287 zbsic2 = 1.e+6 ! " " second " 288 zbsic3 = 0. ! " " thrid " 289 ln_obc_clim= .false. ! climatological obc data files (T) or not (F) 290 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 291 / 292 !----------------------------------------------------------------------- 293 &namagrif ! ("key_agrif") 294 !----------------------------------------------------------------------- 295 nbclineupdate = 3 ! baroclinic update frequency 296 ln_spc_dyn = .true. ! use 0 as special value for dynamics 297 visc_tra = 2880. ! viscosity coeeficient for tracers sponge layer 298 visc_dyn = 2880. ! viscosity coeeficient for dynamics sponge layer 299 / 300 !----------------------------------------------------------------------- 301 &nambdy ! unstructured open boundaries parameters ("key_bdy") 302 !----------------------------------------------------------------------- 303 filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) 304 filbdy_data_T = 'bdydata_grid_T.nc' ! name of data file (T-points) 305 filbdy_data_U = 'bdydata_grid_U.nc' ! name of data file (U-points) 306 filbdy_data_V = 'bdydata_grid_V.nc' ! name of data file (V-points) 307 ln_bdy_clim = .false. ! contain 1 (T) or 12 (F) time dumps and be cyclic 308 ln_bdy_vol = .true. ! total volume correction (see volbdy parameter) 309 ln_bdy_mask = .false. ! boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F) 310 ln_bdy_tides = .true. ! Apply tidal harmonic forcing with Flather condition 311 ln_bdy_dyn_fla = .true. ! Apply Flather condition to velocities 312 ln_bdy_tra_frs = .false. ! Apply FRS condition to temperature and salinity 313 ln_bdy_dyn_frs = .false. ! Apply FRS condition to velocities 314 nbdy_dta = 1 ! = 0, bdy data are equal to the initial state 315 ! = 1, bdy data are read in 'bdydata .nc' files 316 nb_rimwidth = 9 ! width of the relaxation zone 317 volbdy = 0 ! = 0, the total water flux across open boundaries is zero 318 ! = 1, the total volume of the system is conserved 319 / 320 !----------------------------------------------------------------------- 321 &namtide ! tidal forcing at unstructured boundaries 322 !----------------------------------------------------------------------- 323 filtide = 'bdytide_' ! file name root of tidal forcing files 324 tide_cpt = 'M2','S1' ! names of tidal components used 325 tide_speed = 28.984106, 15.000001 ! phase speeds of tidal components (deg/hour) 326 ln_tide_date = .false. ! adjust tidal harmonics for start date of run 327 / 449 328 !!====================================================================== 450 329 !! *** Bottom boundary condition *** … … 454 333 !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") 455 334 !!====================================================================== 456 !----------------------------------------------------------------------- 457 ! nambfr bottom friction 458 !----------------------------------------------------------------------- 459 ! nbotfr type of bottom friction : = 0 : no slip ; = 1 : linear friction 460 ! = 3 : free slip ; = 2 : non linear friction 461 ! bfri1 bottom drag coefficient (linear case) 462 ! bfri2 bottom drag coefficient (non linear case) 463 ! bfeb2 bottom turbulent kinetic energy (m2/s2) 464 &nambfr 465 nbotfr = 1 466 bfri1 = 4.e-4 467 bfri2 = 1.e-3 468 bfeb2 = 2.5e-3 469 / 470 !----------------------------------------------------------------------- 471 ! nambbc bottom temperature boundary condition 472 !----------------------------------------------------------------------- 473 ! ngeo_flux = 0 no geothermal heat flux 474 ! = 1 constant geothermal heat flux 475 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 476 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 477 &nambbc 478 ngeo_flux = 2 479 ngeo_flux_const = 86.4e-3 480 / 481 !----------------------------------------------------------------------- 482 ! nambbl bottom boundary layer scheme 483 !----------------------------------------------------------------------- 484 ! ! diffusive bbl ("key_trabbl") 485 ! ! advective bbl ("key_trabbl_adv") 486 ! atrbbl mixing coefficient of the bottom boundary layer scheme (m2/s) 487 &nambbl 488 atrbbl = 10000. 489 / 490 335 336 !----------------------------------------------------------------------- 337 &nambfr ! bottom friction 338 !----------------------------------------------------------------------- 339 nbotfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction 340 ! = 3 : free slip, = 1 : linear friction 341 bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 342 bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 343 bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 344 / 345 !----------------------------------------------------------------------- 346 &nambbc ! bottom temperature boundary condition 347 !----------------------------------------------------------------------- 348 ngeo_flux = 2 ! geothermal heat flux = 0 no flux considered 349 ! = 1 constant flux 350 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 351 ngeo_flux_const = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 352 / 353 !----------------------------------------------------------------------- 354 &nambbl ! bottom boundary layer scheme 355 !----------------------------------------------------------------------- 356 ! ! diffusive bbl ("key_trabbl") 357 ! ! advective bbl ("key_trabbl_adv") 358 atrbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 359 / 491 360 !!====================================================================== 492 361 !! Tracer (T & S ) namelists … … 499 368 500 369 !----------------------------------------------------------------------- 501 ! nameos ocean physical parameters 502 !----------------------------------------------------------------------- 503 ! neos type of equation of state and Brunt-Vaisala frequency 504 ! ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 505 ! ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 506 ! ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 507 ! ralpha thermal expension coefficient (neos= 1 or 2) 508 ! rbeta saline expension coefficient (neos= 2) 509 &nameos 510 neos = 0 511 ralpha = 2.e-4 512 rbeta = 0.001 513 / 514 !----------------------------------------------------------------------- 515 ! nam_traadv advection scheme for tracer 516 !----------------------------------------------------------------------- 517 ! ln_traadv_cen2 2nd order centered scheme (default T) 518 ! ln_traadv_tvd TVD scheme (default F) 519 ! ln_traadv_muscl MUSCL scheme (default F) 520 ! ln_traadv_muscl2 MUSCL scheme + cen2 at boundaries (default F) 521 ! ln_traadv_ubs UBS scheme (default F) 522 ! ln_traadv_qck QUICKEST scheme (default F) 523 &nam_traadv 524 ln_traadv_cen2 = .true. 525 ln_traadv_tvd = .false. 526 ln_traadv_muscl = .false. 527 ln_traadv_muscl2 = .false. 528 ln_traadv_ubs = .false. 529 ln_traadv_qck = .false. 530 / 531 !----------------------------------------------------------------------- 532 ! nam_traldf lateral diffusion scheme for tracer 533 !----------------------------------------------------------------------- 534 ! Type of the operator : 535 ! ln_traldf_lap laplacian operator (default T) 536 ! ln_traldf_bilap bilaplacian operator (default F) 537 ! Direction of action : 538 ! ln_traldf_level iso-level (default F) 539 ! ln_traldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" when ln_sco=T) 540 ! ln_traldf_iso iso-neutral (default T) (require "key_ldfslp") 541 ! Coefficient 542 ! aht0 horizontal eddy diffusivity for tracers (m2/s) 543 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s) 544 ! aeiv0 eddy induced velocity coefficient (m2/s) (require "key_traldf_eiv") 545 &nam_traldf 546 ln_traldf_lap = .true. 547 ln_traldf_bilap = .false. 548 ln_traldf_level = .false. 549 ln_traldf_hor = .false. 550 ln_traldf_iso = .true. 551 aht0 = 2000. 552 ahtb0 = 0. 553 aeiv0 = 2000. 554 / 555 !----------------------------------------------------------------------- 556 ! namtdp tracer newtonian damping ('key_tradmp') 557 !----------------------------------------------------------------------- 558 ! ndmp type of damping in temperature and salinity 559 ! ! ='latitude', damping poleward of 'ndmp' degrees and function 560 ! ! of the distance-to-coast. Red and Med Seas as ndmp=-1 561 ! ! =-1 damping only in Med and Red Seas 562 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 563 ! nmldmp type of damping in the mixed layer 564 ! ! =0 damping throughout the water column 565 ! ! =1 no damping in the mixed layer defined by avt >5cm2/s 566 ! ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 567 ! sdmp surface time scale for internal damping (days) 568 ! bdmp bottom time scale for internal damping (days) 569 ! hdmp depth of transition between sdmp and bdmp (meters) 570 &namtdp 571 ndmp = -1 572 ndmpf = 1 573 nmldmp = 1 574 sdmp = 50. 575 bdmp = 360. 576 hdmp = 800. 577 / 578 370 &nameos ! ocean physical parameters 371 !----------------------------------------------------------------------- 372 neos = 0 ! type of equation of state and Brunt-Vaisala frequency 373 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 374 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 375 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 376 ralpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 377 rbeta = 0.001 ! saline expension coefficient (neos= 2) 378 / 379 !----------------------------------------------------------------------- 380 &nam_traadv ! advection scheme for tracer 381 !----------------------------------------------------------------------- 382 ln_traadv_cen2 = .true. ! 2nd order centered scheme 383 ln_traadv_tvd = .false. ! TVD scheme 384 ln_traadv_muscl = .false. ! MUSCL scheme 385 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 386 ln_traadv_ubs = .false. ! UBS scheme 387 / 388 !----------------------------------------------------------------------- 389 &nam_traldf ! lateral diffusion scheme for tracer 390 !----------------------------------------------------------------------- 391 ! ! Type of the operator : 392 ln_traldf_lap = .true. ! laplacian operator 393 ln_traldf_bilap = .false. ! bilaplacian operator 394 ! Direction of action : 395 ln_traldf_level = .false. ! iso-level 396 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 397 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 398 ! ! Coefficient 399 aht0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 400 ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 401 aeiv0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 402 / 403 !----------------------------------------------------------------------- 404 &namtdp ! tracer newtonian damping ('key_tradmp') 405 !----------------------------------------------------------------------- 406 ndmp = -1 ! type of damping in temperature and salinity 407 ! ='latitude', damping poleward of 'ndmp' degrees and function 408 ! of the distance-to-coast. Red and Med Seas as ndmp=-1 409 ! =-1 damping only in Med and Red Seas 410 ndmpf = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 411 nmldmp = 1 ! type of damping: =0 damping throughout the water column 412 ! =1 no damping in the mixed layer defined by avt >5cm2/s ) 413 ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 414 sdmp = 50. ! surface time scale for internal damping (days) 415 bdmp = 360. ! bottom time scale for internal damping (days) 416 hdmp = 800. ! depth of transition between sdmp and bdmp (meters) 417 / 579 418 !!====================================================================== 580 419 !! *** Dynamics namelists *** … … 589 428 590 429 !----------------------------------------------------------------------- 591 ! nam_dynadv formulation of the momentum advection 592 !----------------------------------------------------------------------- 593 ! ln_dynadv_vec vector form (T) or flux form (F) (default T) 594 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T) 595 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F) 596 &nam_dynadv 597 ln_dynadv_vec = .true. 598 ln_dynadv_cen2 = .false. 599 ln_dynadv_ubs = .false. 430 &nam_dynadv ! formulation of the momentum advection 431 !----------------------------------------------------------------------- 432 ln_dynadv_vec = .false. ! vector form (T) or flux form (F) 433 ln_dynadv_cen2= .true. ! flux form - 2nd order centered scheme 434 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 600 435 / 601 436 !----------------------------------------------------------------------- 602 ! nam_dynvor option of physics/algorithm (not control by CPP keys) 603 !----------------------------------------------------------------------- 604 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T) 605 ! ln_dynvor_ene " " : energy conserving scheme (default F) 606 ! ln_dynvor_mix " " : mixed scheme (default F) 607 ! ln_dynvor_een " " : energy & enstrophy scheme (default F) 608 &nam_dynvor 609 ln_dynvor_ene = .false. 610 ln_dynvor_ens = .true. 611 ln_dynvor_mix = .false. 612 ln_dynvor_een = .false. 613 / 614 !----------------------------------------------------------------------- 615 ! nam_dynhpg Hydrostatic pressure gradient option 616 !----------------------------------------------------------------------- 617 ! type of pressure gradient scheme (choose one only!) 618 ! ln_hpg_zco z-coordinate - full steps (default T) 619 ! ln_hpg_zps z-coordinate - partial steps (interpolation) 620 ! ln_hpg_sco s-coordinate (standard jacobian formulation) 621 ! ln_hpg_hel s-coordinate (helsinki modification) 622 ! ln_hpg_wdj s-coordinate (weighted density jacobian) 623 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial) 624 ! ln_hpg_rot s-coordinate (ROTated axes scheme) 625 ! parameters 626 ! gamm weighting coefficient (wdj scheme) 627 &nam_dynhpg 628 ln_hpg_zco = .false. 629 ln_hpg_zps = .true. 630 ln_hpg_sco = .false. 631 ln_hpg_hel = .false. 632 ln_hpg_wdj = .false. 633 ln_hpg_djc = .false. 634 ln_hpg_rot = .false. 635 gamm = 0.e0 636 / 637 !----------------------------------------------------------------------- 638 ! namflg algorithm flags (algorithm not control by CPP keys) 639 !----------------------------------------------------------------------- 640 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T) 641 ! centered time scheme (F) 642 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0) 643 &namflg 644 ln_dynhpg_imp = .false. 645 nn_dynhpg_rst = 0 646 / 647 !----------------------------------------------------------------------- 648 ! nam_dynspg surface pressure gradient (CPP key only) 649 !----------------------------------------------------------------------- 650 ! ! explicit free surface ("key_dynspg_exp") 651 ! ! filtered free surface ("key_dynspg_flt") 652 ! ! split-explicit free surface ("key_dynspg_ts") 653 ! ! rigid-lid ("key_dynspg_rl") 654 655 !----------------------------------------------------------------------- 656 ! nam_dynldf lateral diffusion on momentum 657 !----------------------------------------------------------------------- 658 ! Type of the operator : 659 ! ln_dynldf_lap laplacian operator (default T) 660 ! ln_dynldf_bilap bilaplacian operator (default F) 661 ! Direction of action : 662 ! ln_dynldf_level iso-level (default F) 663 ! ln_dynldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" if ln_sco=T) 664 ! ln_dynldf_iso iso-neutral (default T) (require "key_ldfslp") 665 ! Coefficient 666 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s) 667 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s) 668 &nam_dynldf 669 ln_dynldf_lap = .true. 670 ln_dynldf_bilap = .false. 671 ln_dynldf_level = .false. 672 ln_dynldf_hor = .true. 673 ln_dynldf_iso = .false. 674 ahm0 = 40000. 675 ahmb0 = 0. 676 / 677 437 &nam_dynvor ! option of physics/algorithm (not control by CPP keys) 438 !----------------------------------------------------------------------- 439 ln_dynvor_ene = .false. ! enstrophy conserving scheme 440 ln_dynvor_ens = .true. ! energy conserving scheme 441 ln_dynvor_mix = .false. ! mixed scheme 442 ln_dynvor_een = .false. ! energy & enstrophy scheme 443 / 444 !----------------------------------------------------------------------- 445 &nam_dynhpg ! Hydrostatic pressure gradient option 446 !----------------------------------------------------------------------- 447 ln_hpg_zco = .true. ! z-coordinate - full steps 448 ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) 449 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 450 ln_hpg_hel = .false. ! s-coordinate (helsinki modification) 451 ln_hpg_wdj = .false. ! s-coordinate (weighted density jacobian) 452 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 453 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 454 gamm = 0.e0 ! weighting coefficient (wdj scheme) 455 / 456 !----------------------------------------------------------------------- 457 &namflg ! algorithm flags (algorithm not control by CPP keys) 458 !----------------------------------------------------------------------- 459 ln_dynhpg_imp = .false. ! hydrostatic pressure gradient: semi-implicit time scheme (T) 460 ! centered time scheme (F) 461 nn_dynhpg_rst = 0 ! add dynhpg implicit variables in restart ot not (1/0) 462 / 463 !----------------------------------------------------------------------- 464 !nam_dynspg ! surface pressure gradient (CPP key only) 465 !----------------------------------------------------------------------- 466 ! ! explicit free surface ("key_dynspg_exp") 467 ! ! filtered free surface ("key_dynspg_flt") 468 ! ! split-explicit free surface ("key_dynspg_ts") 469 ! ! rigid-lid ("key_dynspg_rl") 470 471 !----------------------------------------------------------------------- 472 &nam_dynldf ! lateral diffusion on momentum 473 !----------------------------------------------------------------------- 474 ! ! Type of the operator : 475 ln_dynldf_lap = .true. ! laplacian operator 476 ln_dynldf_bilap = .false. ! bilaplacian operator 477 ! ! Direction of action : 478 ln_dynldf_level = .false. ! iso-level 479 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 480 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 481 ! Coefficient 482 ahm0 = 40000. ! horizontal eddy viscosity [m2/s] 483 ahmb0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 484 / 678 485 !!====================================================================== 679 486 !! Tracers & Dynamics vertical physics namelists … … 688 495 689 496 !----------------------------------------------------------------------- 690 ! namzdf vertical physics 691 !----------------------------------------------------------------------- 692 ! ln_zdfevd enhanced vertical diffusion (default T) 693 ! ln_zdfnpc Non-Penetrative Convection (default T) 694 ! avm0 vertical eddy viscosity for the dynamic (m2/s) (also background Kz if not "key_zdfcst") 695 ! avt0 vertical eddy diffusivity for tracers (m2/s) (also background Kz if not "key_zdfcst") 696 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 697 ! n_evdm = 0 apply enhanced mixing on tracer only 698 ! = 1 apply enhanced mixing on both tracer and momentum 699 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 700 ! (=F) euler backward (F) 701 ! n_zdfexp number of sub-timestep for time splitting scheme 702 &namzdf 703 ln_zdfevd = .true. 704 ln_zdfnpc = .false. 705 avm0 = 1.2e-4 706 avt0 = 1.2e-5 707 avevd = 100. 708 n_evdm = 1 709 ln_zdfexp = .false. 710 n_zdfexp = 3 711 / 712 !----------------------------------------------------------------------- 713 ! namnpc non penetrative convection 714 !----------------------------------------------------------------------- 715 ! nnpc1 non penetrative convective scheme frequency 716 ! nnpc2 non penetrative convective scheme print frequency 717 &namnpc 718 nnpc1 = 1 719 nnpc2 = 365 720 / 721 !----------------------------------------------------------------------- 722 ! namric richardson number dependent vertical diffusion ("key_zdfric" ) 723 !----------------------------------------------------------------------- 724 ! avmri maximum value of the vertical viscosity 725 ! alp coefficient of the parameterization 726 ! nric coefficient of the parameterization 727 &namric 728 avmri = 100.e-4 729 alp = 5. 730 nric = 2 731 / 732 !----------------------------------------------------------------------- 733 ! namtke turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 734 !----------------------------------------------------------------------- 735 ! ln_rstke flag to restart with tke from a run without tke (default F) 736 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 737 ! ediss coef. of the Kolmogoroff dissipation 738 ! ebb coef. of the surface input of tke 739 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 740 ! emin minimum value of tke (m^2/s^2) 741 ! emin0 surface minimum value of tke (m^2/s^2) 742 ! nitke number of restart iterative loops 743 ! ri_c critic richardson number = 2/9 = 0.22222222 (hard coded) 744 ! nmxl length used = 0 bounded by the distance to surface and bottom 745 ! ! = 1 bounded by the local vertical scale factor 746 ! ! = 2 first vertical derivative of mixing length bounded by 1 747 ! npdl Prandtl number = 0 no vertical prandtl number (avt=avm) 748 ! ! = 1 prandtl number function of richarson number (avt=pdl*avm) 749 ! ! = 2 same as = 1 but a shapiro filter is applied on pdl 750 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 751 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 752 &namtke 753 ln_rstke = .false. 754 ediff = 0.1 755 ediss = 0.7 756 ebb = 3.75 757 efave = 1. 758 emin = 1.e-6 759 emin0 = 1.e-4 760 nitke = 50 761 nmxl = 2 762 npdl = 1 763 navb = 0 764 / 765 !----------------------------------------------------------------------- 766 ! namkpp K-Profile Parameterization dependent vertical mixing ("key_zdfkpp" ) 767 !----------------------------------------------------------------------- 768 ! ln_kpprimix shear instability mixing (default T) 769 ! difmiw constant internal wave viscosity (m2/s) 770 ! difsiw constant internal wave diffusivity (m2/s) 771 ! Riinfty local Richardson Number limit for shear instability 772 ! difri maximum shear mixing at Rig = 0 (m2/s) 773 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 774 ! difcon maximum mixing in interior convection (m2/s) 775 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 776 ! navb = 0/1 flag for constant or profile background avt 777 &namkpp 778 ln_kpprimix = .true. 779 difmiw = 1.e-04 780 difsiw = 0.1e-04 781 Riinfty = 0.8 782 difri = 0.0050 783 bvsqcon = -0.01e-07 784 difcon = 1. 785 navb = 0 786 nave = 1 787 / 788 !----------------------------------------------------------------------- 789 ! namddm double diffusive mixing parameterization ("key_zdfddm") 790 !----------------------------------------------------------------------- 791 ! avts maximum avs for dd mixing 792 ! hsbfr heat/salt buoyancy flux ratio 793 &namddm 794 avts = 1.e-4 795 hsbfr = 1.6 796 / 797 497 &namzdf ! vertical physics 498 !----------------------------------------------------------------------- 499 avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 500 avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 501 ln_zdfnpc = .false. ! convection: Non-Penetrative algorithm (T) or not (F) 502 ln_zdfevd = .true. ! convection: enhanced vertical diffusion (T) or not (F) 503 avevd = 100. ! vertical coefficient for enhanced diffusion scheme [m2/s] 504 n_evdm = 1 ! enhanced mixing apply on tracer (=0) or on tracer and momentum (=1) 505 ln_zdfexp = .false. ! split explicit (T) or implicit (F) time stepping 506 n_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 507 / 508 !----------------------------------------------------------------------- 509 &namnpc ! non penetrative convection 510 !----------------------------------------------------------------------- 511 nnpc1 = 1 ! non penetrative convective scheme computation frequency 512 nnpc2 = 365 ! non penetrative convective scheme print frequency 513 / 514 !----------------------------------------------------------------------- 515 &namric ! richardson number dependent vertical diffusion ("key_zdfric" ) 516 !----------------------------------------------------------------------- 517 avmri = 100.e-4 ! maximum value of the vertical viscosity 518 alp = 5. ! coefficient of the parameterization 519 nric = 2 ! coefficient of the parameterization 520 / 521 !----------------------------------------------------------------------- 522 &namtke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 523 !----------------------------------------------------------------------- 524 ln_rstke = .false. ! restart with tke from a run without tke (T) or not (F) 525 nn_itke = 50 ! number of iterative loops if ln_rstke=T 526 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 527 rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation 528 rn_ebb = 3.75 ! coef. of the surface input of tke 529 rn_efave = 1. ! boost of the tke diffusion ( avtke=rn_efave*avm ) 530 rn_emin = 1.e-6 ! minimum value of tke [m2/s2] 531 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] 532 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom 533 ! = 1 bounded by the local vertical scale factor 534 ! = 2 first vertical derivative of mixing length bounded by 1 535 ! = 3 same criteria as case 2 but applied in a different way 536 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 537 nn_avb = 0 ! profile for constant background used on avt & avm (=1) or not (=0) 538 nn_ave = 1 ! horizontal averaged on avt (=1) or not (=0) 539 ln_mxl0 = .false. ! mixing length scale surface value as function of wind stress (T) or not (F) 540 rn_lmin = 0.4 ! interior buoyancy lenght scale minimum value 541 rn_lmin0 = 0.4 ! surface buoyancy lenght scale minimum value 542 nn_etau = 0 ! exponentially deceasing penetration of tke due to internal & intertial waves 543 ! = 0 no penetration ( O(2 km) resolution) 544 ! = 1 additional tke source 545 ! = 2 additional tke source applied only at the base of the mixed layer 546 nn_htau = 2 ! type of exponential decrease of tke penetration 547 ! = 0 constant 10 m length scale 548 ! = 1 ??? 549 ! = 2 ??? 550 rn_efr = 0.05 ! fraction of surface tke value which penetrates inside the ocean 551 ln_lc = .false. ! Langmuir cell effect 552 rn_lc = 0.15 ! coef. associated to Langmuir cells 553 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 554 / 555 !------------------------------------------------------------------------ 556 &namkpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally: 557 ! "key_kppcustom" or "key_kpplktb") 558 !------------------------------------------------------------------------ 559 ln_kpprimix = .true. ! shear instability mixing 560 difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 561 difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 562 Riinfty = 0.8 ! local Richardson Number limit for shear instability 563 difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 564 bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 565 difcon = 1. ! maximum mixing in interior convection [m2/s] 566 navb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 567 nave = 1 ! constant (=0) or profile (=1) background on avt 568 / 569 !----------------------------------------------------------------------- 570 &namddm ! double diffusive mixing parameterization ("key_zdfddm") 571 !----------------------------------------------------------------------- 572 avts = 1.e-4 ! maximum avs (vertical mixing on salinity) 573 hsbfr = 1.6 ! heat/salt buoyancy flux ratio 574 / 798 575 !!====================================================================== 799 576 !! *** Miscelaneous namelists *** 800 577 !!====================================================================== 801 !! namsol elliptic solver / island / free surface 802 !!====================================================================== 803 804 !----------------------------------------------------------------------- 805 ! namsol elliptic solver / island / free surface 806 !----------------------------------------------------------------------- 807 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 808 ! (=2 successive-over-relaxation: sor) 809 ! (=3 FETI currently it does not work! ("key_feti") 810 ! nsol_arp absolute/relative (0/1) precision convergence test 811 ! nmin minimum of iterations for the SOR solver 812 ! nmax maximum of iterations for the SOR solver 813 ! nmod frequency of test for the SOR solver 814 ! eps absolute precision of the solver 815 ! resmax absolute precision for the SOR solver 816 ! sor optimal coefficient for SOR solver 817 ! epsisl absolute precision on stream function solver 818 ! nmisl maximum pcg iterations for island 819 ! rnu strength of the additional force used in free surface b.c. 820 &namsol 821 nsolv = 1 822 nsol_arp = 0 823 nmin = 300 824 nmax = 800 825 nmod = 10 826 eps = 1.E-6 827 resmax = 1.E-10 828 sor = 1.92 829 epsisl = 1.e-10 830 nmisl = 4000 831 rnu = 1. 578 !! nam_mpp Massively Parallel Processing ("key_mpp_mpi) 579 !! nam_mpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist") 580 !! namctl Control prints & Benchmark 581 !! namsol elliptic solver / island / free surface 582 !!====================================================================== 583 584 !----------------------------------------------------------------------- 585 &namsol ! elliptic solver / island / free surface 586 !----------------------------------------------------------------------- 587 nsolv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 588 ! =2 successive-over-relaxation (sor) 589 ! =3 FETI (fet) ("key_feti") 590 ! =4 sor with extra outer halo 591 nsol_arp = 0 ! absolute/relative (0/1) precision convergence test 592 nmin = 300 ! minimum of iterations for the SOR solver 593 nmax = 800 ! maximum of iterations for the SOR solver 594 nmod = 10 ! frequency of test for the SOR solver 595 eps = 1.e-6 ! absolute precision of the solver 596 resmax = 1.e-10 ! absolute precision for the SOR solver 597 sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 598 epsisl = 1.e-10 ! absolute precision on stream function solver 599 nmisl = 4000 ! maximum pcg iterations for island ("key_islands") 600 rnu = 1. ! strength of the additional force used in filtered free surface 601 / 602 !----------------------------------------------------------------------- 603 &nam_mpp ! Massively Parallel Processing ("key_mpp_mpi) 604 !----------------------------------------------------------------------- 605 c_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 606 ! buffer blocking send or immediate non-blocking sends, resp. 607 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 608 / 609 !----------------------------------------------------------------------- 610 &nam_mpp_dyndist ! Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 611 !----------------------------------------------------------------------- 612 jpni = 1 ! jpni number of processors following i 613 jpnj = 1 ! jpnj number of processors following j 614 jpnij = 1 ! jpnij number of local domains 615 / 616 !----------------------------------------------------------------------- 617 &namctl ! Control prints & Benchmark 618 !----------------------------------------------------------------------- 619 ln_ctl = .false. ! trends control print (expensive!) 620 nprint = 0 ! level of print (0 no extra print) 621 nictls = 0 ! start i indice of control sum (use to compare mono versus 622 nictle = 0 ! end i indice of control sum multi processor runs 623 njctls = 0 ! start j indice of control over a subdomain) 624 njctle = 0 ! end j indice of control 625 isplt = 1 ! number of processors in i-direction 626 jsplt = 1 ! number of processors in j-direction 627 nbench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 628 ! (no physical validity of the results) 629 nbit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test 630 ! of comparison between single and multiple processor runs 832 631 / 833 632 !!====================================================================== … … 842 641 843 642 !----------------------------------------------------------------------- 844 ! namtrd diagnostics on dynamics and/or tracer trends("key_trddyn" and/or "key_trdtra")643 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 845 644 ! or mixed-layer trends ('key_trdmld') 846 645 ! or barotropic vorticity ("key_trdvor") 847 646 !----------------------------------------------------------------------- 848 ! ntrd time step frequency dynamics and tracers trends 849 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 850 ! ln_trdmld_restart restart for ML diagnostics 851 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 852 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 853 &namtrd 854 ntrd = 365 855 nctls = 0 856 ln_trdmld_restart = .false. 857 ucf = 1. 858 ln_trdmld_instant = .false. 859 / 860 !----------------------------------------------------------------------- 861 ! namgap level mean model-data gap ('key_diagap') 862 !----------------------------------------------------------------------- 863 ! ngap time-step frequency of model-data gap computation 864 ! nprg time-step frequency of gap print in model output 865 &namgap 866 ngap = 15 867 nprg = 10 868 / 869 !----------------------------------------------------------------------- 870 ! namspr surface pressure diagnostic 871 !----------------------------------------------------------------------- 872 ! nmaxp maximum of iterations for the solver 873 ! epsp absolute precision of the solver 874 ! niterp number of iteration done by the solver 875 &namspr 876 nmaxp = 1000 877 epsp = 1.e-3 878 niterp = 400 879 / 880 !----------------------------------------------------------------------- 881 ! namflo float parameters ("key_float") 882 !----------------------------------------------------------------------- 883 ! ln_rstflo boolean term for float restart (true or false) 884 ! nwritefl frequency of float output file 885 ! nstockfl frequency of float restart file 886 ! ln_argo Argo type floats (stay at the surface each 10 days) 887 ! ln_flork4 = T trajectories computed with a 4th order Runge-Kutta 888 ! = F (default) computed with Blanke' scheme 889 &namflo 890 ln_rstflo = .false. 891 nwritefl = 75 892 nstockfl = 5475 893 ln_argo = .false. 894 ln_flork4 = .false. 895 / 896 !----------------------------------------------------------------------- 897 ! namptr Poleward Transport Diagnostic 898 !----------------------------------------------------------------------- 899 ! ln_diaptr logical flag for Poleward transport computation 900 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 901 ! need input basins mask file named "subbasins.nc" 902 ! nf_ptr Frequency of computation 903 &namptr 904 ln_diaptr = .false. 905 ln_subbas = .false. 906 nf_ptr = 15 907 / 647 ntrd = 365 ! time step frequency dynamics and tracers trends 648 nctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 649 ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 650 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 651 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 652 ln_trdmld_restart = .false. ! restart for ML diagnostics 653 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 654 / 655 !----------------------------------------------------------------------- 656 &namgap ! level mean model-data gap ('key_diagap') 657 !----------------------------------------------------------------------- 658 ngap = 15 ! time-step frequency of model-data gap computation 659 nprg = 10 ! time-step frequency of gap print in model output 660 / 661 !----------------------------------------------------------------------- 662 &namspr ! surface pressure diagnostic 663 !----------------------------------------------------------------------- 664 nmaxp = 1000 ! maximum of iterations for the solver 665 epsp = 1.e-3 ! absolute precision of the solver 666 niterp = 400 ! number of iteration done by the solver 667 / 668 !----------------------------------------------------------------------- 669 &namflo ! float parameters ("key_float") 670 !----------------------------------------------------------------------- 671 ln_rstflo = .false. ! float restart (T) or not (F) 672 nwritefl = 75 ! frequency of writing in float output file 673 nstockfl = 5475 ! frequency of creation of the float restart file 674 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 675 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 676 ! or computed with Blanke' scheme (F) 677 / 678 !----------------------------------------------------------------------- 679 &namptr ! Poleward Transport Diagnostic 680 !----------------------------------------------------------------------- 681 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 682 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 683 ! (orca configuration only, need input basins mask file named "subbasins.nc" 684 nf_ptr = 15 ! Frequency of ptr computation [time step] 685 / -
CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namelist_ice_ORCA2
r396 r561 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/LIM2 : 1 - dynamics/advection/thermo (namicerun) 3 !! namelists 2 - ice intialisation (namiceini) 4 !! 3 - ice dynamic (namicedyn) 5 !! 5 - ice advection (namicetrp) 6 !! 6 - thermodynamic (namicethd) 7 !! 7 - ice salinity (namicesal) 8 !! 8 - mechanical redistribution of ice (namiceitdme) 9 !! 3 - ice diagnostics (namicedia) 10 !! 9 - ice outputs (namiceout) 1 11 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 ! 3 ! ICE namelist : ice model option and parameter input 4 ! ------------- 5 ! 6 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 7 ! 8 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 9 ! namicerun Share parameters for dynamics/advection/thermo 10 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 11 ! 12 ! ln_limdyn : switch for ice dynamics (true) or not (false) 13 ! acrit(1/2) : minimum fraction for leads in the Northern (Southern) Hemisphere 14 ! hnsdif : computation of temperature in snow (=0.0) or not (=9999.0) 15 ! hicdif : computation of temperature in ice (=0.0) or not (=9999.0) 16 ! 17 &namicerun 18 ln_limdyn = .true. 19 acrit = 1.0e-06 , 1.0e-06 20 hsndif = 0.0 21 hicdif = 0.0 12 13 !----------------------------------------------------------------------- 14 &namicerun ! Share parameters for dynamics/advection/thermo 15 !----------------------------------------------------------------------- 16 cn_icerst_in = "restart_ice_in" ! suffix of ice restart name (input) 17 cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) 18 ln_limdyn = .true. ! ice dynamics (T) or thermodynamics only (F) 19 ln_limdmp = .false. ! restoring ice thickness and fraction leads (T) or not (F) 20 acrit = 1.0e-06 , 1.0e-06 ! minimum fraction for leads in the Northern (Southern) Hemisphere 21 hsndif = 0.0 ! computation of temperature in snow (=0.0) or not 22 hicdif = 0.0 ! computation of temperature in ice (=0.0) or not (=9999.0) 22 23 / 23 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 24 ! namiceini parameters for ice initialisation 25 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 26 ! 27 ! ttest : threshold water temperature for initial sea ice 28 ! hninn : initial snow thickness in the north 29 ! hginn : initial ice thickness in the north 30 ! alinn : initial leads area in the north 31 ! hnins : initial snow thickness in the south 32 ! hgins : initial ice thickness in the south 33 ! alins : initial leads area in the south 34 ! 35 &namiceini 36 ttest = 2.0 37 hninn = 0.5 38 hginn = 3.0 39 alinn = 0.05 40 hnins = 0.1 41 hgins = 1.0 42 alins = 0.1 24 !----------------------------------------------------------------------- 25 &namiceini ! ice initialisation 26 !----------------------------------------------------------------------- 27 ln_limini = .false. ! read the ice initial state in the file 'Ice_initialization.nc' (T) or not (F) 28 ttest = 2.0 ! threshold water temperature for initial sea ice 29 hninn = 0.5 ! initial snow thickness in the north 30 hginn = 3.0 ! initial ice thickness in the north 31 alinn = 0.05 ! initial leads area in the north 32 hnins = 0.1 ! same three parameter in the south 33 hgins = 1.0 ! " " south 34 alins = 0.1 ! " " south 43 35 / 44 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 45 ! namicedia parameters for ice diagnostics 46 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 47 ! 48 ! fmtinf : format of the output values 49 ! nfrinf : number of variables written in one line 50 ! ntmoy : instantaneous values of ice evolution or averaging 51 ! ninfo : frequency of ouputs on file ice_evolu in case of averaging 52 ! 53 &namicedia 54 fmtinf = '1PE13.5 ' 55 nfrinf = 4 56 ntmoy = 1 57 ninfo = 1 36 !----------------------------------------------------------------------- 37 &namicedyn ! ice dynamic 38 !----------------------------------------------------------------------- 39 epsd = 1.0e-20 ! tolerance parameter 40 alpha = 0.5 ! coefficient for semi-implicit coriolis 41 dm = 0.6e+03 ! diffusion constant for dynamics 42 nbiter = 1 ! number of sub-time steps for relaxation 43 nbitdr = 100 ! maximum number of iterations for relaxation 44 om = 0.5 ! relaxation constant 45 resl = 5.0e-05 ! maximum value for the residual of relaxation 46 cw = 5.0e-03 ! drag coefficient for oceanic stress 47 angvg = 0.0 ! turning angle for oceanic stress 48 pstar = 1.0e+04 ! 1st bulk-rheology parameter 49 c_rhg = 20.0 ! 2nd bulk-rhelogy parameter 50 etamn = 0.0e+07 ! minimun value for viscosity 51 creepl = 2.0e-08 ! creep limit 52 ecc = 2.0 ! eccentricity of the elliptical yield curve 53 ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice [m2/s] 58 54 / 59 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 60 ! namicedyn parameters for ice dynamic 61 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 62 ! 63 ! epsd : tolerance parameter 64 ! alpha : coefficient for semi-implicit coriolis 65 ! bound : boundary conditions (=0.0 no-slip, =1.0 free-slip) 66 ! dm : diffusion constant for dynamics. 67 ! nbiter : number of sub-time steps for relaxation 68 ! nbitdr : maximum number of iterations for relaxation 69 ! om : relaxation constant 70 ! resl : maximum value for the residual of relaxation 71 ! cw : drag coefficient for oceanic stress 72 ! angvg : turning angle for oceanic stress 73 ! pstar : first bulk-rheology parameter 74 ! c_rhg : second bulk-rhelogy parameter 75 ! etamn : minimun value for viscosity 76 ! creepl : creep limit 77 ! ecc : eccentricity of the elliptical yield curve 78 ! ahi0 : horizontal eddy diffusivity coefficient for sea-ice (m2/s) 79 ! 80 &namicedyn 81 epsd = 1.0e-20 82 alpha = 0.5 83 dm = 0.6e+03 84 nbiter = 1 85 nbitdr = 100 86 om = 0.5 87 resl = 5.0e-05 88 cw = 5.0e-03 89 angvg = 0.0 90 pstar = 1.0e+04 91 c_rhg = 20.0 92 etamn = 0.0e+07 93 creepl = 2.0e-08 94 ecc = 2.0 95 ahi0 = 350.e0 55 !----------------------------------------------------------------------- 56 &namicetrp ! ice transport 57 !----------------------------------------------------------------------- 58 bound = 0. ! boundary conditions (=0.0 no-slip, =1.0 free-slip) 96 59 / 97 ! 98 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 99 ! namicetrp parameters for ice advection 100 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 101 ! bound : boundary conditions (=0.0 no-slip, =1.0 free-slip) 102 &namicetrp 103 bound = 0. 104 / 105 106 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 107 ! namicethd parameters for thermodynamic computation 108 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 109 ! 60 !----------------------------------------------------------------------- 61 &namicethd ! ice thermodynamic 62 !----------------------------------------------------------------------- 110 63 ! hmelt : maximum melting at the bottom 111 64 ! hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere … … 127 80 ! alphs : coefficient for snow density when snow ice formation 128 81 ! 129 &namicethd 130 hmelt = -0.15 131 hiccrit = 0.3 , 0.3 132 hicmin = 0.2 133 hiclim = 0.05 134 amax = 0.999 135 swiqst = 1. 136 sbeta = 1. 137 parlat = 0.0 138 hakspl = 0.5 139 hibspl = 0.5 140 exld = 2.0 141 hakdif = 1.0 142 thth = 0.2 143 hnzst = 0.1 144 parsub = 0.0 145 alphs = 1.0 82 hmelt = -0.15 ! maximum melting at the bottom 83 hiccrit = 0.3 , 0.3 ! ice thickness for lateral accretion in the Northern (Southern) Hemisphere 84 hicmin = 0.2 ! ice thickness corr. to max. energy stored in brine pocket 85 hiclim = 0.05 ! minimum ice thickness 86 amax = 0.999 ! maximum lead fraction 87 swiqst = 1. ! energy stored in brine pocket (=1) or not (=0) 88 sbeta = 1. ! numerical caracteritic of the scheme for diffusion in ice 89 parlat = 0.0 ! percentage of energy used for lateral ablation 90 hakspl = 0.5 ! slope of distr. for Hakkinen-Mellor's lateral melting 91 hibspl = 0.5 ! slope of distribution for Hibler's lateral melting 92 exld = 2.0 ! exponent for leads-closure rate 93 hakdif = 1.0 ! coefficient for diffusions of ice and snow 94 thth = 0.2 ! threshold thickness for comp. of eq. thermal conductivity 95 hnzst = 0.1 ! thickness of the surf. layer in temp. computation 96 parsub = 0.0 ! switch for snow sublimation or not 97 alphs = 1.0 ! coefficient for snow density when snow ice formation 146 98 / 99 !----------------------------------------------------------------------- 100 &namicedia ! ice diagnostics 101 !----------------------------------------------------------------------- 102 fmtinf ='1PE13.5 ' ! format of the output values 103 nfrinf = 4 ! number of variables written in one line 104 ntmoy = 1 ! instantaneous values of ice evolution or averaging 105 ninfo = 1 ! frequency of ouputs on file ice_evolu in case of averaging 106 / 107 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 108 &namiceout ! parameters for outputs 109 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 110 noumef = 19 ! number of fields 147 111 ! 148 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 149 ! namiceout parameters for outputs 150 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 151 ! 152 ! noumef : number of fields 153 ! field 154 ! title : title of the field 155 ! name : name of the field 156 ! unit : unit of the field 157 ! nc : switch for saving field ( = 1 ) or not ( = 0 ) 158 ! cmulti : multiplicative constant 159 ! cadd : additive constant 160 ! 161 &namiceout 162 noumef = 19 163 field_1 = 'Snow thickness ', 'isnowthi', 'm ', 1, 1.0, 0.0 164 field_2 = 'Ice thickness ', 'iicethic', 'm ', 1, 1.0, 0.0 165 field_3 = 'Ice produced ', 'iiceprod', 'm/kt ', 1, 1.0, 0.0 166 field_4 = 'Ice concentration ', 'ileadfra', '% ', 1, -1.0, 1.0 167 field_5 = 'Ice temperature ', 'iicetemp', 'C ', 1, 1.0, -273.15 168 field_6 = 'Oceanic flux at the ice base ', 'ioceflxb', 'w/m2 ', 1, 1.0, 0.0 169 field_7 = 'Ice velocity u ', 'iicevelu', 'm/s ', 1, 1.0, 0.0 170 field_8 = 'Ice velocity v ', 'iicevelv', 'm/s ', 1, 1.0, 0.0 171 field_9 = 'Sea surface temperature ', 'isstempe', 'C ', 1, 1.0, -273.15 172 field_10 = 'Sea surface salinity ', 'isssalin', 'PSU ', 1, 1.0, 0.0 173 field_11 = 'Total flux at ocean surface ', 'iocetflx', 'w/m2 ', 1, 1.0, 0.0 174 field_12 = 'Solar flux at ocean surface ', 'iocesflx', 'w/m2 ', 1, 1.0, 0.0 175 field_13 = 'Non-solar flux at ocean surface ', 'iocwnsfl', 'w/m2 ', 1, 1.0, 0.0 176 field_14 = 'Salt flux at ocean surface ', 'iocesafl', 'kg/m2/kt', 1, 1.0, 0.0 177 field_15 = 'Wind stress u ', 'iocestru', 'Pa ', 1, 1.0, 0.0 178 field_16 = 'Wind stress v ', 'iocestrv', 'Pa ', 1, 1.0, 0.0 179 field_17 = 'Solar flux at ice/ocean surface ', 'iicesflx', 'w/m2 ', 1, 1.0, 0.0 180 field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2 ', 1, 1.0, 0.0 181 field_19 = 'Snow precipitation ', 'isnowpre', 'kg/day ', 1, 1.0, 0.0 112 ! ! title of the field ! name ! units ! save ! multipl. ! additive ! 113 ! ! ! ! ! or not ! factor ! factor ! 114 field_1 = 'Snow thickness ', 'isnowthi', 'm ', 1 , 1.0 , 0.0 115 field_2 = 'Ice thickness ', 'iicethic', 'm ', 1 , 1.0 , 0.0 116 field_3 = 'Ice produced ', 'iiceprod', 'm/kt ', 1 , 1.0 , 0.0 117 field_4 = 'Ice concentration ', 'ileadfra', '% ', 1 , -1.0 , 1.0 118 field_5 = 'Ice temperature ', 'iicetemp', 'C ', 1 , 1.0 , -273.15 119 field_6 = 'Oceanic flux at the ice base ', 'ioceflxb', 'w/m2 ', 1 , 1.0 , 0.0 120 field_7 = 'Ice velocity u ', 'iicevelu', 'm/s ', 1 , 1.0 , 0.0 121 field_8 = 'Ice velocity v ', 'iicevelv', 'm/s ', 1 , 1.0 , 0.0 122 field_9 = 'Sea surface temperature ', 'isstempe', 'C ', 1 , 1.0 , -273.15 123 field_10 = 'Sea surface salinity ', 'isssalin', 'PSU ', 1 , 1.0 , 0.0 124 field_11 = 'Total flux at ocean surface ', 'iocetflx', 'w/m2 ', 1 , 1.0 , 0.0 125 field_12 = 'Solar flux at ocean surface ', 'iocesflx', 'w/m2 ', 1 , 1.0 , 0.0 126 field_13 = 'Non-solar flux at ocean surface ', 'iocwnsfl', 'w/m2 ', 1 , 1.0 , 0.0 127 field_14 = 'Salt flux at ocean surface ', 'iocesafl', 'kg/m2/kt', 1 , 1.0 , 0.0 128 field_15 = 'Wind stress u ', 'iocestru', 'Pa ', 1 , 1.0 , 0.0 129 field_16 = 'Wind stress v ', 'iocestrv', 'Pa ', 1 , 1.0 , 0.0 130 field_17 = 'Solar flux at ice/ocean surface ', 'iicesflx', 'w/m2 ', 1 , 1.0 , 0.0 131 field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2 ', 1 , 1.0 , 0.0 132 field_19 = 'Snow precipitation ', 'isnowpre', 'kg/day ', 1 , 1.0 , 0.0 182 133 /
Note: See TracChangeset
for help on using the changeset viewer.
