wiki:CarbonCyclePage

Version 17 (modified by nbouttes, 6 years ago) (diff)

--

Description of the Carbon Cycle procedures to use the iLOVECLIM carbon cycle version


I. Preindustrial

To run a pre-industrial (PI) simulation, not specific to the carbon cycle:

  • Modify the file: lbm/sources/choixcomposantes.h to have:

IMSK 1
COMATM 1
ROUTEAU 1 (previous version 0)
EVAPTRS 1
EVAPSI 1 (previous version 0)
F_PALAEO 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!!


To start a (preindustrial) run with the carbon cycle (PI-CC):

1) In general the run is started from the pre-indutrial control simulation without carbon, it is thus a restart simulation.

In ocycc/sources/loveclim_transfer_mod.f90 : change KLSR to 1 to have a restart (KLRS=0 for a fresh start)

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:

IMSK 1
COMATM 1
ROUTEAU 1
EVAPTRS 1
EVAPSI 1
F_PALAEO 0

CYCC 2
OCYCC 1
OLDC14 0
KC14 1
KC14P 0
OXNITREUX 0
WINDINCC 0
O2ATM 0
N2OATM 0
MEDUSA 0
PATH 0

LGMSWITCH 0

4) Initial conditions in data/default (cf. "README.txt").

To start from the current PI-CC standard conditions, add: -s 65800 and -S /home/acclimate/smoreira/_PROJS_/004_iLOVECLIM_MED/20171211_iloveclim_parallel/trunk/wkdir/op02_NOMED_5ky_ecbilt_clio/run/restartdata/res065800

The previous initial state was: data/default/ic50800

II. LGM


To start a LGM run (LGM):

LGM Method 1 with choixcomposantes.h:

  • In lbm/sources/choixcomposantes.h:

F_PALAEO 1
LGMSWITCH 1

  • Set the palaeo_year to 21000 in ludus-code/ecbilt/parameters/palaeo_transient_parameters.dat
  • Start run with something like:

./bin/loveclim-new-run -l lgm_nath -s 10000 -n 3000 -r 1000 –P scenario/lgm-pmip2-reb/parameters -I scenario/lgm-pmip2-reb/inputdata –S data/restart_lgm_Nath/ic010000

Note: at the moment restart from PI run instead of the given (old) LGM simulation

LGM Method 2:

  • change the GHG in ghg.dat (ludus-code/ecbilt/inputdata/GHG.dat)
  • change the orbital parameters in namelist file

and do as for a PI restart run


To start a LGM run with the carbon cycle (LGM-CC):

./bin/loveclim-new-run -l mon_joli_lgm_CC -s 65800 -n 1000 -r 250 –P scenario/lgm-pmip3-cc1/parameters -I scenario/lgm-pmip3-cc1/inputdata -S /home/acclimate/smoreira/_PROJS_/004_iLOVECLIM_MED/20171211_iloveclim_parallel/trunk/wkdir/op02_NOMED_5ky_ecbilt_clio/run/restartdata/res065800

To start a LGM run with the carbon cycle and Pa-Th (LGM-CC-PaTh?):

  • In lbm/sources/choixcomposantes.h:

PATH 1

  • Set restart_path = .TRUE. in iloveclim/ocycc/sources/path_mod.f90
  • Check that the four Pa-Th outputs are set to T in ludus-code/clio/parameters/netcdfout.param-scal17

./bin/loveclim-new-run -l testLGM_CC_PaTh -s 67800 -n 1000 -r 250 –P scenario/lgm-pmip2-reb-CC/parameters -I scenario/lgm-pmip2-reb-CC/inputdata -S /home/acclimate/lmissi/iloveclim_lr/wkdir/preind_eq_2ky_v1_15Mar18_ecbilt_clio/run/restartdata/res067800


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


To force the ice caps topography in LGM run using GRISLI:

  • In lbm/sources/choixcomposantes.h:

F_PALAEO 2
ISM 3
LGMSWITCH 1

  • Check the choice of reference topography in /grisli/SOURCES/Fichiers-parametres/hemin40_param_list.dat
  • Check the path in iloveclim/config/config.obelix
  • To impose the LGM bathymetry, fractoc.dat, bath.om and mozaic.w are needed in inputdata/ directory
  • GRISLI needs to know the LGM bathymetry as well. Create a LGM.dat file in /grisli/INPUT/Forcage/ containing it

Beware! The option ISM 3 is not compatible with cold conditions in the model for now - unstable


To fix the vegetation to the PI in emic.param set VECODE option to F instead of T

III. 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