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 500 - 1000 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(500,1000)

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

plt.ylim(-0.6, 0.4)

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)