Version 7 (modified by roche, 6 years ago) (diff) |
---|
Description of the Carbon Cycle procedures to use the iLOVECLIM carbon cycle version
Not specific to the carbon cycle:
- Modify the file: lbm/sources/choixcomposantes.h to have:
IMSK 1
COMATM 1
ROUTEAU 0
EVAPTRS 1
EVAPSI 0
- To start a run:
nohup ./bin/loveclim-new-run -l LH -n 500 -r 250 > run_LH &
starts a simulation called LH which lasts 500 years and writes restart files every 250 years, starting from the default initial conditions.
- To see more options:
./bin/loveclim-new-run -h
- To modify the ocean outputs: in clio/parameters/netcdfout.param
Beware'' for the carbon cycle: output set to 100 per file![[BR]]
To start a run with the carbon cycle:
1) In ocycc/sources/loveclim_transfer_mod.f90 : change KLSR to 1 to have a restart
2) In veccarb/sources/carbone_co2.f90: new_run_c should be set to 0 (if run without c13 values for the vegetation-> set to 1)
3) In lbm/sources/choixcomposantes.h:
CYCC 2
OCYCC 1
KC14 1
4) Initial conditions in data/default (cf. "README.txt"). To start from the initial conditions, add: -s 50800 and -S data/default/ic50800
To start an LGM run:
4) Set gammat to 0 in clio/sources/thersf.f (Suppression of fluxes from ice-shelves, not coherent with LGM topography in Antarctica)
5) Change initial conditions
Start run with something like:
For LGM:
./bin/loveclim-new-run -l lgm_nath -s 10000 -n 3000 -r 1000 –P scenario/lgm-pmip3-v01/parameters -I scenario/lgm-pmip3-v01/inputdata –S data/restart_lgm_Nath/ic010000
For LGM with carbon cycle:
./bin/loveclim-new-run -l mon_joli_lgm_CC -s 17000 -n 1000 -r 250 –P scenario/lgm-pmip3-cc1/parameters -I scenario/lgm-pmip3-cc1/inputdata –S data/default/lgm_CC_v01/res017000
Chercher dans ocycc/sources/init_mb.f ET ocycc/sources/initmbiopar.f le drapeau LGMSWITCH : le définir à 1 pour le LGM (depuis la révision 220) (sinon 0 ?)
To simulate the sinking of brines:
- Set frac to a value between 0 (no brines) and 1 (maximum effect) in brines.f.
- Set use_brines to 1 in clio/sources/clio.f
The different atmospheric CO2 concentrations
- PGACO2 = Radiative pCO2 for the atmosphere. In emic.f :
! Radiative: IF(.NOT.lradCO2) PGACO2=PCO2ref
- PCO2ref = constant to pre-industrial pCO2 levels ...
in emic.f :
PCO2ref=277.4D0
- patmCO2 is used to set the CO2 seen by the vegetation. In emic.f:
! Fertilization: IF(lferCO2) then patmCO2=PGACO2 ELSE patmCO2=PCO2ref ENDIF ! lferCO2 !
- And used in VECODE as: (ccparam.f)
nppt=npp*(1.0+(betat*LOG(co2ghg/280.))) nppg=npp*(1.0+(betag*LOG(co2ghg/280.)))
- PA_C is the current CO2 in the carbon cycle for the atmosphere
- PA0_C is the PA_C at restart
- Link between CO2s: PGACO2 is reset in atmphy0.f from ghg(1), read in GHG.dat input file