| 366 | |
| 367 | ======================================================== |
| 368 | THE FOLLOWING IS THE SETUP NV USES: |
| 369 | ======================================================== |
| 370 | === 3.2 ENSEMBLE: fluxnet sites [http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation] === |
| 371 | ==== 3.2.1 create new experiment folder ==== |
| 372 | {{{ |
| 373 | cd ../../config/ORCHIDEE_OL/ |
| 374 | cp -fr ENSEMBLE/ ENSEMBLE_CNP/ |
| 375 | cd ENSEMBLE_CNP/ |
| 376 | }}} |
| 377 | |
| 378 | ==== 3.2.2 tell libIGCM about the new input files (write me an email if you need them) ==== |
| 379 | The configuration which are used are stored in |
| 380 | ../SPINUP/SUBJOB/OOL_SEC_STO/COMP/ |
| 381 | so, do as in 4.1.2: change the path of soil_param.nc and add the new files |
| 382 | |
| 383 | sechiba.card: |
| 384 | {{{ |
| 385 | (/home/scratch01/dgoll/ORCHIDEE/ORC-CNP/simulations/simple_forcing_DSG/soils_param.nc, .), \ |
| 386 | }}} |
| 387 | stomate.card: |
| 388 | {{{ |
| 389 | ListNonDel= (${R_IN}/SRF/reftemp.nc, .) \ |
| 390 | (/home/scratch01/dgoll/ORCHIDEE/ORC-CNP/simulations/simple_forcing_DSG/USDA_SoilSuborder.nc, .), \ |
| 391 | (/home/scratch01/dgoll/ORCHIDEE/ORC-CNP/simulations/simple_forcing_DSG/lithology.nc, .), \ |
| 392 | (/home/scratch01/dgoll/ORCHIDEE/ORC-CNP/simulations/simple_forcing_DSG/deposition.nc, .) |
| 393 | }}} |
| 394 | WARNING: the deposition file is a static field at the moment; it does not include N deposition either. At the moment N deposition is set to P deposition in the ORCHIDEE code. Anyway N deposition is not accounted for in the calculations for the N & P cycle for now. |
| 395 | |
| 396 | copy a run.def with all parameters needed to PARAM, for example this one: |
| 397 | {{{ |
| 398 | cp ~/ORCHIDEE/DEFS_n_CARDS/run.def_NV PARAM/run.def |
| 399 | }}} |
| 400 | |
| 401 | ==== 3.2.3 adjust the libIGCM ==== |
| 402 | modify the configuration of the SPINUP. These variables were originally used for the outdated spinup method using FORCESOIL. |
| 403 | Thus the variable names are misleading. |
| 404 | |
| 405 | modify them in ensemble.card and fluxnet.card: |
| 406 | {{{ |
| 407 | [SPINUP] |
| 408 | # Initialisation for spin-up : |
| 409 | # orchidee with sechiba alone (!!! if ok_stomate == n !!!) |
| 410 | duree_nostomate=0 |
| 411 | # orchidee with stomate |
| 412 | duree_inistomate=1 |
| 413 | # teststomate (only if duree_inistomate > 0) |
| 414 | duree_offlineini=0 |
| 415 | |
| 416 | # Loop over ORCHIDEE runs (used for spin-up) |
| 417 | # The whole job is restarted n_iter times |
| 418 | n_iter=1 |
| 419 | # orchidee with sechiba (and stomate if ok_stomate=y below) |
| 420 | duree_sechiba=500 |
| 421 | # teststomate |
| 422 | duree_stomate=0 |
| 423 | # forcesoil |
| 424 | duree_carbonsol=0 |
| 425 | |
| 426 | # Final run (full ORCHIDEE) |
| 427 | # This last parameter must be non-zero. |
| 428 | duree_final=1 |
| 429 | }}} |
| 430 | |
| 431 | |
| 432 | In case you want to check if the spinup works well, increase output level and precent libIGCM from removing the files. |
| 433 | modify in fluxnet.card and ensemble.card: |
| 434 | {{{ |
| 435 | # level of the history files for each iteration of the SPINUP |
| 436 | level_hist=10 |
| 437 | # level of the history files for the final iteration of the SPINUP |
| 438 | level_hist_final=10 |
| 439 | |
| 440 | # DEBUG mode for SPINUP |
| 441 | # |
| 442 | # This mode keep all SPINUP directory in ARCHIVE |
| 443 | # If disable, all ARCHIVE is automaticly cleaned. |
| 444 | DEBUG_SPIN=y |
| 445 | # If you don't want to keep old spinup steps, but last one |
| 446 | CONSERVE=y |
| 447 | }}} |
| 448 | |
| 449 | add the variables to turn on/off the nutrients cycles to PARAM/run.def: |
| 450 | {{{ |
| 451 | # switch for N cycle |
| 452 | STOMATE_OK_NCYCLE=y |
| 453 | # switch for P cycle |
| 454 | STOMATE_OK_PCYCLE=y |
| 455 | # switch to either impose N content of vegetation (to surpass N limitation): |
| 456 | STOMATE_IMPOSE_CN=n |
| 457 | # switch to either impose N content of vegetation (to surpass P limitation): |
| 458 | STOMATE_IMPOSE_NP=n |
| 459 | # do not fake the initial height (DOFOCO parameter) |
| 460 | LD_FAKE_HEIGHT=n |
| 461 | |
| 462 | Nammonium_FILE =NONE |
| 463 | Nnitrate_FILE =NONE |
| 464 | Nfert_FILE=NONE |
| 465 | Nbnf_FILE=NONE |
| 466 | NINPUT_UPDATE=-1 |
| 467 | }}} |
| 468 | if you want to debug the model, set printlev to 4 for full debug information in PARAM/run.def: |
| 469 | {{{ |
| 470 | PRINTLEV=4 |
| 471 | }}} |
| 472 | |
| 473 | you can add defaults for the nutrient switches into PARAM/orchidee.default |
| 474 | {{{ |
| 475 | # STOMATE_OK_NCYLE ([FLAG]) : Activate the nitrogen cycle {OK_STOMATE} |
| 476 | STOMATE_OK_NCYCLE = n |
| 477 | |
| 478 | # STOMATE_OK_PCYLE ([FLAG]) : Activate the phosphorus cycle {OK_STOMATE} |
| 479 | STOMATE_OK_PCYCLE = n |
| 480 | |
| 481 | # STOMATE_MASS_CONSERVATION ([FLAG]) : activate mass conservation checks which force stop if violated {OK_STOMATE} |
| 482 | STOMATE_MASS_CONSERVATION = n |
| 483 | |
| 484 | # STOMATE_DSG_DEBUG ([FLAG]) : activate checks for stoichiometry and negative pools {OK_STOMATE} |
| 485 | STOMATE_DSG_DEBUG = n |
| 486 | |
| 487 | # LD_FAKE_HEIGHT ([FLAG]) : fakes height of vegetation (ORCHIDEE-CAN related) {OK_STOMATE} |
| 488 | LD_FAKE_HEIGHT = n |
| 489 | |
| 490 | }}} |
| 491 | |
| 492 | REMARKS: |
| 493 | Do not enable XIOS, it is per default disabled in "PARAM/orchidee.default" |
| 494 | Do not enable river routing in site scale simulation, is per default disabled in "fluxnet.card" |
| 495 | |
| 496 | ==== 3.2.4 adjust the job as usual ==== |
| 497 | modify config.card as usual ([http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation]) |
| 498 | |
| 499 | modify run.def as usual ([http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation]) |
| 500 | plus |
| 501 | if you want to archive, for example in the IMBALANCE-P storage; add to config.card: |
| 502 | {{{ |
| 503 | ARCHIVE=/home/surface3/dgoll |
| 504 | }}} |
| 505 | |
| 506 | modify fluxnet.card as usual (no information): |
| 507 | - if you want to keep the output files of the spinup you need to set the following variables |
| 508 | {{{ |
| 509 | # DEBUG mode for SPINUP |
| 510 | # |
| 511 | # This mode keep all SPINUP directory in ARCHIVE |
| 512 | # If disable, all ARCHIVE is automaticly cleaned. |
| 513 | DEBUG_SPIN=y |
| 514 | # If you don't want to keep old spinup steps, but last one |
| 515 | CONSERVE=y |
| 516 | }}} |
| 517 | |
| 518 | ==== 3.2.5 create job ==== |
| 519 | {{{ |
| 520 | ../../../libIGCM/ins_job |
| 521 | }}} |
| 522 | |
| 523 | ==== 3.2.6 adjust job (optional) ==== |
| 524 | |
| 525 | ==== 3.2.7 start simulations ==== |
| 526 | {{{ |
| 527 | my_prompt>> ksh |
| 528 | my_prompt>> nohup ./Job_ENSEMBLE fluxnet > out_Job 2>&1 & |
| 529 | }}} |
| 530 | ==== 3.2.8 check simulations ==== |
| 531 | to see the decomposition of the spatial domain among the processors check in the run dir "Load_balance_orchidee.dat" |
| 532 | |
| 533 | |
| 534 | on curie |
| 535 | {{{ |
| 536 | ccc_mstat -u userid |
| 537 | }}} |
| 538 | {{{ |
| 539 | ccc_mdel jobid |
| 540 | }}} |
| 541 | |
| 542 | ==== 3.2.8 post processing ==== |
| 543 | ensure the Jobname and the path to the reference simulation is set in fluxnet.card; in case there is no reference set it to NONE |
| 544 | {{{ |
| 545 | # History file of former ORCHIDEE runs (Reference) to compare with the current simulations |
| 546 | reference_base_path=/home/scratch01/dgoll/IGCM_OUT/OL2/DEVT/ensemble |
| 547 | reference_config=DSG04EnsCN |
| 548 | }}} |