Yhteenveto

Source of data is

ftp://ftp.ncdc.noaa.gov/pub/data/paleo/speleothem/asia/india/mawmluh2012.txt

^[1]

^[2]

Python src to produce image In python 3.8

^[3]

Lisenssi

Minä, tämän teoksen tekijänoikeuksien haltija, julkaisen täten tämän teoksen seuraavalla lisenssillä:

Tämä tiedosto on lisensoitu Creative Commons Nimeä-JaaSamoin 4.0 Kansainvälinen -lisenssillä.

Voit:

• jakaa – kopioida, levittää ja esittää teosta
• remiksata – valmistaa muutettuja teoksia

Seuraavilla ehdoilla:

• nimeäminen – Sinun on mainittava lähde asianmukaisesti, tarjottava linkki lisenssiin sekä merkittävä, mikäli olet tehnyt muutoksia. Voit tehdä yllä olevan millä tahansa kohtuullisella tavalla, mutta et siten, että annat ymmärtää lisenssinantajan suosittelevan sinua tai teoksen käyttöäsi.
• jaa samoin – Jos muutat tai perustat tähän työhön, voit jakaa tuloksena syntyvää työtä vain tällä tai tämän kaltaisella lisenssillä.

https://creativecommons.org/licenses/by-sa/4.0CC BY-SA 4.0 Creative Commons Attribution-Share Alike 4.0 truetrue

1. ↑ Berkelhammer, M., et al. 2012. Mawmluh Cave, India Holocene Stalagmite d18O Data. IGBP PAGES/World Data Center for Paleoclimatology Data Contribution Series # 2012-059. NOAA/NCDC Paleoclimatology Program, Boulder CO, USA. ORIGINAL REFERENCE: Berkelhammer, M., A. Sinha, L. Stott, H. Cheng, F. Pausata, and K. Yoshimura. In press 2012. An Abrupt shift in the Indian Monsoon 4,000 years ago. Geophysical Monograph, American Geophysical Union. doi: 10.1029/2012GM001207.
2. ↑ Year_BP;delta_18O_Calcite_Permil_VPDB 3653;-5.98978 3659;-6.34279 3665;-6.371348 3671;-6.156792 3677;-6.107148 3683;-6.138709 3689;-6.265227 3695;-6.420271 3701;-6.45645 3707;-6.474284 3713;-6.600483 3719;-6.737324 3725;-6.996647 3731;-7.168338 3737;-6.738413 3743;-6.388649 3749;-6.366517 3755;-6.401711 3761;-6.605557 3767;-6.723481 3773;-6.637361 3779;-6.626879 3785;-6.775418 3791;-6.808547 3797;-6.606922 3803;-6.506546 3809;-6.625383 3815;-6.782594 3821;-7.038703 3827;-7.190117 3833;-6.974262 3839;-6.781875 3845;-6.726361 3851;-6.672892 3857;-6.653369 3863;-6.633846 3869;-6.663875 3875;-6.696148 3881;-6.680064 3887;-6.655423 3893;-6.168229 3899;-5.517952 3905;-4.867675 3911;-4.627194 3917;-4.593027 3923;-4.704356 3929;-4.723486 3935;-5.45048 3941;-5.786952 3947;-4.835108 3953;-4.558823 3959;-4.756806 3965;-5.083706 3971;-4.68298 3977;-4.614502 3983;-4.710526 3989;-4.850247 3995;-4.693171 4001;-4.707682 4007;-4.70741 4013;-4.707863 4019;-4.842658 4025;-4.886334 4031;-4.65784 4037;-4.456344 4043;-4.484125 4049;-4.511906 4055;-4.897943 4061;-5.302656 4067;-5.692039 4073;-5.890646 4079;-5.67964 4085;-5.468633 4091;-5.669297 4097;-5.867621 4103;-5.775672 4109;-5.428576 4115;-5.333582 4121;-5.709503 4127;-6.113525 4133;-6.010202 4139;-5.441222 4145;-5.495727 4151;-5.702946 4157;-6.183804 4163;-6.113916 4169;-5.838163 4175;-5.725463 4181;-5.909567 4187;-5.849016 4193;-5.712925 4199;-5.68069 4205;-5.706087 4211;-5.178929 4217;-5.584249 4223;-5.599146 4229;-5.747592 4235;-5.871376 4241;-5.728776 4247;-5.649844 4253;-5.571714 4259;-5.586808 4265;-5.73337 4271;-5.879932 4277;-5.786213 4283;-5.685634 4289;-5.663336 4295;-5.755051 4301;-5.603889 4307;-5.346635 4313;-6.236417 4319;-6.443348 4325;-6.402271 4331;-6.588203 4337;-6.548169 4343;-6.212441 4349;-6.397603 4355;-6.566164 4361;-6.400625 4367;-6.235085 4373;-6.069546 4379;-5.849665 4385;-5.676389 4391;-5.625474 4397;-5.585594 4403;-5.59135
3. ↑
   1. 2. mawmluh cave delta-o18 4200-1800 delta-oi8 csv "spectra drawing"
      3. test v "ks8.py"
   import csv import matplotlib import matplotlib.pyplot as plt from matplotlib.colors import LinearSegmentedColormap import matplotlib.ticker as mtick from matplotlib.ticker import FormatStrFormatter from matplotlib import cm
   1. from colorspacious import cspace_converter
   from collections import OrderedDict import matplotlib.colors as colors import matplotlib.cbook as cbook import numpy as np from scipy.interpolate import interp1d fname1='mawmluh3.csv' beg1=3800 end1=4200 step1=1 SMALL_SIZE = 14 MEDIUM_SIZE = 16 BIGGER_SIZE = 21 numbero=0 numbero2=0 stakkastep1=100 stakkan=[] stakkan2=[] xxx=[] yyy=[] sxx=[] syy=[] ixx=[] iyy=[] sixx=[] siyy=[] pixx=[] piyy=[] plt.rc('font', size=BIGGER_SIZE) # controls default text sizes plt.rc('axes', titlesize=MEDIUM_SIZE) # fontsize of the axes title plt.rc('axes', labelsize=MEDIUM_SIZE) # fontsize of the x and y labels plt.rc('xtick', labelsize=MEDIUM_SIZE) # fontsize of the tick labels plt.rc('ytick', labelsize=SMALL_SIZE) # fontsize of the tick labels plt.rc('legend', fontsize=MEDIUM_SIZE) # legend fontsize plt.rc('figure', titlesize=BIGGER_SIZE) # fontsize of the figure title def load_interpolate(fname1, beg1, end1, step1): global numbero global numbero2 global stakkastep1 global stakkan global stakkan2 global xxx global yyy global sxx global syy global ixx global iyy global sixx global siyy nn=0 with open(fname1, newline=) as csvfile: spamreader = csv.reader(csvfile, delimiter=';', quotechar='|') for row in spamreader: if(nn>0): xx=float(row[0]) yy=float(row[1]) yyy.append(yy) xxx.append(xx) numbero=numbero+1 nn=1 minx=1e8 miny=1e8 maxx=-1e8 maxy=-1e8 for n in range(0, numbero): xx=xxx[n] yy=yyy[n] if xx < minx: minx=xx if yy < miny: miny=yy if xx > maxx: maxx=xx if yy > maxy: maxy=yy #print(yyy) #print (miny) deltax=maxx-minx deltay=maxy-miny for n in range(0, numbero): xx=xxx[n] yy=yyy[n] sx=(xx-minx)/deltax sy=(yy-miny)/deltay syy.append(sy) sxx.append(sx) f1 = interp1d(xxx, yyy, kind='cubic') numnero2=0 for n in range(beg1, end1,step1): xx=float(n) yya=f1(xx) yy=float(yya) iyy.append(yy) ixx.append(xx) numbero2=numbero2+1 numbero2=len(ixx) #print(ixx) #print(iyy) minx2=1e8 miny2=1e8 maxx2=-1e8 maxy2=-1e8 for n in range(0, numbero2): xx=ixx[n] yy=iyy[n] if xx < minx2: minx2=xx if yy < miny2: miny2=yy if xx > maxx2: maxx2=xx if yy > maxy2: maxy2=yy dim2=end1-beg1 dx2=dim2/step1 #for n in range(beg1, end1,step1): for n in range(0, numbero2): yy=iyy[n] xx=ixx[n] sx=(xx-beg1)/dim2 sy=(yy-miny2)/deltay siyy.append(sy) sixx.append(sx) #print (siyy) for n in range(0, numbero2,stakkastep1): xx=ixx[n] stakkan2.append(int(xx)) xx=xx+stakkastep1 stakkan2.append(int(xx)) return(0) def wanha_my_colordict(y): green = [(0.0, 0.0, 0.0), (0.5, 0, 0), (1.0, 0, 0.0)] red = [(0.0, 0.0, 0.0), (0.5, y, y), (1.0, y, y)] blue = [(0.0, 1, 1), (0.5, 1, 1),(1.0,0.0,0.0)] colordict = dict(red=red, green=green, blue=blue) bluegreenmap = LinearSegmentedColormap('bluegreen', colordict, 10000) return bluegreenmap def my_colordict(y): green = [(0.0, 1.0, 1.0), (0.5, 1-y, 1-y), (1.0, 0.0, 0.0)] red = [(0.0, 0.0, 0.0), (0.5, y, y), (1.0, y, y)] blue = [(0.0, 1, 1), (0.5, 1-y, 1-y),(1.0,0.0,0.0)] colordict = dict(red=red, green=green, blue=blue) bluegreenmap = LinearSegmentedColormap('bluegreen', colordict, 10000) return bluegreenmap def plot_test_one(): global numbero global numbero2 global stakkastep1 global stakkan global stakkan2 global xxx global yyy global sxx global syy global ixx global iyy global sixx global siyy global pixx global piyy piaxx=[] piayy=[] sryy=[] for n in range(0, numbero2): xx=sixx[n] yy=siyy[n] px=(xx)*1200/3 py=(1-yy)*400/3 piaxx.append(px) piayy.append(py) pixx=piaxx piyy=piayy[::-1] sryy=siyy[::-1] ta=xxx sa=yyy t=xxx t2=sixx s=syy s2=siyy ## interp sr=s[::-1] sr2=s2[::-1] stakkanr=stakkan[::-1] stakkan2r=stakkan2[::-1] #print ("stakkan") print (stakkan2r) fig = plt.figure(figsize=(12,3)) ax = fig.add_subplot(111) #fig, ax = plt.subplots() a=np.outer(np.ones(2),sr2) fact = 1.0/255.0 my_cmap2=my_colordict(1) #ax.set(xlabel='Age (BP)', ylabel='(blue moist/ red dry)', title='Mawmluh Cave d-O18') plt.xlabel('Age (yr BP)') plt.ylabel('d-O18') plt.title('Mawmluh Cave') plt.yticks([]) plt.xticks(np.arange(0, numbero2, step=100), stakkan2r) plt.imshow(a,aspect='auto', cmap ='jet') #ax.patch.set_facecolor('white') ax.patch.set_alpha(0.9) ax.plot(pixx, sryy, linewidth=5, color='darkblue',alpha=0.7) plt.savefig("../Downloads/haista.svg", transparent=True) plt.show() return load_interpolate(fname1, beg1, end1, step1)
   1. print ("pppp")
   2. print(stakkan2)
   plot_test_one()