File:Ngrip 1yr data do18 ca 1 2000 ad 4.svg
Original file (SVG file, nominally 1,471 × 601 pixels, file size: 132 KB)
Captions
Captions
Summary[edit]
| Description |
English: NGRIP oxygen isotope 18 data ca. last 2000 years |
||
| Date | |||
| Source | Own work | ||
| Author | Merikanto | ||
| Permission (Reusing this file) |
I, the copyright holder of this work, hereby publish it under the following licenses:
   This file is licensed under the Creative Commons Attribution-Share Alike 4.0 International, 3.0 Unported, 2.5 Generic, 2.0 Generic and 1.0 Generic license.
You may select the license of your choice. |
PDF of data is downloaded from dk ice and climate site, and further processed to .csv file to copy paste to text editor and finally done .Scidavis visualization.
Data is Kaufman et al 2009, 1 yr samples, NGRIP
High resolution δ18O data from Agassiz, DYE-3, NGRIP and Renland used in Kaufman et al., Science, 2009.
https://www.iceandclimate.nbi.ku.dk/data/
https://www.iceandclimate.nbi.ku.dk/data/Kaufman_etal_2009_data_29sep2009.pdf
Citation of data
MAIN REFERENCES The dating of the NGRIP and DYE-3 cores is described in the paper: B.M. Vinther, H.B. Clausen, S.J. Johnsen, S.O. Rasmussen, K.K. Andersen, S.L. Buchardt, D. Dahl-Jensen, I.K. Seierstad, M.-L. Siggaard-Andersen, J.P. Steffensen, A.M. Svensson, J. Olsen, and J. Heinemeier A synchronized dating of three Greenland ice cores throughout the Holocene J. Geophys. Res., vol. 111, D13102, 2006, doi:10.1029/2005JD006921 The dating of the A79 core is described in the paper: B.M. Vinther, H.B. Clausen, D.A. Fisher, R.M. Koerner, S.J. Johnsen, K.K. Andersen, D. Dahl-Jensen, S.O. Rasmussen, J.P. Steffensen, and A.M. Svensson Synchronizing Ice Cores from the Renland and Agassiz Ice Caps to the Greenland Ice Core Chronologyy g g p gy J. Geophys. Res., 113, D08115, 2008, doi:10.1029/2007JD009143. DATA USED FOR FORMING DECADAL MEANS IN: D.S. Kaufman, D.P. Schneider, N.P. McKay, C.M. Ammann, R.S. Bradley, K.R. Briffa, G.H. Miller, B.L. Otto-Bliesner, J.T. Overpeck, B.M. Vinther, Arctic Lakes 2k Project Members Recent warming reverses long-term Arctic cooling. Science 325, 1236-1239, 2009.
Python code
- drawing climate data diagram in python 3
- version 4000.0003
- 04.12.2022
-
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
from scipy import interpolate
from matplotlib.ticker import (MultipleLocator, AutoMinorLocator)
import scipy.signal
import matplotlib.colors as colors
import math as math
datafilename="data1.csv"
captioni="NGRIP delta-O18"
savename="ngrip_ca_20_1970_ad_4.svg"
beginx=50
endx=1950
- beginx=300
- endx=700
- beginx=1800
- endx=1950
beginy=-30
endy=-41
figsizex=18
figsizey=7
yearspan=int(math.fabs(endx-beginx))
dfin0=pd.read_csv(datafilename, sep=";")
- lst1=['Year_AD','NGRIP']
- dfin1 = dfin0[dfin0.columns.intersection(lst1)]
- x0=dfin1['Year_AD']
- y0=dfin1['NGRIP']
s1='Year_AD'
s2='NGRIP'
lst1=[s1, s2]
dfin1 = dfin0[dfin0.columns.intersection(lst1)]
x0=dfin1[s1]
y0=dfin1[s2]
x=np.array(x0)
y=np.array(y0)
- print (y)
- quit(-1)
size0=14
size1=16
size2=18
size3=22
x_sm = np.array(x)
y_sm = np.array(y)
x_smooth = np.linspace(beginx,endx, yearspan*10)
funk1 = interpolate.interp1d(x_sm, y_sm, kind="cubic")
y_smooth = funk1(x_smooth)
y_savgol1 = scipy.signal.savgol_filter(y_smooth,11, 3)
y_savgol2 = scipy.signal.savgol_filter(y_smooth,501, 3)
- a=np.outer(np.ones(2),y_smooth)
a=np.outer(np.ones(2),y_savgol2)
cmap1="RdBu"
cmap1="hsv"
cmap1="seismic"
cmap1="bwr"
clam1="coolwarm1"
clam1="RdBu_r"
print (type(cmap1))
- fig=plt.gca()
- plt.figure(figsize=(figsizex, figsizey))
plt.xlim([beginx, endx])
plt.rcParams["figure.figsize"] = [figsizex, figsizey]
- plt.imshow(a,aspect='auto', cmap=cmap1, extent=[beginx, endx,-50,-30], alpha=0.5)
- plt.imshow(a,aspect='auto', cmap=cmap1, extent=[beginx, endx,endy,beginy], alpha=0.5)
plt.imshow(a,aspect='auto', cmap=clam1, vmin=endy-0, vmax=beginy+0, extent=[beginx, endx,endy,beginy], alpha=1.0)
- plt.xlim(beginx, endx)
- plt.ylim(endy, beginy)
- plt.plot(x, y, color="#7f0000", alpha=1.0, linewidth=1)
plt.plot(x, y, color="#0000ff", alpha=0.7, linewidth=1)
- plt.plot(x_smooth, y_savgol1, color="#0000ff", alpha=0.33, linewidth=1)
plt.plot(x_smooth, y_savgol2, color="#0000ff", linewidth=2)
plt.title(captioni, fontsize=size3, color="#0000af")
plt.xlabel('Year AD', color="darkgreen", fontsize=21)
plt.ylabel('delta-O18', color="darkgreen", fontsize=21)
plt.xticks(fontsize=17)
plt.yticks(fontsize=17)
- plt.gca().invert_xaxis()
- mng = plt.get_current_fig_manager()
- mng.window.state('zoomed')
plt.savefig('figure.svg')
plt.show()
Licensing[edit]
- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
File history
Click on a date/time to view the file as it appeared at that time.
| Date/Time | Thumbnail | Dimensions | User | Comment | |
|---|---|---|---|---|---|
| current | 07:49, 4 December 2022 | | 1,471 × 601 (132 KB) | Merikanto (talk | contribs) | Update |
| 17:54, 26 November 2022 |  | 1,516 × 630 (137 KB) | Merikanto (talk | contribs) | Uploaded own work with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
