Changes between Version 66 and Version 67 of CarbonCyclePage

Timestamp:

02/12/21 11:47:45 (3 years ago)

Author:

nbouttes

Comment:

--

CarbonCyclePage

@@ -367 +367 @@
 
 '''3.9. Transient CO2 concentration historic and futur (with corals)'''[[BR]]
+This is to run simulations wih forced CO2 concentration in the radiative code. It uses scenarios from CMIP6, combining historic forcings + ssp126 or ssp585. It is currenlty set with carbon cycle and corals, but if the corresponding flags are changed in choixcomposantes.h it should still work. [[BR]]
 
 '''Pratical:''' [[BR]]
-1) set date in ludus-code/ecbilt/parameters/palaeo_transient_parameter.dat [[BR]]
-For 1850 AD: 100 BP [[BR]]
-
-2) run  with: [[BR]]
-./bin/run-iloveclim -l PI_ssp585_concentration_2 -n 2 -r 2 -k -z -F scenario/ssp585/sources_concentration -I scenario/ssp585/inputdata_concentration -s 23000 -S /home/clim01/nbouttes/iloveclim_coastal/data/PI-CC_r1309_2_ecbilt_clio/ic023000 -v 1
+1) set values in ludus-code/ecbilt/parameters/palaeo_transient_parameter.dat: [[BR]]
+Note that the date is in BP and 1850 AD= 100 BP [[BR]]
+palaeo_year= 100 [[BR]]
+ref_irunlabel=23000 [[BR]]
+accel_rt=1 [[BR]]
+equil_flag=0 [[BR]]
+
+
+2) run  with (using scenario): [[BR]]
+./bin/run-iloveclim -l PI_ssp585_concentration_2 -n 2 -r 2 -k -F scenario/ssp585/sources_concentration -I scenario/ssp585/inputdata_concentration -s 23000 -S /home/clim01/nbouttes/iloveclim_coastal/data/PI-CC_r1309_2_ecbilt_clio/ic023000 -v 1
 
 
 '''Notes:''' [[BR]]
-This runs with carbon cycle flages + coral flag + F_PALAEO=2 [[BR]]
+This runs with carbon cycle flags + coral flag + F_PALAEO=2 [[BR]]
 
 Modifications have been made in: [[BR]]
 
-- comphys.f90: change ghgscen to the number of lines with values in GHG.dat
+- comphys.f90: change ghgscen to the number of lines with values in GHG.dat (ghgscen=2501 for ssp)
 - paleo_timer_mod.f90: 
-change start_year_ghg with more recent (or future) value at top of GHG.dat [[BR]]
-change nb_year_ghg with oldest value at last line of GHG.dat [[BR]]
+change start_year_ghg with more recent (or future) value at top of GHG.dat (start_year_ghg=-550 for ssp; -550 BP = 2500 AD, by definition 0 BP=1950 AD)) [[BR]]
+change nb_year_ghg with number of values in GHG.dat (2501 for ssp, should be the same as ghgscen) ([[BR]]
 
 - GHG.dat: [[BR]]
 starts with more recent value at line 2, and finishes with oldest value at last line [[BR]]
-value on first line is the number of lines with data [[BR]]
+value on first line is the number of lines with data (same as ghgscen and nb_year_ghg) [[BR]]
 
 [[BR]]