Changes between Version 3 and Version 4 of DevelopmentActivities/ORCHIDEE-CNP/howtoUse


Ignore:
Timestamp:
2016-03-16T14:45:55+01:00 (8 years ago)
Author:
dgoll
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/ORCHIDEE-CNP/howtoUse

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