369 | | 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]] |
370 | | |
371 | | '''Pratical:''' [[BR]] |
372 | | 1) set values in ludus-code/ecbilt/parameters/palaeo_transient_parameter.dat: [[BR]] |
373 | | Note that the date is in BP and 1850 AD= 100 BP [[BR]] |
374 | | palaeo_year= 100 [[BR]] |
375 | | ref_irunlabel=23000 [[BR]] |
376 | | accel_rt=1 [[BR]] |
377 | | equil_flag=0 [[BR]] |
378 | | |
379 | | |
380 | | 2) run with (using scenario): [[BR]] |
381 | | ./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 |
382 | | |
383 | | |
384 | | '''Notes:''' [[BR]] |
385 | | This runs with carbon cycle flags + coral flag + F_PALAEO=2 [[BR]] |
386 | | |
387 | | Modifications have been made in: [[BR]] |
388 | | |
389 | | - comphys.f90: change ghgscen to the number of lines with values in GHG.dat (ghgscen=2501 for ssp) |
390 | | - paleo_timer_mod.f90: |
391 | | 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]] |
392 | | change nb_year_ghg with number of values in GHG.dat (2501 for ssp, should be the same as ghgscen) ([[BR]] |
393 | | |
394 | | - GHG.dat: [[BR]] |
395 | | starts with more recent value at line 2, and finishes with oldest value at last line [[BR]] |
396 | | value on first line is the number of lines with data (same as ghgscen and nb_year_ghg) [[BR]] |
| 369 | see: |
| 370 | |
| 371 | http://forge.ipsl.jussieu.fr/ludus/wiki/palaeo_transient |
| 372 | |
| 373 | |