Changes between Version 3 and Version 4 of palaeo_transient

Timestamp:

02/19/21 13:28:17 (3 years ago)

Author:

nbouttes

Comment:

--

palaeo_transient

@@ -83 +83 @@
 == Happy transient palaeo simulations! == 
 ** // problems, bugs etc. report to dmr // **
+
+----
+
+
+== Transient runs with CO2 concentration (historic and futur) ==
+This is to run simulations wih forced CO2 concentration in the radiative code (with F_PALEO=2). It uses scenarios from CMIP6, combining historic forcings + ssp126 or ssp585. With this following configuration the orbital parameters also change (although not important for a few hundred years) but the ice sheets do not evolve. [[BR]]
+
+BEWARE: these scenarios have been set to the carbon cycle, need to remove the carbon cycle flags (inlcuding coral) in choixcomposantes.h if without. [[BR]]
+
+'''Pratical:''' [[BR]]
+1) set values in ludus-code/ecbilt/parameters/palaeo_transient_parameter.dat: [[BR]]
+Note that the date is in BP (1850 AD= 100 BP) [[BR]]
+palaeo_year= 100 [[BR]]
+ref_irunlabel=23000 [[BR]]
+accel_rt=1 [[BR]]
+equil_flag=0 [[BR]]
+
+
+2) use scenarios: [[BR]]
+- For SSP126:
+{{{
+-F scenario/ssp126/sources_concentration -I scenario/ssp126/inputdata_concentration
+}}}
+
+- For SSP585:
+{{{
+-F scenario/ssp585/sources_concentration -I scenario/ssp585/inputdata_concentration
+}}}
+
+exemple (with carbon cycle restart):
+{{{
+./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 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 (ghgscen=2501 for ssp)
+- paleo_timer_mod.f90: 
+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 (same as ghgscen and nb_year_ghg) [[BR]]
+
+[[BR]]