File:Changes of temperature of world 1 600 ad 01 01.png

2. 1. plot LMR data
3. https://atmos.washington.edu/~hakim/lmr/LMRv2/
4. Last Millennium Reanalysis (LMR) Project Data Page - version 2
   1. LMR v2.1
   2. python xarray
6. 2023.05.13 0000.0000

import xarray as xr import pandas as pd import matplotlib.pyplot as plt import matplotlib.colors as colors from matplotlib.ticker import (MultipleLocator, AutoMinorLocator) import matplotlib.ticker as tck import datetime import numpy as np from scipy import interpolate import scipy.signal import math as math

captioni="Mean temperature changes of world between 1 - 600 AD"

1. subtitle="deg C"

file_path = './data/gmt_MCruns_ensemble_full_LMRv2.1.nc'

1. file_path = './data/nhmt_MCruns_ensemble_full_LMRv2.1.nc'
2. file_path = './data/shmt_MCruns_ensemble_full_LMRv2.1.nc'

ds = xr.open_dataset(file_path)

k1=ds.gmt[:]

k2=np.array(k1)

1. y=k2[:,1,1]

1. y=np.mean(axis=1)

y=np.mean(k2, axis=(1, 2))

print(np.shape(y))

1. quit(-1)

1. for iy in range(100):
2. for ix in range(20):
3. y2=k2[:,ix,iy]
4. plt.plot(y2,color="#000070", alpha=0.01, linewidth=1 )

yys1=[] yys2=[]

for iy in range(100): yys1.append(np.mean(k2, axis=2))

for iy in range(20): yys2.append(np.mean(k2, axis=1))

print (np.shape(yys1[1]))

1. quit()

1. for iy in range(20):
2. plt.plot(yys2[iy], color="#0000ff", alpha=0.01, linewidth=1 )

for iy in range(100): plt.plot(yys1[iy], color="#8080ff", alpha=0.005, linewidth=1 )

plt.plot(y, color="#7f0000", alpha=0.8, linewidth=2)

1. plt.grid(True)

plt.xlim(00,600)

1. plt.ylim(-0.5, 0.2)
   1. plt.ylim(-1, 1)

plt.ylim(-0.6, 0.5)

1. plt.axhline(y = ylimit1, color = 'g', linestyle = ':', linewidth=2)

1. plt.locator_params(axis='x', nbins=10)
2. plt.axis.xaxis.set_minor_locator(MultipleLocator(5))
3. plt.axes().xaxis.set_minor_locator(MultipleLocator(5))
4. plt.axes().yaxis.set_minor_locator(tck.AutoMinorLocator())

plt.gca().xaxis.set_major_locator(plt.MultipleLocator(100)) plt.gca().xaxis.set_minor_locator(plt.MultipleLocator(20))

1. plt.gca().yaxis.set_major_locator(plt.MultipleLocator(100))
2. plt.gca().yaxis.set_minor_locator(plt.MultipleLocator(10))

plt.grid(which='major', linestyle=':', linewidth='0.5', color='blue', alpha=.5)

1. Customize the minor grid

plt.grid(which='minor', linestyle=':', linewidth='0.5', color='green', alpha=.5)

plt.title(captioni, fontsize=18, color="#0000af")

1. plt.suptitle(subtitle,fontsize=14)

plt.xlabel('Year AD', color="darkgreen", fontsize=21) plt.ylabel("Delta T deg Celsius", color="darkgreen", fontsize=21) plt.xticks(fontsize=17) plt.yticks(fontsize=17)

1. plt.gca().invert_yaxis()

1. mng = plt.get_current_fig_manager()
2. mng.window.state('zoomed')

plt.savefig('figure.svg')

plt.show()

quit(-1)