File:Changes of temperature of world 1000 1500 ad 01 01.png

###########################################
#
## plot LMR data 
# https://atmos.washington.edu/~hakim/lmr/LMRv2/
# Last Millennium Reanalysis (LMR) Project Data Page - version 2
## LMR v2.1
## python xarray
##
#
# 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 1000 - 1500 AD"
#subtitle="deg C"

file_path = './data/gmt_MCruns_ensemble_full_LMRv2.1.nc'
#file_path = './data/nhmt_MCruns_ensemble_full_LMRv2.1.nc'
#file_path = './data/shmt_MCruns_ensemble_full_LMRv2.1.nc'

ds = xr.open_dataset(file_path)

k1=ds.gmt[:]

k2=np.array(k1)

#y=k2[:,1,1]

#y=np.mean(axis=1)

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

print(np.shape(y))

#quit(-1)

#for iy in range(100):
#	for ix in range(20):
#		y2=k2[:,ix,iy]
#		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]))

#quit()

#for iy in range(20):
#	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)
#plt.grid(True)
 
plt.xlim(1000,1500)
#plt.ylim(-0.5, 0.2)
##plt.ylim(-1, 1)
plt.ylim(-0.5, 0.2)

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

#plt.locator_params(axis='x', nbins=10) 
#plt.axis.xaxis.set_minor_locator(MultipleLocator(5))
#plt.axes().xaxis.set_minor_locator(MultipleLocator(5))
#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))

#plt.gca().yaxis.set_major_locator(plt.MultipleLocator(100))
#plt.gca().yaxis.set_minor_locator(plt.MultipleLocator(10))

plt.grid(which='major', linestyle=':', linewidth='0.5', color='blue', alpha=.5)
# Customize the minor grid
plt.grid(which='minor', linestyle=':', linewidth='0.5', color='green', alpha=.5)

plt.title(captioni, fontsize=18, color="#0000af")
#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)

#plt.gca().invert_yaxis()

#mng = plt.get_current_fig_manager()
#mng.window.state('zoomed')

plt.savefig('figure.svg')

plt.show()

quit(-1)