| 1 | = ENSEMBLE: fluxnet sites [http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation] = |
| 2 | [[Image(fluxnet_sequence.png)]] |
| 3 | Figure shows the sequence of simulations: the labels(TOP), the time length(MIDDLE) and the spinup_analytic(BOTTOM). Adopted from NV scheme. |
| 4 | |
| 5 | |
| 6 | == 1.1 create new experiment folder == |
| 7 | {{{ |
| 8 | cd ../../config/ORCHIDEE_OL/ |
| 9 | cp -fr ENSEMBLE/ ENSEMBLE_CNP/ |
| 10 | cd ENSEMBLE_CNP/ |
| 11 | }}} |
| 12 | |
| 13 | == 1.2 tell libIGCM about the new input files == |
| 14 | For now, I store all files needed on the IMBALANCE-P disk on obelix (you find the path below). Please ensure you are using the newest revision (rev00 is not fully tested). |
| 15 | |
| 16 | The configuration which are used are stored in |
| 17 | ../SPINUP/SUBJOB/OOL_SEC_STO/COMP/ |
| 18 | so, do as in 4.1.2: change the path of soil_param.nc and add the new files |
| 19 | |
| 20 | sechiba.card: |
| 21 | {{{ |
| 22 | (/home/surface3/common/CNP_files/rev00/soils_param.nc, .), \ |
| 23 | }}} |
| 24 | stomate.card: |
| 25 | {{{ |
| 26 | ListNonDel= (${R_IN}/SRF/reftemp.nc, .) \ |
| 27 | (/home/surface3/common/CNP_files/rev00/USDA_SoilSuborder.nc, .), \ |
| 28 | (/home/surface3/common/CNP_files/rev00/lithology.nc, .), \ |
| 29 | (/home/surface3/common/CNP_files/rev00/deposition.nc, .) |
| 30 | }}} |
| 31 | 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. |
| 32 | |
| 33 | exchange the spinup.driver file in ../SPINUP/COMP/ with the new modified driver from NV. |
| 34 | {{{ |
| 35 | cp /home/surface3/common/CNP_files/rev00/spinup.driver ../SPINUP/COMP/ |
| 36 | }}} |
| 37 | |
| 38 | |
| 39 | copy a run.def with all parameters needed for the ORCHIDEE-CAN parts of the model to PARAM, for example this one: |
| 40 | {{{ |
| 41 | cp /home/surface3/common/CNP_files/rev00/run.def_NV2 PARAM/run.def |
| 42 | }}} |
| 43 | |
| 44 | == 1.3 adjust the libIGCM == |
| 45 | modify the configuration of the SPINUP. These variables were originally used for the outdated spinup method using FORCESOIL. |
| 46 | Thus the variable names are misleading. |
| 47 | |
| 48 | modify them in ensemble.card and fluxnet.card: |
| 49 | {{{ |
| 50 | [SPINUP] |
| 51 | # Initialisation for spin-up : |
| 52 | # orchidee with sechiba alone (!!! if ok_stomate == n !!!) |
| 53 | duree_nostomate=0 |
| 54 | # normal SEC_STO simulation to fill biomass pools (equilibrate litter fluxes) of duree_inistomate years |
| 55 | duree_inistomate=1000 |
| 56 | # teststomate (only if duree_inistomate > 0); DSG: don't touch this one |
| 57 | duree_offlineini=0 |
| 58 | |
| 59 | # Loop over ORCHIDEE runs (used for spin-up) |
| 60 | # The analytical spinup is restarted n_iter times ... |
| 61 | n_iter=3 |
| 62 | # ... for duree_sechiba years |
| 63 | duree_sechiba=50 |
| 64 | # teststomate; DSG: don't touch this one |
| 65 | duree_stomate=0 |
| 66 | # forcesoil; DSG: don't touch this one |
| 67 | duree_carbonsol=1 |
| 68 | |
| 69 | # Final run after the analytical spinup |
| 70 | # This last parameter must be non-zero. |
| 71 | duree_final=100 |
| 72 | |
| 73 | |
| 74 | }}} |
| 75 | |
| 76 | make sure that the variable SPINUP_ANALYTIC is set to "n": |
| 77 | {{{ |
| 78 | SPINUP_ANALYTIC=n |
| 79 | }}} |
| 80 | add the variables to turn on/off the nutrients cycles to PARAM/run.def: |
| 81 | {{{ |
| 82 | # switch for N cycle |
| 83 | STOMATE_OK_NCYCLE=y |
| 84 | # switch for P cycle |
| 85 | STOMATE_OK_PCYCLE=y |
| 86 | # switch to either impose N content of vegetation (to surpass N limitation): |
| 87 | STOMATE_IMPOSE_CN=n |
| 88 | # switch to either impose N content of vegetation (to surpass P limitation): |
| 89 | STOMATE_IMPOSE_NP=n |
| 90 | # do not fake the initial height (DOFOCO parameter) |
| 91 | LD_FAKE_HEIGHT=n |
| 92 | |
| 93 | Nammonium_FILE =NONE |
| 94 | Nnitrate_FILE =NONE |
| 95 | Nfert_FILE=NONE |
| 96 | Nbnf_FILE=NONE |
| 97 | NINPUT_UPDATE=-1 |
| 98 | }}} |
| 99 | if you want to debug the model, set printlev to 4 for full debug information in PARAM/run.def: |
| 100 | {{{ |
| 101 | PRINTLEV=4 |
| 102 | }}} |
| 103 | |
| 104 | you can add defaults for the nutrient switches into PARAM/orchidee.default |
| 105 | {{{ |
| 106 | # STOMATE_OK_NCYLE ([FLAG]) : Activate the nitrogen cycle {OK_STOMATE} |
| 107 | STOMATE_OK_NCYCLE = n |
| 108 | |
| 109 | # STOMATE_OK_PCYLE ([FLAG]) : Activate the phosphorus cycle {OK_STOMATE} |
| 110 | STOMATE_OK_PCYCLE = n |
| 111 | |
| 112 | # STOMATE_MASS_CONSERVATION ([FLAG]) : activate mass conservation checks which force stop if violated {OK_STOMATE} |
| 113 | STOMATE_MASS_CONSERVATION = n |
| 114 | |
| 115 | # STOMATE_DSG_DEBUG ([FLAG]) : activate checks for stoichiometry and negative pools {OK_STOMATE} |
| 116 | STOMATE_DSG_DEBUG = n |
| 117 | |
| 118 | # LD_FAKE_HEIGHT ([FLAG]) : fakes height of vegetation (ORCHIDEE-CAN related) {OK_STOMATE} |
| 119 | LD_FAKE_HEIGHT = n |
| 120 | |
| 121 | }}} |
| 122 | |
| 123 | REMARKS: |
| 124 | Do not enable XIOS, it is per default disabled in "PARAM/orchidee.default". |
| 125 | Do not enable river routing in site scale simulation, is per default disabled in "fluxnet.card" |
| 126 | |
| 127 | == 1.4 adjust the job as usual to your specific needs == |
| 128 | 1. modify config.card following [http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation] |
| 129 | give information on Job Name, Tag Name, Dates and Restart. Make sure the path ARCHIVE is correct for the machine you are using. You must set the ARCHIVE on curie or your simulation will crash. |
| 130 | for example in the IMBALANCE-P storage; add to config.card: |
| 131 | {{{ |
| 132 | ARCHIVE=/home/surface3/dgoll |
| 133 | }}} |
| 134 | 2. modify run.def following [http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation] |
| 135 | |
| 136 | 3. modify fluxnet.card following [http://forge.ipsl.jussieu.fr/orchidee/wiki/Scripts/FluxnetValidation] |
| 137 | |
| 138 | In case you want to check if the spinup works well, increase output level of the spinup simulations & prevent libIGCM from removing these files. |
| 139 | The files can be found in the respective site folders in the experiment folder. |
| 140 | |
| 141 | modify in fluxnet.card and ensemble.card: |
| 142 | {{{ |
| 143 | # level of the history files for each iteration of the SPINUP |
| 144 | level_hist=10 |
| 145 | # level of the history files for the final iteration of the SPINUP |
| 146 | level_hist_final=10 |
| 147 | |
| 148 | # DEBUG mode for SPINUP |
| 149 | # |
| 150 | # This mode keep all SPINUP directory in ARCHIVE |
| 151 | # If disable, all ARCHIVE is automaticly cleaned. |
| 152 | DEBUG_SPIN=y |
| 153 | # If you don't want to keep old spinup steps, but last one |
| 154 | CONSERVE=y |
| 155 | }}} |
| 156 | You find the files in the folder "JobName"STOI in the IGCM_OUT path for the STOI simulation. Exchange "STOI" with "ORC-1", "FIN" etc for the simulations you want to check. |
| 157 | |
| 158 | == 1.5 create job == |
| 159 | {{{ |
| 160 | ../../../libIGCM/ins_job |
| 161 | }}} |
| 162 | |
| 163 | == 1.6 adjust job == |
| 164 | '''on obelix''' |
| 165 | change the queue to longp to increase the allocated time (the default queue mediump is too short and you will exceed wall clock time). |
| 166 | {{{ |
| 167 | #PBS -q longp |
| 168 | }}} |
| 169 | |
| 170 | If you want to debug on obelix: set the variable RUN_DIR_PATH to disk with enough space (the default path is not accessible for the user) |
| 171 | {{{ |
| 172 | # for example on scratch01 |
| 173 | RUN_DIR_PATH=/home/scratch01/dgoll/RUN_DIR |
| 174 | }}} |
| 175 | ACHTUNG: On curie, do not set RUN_DIR_PATH in the Job_*. If you do, the configuration will crash. You find the default path to the RUN_DIR_PATH in the section "3.2.8 check simulations". |
| 176 | |
| 177 | |
| 178 | == 1.7 start simulations == |
| 179 | {{{ |
| 180 | ksh |
| 181 | nohup ./Job_ENSEMBLE fluxnet > out_Job 2>&1 & |
| 182 | }}} |
| 183 | |
| 184 | == 1.8 check simulations == |
| 185 | to see the decomposition of the spatial domain among the processors check in the run dir "Load_balance_orchidee.dat" |
| 186 | '''on curie''' |
| 187 | The RUN_DIR is in |
| 188 | {{{ |
| 189 | /ccc/scratch/cont003/dsm/<userid>/RUN_DIR |
| 190 | }}} |
| 191 | |
| 192 | to check the status of the job; type |
| 193 | {{{ |
| 194 | ccc_mstat -u <userid> |
| 195 | }}} |
| 196 | to delete a job; type |
| 197 | {{{ |
| 198 | ccc_mdel <jobid> |
| 199 | }}} |
| 200 | === 1.8.1 what to do when a simulation crash === |
| 201 | You cannot use clean_month and clean_year. You should remove the sites that crashed from the IGCM_OUT path |
| 202 | and rerun the same experiment with a fluxnet.card containing only the sites which crashed. |
| 203 | |
| 204 | |
| 205 | == 1.9 post processing == |
| 206 | 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 |
| 207 | {{{ |
| 208 | # History file of former ORCHIDEE runs (Reference) to compare with the current simulations |
| 209 | reference_base_path=/home/scratch01/dgoll/IGCM_OUT/OL2/DEVT/ensemble |
| 210 | reference_config=DSG04EnsCN |
| 211 | }}} |
| 212 | start the skript |
| 213 | {{{ |
| 214 | ksh |
| 215 | nohup ./Job_Post_FLUXNET fluxnet > out_Post 2>&1 & |
| 216 | }}} |
| 217 | you will find the graphics in the experiment folder in the IGCM_OUT path. There is a pdf in the folder starting with "fluxnet_taylor_diff". |