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 matplotlib.ticker import (MultipleLocator, AutoMinorLocator)

from scipy import interpolate import scipy.signal

from scipy.interpolate import CubicSpline

def running_mean(x, N):

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

minx1=8000 maxx1=8400 miny1=-33.4 maxy1=-36.2

savsize1=5 # savgol filtering buffer size

datafilename1="ngrip_1_ngrip_2_grip_holocene1.csv" captioni="8.2 ka event on NGRIP" savename="neem_ch4_yd_pb_1.svg"

1. age_b2k;dO18_ngrip1;dO18_ngrip2;dO18_grip

column01="age_b2k" column02="dO18_ngrip1"

figsizex=16 figsizey=8

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

dfin00=pd.read_csv(datafilename1, sep=";")

lst01=[column01, column02]

dfin01 = dfin00[dfin00.columns.intersection(lst01)]

x0=dfin00[column01] y0=dfin01[column02]

x=x0[::2] y=y0[::2]

dfin01[column02].rolling(11).mean()

1. x=x0
2. y=y0

1. points = zip(x0, y0)

1. points = sorted(points, key=lambda point: point[0])
2. x, y = zip(*points)

1. x=np.array(x).astype(float)
2. y=np.array(y).astype(float)

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

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

1. print(x)

x_sm = x y_sm = y

1. print (x_sm)

x_smooth = np.linspace(x_sm.min(), x_sm.max(), 2000) x_smooth= np.linspace(minx1,maxx1, 2000)

funk1 = interpolate.interp1d(x_sm, y_sm, kind="cubic") y_smooth = funk1(x_smooth)

xs = np.arange(minx1,maxx1,1)

1. lenu=len(x)
2. xn=np.arange(1,lenu,1)

1. cs = CubicSpline(xn, y)

fig, ax1 = plt.subplots()

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

ax1.set_xlim(minx1,maxx1) ax1.set_ylim(miny1, maxy1)

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

plt.gca().invert_xaxis()

plt.gca().invert_yaxis()

ax1.set_ylabel('CH4 (ppbv)', color='#0000ff', fontsize=size2+2)

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

ax1.plot(x_smooth,y_smooth, color="#0000ff", linewidth=4,label="NGRIP delta-O18") ax1.plot(x,y_savgol, color="#5fA05f", linewidth=3, label="SavGol filter,"+str(71)+" and 3")

1. ax1.plot(x,cs(xs), color="#007f00", linewidth=4,label="Smoothed")

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

ax1.xaxis.set_minor_locator(MultipleLocator(100)) ax1.xaxis.set_minor_locator(MultipleLocator(20))

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)

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

plt.show()