1. 3. drawing climate diagram in python 3
   4. version 2.11
   5. 11.9.2020

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

def running_mean(x, N):

   cumsum = np.cumsum(np.insert(x, 0, 0)) 
   return (cumsum[N:] - cumsum[:-N]) / float(N)

datafilename="ngrip1.csv" captioni="Bölling-Alleröd period in NGRIP ice core" savename="ngrip_dryas.svg"

figsizex=16 figsizey=8

1. x = []
2. y = []
3. y2= []

dfin0=pd.read_csv(datafilename, sep=";")

lst1=['gicc05_age','delta_O18']

dfin1 = dfin0[dfin0.columns.intersection(lst1)]

x0=dfin1['gicc05_age'] y0=dfin1['delta_O18']

1. y20=dfin1['GISP_dO18']
2. y30=dfin1['GISP2_dO18']

x=np.array(x0) y=np.array(y0)

size0=14 size1=16 size2=18 size3=24

1. y_savgol = scipy.signal.savgol_filter(y,31, 3)

y_savgol = scipy.signal.savgol_filter(y,71, 3)

1. y_running = running_mean(y, 31)

x_sm = np.array(x) y_sm = np.array(y) x_smooth = np.linspace(x_sm.min(), x_sm.max(), 20000) funk1 = interpolate.interp1d(x_sm, y_sm, kind="cubic") y_smooth = funk1(x_smooth)

fig, ax1 = plt.subplots()

1. ax1.axis((11600,14000,0,ymax1))

ax1.set_xlim(12000,16000) ax1.set_ylim(-33.0, -46.0)

1. ax1.set_ylim(-35.0, -42.0)

plt.gca().invert_xaxis() plt.gca().invert_yaxis()

ax1.set_ylabel('delta-O18', color='#0000ff', fontsize=size2+2)

ax1.plot(x,y, color="#b0b0ff", linewidth=1,label="NGRIP delta-O18")

1. ax1.plot(x_smooth,y_smooth, color="#0000ff", linewidth=3,label="NGRIP delta-O18")

ax1.plot(x,y_savgol, color="#0000FF", linewidth=4, label="SavGol filter, 71 and 3")

1. ax1.plot(x,y_running, color="#FF0000", linewidth=3)

1. ax1.plot(x,data_avg1, color="#ff0000", linewidth=2, linestyle=":", label="Average of NGRIP, GISP, GISP2 delta-O18")

ax1.tick_params(axis='both', which='major', labelsize=size2)

ax1.xaxis.set_minor_locator(MultipleLocator(1000)) ax1.xaxis.set_minor_locator(MultipleLocator(200))

ax1.yaxis.set_minor_locator(MultipleLocator(1.0)) ax1.yaxis.set_minor_locator(MultipleLocator(0.1))

ax1.grid(which='major', linestyle='-', linewidth='0.1', color='black') ax1.grid(which='minor', linestyle=':', linewidth='0.1', color='black')

ax1.set_xlabel('Age BP', color="darkgreen", fontsize=size2)

ax1.set_title(captioni, fontsize=size3, color="#0000af")

plt.legend(fontsize=size0)

fig = plt.gcf() fig.set_size_inches(figsizex, figsizey, forward=True)

plt.savefig(savename, format="svg", dpi = 100)

plt.show()