Changeset 445 for CONFIG/trunk/IPSLCM5/EXP00
- Timestamp:
- 12/04/08 17:27:19 (16 years ago)
- Location:
- CONFIG/trunk/IPSLCM5/EXP00
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
CONFIG/trunk/IPSLCM5/EXP00/COMP/nemo.driver
r387 r445 14 14 fi 15 15 16 ##--Transformation namelist fortran 77 --> fortran 90 17 # ? ? Is it really necessary ? ? 18 #sed -e /:/d -e s/" \&END"/"\/"/ -e s/" \&"/"\&"/ ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} > ${SUBMIT_DIR}/PARAM/namelist 19 16 # Local function to find namelists parameters 17 supergrep () { 18 grep "^ *$1 *=" $2 | sed -e "s% *\!.*%%" 19 } 20 20 21 ##--Variables used by OPA -- 21 22 … … 32 33 # nmsh =1 create a mesh file (coordinates, scale factors, masks) 33 34 34 PAT_NO=$( grep 'no *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v off | grep -v \!)35 PAT_CEXPER=$( grep 'cexper *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)36 PAT_NIT000=$( grep 'nit000 *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)37 PAT_NITEND=$( grep 'nitend *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)38 PAT_NBISEX=$( grep 'nbisex *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE}| sed 's/ nbisex =//' |sed '/^[ ]$/d'|sed 's/,//' | grep -v \!)39 PAT_NWRITE=$( grep 'nwrite *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)40 PAT_RESTAR=$( grep 'ln_rstart *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)41 PAT_NSTOCK=$( grep 'nstock *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)42 PAT_NRSTAR=$( grep 'nrstdt *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)43 PAT_NPRG=$( grep 'nprg *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)44 PAT_NDATE0=$( grep 'ndate0 *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)45 PAT_NMSH=$( grep 'nmsh *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)46 ####### PAT_NWRIHF=$( grep 'nwrihf' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 47 48 ORCA_NSTOCK=$( echo $PAT_NSTOCK | sed "s/[a-z,A-Z,=]//g" )49 ORCA_RDT=$( egrep 'rdt *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | sed 's/ *rdt *=//' | sed 's/\.//' | grep -v \!)35 PAT_CEXPER=$( supergrep cexper ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 36 PAT_NIT000=$( supergrep nit000 ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 37 PAT_NITEND=$( supergrep nitend ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 38 PAT_NWRITE=$( supergrep nwrite ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 39 PAT_RESTAR=$( supergrep ln_rstart ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 40 PAT_NSTOCK=$( supergrep nstock ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 41 PAT_NRSTAR=$( supergrep nrstdt ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 42 PAT_NDATE0=$( supergrep ndate0 ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 43 PAT_NMSH=$( supergrep nmsh ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 44 45 ORCA_NSTOCK=$( echo $PAT_NSTOCK | sed "s/[a-z,A-Z,=]//g" ) 46 ORCA_RDT=$( supergrep rdt ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | sed 's/ *rdt *=//' | sed 's/\. *,//' ) 47 48 ##--Transformation namelist fortran 77 --> fortran 90 49 # ? ? Is it really necessary ? ? 50 #sed -e /:/d -e s/" \&END"/"\/"/ -e s/" \&"/"\&"/ ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} > ${SUBMIT_DIR}/PARAM/namelist 50 51 51 52 # Period Length In Days between DateBegin and first day of calendar 0001 01 01 … … 88 89 (( ORCA_NWRITE = 0 )) ;; 89 90 esac 90 ### (( ORCA_NWRITE = 1 )) ;91 ### PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}d ;92 ### R_OUT_OCE_NWRITE=${R_OUT_OCE_O_D} ;93 ### R_OUT_ICE_NWRITE=${R_OUT_ICE_O_D} ;94 ### WF1=${factor}D ;95 96 91 97 92 ##--Write Frequency Purpose .... concerning the second opa file type. -
CONFIG/trunk/IPSLCM5/EXP00/PARAM/namcouple_ORCA2xLMD144142
r408 r445 29 29 # 30 30 $CHANNEL 31 31 <lib_mpi> 32 32 <nproc_atm> <nproc_atm> 33 33 <nproc_oce> <nproc_oce> … … 38 38 # 39 39 $NFIELDS 40 2 140 22 41 41 $END 42 42 ############################################################################### … … 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 ############################################################################### … … 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. … … 320 321 # 321 322 ######################################################################### 322 # Field 12 : tprecip = ( rain_ice + snow_ice ) - evap_ice(a->o flx 7)323 # 324 COPEFWAT O IPr-Sub 25 86400 5flxat.nc EXPORTED323 # Field 12 : emp_ice = evap_ice - ( rain_ice + snow_ice ) (a->o flx 7) 324 # 325 COPEFWAT OPre-Sub 25 86400 4 flxat.nc EXPORTED 325 326 tlmd torc LAG=1800 326 327 P 0 P 2 327 INVERT CHECKIN MOZAIC BLASNEWCHECKOUT328 INVERT CHECKIN MOZAIC CHECKOUT 328 329 # Reverse 329 330 NORSUD WSTEST … … 332 333 # Interpolation method ou parametres mozaic 333 334 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 335 # CHECKOUT: indicate computation of global, land and sea field integrals. 336 INT=1 337 # 338 ########################################################################## 339 # Field 13 : sprecip = evap_ice - snow_ice (o->a flx 6) 340 # 341 COPEFICE OISub-Sn 25 86400 4 flxat.nc EXPORTED 343 342 tlmd torc LAG=1800 344 343 P 0 P 2 345 INVERT CHECKIN MOZAIC BLASNEWCHECKOUT344 INVERT CHECKIN MOZAIC CHECKOUT 346 345 # Reverse 347 346 NORSUD WSTEST … … 350 349 # Interpolation method ou parametres mozaic 351 350 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 351 # CHECKOUT: indicate computation of global, land and sea field integrals. 352 INT=1 353 # 354 ########################################################################## 355 # Field 14 : emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 356 # 357 COTOSPSU OOEva-Pr 28 86400 4 flxat.nc EXPORTED 361 358 tlmd torc LAG=1800 362 359 P 0 P 2 … … 368 365 # Interpolation method ou parametres mozaic 369 366 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 374 ########################################################################## 375 # Field 15 : solar heat flux on ocean (a->o flx 2) 367 # CHECKOUT: indicate computation of global, land and sea field integrals. 368 INT=1 369 # 370 ########################################################################## 371 # Field 15 : Evaporation Ice (a->o flx 8) 372 # 373 COICEVAP OIceEvap 27 86400 4 flxat.nc EXPORTED 374 tlmd torc LAG=1800 375 P 0 P 2 376 INVERT CHECKIN MOZAIC CHECKOUT 377 # Reverse 378 NORSUD WSTEST 379 # CHECKIN: indicate computation of global, land and sea field integrals. 380 INT=1 381 # Interpolation method ou parametres mozaic 382 wa2o.flx 92 2 32 383 # CHECKOUT: indicate computation of global, land and sea field integrals. 384 INT=1 385 # 386 ########################################################################## 387 # Field 16 : solar heat flux on ocean (a->o flx 2) 376 388 # 377 389 COSWFLDO O_QsrOce 7 86400 4 flxat.nc EXPORTED … … 387 399 # CHECKOUT: indicate computation of global, land and sea field integrals. 388 400 INT=1 401 # 389 402 ######################################################################### 390 # Field 1 6: Non solar heat flux on ocean (a->o flx 4)403 # Field 17 : Non solar heat flux on ocean (a->o flx 4) 391 404 # 392 405 CONSFLDO O_QnsOce 6 86400 4 flxat.nc EXPORTED … … 402 415 # CHECKOUT: indicate computation of global, land and sea field integrals. 403 416 INT=1 404 ########################################################################## 405 # Field 17 : solar heat flux on ice (a->o flx 1) 417 # 418 ########################################################################## 419 # Field 18 : solar heat flux on ice (a->o flx 1) 406 420 # 407 421 COSHFLIC O_QsrIce 7 86400 4 flxat.nc EXPORTED … … 419 433 # 420 434 ########################################################################## 421 # Field 1 8: Non solar heat flux on ice (a->o flx 3)435 # Field 19 : Non solar heat flux on ice (a->o flx 3) 422 436 # 423 437 CONSFLIC O_QnsIce 6 86400 4 flxat.nc EXPORTED … … 435 449 # 436 450 ############################################################################ 437 # Field 19: Non solar heat flux derivative (o->a flx 5)451 # Field 20 : Non solar heat flux derivative (o->a flx 5) 438 452 # 439 453 CODFLXDT O_dQnsdT 35 86400 4 flxat.nc EXPORTED … … 449 463 # CHECKOUT: indicate computation of global, land and sea field integrals. 450 464 INT=1 465 # 451 466 ########################################################################### 452 # Field 2 0: iceberg calving (a->o flx 12)453 # 454 CRWOCEIS OCalving 32 86400 5flxat.nc EXPORTED467 # Field 21 : iceberg calving (a->o flx 12) 468 # 469 CRWOCEIS OCalving 32 86400 4 flxat.nc EXPORTED 455 470 tlmd torc LAG=1800 456 471 P 0 P 2 457 INVERT CHECKIN MOZAIC BLASNEWCHECKOUT472 INVERT CHECKIN MOZAIC CHECKOUT 458 473 # Reverse 459 474 NORSUD WSTEST … … 462 477 # Interpolation method ou parametres mozaic 463 478 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 479 # CHECKOUT: indicate computation of global, land and sea field integrals. 480 INT=1 481 # 468 482 ############################################################################ 469 # Field 2 1: run-off (a->o flx 11)483 # Field 22 : run-off (a->o flx 11) 470 484 # 471 485 CRWOCERD O_Runoff 32 86400 5 flxat.nc EXPORTED … … 481 495 # change units from m/s to kg/s, change of sign 482 496 # weights convert from kg/s to kg/m^2/s 483 -1000.0 0497 1000.0 0 484 498 # CHECKOUT: indicate computation of global, land and sea field integrals. 485 499 INT=1 -
CONFIG/trunk/IPSLCM5/EXP00/PARAM/namcouple_ORCA2xLMD9671
r375 r445 38 38 # 39 39 $NFIELDS 40 2 140 22 41 41 $END 42 42 ############################################################################### … … 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 … … 106 106 # -------------------- 107 107 # Field 1 : Sea surface temperature (o->a 1) 108 # 108 109 O_SSTSST SISUTESW 1 86400 5 sstoc.nc EXPORTED 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) … … 182 181 # ------------------------- 183 182 ############################################################################ 184 # Field 17: wind stress along X axis (a->o tau 1)183 # Field 5 : wind stress along X axis (a->o tau 1) 185 184 # 186 185 COTAUXXU O_OTaux1 23 86400 6 flxat.nc EXPORTED … … 202 201 # 203 202 ########################################################################### 204 # 205 # 206 # Field 18 : stress along Y axis (a->o tau 1) 203 # Field 6 : stress along Y axis (a->o tau 1) 207 204 # 208 205 COTAUYYU O_OTauy1 23 86400 6 flxat.nc EXPORTED … … 224 221 # 225 222 ########################################################################### 226 # 227 # 228 # Field 19 :wind stress along Z axis (a->o tau 1) 223 # Field 7 :wind stress along Z axis (a->o tau 1) 229 224 # 230 225 COTAUZZU O_OTauz1 23 86400 6 flxat.nc EXPORTED … … 246 241 # 247 242 ########################################################################## 248 # 249 # Field 20 : wind stress along X axis 2 (a->o tau 2) 243 # Field 8 : wind stress along X axis 2 (a->o tau 2) 250 244 # 251 245 COTAUXXV O_OTaux2 24 86400 6 flxat.nc EXPORTED … … 267 261 # 268 262 ########################################################################## 269 # 270 # Field 21 : wind stress along Y axis 2 (a->o tau 2) 263 # Field 9 : wind stress along Y axis 2 (a->o tau 2) 271 264 # 272 265 COTAUYYV O_OTauy2 24 86400 6 flxat.nc EXPORTED … … 288 281 # 289 282 ########################################################################## 290 # 291 # Field 22 : wind stress along Z axis 2 (a->o tau 2) 283 # Field 10 : wind stress along Z axis 2 (a->o tau 2) 292 284 # 293 285 COTAUZZV O_OTauz2 24 86400 6 flxat.nc EXPORTED … … 309 301 ## 310 302 ########################################################################## 311 #312 303 # Field 11 : 10-m wind speed (a->o tau 2) 313 304 # … … 330 321 # 331 322 ######################################################################### 332 # 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 323 # Field 12 : emp_ice = evap_ice - ( rain_ice + snow_ice ) (a->o flx 7) 324 # 325 COPEFWAT OPre-Sub 25 86400 4 flxat.nc EXPORTED 337 326 tlmd torc LAG=1800 338 327 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 328 INVERT CHECKIN MOZAIC CHECKOUT 329 # Reverse 330 NORSUD WSTEST 331 # CHECKIN: indicate computation of global, land and sea field integrals. 332 INT=1 333 # Interpolation method ou parametres mozaic 334 wa2o.flx 92 2 14 348 335 # CHECKOUT: indicate computation of global, land and sea field integrals. 349 336 INT=1 350 337 # 351 338 ########################################################################## 352 # 353 # Field 11 : sprecip = snow_ice - evap_ice (o->a flx 6) 354 # 355 COPEFICE OISn-Sub 25 86400 5 flxat.nc EXPORTED 339 # Field 13 : sprecip = evap_ice - snow_ice (o->a flx 6) 340 # 341 COPEFICE OISub-Sn 25 86400 4 flxat.nc EXPORTED 356 342 tlmd torc LAG=1800 357 343 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 344 INVERT CHECKIN MOZAIC CHECKOUT 345 # Reverse 346 NORSUD WSTEST 347 # CHECKIN: indicate computation of global, land and sea field integrals. 348 INT=1 349 # Interpolation method ou parametres mozaic 350 wa2o.flx 92 2 14 367 351 # CHECKOUT: indicate computation of global, land and sea field integrals. 368 352 INT=1 369 353 ######################################################################### 370 # 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 354 # Field 14 : emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 355 # 356 COTOSPSU OOEva-Pr 28 86400 4 flxat.nc EXPORTED 375 357 tlmd torc LAG=1800 376 358 P 0 P 2 … … 382 364 # Interpolation method ou parametres mozaic 383 365 wa2o.flx 92 2 14 384 # change of sign385 ###-1.0 0386 366 # CHECKOUT: indicate computation of global, land and sea field integrals. 387 367 INT=1 388 368 ########################################################################## 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 369 # Field 15 : Evaporation Ice (a->o flx 8) 370 # 371 COICEVAP OIceEvap 27 86400 4 flxat.nc EXPORTED 372 tlmd torc LAG=1800 373 P 0 P 2 374 INVERT CHECKIN MOZAIC CHECKOUT 375 # Reverse 376 NORSUD WSTEST 377 # CHECKIN: indicate computation of global, land and sea field integrals. 378 INT=1 379 # Interpolation method ou parametres mozaic 380 wa2o.flx 92 2 14 381 # CHECKOUT: indicate computation of global, land and sea field integrals. 382 INT=1 407 383 # 408 384 ########################################################################### 409 # Field 6 : solar heat flux on ocean (a->o flx 2)385 # Field 16 : solar heat flux on ocean (a->o flx 2) 410 386 # 411 387 COSWFLDO O_QsrOce 7 86400 4 flxat.nc EXPORTED … … 421 397 # CHECKOUT: indicate computation of global, land and sea field integrals. 422 398 INT=1 399 # 423 400 ######################################################################### 424 401 # 425 # Field 8: Non solar heat flux on ocean (a->o flx 4)402 # Field 17 : Non solar heat flux on ocean (a->o flx 4) 426 403 # 427 404 CONSFLDO O_QnsOce 6 86400 4 flxat.nc EXPORTED … … 439 416 # 440 417 ########################################################################### 441 # Field 5: solar heat flux on ice (a->o flx 1)418 # Field 18 : solar heat flux on ice (a->o flx 1) 442 419 # 443 420 COSHFLIC O_QsrIce 7 86400 4 flxat.nc EXPORTED … … 454 431 INT=1 455 432 # 456 #457 433 ############################################################################ 458 # 459 # Field 7 : Non solar heat flux on ice (a->o flx 3) 434 # Field 19 : Non solar heat flux on ice (a->o flx 3) 460 435 # 461 436 CONSFLIC O_QnsIce 6 86400 4 flxat.nc EXPORTED … … 473 448 # 474 449 ############################################################################ 475 # 476 # Field 9 : Non solar heat flux derivative (o->a flx 5) 450 # Field 20 : Non solar heat flux derivative (o->a flx 5) 477 451 # 478 452 CODFLXDT O_dQnsdT 35 86400 4 flxat.nc EXPORTED … … 490 464 # 491 465 ######################################################################### 492 # 493 # Field 16 : iceberg calving (a->o flx 12) 494 # 495 CRWOCEIS OCalving 32 86400 5 flxat.nc EXPORTED 466 # Field 21 : iceberg calving (a->o flx 12) 467 # 468 CRWOCEIS OCalving 32 86400 4 flxat.nc EXPORTED 496 469 tlmd torc LAG=1800 497 470 P 0 P 2 498 INVERT CHECKIN MOZAIC BLASNEWCHECKOUT471 INVERT CHECKIN MOZAIC CHECKOUT 499 472 # Reverse 500 473 NORSUD WSTEST … … 503 476 # Interpolation method ou parametres mozaic 504 477 wa2o.cal 98 5 36 505 # change of sign 506 -1.0 0 507 # CHECKOUT: indicate computation of global, land and sea field integrals. 508 INT=1 478 # CHECKOUT: indicate computation of global, land and sea field integrals. 479 INT=1 480 # 509 481 ######################################################################### 510 # 511 # Field 15 : river run-off (a->o flx 11) 482 # Field 22 : river run-off (a->o flx 11) 512 483 # 513 484 CRWOCERD O_Runoff 32 86400 5 flxat.nc EXPORTED … … 523 494 # change units from m/s to kg/s, change of sign 524 495 # 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 ########################################################################## 496 1000.0 0 497 # CHECKOUT: indicate computation of global, land and sea field integrals. 498 INT=1 499 # 550 500 ########################################################################### 551 501 # -
CONFIG/trunk/IPSLCM5/EXP00/PARAM/namelist_ORCA2
r385 r445 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 rdtbt = 90. ! 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 190 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'191 cn_rcv_tau_nature = 'oce only' ! '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' … … 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.06 ! 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/trunk/IPSLCM5/EXP00/PARAM/namelist_ice_ORCA2
r375 r445 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.