When using the data please cite:
Brugger, J.; Hofmann, M.; Petri, S.; Feulner, G. (2019): Model output for the publication: "On the
sensitivity of the Devonian climate to continental configuration, vegetation cover, orbital
configuration, CO_2 concentration and insolation". GFZ Data Services.
http://doi.org/10.5880/PIK.2019.002
The data are supplementary material to:
Brugger, J.; Hofmann, M.; Petri, S.; Feulner, G. (2019): "On the sensitivity of the Devonian climate
to continental configuration, vegetation cover, orbital configuration, CO_2 concentration and
insolation". Paleoceanography and Paleoclimatology.
https://doi.org/10.1029/2019PA003562
92M Jul 18 14:01 Figures.tar.xz
devonian_topography_vegetation.i
function func devn_proc_topo( timeslice=, manualmode=, test1=, test2= )
convert_to_ascii_frlnd_cont.py
writes land fraction for each grid cell of the 3 different
continental configurations in ModelOutput_p2s_frlnd_.….dat
file
readoutfile…
calculates surface air temperature over continents (annual global mean, averaged
over 1000 years)
readoutfileps_tsland | continents.py | continental sensitivity |
vegetation.py |
Early, Middle, Late Devonian vegetation and bare, shrubcover, treecover extreme cases | |
LeHiralbedo.py | treecover and bare land with albedo values from Le Hir 2011 | |
bestguess.py | bestguess.py |
CO2_delta18O_Devonian_corrected.pdf
)
delta18ODevonian_correcteddata.py
CO2_delta18O_Devonian_corrected.pdf
)co2fit.xls
datapointsco2.xlsx
Joachimski18oDevonian.xlsx
AllMaps.pdf
)
convert_to_ascii_…
files write relevant output variables in .dat files
readoutfilecontinents.py
generates maps using the model output written in the .dat files
DevonianVegetation.pdf
)
Schematic vegetation sketches generated using inkscape (svg files)
Vegetation maps generated by devonian_fig3.i
, function devn_plot_veg( timeslice= )
Vegetation_Devoniantimes_plus_extremescenarios_1000years.pdf
)
readoutfile_plot_vegminbare_4x5_1000years.py
reads relevant output variables for
different vegetation runs from model output, based on 1000 years snapshot data (yearly
snapshots) saved in snapshots_potsdam2-1000years-yearmean_as_e_qi.nc
as well as
surface air temperature from 1000 years average in history_p2.nc
, and generates
the figure based on the data in the .dat files
devn_temp_orbital.pdf
)
generated by devonian_fig5.i
, function devn_plot_orbital(season=)
Modes.pdf
)
convert_to_ascii.py
writes relevant output variables in .dat files; the run used has a short
time step (1 hour) and does not split ocean timesteps
readoutfilemodesplot.py
generates the figure (Modes.pdf
) based on the data in the .dat
files
modestransition.pdf
)
plot inside is generated in the same way as Fig.6 (Modes_inplot.pdf
)
4modes.pdf
)
convert_to_ascii.py
writes relevant output variables in .dat files; the run used has a short
time step (1 hour) and does not split ocean timesteps
readoutfilemom.py
generates the figure (4modes.pdf
) based on the data in the .dat
files
GLB_OVER.pdf
)
convert_to_ascii.py
writes global overturning in ModelOutput_momh_24p5_3600s.dat
files;
the run used has a short time step (1 hour) and does not split ocean timesteps
readoutfile_streamfunction.py
generates the figure (GLB_OVER.pdf
) based on the data in
the .dat files
Bestguess_maps_withsst.eps
)
convert_to_ascii_p2history.py
writes annual mean surface air temperature in .dat files
convert_to_ascii_isis.py
writes annual mean ice fraction in .dat files
convert_to_ascii_temp.py
writes ocean surface temperature and surface velocity in .dat
files
readoutfile_withtemp.py
generates the figure (Bestguess_maps_withsst.eps
) based on the
data in the .dat files
devonian_fig10.i
, func devn_plot_sens
SST_Devonian_correcteddata_withfit.pdf
)
tempDevonianwithfit_correcteddata.py
black dots: modeled temperatures, annual average between 30 and 10°S; error bars
indicate temperature range from 30 to 10°S
red and blue crosses: temperature estimates from δ18O data (Joachimski et al. 2009),
using different calibration standards
tempJoachimski2009withfit.xlsx
temperature data from Joachimski et al. 2009
Column 1: age,
2: δ18O,
3: temperature;
other columns are Loess fit with uncertainties which we do not use here
tempJoachimskicorrected.xlsx
temperature data from Joachimski et al. 2009, personal communication, using a different
calibration standard than in 2009 (Lecuyer et al. 2003) and a different temperature
equation (Lecuyer et al. 2013) in column 3, different one (Puc&eacut;at et al. 2010) which is not
shown in our plot in column 4.
Column 1: time
Column 2: δ18O
The LOESS fits shown in the plot are created with Python and therefore differ slightly from
the one in Joachimski et al. 2009
c3beta_devn_380Ma_1500ppm_1319Wm2_O23p5_E0p000_P000_bare
c3beta | model version |
devn | period in Earth history (Devonian) |
380Ma | continental configuration |
1500ppm | CO2 concentration |
1319Wm2 | solar constant |
O23p5 | obliquity in degree |
E0p000 | eccentricity |
P000 | precession in degree |
bare | vegetation cover |
The sensitivity runs for the different variables are in the corresponding folder:
bestguess
co2
continents
flips
orbital
solarconstant
vegetation
c3beta_devn_380Ma_1500ppm_1319Wm2_O23p5_E0p000_P000_bare
is the Middle Devonian standard run which is used for comparison in most of the sensitivity runs.
11G Jul 18 16:23 c3beta_devn_380Ma_1500ppm_1319Wm2_O23p5_E0p000_P000_bare.tar.xz
vegetation
folder:c3beta_devn_380Ma_1500ppm_1319Wm2_O23p5_E0p000_P000_…
different Devonian vegetation distributions:earlydevveg
middledevveg
latedevveg
extreme scenarios: covering all continents with bare land, shrub, trees:
bare
shrubcover
treecover
using the albedo values of Le Hir et al. 2011 for bare land and tree-covered continents:
LeHir_albedo_bare
LeHir_albedo_treecover
flips
folder:c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_timestep0p25
test influence of 6 hour timestep
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_3600s
shorter time step 1 hour, this is the run our evaluation in Chapter 4 is based on
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_ifort15_rconvect_restartmittel3502
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_ifort15_rconvect_restartmittel3602
test of influence of different convection scheme and restart from different times of standard run
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_ifort17_newstart
new start with different compiler: ifort17
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_ifort17_restartgross3402
restart from standard run, changing the compiler at t=3402 on a large fluctuation
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_ifort17_restartklein3202
restart from standard run, changing the compiler at t=3202 on a small fluctuation
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_ifort17_restartmittel4602
restart from standard run, changing the compiler at t=3202 on a medium fluctuation
c3beta_devn_380Ma_1500ppm_1319Wm2_O24p5_E0p069_P315_bare_snapshots
standard run with snapshots saved in a 10-year interval starting after 4000 years to 6000
years
snapshots
contain monthly values of the variables for the model year given in the name
Ocean model data: | snapshots.00….01.01.dta.nc overturn.00….01.01.dta.nc |
Atmosphere model data: | snapshots_potsdam2.00….01.01.dta.nc |
Ice model data: | snapshots_isis.00….01.01.dta.nc |
snapshots_potsdam2-1000years-yearmean_….nc
(atmosphere) andsnapshots-1000years-yearmean_….nc
(ocean)
history
files contain yearly values of the variables for each model year the simulation was run
Ocean model data: | history.nc |
Atmosphere model data: | history_p2.nc |
Ice model data: | history_isis.nc |
topog.dta.nc
contains information about topography and cells
Continent_maps_diff.pdf
)
convert_to_ascii_p2history.py file writes tsann in ModelOutput_p2h_.….dat
files for each
continental configuration
readoutfilediff.py
generates surface air temperature difference maps using the model
output written in the .dat files
Continents_maps_sstuv
)
convert_to_ascii_temp.py
file writes relevant output variables in ModelOutput_momh_…_continents.dat
files for each continental configuration
readoutfile_withtemp.py
generates maps using the model output written in the .dat files
Solar_maps_diff.pdf
)
convert_to_ascii_p2history.py
file writes tsann in ModelOutput_p2h_.….dat
files for each
solar constant
readoutfilediff.py
generates surface air temperature difference maps using the model
output written in the .dat files
CO2_maps_diff.pdf
)
convert_to_ascii_p2history.py
file writes tsann in ModelOutput_p2h_.….dat
files for each
CO2 concentration
readoutfilediff.py
generates surface air temperature difference maps using the model
output written in the .dat files
Vegetation_Devoniantimes_clouds_1000years.pdf
)
readoutfileps_plot_vegminbare_4x5_1000years.py
uses model output for total cloud
cover (1000 years from yearly snapshots) to generate difference maps of total cloud cover
for different vegetation cover
Albedo_maps_diff.pdf
)
convert_to_ascii_p2history.py
file writes tsann in ModelOutput_p2h_.….dat
files for each
albedo used constant
readoutfilediff.py
generates surface air temperature difference maps using the model
output written in the .dat files
eminp_maps_diff.pdf
)
convert_to_ascii_p2history.py
file writes annual evaportion and precipitation in
ModelOutput_p2h_.dat
file
readoutfilediff.py
generates difference map using the model output written in the .dat files