diff --git a/Python code/Difference between myEels2.f90 and MyEels3.f90 .txt b/Python code/Difference between myEels2.f90 and MyEels3.f90 .txt new file mode 100644 index 0000000000000000000000000000000000000000..1cd301a7a2dde5798157f2a0668638a49aa83fb3 --- /dev/null +++ b/Python code/Difference between myEels2.f90 and MyEels3.f90 .txt @@ -0,0 +1,28 @@ +% Date: 2024-04-18 15:37:15 - 3 differences +% /Users/michael/Developer/EELS2 git-svn/EELS2/Python code/myEels2.f90 +% /Users/michael/Developer/EELS2 git-svn/EELS2/Python code/myEels3.f90 +% Use up and down arrows to move through differences, ⌃⌥-left and ⌃⌥-right to patch changes from one window to the other. +978c978,979 +< if (osc(1,j) > 0.) then ! positive TO mode: 'Kurosa' form +--- +> +> if (osc(1,j) > 0.) then ! positive TO mode: 'Kurosawa' form: _Multiplicative_ phonon mode + +982,987c983,988 +< else ! Negative TO mode means: treat as term added to epsilon +< if (osc(3,j) > 0) then ! it is a Lorentz oscillator +< addeps = addeps + osc(1,j)**2 * osc(2,j) /dcmplx(osc(1,j)**2 + wn2, wn*osc(3,j)) +< else ! it is a Drude term +< addeps = addeps - osc(1,j)**2/dcmplx(wn2, -1*wn*osc(3,j)) +< end if +--- +> +> else if (osc(1,j) < 0.) then! Negative TO mode means: _Additive_ Lorentz oscillator with Q +> addeps = addeps + osc(1,j)**2 * osc(2,j) /dcmplx(osc(1,j)**2 + wn2, wn*osc(3,j)) +> +> else ! osc(1,j) = 0 -> it is a Drude term +> addeps = addeps - dcmplx(osc(1,j+m)**2, wn*(osc(3,j)-osc(3,j+m))) /dcmplx(wn2, wn*osc(3,j)) + +988a990 +> + diff --git a/Python code/Difference between myEels3.f90 and myEels4.f90 .txt b/Python code/Difference between myEels3.f90 and myEels4.f90 .txt new file mode 100644 index 0000000000000000000000000000000000000000..cb7a99c04f97a5f05794b8b443c78d1597f0ed9e --- /dev/null +++ b/Python code/Difference between myEels3.f90 and myEels4.f90 .txt @@ -0,0 +1,9 @@ +% Date: 2024-04-18 15:39:59 - 1 difference +% /Users/michael/Developer/EELS2 git-svn/EELS2/Python code/myEels3.f90 +% /Users/michael/Developer/EELS2 git-svn/EELS2/Python code/myEels4.f90 +% Use up and down arrows to move through differences, ⌃⌥-left and ⌃⌥-right to patch changes from one window to the other. +993c993 +< eps(l) = epsinf(l) * nomi / deno + addeps +--- +> eps(l) = epsinf(l) * (nomi / deno + addeps) ! brackets changed by HHe 230915 + diff --git a/Python code/HREELS.py b/Python code/HREELS.py index 75093c662bac69e6f7dab6bdba21b052e8794de1..6663755c9a75a67baf0705c2b1824c98de6ab64c 100755 --- a/Python code/HREELS.py +++ b/Python code/HREELS.py @@ -13,13 +13,6 @@ from matplotlib.transforms import offset_copy # from pylab import figure from matplotlib.pyplot import figure -if os.name == 'nt': - diskd = 'D:\\' -elif os.name == 'posix': - diskd = '/mnt/d/' -else: - print('Define diskd! System not known ...') - def myPath(path): '''Convert any Windows or Linux path to a Linux path. Adds a slash at the end, if missing.''' if path == None: @@ -481,7 +474,7 @@ class HREELS: def myMain(): fig = None - datapath = diskd+ "/Data/Python/HREELS/expHREELS/data" + datapath = "data/gph_dat/" data1 = HREELS("I2L02",datapath) if data1.valid: data1.info() @@ -499,4 +492,4 @@ def myMain(): if __name__ == '__main__': - myMain() \ No newline at end of file + myMain() diff --git a/Python code/calcHREELS.py b/Python code/calcHREELS kms?.py similarity index 100% rename from Python code/calcHREELS.py rename to Python code/calcHREELS kms?.py diff --git a/Python code/calcHREELS20.py b/Python code/calcHREELS20.py new file mode 100755 index 0000000000000000000000000000000000000000..eccf213666d30a3b3dcec08890c6acb79e484a1e --- /dev/null +++ b/Python code/calcHREELS20.py @@ -0,0 +1,301 @@ +#!/usr/bin/env python3 +import numpy as np +import json +import sys, re, os +from libhreels.HREELS import myPath +from copy import deepcopy + +import scipy.integrate as integrate + +libDir = os.path.dirname(os.path.realpath(__file__)) + +try: + from libhreels import myEels20 as LambinEELS # wrapper for myEels20.f90 + from libhreels import myBoson as LambinBoson # wrapper for myBoson.f90 +except: + print('myEels20 and MyBoson are not available here (Check your version)') + print('''Make sure the Fortran routines 'myEels20' and 'myBoson' + have been complied with the proper f2py for the right python version!!''') + print('\n\n\n') + +# Experimental setup as dictionary: +setup = { + "e0": 4.0, + "theta": 60., + "phia": 0.33, + "phib": 2.0, + "temperature": 298., + "debug": False +} +# Instrumental function describing elastic peak shape: +instrument = { + "width": 18., + "intensity": 100000., + "asym": 0.01, + "gauss": 0.88 +} + + +def importMaterials(string='', path=libDir): + ''' Returns a dictionary with all phonon parameters for the material provided + as string argument. If the string is empty or not matching, a list of all available + materials is printed. + ''' + file = os.path.join(myPath(path),'materials20.json') + with open(file) as json_file: + materials = json.load(json_file) + try: + mat = materials[string] + + #Following if-case for preventing old material data to be used in wrong terminology + if mat["wLO"][0]==1 and mat["gLO"][0]==1 and mat["wTO"][0]<0 and mat["gTO"][0]<0: + print('It seems you are trying to load a Material from an older version. Parameter will be altered to fit the current Version.') + wTO = 0 + gTO = -mat["gTO"][0] + wLO = -mat["wTO"][0] + gLO = -mat["gTO"][0] + mat["wTO"][0] = wTO + mat["gTO"][0] = gTO + mat["wLO"][0] = wLO + mat["gLO"][0] = gLO + except: + print('No data for material >>{}<< found in {} materials.json!!'.format(string, path)) + print('Available materials:\n{}\n'.format(materials.keys())) + mat = 'None' + return mat + +def addDrude(wLOPlasma, gLOPlasma, material, gTOPlasma='None'): + ''' Adds a generalized Drude response to the materials properties (which are provided + as last argument) and returns a new materials dictionary with all phonon parameters. Note + that at least the eps_infinity has to given before. + ''' + if gTOPlasma == 'none': + gTOPlasma = gLOPlasma + newMaterial = deepcopy(material) + try: + if len(newMaterial['wTO']) > 0: + newMaterial['wTO'] += [0.] + newMaterial['gTO'] += [gTOPlasma] + newMaterial['wLO'] += [wLOPlasma] + newMaterial['gLO'] += [gLOPlasma] + if newMaterial.get('Q'): + newMaterial['Q'] += [0] + return newMaterial + except: + print('Cannot add Drude to material',material) + return material + +################################################################################ +################################################################################ +class lambin: + def __init__(self, film, setup=setup, instrument=instrument): + self.e0 = setup['e0'] + self.theta = setup['theta'] + self.phia = setup['phia'] + self.phib = setup['phib'] + self.temperature = setup['temperature'] + self.debug = setup['debug'] + self.width = instrument['width'] + self.gauss = instrument['gauss'] + self.intensity = instrument['intensity'] + self.asym = instrument['asym'] + self.layers = len(film) # number of layers + self.neps = self.layers + # name_size = self.layers + self.name = []; self.thick=[]; self.listNOsci=[]; self.epsinf =[]; Q = [] + allTO=[]; allgTO=[]; allgLO=[]; nDrude=0; Qdummy = [] + name2 = [] + for layer in film: + try: + a = layer[0]['name'] + except: + a = 'None' + self.name.append('{:<10}'.format(a[:10])) # film name and material + name2.append(a) + try: + a = layer[1] + except: + a = 10000. + self.thick.append(a) + self.epsinf.append(layer[0]['eps']) + nTO = 2 * len(layer[0]['wTO']) + allTO.extend(layer[0]['wTO']) + allgTO.extend(layer[0]['gTO']) + allTO.extend(layer[0]['wLO']) + allgTO.extend(layer[0]['gLO']) + qList = layer[0].get('Q') + if qList: + Q.extend(layer[0]['Q']) + Q.extend(len(layer[0]['Q'])*[0.]) + else: + Q.extend(2* len(layer[0]['wTO'])*[0.]) + self.listNOsci.append(nTO) + + if len(allTO)!=sum(self.listNOsci) or len(allgTO)!=sum(self.listNOsci): + print('Error in materials: ', layer[0]) + if len(Q)!=sum(self.listNOsci) : + print('Error in materials (Check Q): ', layer[0]) + self.wOsc = np.array(allTO) + self.gOsc = np.array(allgTO) + self.osc = np.array([self.wOsc, np.array(Q), self.gOsc]) + # print('[self.wOsc, np.array(Q), self.gOsc]: \n',self.osc) + return + + def calcSurfaceLoss(self,x): + ''' Calculate the surface loss spectrum for the array of x, which needs to be an equidistant array. + All parameters are defined in the class __init__() call.''' + wmin = min(x) + wmax = max(x)-0.001 + dw = (wmax-wmin)/(len(x)-1) # assumes that x is an equidistant array + wn_array_size = len(x) # size of array for x and epsilon (wn_array, loss_array) + nper = 1. + contrl = '{:<10}'.format('None'[:10]) # Can be 'image' to include image charge + mode = '{:<10}'.format('kurosawa'[:10]) + wn_array,loss_array = LambinEELS.mod_doeels.doeels(self.e0,self.theta,self.phia,self.phib, + wmin,wmax,dw,self.layers,self.neps,nper,self.name, + self.thick,self.epsinf,self.listNOsci,self.osc,contrl,mode,wn_array_size) + i=0 + for item in wn_array: + if item > 0: break + i += 1 + return wn_array[i-1:], loss_array[i-1:] + + def calcHREELS(self,x, normalized=True, areanormalized=False): + emin = min(x) + emax = max(x)-0.001 + norm = 1 + xLoss,loss_array = self.calcSurfaceLoss(x) + wmin = min(xLoss) + wmax = max(xLoss) + xOut,spectrum,n = LambinBoson.doboson3(self.temperature,self.width,self.gauss,self.asym, + emin,emax,wmin,wmax,loss_array,self.debug,len(loss_array)) + if normalized: + norm = max(spectrum[:n]) + if areanormalized: #edit by HHE + try: + areanormalize_xstart = np.argmin(abs(x+100.)) #seems to be oddly complicated, but is way more stable than x.index(-100.) or where() + except: + areanormalize_xstart = 0 + try: + areanormalize_xend = np.argmin(abs(x-1000.)) + except: + areanormalize_xend = len(x) + cropped_spectra=spectrum[areanormalize_xstart:areanormalize_xend] + cropped_x=x[areanormalize_xstart:areanormalize_xend] + + norm=integrate.simps(cropped_spectra, dx=x[areanormalize_xstart+1]-x[areanormalize_xstart]) + + # else: + # print("not normalized") + return xOut[:len(x)], spectrum[:len(x)]/norm + + def calcEps(self, x): + epsArray = [] + nOsci = len(self.wOsc) + for wn in x: + yn = LambinEELS.mod_doeels.seteps(self.listNOsci,nOsci,self.osc,self.epsinf,wn,self.layers) + epsArray.append(yn) + return np.transpose(np.array(epsArray)) + +#################################################################################### +def myMain(): + import matplotlib.pyplot as plt + import numpy as np + import os + from libhreels import dielectrics20 as dielectrics + + x = np.linspace(-100.,1000,400) + material = {'eps': 1., + 'wTO': [-200, -750], + 'gTO': [ 12, 8], + 'wLO': [ 1, 1], # this parameter is irrelvant if wTO is negativ + 'gLO': [ 1, 1], # this parameter is irrelvant if wTO is negativ + 'Q' : [ 10, 15]} + material2 = {'eps': 1., + 'wTO': [-200, -750], + 'gTO': [ 12, 8], + 'wLO': [ 1, 1], # this parameter is irrelvant if wTO is negativ + 'gLO': [ 1, 1], # this parameter is irrelvant if wTO is negativ + 'Q' : [ 0.1, 15]} + material3 = {'eps': 1., + 'wTO': [-200, -750, 0], + 'gTO': [ 12, 8, 50], + 'wLO': [ 1, 1, 2000], # this parameter is irrelvant if wTO is negativ + 'gLO': [ 1, 1, 50], # this parameter is irrelvant if wTO is negativ + 'Q' : [ 10, 15, 0]} + + film1 = lambin(film=[[material,10000.]]) + film2 = lambin(film=[[material2,10000.]]) + film3 = lambin(film=[[material3,10000.]]) + eps3= film3.calcEps(x)[0] + eps2= film2.calcEps(x)[0] + eps = film1.calcEps(x)[0] + plt.plot(x,np.real(eps)) + plt.plot(x,np.real(eps2)) + plt.plot(x,np.real(eps3)) + + ############# Comparison with dielectrics ################ + osci1 = dielectrics.simpleOscillator(material3['wTO'][0], + material3['Q' ][0], + gTO=material3['gTO'][0]) + osci2 = dielectrics.simpleOscillator(material3['wTO'][1], + material3['Q' ][1], + gTO=material3['gTO'][1]) + drude = dielectrics.drude(material3['wLO'][2],material3['gTO'][2],material3['gLO'][2]) + epsInfinity = material3['eps'] + eps_dielectrics = epsInfinity * (osci1(x) + osci2(x) + (1 + drude(x))) + plt.plot(x,np.real(eps_dielectrics),linestyle='dotted') + + + + plt.ylabel(r'$Re(\epsilon)$') + plt.xlabel('Energy Loss (cm$^{-1}$)') + plt.xlim(left=5) + # plt.ylim(-6000,6000) + + plt.text(0.99, 0.01,os.path.basename(__file__), fontsize=10, ha='right', va='bottom', transform=plt.gcf().transFigure) + output_filename = os.path.splitext(__file__)[0] + '.png' + plt.savefig(output_filename) + + plt.show() + + plt.plot(x,np.imag(eps)) + plt.plot(x,np.imag(eps2)) + plt.plot(x,np.imag(eps3)) + plt.plot(x,np.imag(eps_dielectrics)) + plt.xlim(left=0) + plt.ylim(-1500,1500) + plt.show() + + # plt.plot(x,np.imag(dielectrics.sigma(eps,x))) + # plt.plot(x,np.imag(dielectrics.sigma(eps2,x))) + # plt.plot(x,np.imag(dielectrics.sigma(eps3,x))) + # plt.xlim(left=0) + # plt.ylim(-1500,1500) + # plt.show() + + xs, spectrum = film3.calcHREELS(x,normalized=True,areanormalized=False) + plt.plot(xs[:-1],spectrum[:-1], label='normalized=False') + plt.show() + + material4 = {'eps': 1., + 'wTO': [ 200, 750, 0], + 'gTO': [ 12, 8, 50], + 'wLO': [ 600, 950, 2000], + 'gLO': [ 10, 10, 50], + 'Q' : [ 10, 15, 0]} + + # material4 = {'eps': 1., + # 'wTO': [ 200, 750, 0], + # 'gTO': [ 12, 8, 50], + # 'wLO': [ 300, 950, 2000], + # 'gLO': [ 10, 10, 50]} + + + film4 = lambin(film=[[material4,10000.]]) + xs, spectrum = film4.calcHREELS(x,normalized=True,areanormalized=False) + plt.plot(xs[:-1],spectrum[:-1], label='normalized=False') + plt.show() + +if __name__ == '__main__': + myMain() \ No newline at end of file diff --git a/Python code/data/gph_dat/I2L02.GPH b/Python code/data/gph_dat/I2L02.GPH new file mode 100755 index 0000000000000000000000000000000000000000..bb6c7d31458ad2a61347fb9da04a3ab9b14a5f40 --- /dev/null +++ b/Python code/data/gph_dat/I2L02.GPH @@ -0,0 +1,2623 @@ +%%%%% 001 # General data + 6.10 2.00 DELTA FFFFFDE0 # P-Version, C-Version, Type +22.02.2019 11:25 # Date, Time +%%%%% 010 # Comment + + + + + + + + + + + + + + + + + + + + +%%%%% 018 # +0 # 50V Scan state +%%%%% 002 # Modul data +Kat. Ver. 2.3399999142 2.3399999142 0.0000000000 5.0000000000 0.0000760000 80 0 0 1 0 +PK -2.1926999092 -2.1926999092 -10.9141559601 11.1900310516 0.0003372831 21 0 0 1 FFFFFFBC +Rep -1.2324999571 -1.2324999571 -10.9234685898 11.0292186737 0.0003349714 22 0 0 1 FFFFFF78 +A1 links 18.0060005188 18.0060005188 -5.9583468437 66.3768157959 0.0011037469 24 0 0 1 FFFFFF34 +A1 rechts 22.5550003052 22.5550003052 -6.0046720505 66.3250579834 0.0011036640 26 0 0 1 FFFFFEF0 +A2 oben -8.0895996094 -8.0895996094 -10.9080619812 11.0086555481 0.0003344226 27 0 0 1 FFFFFEAC +A2 unten -7.1054000854 -7.1054000854 -10.9495010376 11.0446872711 0.0003356047 28 0 0 1 FFFFFE68 +A3 links 1.4486000538 1.4486000538 -10.9238748550 11.0282192230 0.0003349624 29 0 0 1 FFFFFE24 +A3 rechts 0.8749700189 0.8749700189 -10.9453125000 11.0509996414 0.0003356371 2A 0 0 1 FFFFFDE0 +EVM 0.1369300038 0.1369300038 -2.5876874924 4.6144781113 0.0001098963 31 0 0 1 FFFFFD9C +VM innen 0.6869300008 0.6869300008 -2.5837626457 4.6172842979 0.0001098793 32 0 0 1 FFFFFD58 +VM au�en -0.8130699992 -0.8130699992 -2.1036343575 0.1207878143 0.0000339420 33 0 0 1 FFFFFD14 +VM D1 -0.7231900096 -0.7231900096 -2.1016812325 0.1209521890 0.0000339147 34 0 0 1 FFFFFCD0 +VM D2 -0.7231900096 -0.7231900096 -2.0959563255 0.1164074987 0.0000337580 35 0 0 1 FFFFFC8C +HME -0.7525577545 -0.7525577545 -2.5943999290 4.6251878738 0.0001101622 36 0 0 1 FFFFFC48 +HM innen -0.5025477409 -0.5025477409 -2.5880718231 4.6075630188 0.0001097967 37 0 0 1 FFFFFC04 +HM au�en -1.0025677681 -1.0025677681 -2.1031875610 0.1154775023 0.0000338541 38 0 0 1 FFFFFBC0 +HM D1 -1.0006676912 -1.0006676912 -2.1032655239 0.1165978163 0.0000338724 39 0 0 1 FFFFFB7C +HM D2 -1.0000277758 -1.0000277758 -2.1007342339 0.1186871827 0.0000338657 3A 0 0 1 FFFFFB38 +HMA -0.7525476813 -0.7525476813 -2.5927872658 4.6134500504 0.0001099585 3B 0 0 1 FFFFFAF4 +B1 oben -0.8422180414 -0.8423435688 -16.3437500000 16.5563125610 0.0001255038 2B 1 0 1 FFFFFAB0 +B1 unten -0.8301697969 -0.8300442696 -16.3627204895 16.6008110046 0.0001257459 2C 1 0 1 FFFFFA6C +B2 1.1802383661 1.1821393967 -5.1524972916 119.4303131104 0.0019009828 2E 0 0 1 FFFFFA28 +B3 0.9902987480 0.9712957144 -5.3829998970 119.1553115845 0.0019003039 44 0 0 1 FFFFF9E4 +B4 oben -0.8142320514 -0.8154810667 -16.2590332031 16.4820003510 0.0001248971 45 1 0 1 FFFFF9A0 +B4 unten -0.8162304759 -0.8149814606 -16.2824687958 16.4784374237 0.0001249729 46 1 0 1 FFFFF95C +AM -0.7508257627 -0.7508257627 -2.6005969048 4.6468062401 0.0001105866 47 0 0 1 FFFFF918 +A innen -0.5008258820 -0.5008258820 -2.5957968235 4.6316623688 0.0001102823 48 0 0 1 FFFFF8D4 +A au�en -1.0008256435 -1.0008256435 -2.1033282280 0.1158340648 0.0000338617 49 0 0 1 FFFFF890 +A D1 -0.9750257730 -0.9750257730 -2.1015281677 0.1148478091 0.0000338192 4A 0 0 1 FFFFF84C +A D2 -1.0065255165 -1.0065255165 -2.0990469456 0.1100259423 0.0000337078 4B 0 0 1 FFFFF808 +Cht Eing. -5.0223627090 -5.0223627090 -5.2474431992 5.3180408478 0.0001612165 51 0 0 1 FFFFF7C4 +C1 oben 0.5900969505 0.5900969505 -5.2440004349 5.3125343323 0.0001610799 52 0 0 1 FFFFF780 +C1 unten 0.5900969505 0.5900969505 -5.2433624268 5.3155498505 0.0001611162 53 0 0 1 FFFFF73C +C2 2.2117373943 2.2117373943 -5.2518939972 5.3111190796 0.0001611788 54 0 0 1 FFFFF6F8 +C3 links 2.2332372665 2.2332372665 -5.2459716797 5.3085346222 0.0001610490 55 0 0 1 FFFFF6B4 +C3 rechts 5.2665371895 5.2665371895 -5.2543592453 5.3018937111 0.0001610756 56 0 0 1 FFFFF670 +Cht Geh -15.0000000000 -15.0000000000 -16.4459056854 16.6480636597 0.0005049739 61 0 0 1 FFFFF62C +E0 4.0000000000 4.0000000000 -9.2576847076 259.5778198242 0.0041021043 64 0 0 1 FFFFF5E8 +Probe 7.5999999046 7.5999999046 -9.4702816010 258.8974914551 0.0040949672 66 0 0 1 FFFFF5A4 +Scan fine 0.0000000000 0.0000000000 -0.5422697067 1.2692217827 0.0000276412 4C 0 0 1 FFFFF560 +Scan coarse 0.0000000000 0.0000000000 -1.3911311626 51.0064392090 0.0007995235 68 0 0 1 FFFFF51C +50V Rel. 0.2325034440 0.2325034440 0.2325034440 1.9298999310 0.0000259002 6A 0 0 0 FFFFF4D8 +%%%%% 003 # Segment +0 # Flag f�r Segmentauswahl +2 # max. number of scantypes + 0.0000000000 -0.0099508204 0.2232853472 0.0001934882 0.0000000000 # dummy min max inc dummy + 0.0000000000 -0.0127923759 0.0927447230 0.0015990470 0.0000000000 # dummy min max inc dummy +2 # Segmentanzahl + 0.1000000015 -0.0061991666 0.0681908354 0.0001934882 0.5000000000 # wait start stop inc meas + 0.1000000015 0.1983733326 0.2727633417 0.0001934882 0.5000000000 # wait start stop inc meas + 0.2000000030 -0.0619916692 0.1859750003 0.0015990470 2.0000000000 # wait start stop inc meas + 0.2000000030 -0.0619916692 0.1859750003 0.0063961879 2.0000000000 # wait start stop inc meas +%%%%% 004 # Channel data + 0 # main channel +CH--1 V 0.0000000000 1.0000000000 1249404 1 1 FFFFC7C0 +CH--2 V 0.0000000000 33.3300018311 0 0 1 FFFFDE40 +CH--3 V 0.0000000000 1.0000000000 0 0 1 FFFFCB80 +CH--4 V 0.0000000000 1.0000000000 0 0 1 FFFFF4C0 +CH--5 V 0.0000000000 1.0000000000 0 0 1 FFFFE200 +%%%%% 005 # Averagedata +D:\DATEN\HREELS\BATCH\SAV\I2L02 # Filename Measure Values +0031 # Flags +0018 # Anzahl +000 # Start extension +006 # Save extension +D:\DATEN\HREELS\BATCH\SAV\I2L02v # Filename Average Values +D:\DATEN\HREELS\BATCH\SAV\I2L02a # Filename Export Measure Values +D:\DATEN\HREELS\BATCH\SAV\I2L02av # Filename Export Average Values +%%%%% 006 # Options + 8065.6000976562 # xFactor + 0.0000000000 # xOffset +1000000.0000000000 # warning level +%%%%% 008 # Timedata + 0.5000000000 # Meas time + 2.0000000000 # startTime + 0.1000000015 # sleepTime +%%%%% 009 # Graphicoptionen +0.000000 2200.000000# X-Zoom +0.000000 10000.000000# Y-Zoom +%%%%% 012 # Filament + 3.00000 # Alarm level + 2.34000 # setPoint + 0.25000 # rampdeltavoltage + 1.00000 # rampdeltatime + 0.00000 # amperedeltatime +001 # Alarm check +001 # Filament state + 2.34290 # Volt + 1.47880 # Ampere +%%%%% 013 # Virtual moduls +A1 20.2805004120 20.2805004120 -6.0046720505 66.3768157959 0.0011036640 3 4 1 FFFFD780 +dA1 -4.5489997864 -4.5489997864 -72.3814849854 72.3814849854 0.0022073281 3 4 0 FFFFC340 +A2 -7.5974998474 -7.5974998474 -10.9495010376 11.0446872711 0.0003344226 5 6 3 FFFFF700 +dA2 -0.9841995239 -0.9841995239 -21.9941883087 21.9941883087 0.0006688451 5 6 2 FFFFE2C0 +A3 1.1617850065 1.1617850065 -10.9453125000 11.0509996414 0.0003349624 7 8 5 FFFFCE80 +dA3 0.5736300349 0.5736300349 -21.9963111877 21.9963111877 0.0006699248 7 8 4 FFFFBA40 +EVM 0.2000000030 0.2000000030 -10.0000000000 10.0000000000 0.0001098963 -1 -1 -1 FFFFEE00 +UVM -0.0630699992 -0.0630699992 -2.5837626457 4.6172842979 0.0000339420 10 11 8 FFFFD9C0 +dUVM 1.5000000000 1.5000000000 -7.2010469437 7.2010469437 0.0000678840 10 11 7 FFFFC580 +DVM -0.6601200104 -0.6601200104 -2.1016812325 0.1209521890 0.0000337580 -1 -1 -1 FFFFF940 +dDVM 0.0000000000 0.0000000000 -2.2226333618 2.2226333618 0.0000675160 12 13 9 FFFFE500 +UM -0.7525577545 -0.7525577545 -2.5880718231 4.6075630188 0.0000338541 15 16 12 FFFFD0C0 +dUM 0.5000200272 0.5000200272 -7.1956348419 7.1956348419 0.0000677083 15 16 11 FFFFBC80 +DM -0.2477899790 -0.2477899790 -2.1032655239 0.1186871827 0.0000338657 -1 -1 -1 FFFFF040 +dDM -0.0006399155 -0.0006399155 -2.2219526768 2.2219526768 0.0000677314 -1 -1 -1 FFFFDC00 +B1 -0.8361939192 -0.8361939192 -16.3627204895 16.6008110046 0.0001255038 20 21 16 FFFFC7C0 +dB1 -0.0120482445 -0.0120482445 -32.9635314941 32.9635314941 0.0002510076 20 21 15 FFFFFB80 +B2 1.1802383661 1.1802383661 -5.1524972916 119.4303131104 0.0019009828 22 -1 -1 FFFFE740 +B3 0.9902987480 0.9902987480 -5.3829998970 119.1553115845 0.0019003039 23 -1 -1 FFFFD300 +B4 -0.8152312636 -0.8152312636 -16.2824687958 16.4820003510 0.0001248971 24 25 20 FFFFBEC0 +dB4 0.0019984245 0.0019984245 -32.7644691467 32.7644691467 0.0002497943 24 25 19 FFFFF280 +UA -0.7508257627 -0.7508257627 -2.5957968235 4.6316623688 0.0000338617 27 28 22 FFFFDE40 +dUA 0.4999997616 0.4999997616 -7.2274589539 7.2274589539 0.0000677235 27 28 21 FFFFCA00 +DA -0.2399498820 -0.2399498820 -2.1015281677 0.1148478091 0.0000337078 29 30 24 FFFFFDC0 +dDA 0.0314997435 0.0314997435 -2.2163760662 2.2163760662 0.0000674156 -1 -1 -1 FFFFE980 +C1 0.5900969505 0.5900969505 -5.2440004349 5.3155498505 0.0001610799 32 33 26 FFFFD540 +dC1 0.0000000000 0.0000000000 -10.5595502853 10.5595502853 0.0003221599 32 33 25 FFFFC100 +C3 3.7498872280 3.7498872280 -5.2543592453 5.3085346222 0.0001610490 35 36 28 FFFFF4C0 +dC3 -3.0332999229 -3.0332999229 -10.5628938675 10.5628938675 0.0003220980 35 36 27 FFFFE080 +C2 2.2117373943 2.2117373943 -5.2518939972 5.3111190796 0.0001611788 34 -1 -1 FFFFCC40 +E0 4.0000000000 4.0000000000 -9.2576847076 259.5778198242 0.0041021043 38 -1 -1 0 +delta phi 3.5999999046 3.5999999046 -5.0000000000 5.0000000000 0.0040949672 -1 -1 -1 FFFFEBC0 +Scan fine 0.0000000000 0.0000000000 -0.5422697067 1.2692217827 0.0000276412 40 -1 -1 FFFFD780 +AM -0.7508257627 -0.7508257627 -2.6005969048 4.6468062401 0.0001105866 26 -1 -1 FFFFC340 +Cht Eing. -5.0223627090 -5.0223627090 -5.2474431992 5.3180408478 0.0001612165 31 -1 -1 FFFFF700 +Rep -1.2324999571 -1.2324999571 -10.9234685898 11.0292186737 0.0003349714 2 -1 -1 FFFFE2C0 +PK -2.1926999092 -2.1926999092 -10.9141559601 11.1900310516 0.0003372831 1 -1 -1 FFFFCE80 +HME -0.7525577545 -0.7525577545 -2.5943999290 4.6251878738 0.0001101622 14 -1 -1 FFFFBA40 +Cht Geh -15.0000000000 -15.0000000000 -16.4459056854 16.6480636597 0.0005049739 37 -1 -1 FFFFEE00 +PKK 0.0000000000 0.0000000000 -99.0000000000 99.0000000000 0.0010000000 -1 -1 -1 FFFFD9C0 +HMA 0.0000100732 0.0000100732 -5.0000000000 5.0000000000 0.0001099585 -1 -1 -1 FFFFC580 +Mono Mitte 0.0000000000 0.0000000000 0.0000000000 0.0000000000 0.0010000000 -1 -1 -1 FFFFF940 +Scan coarse 0.0000000000 0.0000000000 -1.3911311626 51.0064392090 0.0007995235 41 -1 -1 FFFFE500 +Passenergy 0.2724998593 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFD0C0 +Kinetic En 5.0233259201 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFBC80 +Retard Rat 18.4342327118 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFF040 +Offset -1.0233256817 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFDC00 +B3 norm 0.4008548260 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFC7C0 +B4 norm 0.0414256267 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFFB80 +dB4 norm 0.0003978290 0.0000000000 -500.0000000000 500.0000000000 0.0010000000 -1 -1 -1 FFFFE740 +%%%%% 014 # Time Ramp Parameters +0:10:0 # Duration +0:0:0:600 # Step +%%%%% 015 # Measure Modes +0 # scan type +40 # scan index +0 # ramp type type +%%%%% 016 # Formular constants +F1_CON1 0.59600 +F1_CON2 1.00000 +F1_CON3 -0.38000 +F1_CON4 1.00000 +F2_CON1 0.54500 +F2_CON2 1.00000 +F2_CON3 -0.49800 +F3_CON1 0.54500 +F3_CON2 -0.49800 +F4_CON1 0.01730 +F4_CON2 0.37700 +F4_CON3 0.37700 +F4_CON4 0.01730 +F4_CON5 1.00000 +F4_CON6 1.00000 +%%%%% 017 # Measure Results +FWHM: no results +Peakfinder: no results +%%%%% 019 # + # Lenscurvename +0 # use lenscurve +0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 # Compensation Line lensType X1 Y1 X2 Y2 XOff YOff Gain +1 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 # Compensation Line lensType X1 Y1 X2 Y2 XOff YOff Gain +2 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 # Compensation Line lensType X1 Y1 X2 Y2 XOff YOff Gain +%%%%% 020 # Measure Global +02-22-2019 09:35:30 # Startzeit +02-22-2019 11:13:49 # Endzeit + 18 # Averageanzahl +10000 # Kan�le +%%%%% 021 # Measure Segment +1 # Number of Segments +1206 0.1000000015 -0.0099508204 0.2232853472 0.0001934882 0.5000000000 # Segmentsize wait start stop inc meastime +%%%%% 022 # Measure Data +CH--1 +S 0000 -0.009951 +M 0000 171.166672 +S 0001 -0.009757 +M 0001 174.166672 +S 0002 -0.009564 +M 0002 180.833328 +S 0003 -0.009370 +M 0003 180.166672 +S 0004 -0.009177 +M 0004 207.500000 +S 0005 -0.008983 +M 0005 196.500000 +S 0006 -0.008790 +M 0006 212.500000 +S 0007 -0.008596 +M 0007 229.166672 +S 0008 -0.008403 +M 0008 212.666672 +S 0009 -0.008209 +M 0009 238.000000 +S 0010 -0.008016 +M 0010 239.666672 +S 0011 -0.007822 +M 0011 247.500000 +S 0012 -0.007629 +M 0012 259.833344 +S 0013 -0.007435 +M 0013 277.500000 +S 0014 -0.007242 +M 0014 277.666656 +S 0015 -0.007049 +M 0015 286.333344 +S 0016 -0.006855 +M 0016 309.833344 +S 0017 -0.006662 +M 0017 334.000000 +S 0018 -0.006468 +M 0018 348.000000 +S 0019 -0.006275 +M 0019 376.666656 +S 0020 -0.006081 +M 0020 419.166656 +S 0021 -0.005888 +M 0021 431.166656 +S 0022 -0.005694 +M 0022 464.833344 +S 0023 -0.005501 +M 0023 508.666656 +S 0024 -0.005307 +M 0024 552.833313 +S 0025 -0.005114 +M 0025 568.333313 +S 0026 -0.004920 +M 0026 644.166687 +S 0027 -0.004727 +M 0027 710.833313 +S 0028 -0.004533 +M 0028 796.833313 +S 0029 -0.004340 +M 0029 889.833313 +S 0030 -0.004146 +M 0030 995.166687 +S 0031 -0.003953 +M 0031 1134.000000 +S 0032 -0.003759 +M 0032 1295.500000 +S 0033 -0.003566 +M 0033 1502.333374 +S 0034 -0.003372 +M 0034 1753.666626 +S 0035 -0.003179 +M 0035 2104.166748 +S 0036 -0.002985 +M 0036 2473.333252 +S 0037 -0.002792 +M 0037 3005.166748 +S 0038 -0.002598 +M 0038 3611.500000 +S 0039 -0.002405 +M 0039 4460.833496 +S 0040 -0.002211 +M 0040 5459.166504 +S 0041 -0.002018 +M 0041 6694.000000 +S 0042 -0.001824 +M 0042 8152.000000 +S 0043 -0.001631 +M 0043 9975.166992 +S 0044 -0.001437 +M 0044 12174.500000 +S 0045 -0.001244 +M 0045 14396.333008 +S 0046 -0.001050 +M 0046 16924.333984 +S 0047 -0.000857 +M 0047 19034.000000 +S 0048 -0.000663 +M 0048 21410.500000 +S 0049 -0.000470 +M 0049 23335.500000 +S 0050 -0.000276 +M 0050 25068.166016 +S 0051 -0.000083 +M 0051 25830.666016 +S 0052 0.000111 +M 0052 25772.333984 +S 0053 0.000304 +M 0053 24502.666016 +S 0054 0.000498 +M 0054 22729.166016 +S 0055 0.000691 +M 0055 20179.500000 +S 0056 0.000885 +M 0056 17682.833984 +S 0057 0.001078 +M 0057 14886.333008 +S 0058 0.001271 +M 0058 12403.000000 +S 0059 0.001465 +M 0059 10260.666992 +S 0060 0.001658 +M 0060 8307.666992 +S 0061 0.001852 +M 0061 6737.333496 +S 0062 0.002045 +M 0062 5502.333496 +S 0063 0.002239 +M 0063 4472.333496 +S 0064 0.002432 +M 0064 3686.166748 +S 0065 0.002626 +M 0065 3101.500000 +S 0066 0.002819 +M 0066 2659.166748 +S 0067 0.003013 +M 0067 2307.500000 +S 0068 0.003206 +M 0068 1979.000000 +S 0069 0.003400 +M 0069 1726.833374 +S 0070 0.003593 +M 0070 1558.500000 +S 0071 0.003787 +M 0071 1382.666626 +S 0072 0.003980 +M 0072 1241.000000 +S 0073 0.004174 +M 0073 1122.333374 +S 0074 0.004367 +M 0074 1074.166626 +S 0075 0.004561 +M 0075 961.833313 +S 0076 0.004754 +M 0076 888.166687 +S 0077 0.004948 +M 0077 846.666687 +S 0078 0.005141 +M 0078 779.333313 +S 0079 0.005335 +M 0079 704.833313 +S 0080 0.005528 +M 0080 690.166687 +S 0081 0.005722 +M 0081 633.500000 +S 0082 0.005915 +M 0082 621.000000 +S 0083 0.006109 +M 0083 619.666687 +S 0084 0.006302 +M 0084 561.500000 +S 0085 0.006496 +M 0085 551.500000 +S 0086 0.006689 +M 0086 494.333344 +S 0087 0.006883 +M 0087 486.500000 +S 0088 0.007076 +M 0088 466.333344 +S 0089 0.007270 +M 0089 459.333344 +S 0090 0.007463 +M 0090 437.000000 +S 0091 0.007657 +M 0091 413.833344 +S 0092 0.007850 +M 0092 405.166656 +S 0093 0.008044 +M 0093 384.333344 +S 0094 0.008237 +M 0094 380.833344 +S 0095 0.008431 +M 0095 348.000000 +S 0096 0.008624 +M 0096 343.000000 +S 0097 0.008818 +M 0097 330.000000 +S 0098 0.009011 +M 0098 339.833344 +S 0099 0.009205 +M 0099 313.000000 +S 0100 0.009398 +M 0100 320.000000 +S 0101 0.009591 +M 0101 299.000000 +S 0102 0.009785 +M 0102 279.333344 +S 0103 0.009978 +M 0103 277.833344 +S 0104 0.010172 +M 0104 266.333344 +S 0105 0.010365 +M 0105 257.166656 +S 0106 0.010559 +M 0106 282.666656 +S 0107 0.010752 +M 0107 246.166672 +S 0108 0.010946 +M 0108 241.500000 +S 0109 0.011139 +M 0109 242.833328 +S 0110 0.011333 +M 0110 229.833328 +S 0111 0.011526 +M 0111 217.333328 +S 0112 0.011720 +M 0112 214.500000 +S 0113 0.011913 +M 0113 206.833328 +S 0114 0.012107 +M 0114 195.500000 +S 0115 0.012300 +M 0115 209.166672 +S 0116 0.012494 +M 0116 198.833328 +S 0117 0.012687 +M 0117 206.500000 +S 0118 0.012881 +M 0118 195.166672 +S 0119 0.013074 +M 0119 193.333328 +S 0120 0.013268 +M 0120 191.000000 +S 0121 0.013461 +M 0121 187.000000 +S 0122 0.013655 +M 0122 180.333328 +S 0123 0.013848 +M 0123 186.833328 +S 0124 0.014042 +M 0124 171.166672 +S 0125 0.014235 +M 0125 176.000000 +S 0126 0.014429 +M 0126 164.000000 +S 0127 0.014622 +M 0127 170.000000 +S 0128 0.014816 +M 0128 160.000000 +S 0129 0.015009 +M 0129 161.833328 +S 0130 0.015203 +M 0130 171.666672 +S 0131 0.015396 +M 0131 155.333328 +S 0132 0.015590 +M 0132 173.666672 +S 0133 0.015783 +M 0133 161.666672 +S 0134 0.015977 +M 0134 155.000000 +S 0135 0.016170 +M 0135 164.500000 +S 0136 0.016364 +M 0136 141.166672 +S 0137 0.016557 +M 0137 147.333328 +S 0138 0.016751 +M 0138 154.666672 +S 0139 0.016944 +M 0139 160.500000 +S 0140 0.017138 +M 0140 155.000000 +S 0141 0.017331 +M 0141 159.500000 +S 0142 0.017524 +M 0142 155.500000 +S 0143 0.017718 +M 0143 161.333328 +S 0144 0.017911 +M 0144 159.000000 +S 0145 0.018105 +M 0145 154.333328 +S 0146 0.018298 +M 0146 148.333328 +S 0147 0.018492 +M 0147 147.666672 +S 0148 0.018685 +M 0148 152.500000 +S 0149 0.018879 +M 0149 154.000000 +S 0150 0.019072 +M 0150 164.000000 +S 0151 0.019266 +M 0151 149.833328 +S 0152 0.019459 +M 0152 147.500000 +S 0153 0.019653 +M 0153 160.333328 +S 0154 0.019846 +M 0154 173.333328 +S 0155 0.020040 +M 0155 191.833328 +S 0156 0.020233 +M 0156 168.666672 +S 0157 0.020427 +M 0157 188.500000 +S 0158 0.020620 +M 0158 204.166672 +S 0159 0.020814 +M 0159 188.833328 +S 0160 0.021007 +M 0160 205.166672 +S 0161 0.021201 +M 0161 203.666672 +S 0162 0.021394 +M 0162 183.833328 +S 0163 0.021588 +M 0163 197.000000 +S 0164 0.021781 +M 0164 214.833328 +S 0165 0.021975 +M 0165 198.166672 +S 0166 0.022168 +M 0166 191.500000 +S 0167 0.022362 +M 0167 187.666672 +S 0168 0.022555 +M 0168 169.000000 +S 0169 0.022749 +M 0169 174.500000 +S 0170 0.022942 +M 0170 170.833328 +S 0171 0.023136 +M 0171 169.000000 +S 0172 0.023329 +M 0172 164.666672 +S 0173 0.023523 +M 0173 162.166672 +S 0174 0.023716 +M 0174 170.666672 +S 0175 0.023910 +M 0175 172.166672 +S 0176 0.024103 +M 0176 159.166672 +S 0177 0.024297 +M 0177 158.333328 +S 0178 0.024490 +M 0178 163.500000 +S 0179 0.024684 +M 0179 154.500000 +S 0180 0.024877 +M 0180 148.000000 +S 0181 0.025071 +M 0181 149.166672 +S 0182 0.025264 +M 0182 153.333328 +S 0183 0.025457 +M 0183 148.333328 +S 0184 0.025651 +M 0184 140.000000 +S 0185 0.025844 +M 0185 149.333328 +S 0186 0.026038 +M 0186 149.833328 +S 0187 0.026231 +M 0187 158.333328 +S 0188 0.026425 +M 0188 155.166672 +S 0189 0.026618 +M 0189 151.500000 +S 0190 0.026812 +M 0190 143.666672 +S 0191 0.027005 +M 0191 141.333328 +S 0192 0.027199 +M 0192 146.666672 +S 0193 0.027392 +M 0193 146.500000 +S 0194 0.027586 +M 0194 150.666672 +S 0195 0.027779 +M 0195 138.000000 +S 0196 0.027973 +M 0196 152.833328 +S 0197 0.028166 +M 0197 158.833328 +S 0198 0.028360 +M 0198 149.333328 +S 0199 0.028553 +M 0199 156.500000 +S 0200 0.028747 +M 0200 155.666672 +S 0201 0.028940 +M 0201 155.333328 +S 0202 0.029134 +M 0202 149.333328 +S 0203 0.029327 +M 0203 152.000000 +S 0204 0.029521 +M 0204 157.166672 +S 0205 0.029714 +M 0205 155.000000 +S 0206 0.029908 +M 0206 149.166672 +S 0207 0.030101 +M 0207 161.833328 +S 0208 0.030295 +M 0208 166.500000 +S 0209 0.030488 +M 0209 165.000000 +S 0210 0.030682 +M 0210 173.833328 +S 0211 0.030875 +M 0211 169.666672 +S 0212 0.031069 +M 0212 156.666672 +S 0213 0.031262 +M 0213 178.166672 +S 0214 0.031456 +M 0214 162.666672 +S 0215 0.031649 +M 0215 180.166672 +S 0216 0.031843 +M 0216 173.000000 +S 0217 0.032036 +M 0217 182.166672 +S 0218 0.032230 +M 0218 171.500000 +S 0219 0.032423 +M 0219 173.333328 +S 0220 0.032617 +M 0220 174.500000 +S 0221 0.032810 +M 0221 184.666672 +S 0222 0.033004 +M 0222 185.000000 +S 0223 0.033197 +M 0223 171.833328 +S 0224 0.033390 +M 0224 172.500000 +S 0225 0.033584 +M 0225 183.333328 +S 0226 0.033777 +M 0226 180.666672 +S 0227 0.033971 +M 0227 184.166672 +S 0228 0.034164 +M 0228 186.166672 +S 0229 0.034358 +M 0229 181.000000 +S 0230 0.034551 +M 0230 180.166672 +S 0231 0.034745 +M 0231 188.833328 +S 0232 0.034938 +M 0232 198.166672 +S 0233 0.035132 +M 0233 185.666672 +S 0234 0.035325 +M 0234 185.833328 +S 0235 0.035519 +M 0235 180.166672 +S 0236 0.035712 +M 0236 205.500000 +S 0237 0.035906 +M 0237 177.000000 +S 0238 0.036099 +M 0238 176.000000 +S 0239 0.036293 +M 0239 189.166672 +S 0240 0.036486 +M 0240 180.500000 +S 0241 0.036680 +M 0241 172.833328 +S 0242 0.036873 +M 0242 173.833328 +S 0243 0.037067 +M 0243 183.166672 +S 0244 0.037260 +M 0244 184.333328 +S 0245 0.037454 +M 0245 180.500000 +S 0246 0.037647 +M 0246 172.333328 +S 0247 0.037841 +M 0247 173.833328 +S 0248 0.038034 +M 0248 163.666672 +S 0249 0.038228 +M 0249 168.333328 +S 0250 0.038421 +M 0250 160.166672 +S 0251 0.038615 +M 0251 179.666672 +S 0252 0.038808 +M 0252 166.833328 +S 0253 0.039002 +M 0253 167.666672 +S 0254 0.039195 +M 0254 175.000000 +S 0255 0.039389 +M 0255 171.000000 +S 0256 0.039582 +M 0256 161.666672 +S 0257 0.039776 +M 0257 167.333328 +S 0258 0.039969 +M 0258 160.333328 +S 0259 0.040163 +M 0259 155.166672 +S 0260 0.040356 +M 0260 163.500000 +S 0261 0.040550 +M 0261 173.833328 +S 0262 0.040743 +M 0262 165.166672 +S 0263 0.040937 +M 0263 175.500000 +S 0264 0.041130 +M 0264 161.166672 +S 0265 0.041323 +M 0265 177.666672 +S 0266 0.041517 +M 0266 169.333328 +S 0267 0.041710 +M 0267 170.333328 +S 0268 0.041904 +M 0268 153.500000 +S 0269 0.042097 +M 0269 170.333328 +S 0270 0.042291 +M 0270 161.500000 +S 0271 0.042484 +M 0271 168.333328 +S 0272 0.042678 +M 0272 179.500000 +S 0273 0.042871 +M 0273 163.500000 +S 0274 0.043065 +M 0274 175.333328 +S 0275 0.043258 +M 0275 171.166672 +S 0276 0.043452 +M 0276 175.666672 +S 0277 0.043645 +M 0277 167.833328 +S 0278 0.043839 +M 0278 174.500000 +S 0279 0.044032 +M 0279 184.833328 +S 0280 0.044226 +M 0280 168.333328 +S 0281 0.044419 +M 0281 174.500000 +S 0282 0.044613 +M 0282 182.166672 +S 0283 0.044806 +M 0283 177.166672 +S 0284 0.045000 +M 0284 173.166672 +S 0285 0.045193 +M 0285 185.333328 +S 0286 0.045387 +M 0286 187.333328 +S 0287 0.045580 +M 0287 199.666672 +S 0288 0.045774 +M 0288 193.500000 +S 0289 0.045967 +M 0289 207.166672 +S 0290 0.046161 +M 0290 183.166672 +S 0291 0.046354 +M 0291 207.000000 +S 0292 0.046548 +M 0292 197.666672 +S 0293 0.046741 +M 0293 202.666672 +S 0294 0.046935 +M 0294 208.666672 +S 0295 0.047128 +M 0295 216.500000 +S 0296 0.047322 +M 0296 219.166672 +S 0297 0.047515 +M 0297 211.833328 +S 0298 0.047709 +M 0298 226.500000 +S 0299 0.047902 +M 0299 216.500000 +S 0300 0.048096 +M 0300 211.000000 +S 0301 0.048289 +M 0301 241.166672 +S 0302 0.048483 +M 0302 250.500000 +S 0303 0.048676 +M 0303 242.833328 +S 0304 0.048870 +M 0304 247.333328 +S 0305 0.049063 +M 0305 258.000000 +S 0306 0.049256 +M 0306 275.333344 +S 0307 0.049450 +M 0307 278.166656 +S 0308 0.049643 +M 0308 268.333344 +S 0309 0.049837 +M 0309 279.166656 +S 0310 0.050030 +M 0310 290.666656 +S 0311 0.050224 +M 0311 301.833344 +S 0312 0.050417 +M 0312 299.666656 +S 0313 0.050611 +M 0313 313.500000 +S 0314 0.050804 +M 0314 326.333344 +S 0315 0.050998 +M 0315 346.000000 +S 0316 0.051191 +M 0316 351.666656 +S 0317 0.051385 +M 0317 353.833344 +S 0318 0.051578 +M 0318 347.333344 +S 0319 0.051772 +M 0319 402.500000 +S 0320 0.051965 +M 0320 392.000000 +S 0321 0.052159 +M 0321 423.333344 +S 0322 0.052352 +M 0322 435.833344 +S 0323 0.052546 +M 0323 453.500000 +S 0324 0.052739 +M 0324 465.333344 +S 0325 0.052933 +M 0325 497.500000 +S 0326 0.053126 +M 0326 521.333313 +S 0327 0.053320 +M 0327 551.666687 +S 0328 0.053513 +M 0328 581.333313 +S 0329 0.053707 +M 0329 606.166687 +S 0330 0.053900 +M 0330 647.500000 +S 0331 0.054094 +M 0331 713.333313 +S 0332 0.054287 +M 0332 774.333313 +S 0333 0.054481 +M 0333 796.333313 +S 0334 0.054674 +M 0334 846.666687 +S 0335 0.054868 +M 0335 918.333313 +S 0336 0.055061 +M 0336 988.000000 +S 0337 0.055255 +M 0337 1094.166626 +S 0338 0.055448 +M 0338 1193.333374 +S 0339 0.055642 +M 0339 1251.000000 +S 0340 0.055835 +M 0340 1341.500000 +S 0341 0.056029 +M 0341 1486.333374 +S 0342 0.056222 +M 0342 1573.666626 +S 0343 0.056416 +M 0343 1673.166626 +S 0344 0.056609 +M 0344 1746.000000 +S 0345 0.056803 +M 0345 1895.666626 +S 0346 0.056996 +M 0346 2007.333374 +S 0347 0.057190 +M 0347 2145.500000 +S 0348 0.057383 +M 0348 2199.833252 +S 0349 0.057576 +M 0349 2255.000000 +S 0350 0.057770 +M 0350 2339.000000 +S 0351 0.057963 +M 0351 2317.166748 +S 0352 0.058157 +M 0352 2284.500000 +S 0353 0.058350 +M 0353 2278.166748 +S 0354 0.058544 +M 0354 2179.666748 +S 0355 0.058737 +M 0355 2106.833252 +S 0356 0.058931 +M 0356 2008.333374 +S 0357 0.059124 +M 0357 1873.500000 +S 0358 0.059318 +M 0358 1701.000000 +S 0359 0.059511 +M 0359 1554.500000 +S 0360 0.059705 +M 0360 1447.500000 +S 0361 0.059898 +M 0361 1307.833374 +S 0362 0.060092 +M 0362 1159.500000 +S 0363 0.060285 +M 0363 1083.333374 +S 0364 0.060479 +M 0364 976.833313 +S 0365 0.060672 +M 0365 884.333313 +S 0366 0.060866 +M 0366 761.000000 +S 0367 0.061059 +M 0367 716.000000 +S 0368 0.061253 +M 0368 680.833313 +S 0369 0.061446 +M 0369 612.500000 +S 0370 0.061640 +M 0370 545.166687 +S 0371 0.061833 +M 0371 521.333313 +S 0372 0.062027 +M 0372 472.333344 +S 0373 0.062220 +M 0373 446.500000 +S 0374 0.062414 +M 0374 413.166656 +S 0375 0.062607 +M 0375 391.666656 +S 0376 0.062801 +M 0376 374.000000 +S 0377 0.062994 +M 0377 374.500000 +S 0378 0.063188 +M 0378 338.666656 +S 0379 0.063381 +M 0379 308.000000 +S 0380 0.063575 +M 0380 298.500000 +S 0381 0.063768 +M 0381 286.500000 +S 0382 0.063962 +M 0382 271.833344 +S 0383 0.064155 +M 0383 277.666656 +S 0384 0.064349 +M 0384 252.500000 +S 0385 0.064542 +M 0385 251.500000 +S 0386 0.064736 +M 0386 248.000000 +S 0387 0.064929 +M 0387 237.666672 +S 0388 0.065123 +M 0388 244.833328 +S 0389 0.065316 +M 0389 245.000000 +S 0390 0.065510 +M 0390 224.500000 +S 0391 0.065703 +M 0391 236.333328 +S 0392 0.065897 +M 0392 226.333328 +S 0393 0.066090 +M 0393 215.833328 +S 0394 0.066284 +M 0394 212.000000 +S 0395 0.066477 +M 0395 218.500000 +S 0396 0.066670 +M 0396 210.166672 +S 0397 0.066864 +M 0397 199.833328 +S 0398 0.067057 +M 0398 217.333328 +S 0399 0.067251 +M 0399 221.000000 +S 0400 0.067444 +M 0400 215.666672 +S 0401 0.067638 +M 0401 206.166672 +S 0402 0.067831 +M 0402 210.333328 +S 0403 0.068025 +M 0403 220.666672 +S 0404 0.068218 +M 0404 209.666672 +S 0405 0.068412 +M 0405 206.333328 +S 0406 0.068605 +M 0406 202.166672 +S 0407 0.068799 +M 0407 220.666672 +S 0408 0.068992 +M 0408 210.000000 +S 0409 0.069186 +M 0409 220.166672 +S 0410 0.069379 +M 0410 206.500000 +S 0411 0.069573 +M 0411 216.166672 +S 0412 0.069766 +M 0412 216.833328 +S 0413 0.069960 +M 0413 208.833328 +S 0414 0.070153 +M 0414 220.833328 +S 0415 0.070347 +M 0415 215.666672 +S 0416 0.070540 +M 0416 221.666672 +S 0417 0.070734 +M 0417 226.833328 +S 0418 0.070927 +M 0418 225.666672 +S 0419 0.071121 +M 0419 207.166672 +S 0420 0.071314 +M 0420 211.000000 +S 0421 0.071508 +M 0421 226.500000 +S 0422 0.071701 +M 0422 213.333328 +S 0423 0.071895 +M 0423 235.833328 +S 0424 0.072088 +M 0424 242.666672 +S 0425 0.072282 +M 0425 231.500000 +S 0426 0.072475 +M 0426 222.666672 +S 0427 0.072669 +M 0427 224.333328 +S 0428 0.072862 +M 0428 246.000000 +S 0429 0.073056 +M 0429 223.500000 +S 0430 0.073249 +M 0430 230.000000 +S 0431 0.073443 +M 0431 257.833344 +S 0432 0.073636 +M 0432 233.166672 +S 0433 0.073830 +M 0433 239.000000 +S 0434 0.074023 +M 0434 246.166672 +S 0435 0.074217 +M 0435 255.500000 +S 0436 0.074410 +M 0436 261.166656 +S 0437 0.074604 +M 0437 252.500000 +S 0438 0.074797 +M 0438 251.833328 +S 0439 0.074991 +M 0439 263.833344 +S 0440 0.075184 +M 0440 267.666656 +S 0441 0.075378 +M 0441 245.166672 +S 0442 0.075571 +M 0442 255.666672 +S 0443 0.075765 +M 0443 266.833344 +S 0444 0.075958 +M 0444 273.833344 +S 0445 0.076152 +M 0445 285.500000 +S 0446 0.076345 +M 0446 283.000000 +S 0447 0.076539 +M 0447 294.166656 +S 0448 0.076732 +M 0448 286.000000 +S 0449 0.076926 +M 0449 302.833344 +S 0450 0.077119 +M 0450 308.333344 +S 0451 0.077313 +M 0451 293.166656 +S 0452 0.077506 +M 0452 296.666656 +S 0453 0.077700 +M 0453 302.333344 +S 0454 0.077893 +M 0454 314.500000 +S 0455 0.078086 +M 0455 313.833344 +S 0456 0.078280 +M 0456 310.833344 +S 0457 0.078473 +M 0457 346.333344 +S 0458 0.078667 +M 0458 328.166656 +S 0459 0.078860 +M 0459 340.000000 +S 0460 0.079054 +M 0460 349.000000 +S 0461 0.079247 +M 0461 372.166656 +S 0462 0.079441 +M 0462 365.666656 +S 0463 0.079634 +M 0463 372.833344 +S 0464 0.079828 +M 0464 371.000000 +S 0465 0.080021 +M 0465 365.166656 +S 0466 0.080215 +M 0466 386.833344 +S 0467 0.080408 +M 0467 408.500000 +S 0468 0.080602 +M 0468 393.333344 +S 0469 0.080795 +M 0469 431.000000 +S 0470 0.080989 +M 0470 435.500000 +S 0471 0.081182 +M 0471 448.666656 +S 0472 0.081376 +M 0472 439.000000 +S 0473 0.081569 +M 0473 440.000000 +S 0474 0.081763 +M 0474 469.333344 +S 0475 0.081956 +M 0475 489.000000 +S 0476 0.082150 +M 0476 457.166656 +S 0477 0.082343 +M 0477 538.833313 +S 0478 0.082537 +M 0478 514.500000 +S 0479 0.082730 +M 0479 538.500000 +S 0480 0.082924 +M 0480 538.666687 +S 0481 0.083117 +M 0481 565.666687 +S 0482 0.083311 +M 0482 603.333313 +S 0483 0.083504 +M 0483 627.166687 +S 0484 0.083698 +M 0484 636.000000 +S 0485 0.083891 +M 0485 661.500000 +S 0486 0.084085 +M 0486 685.833313 +S 0487 0.084278 +M 0487 728.500000 +S 0488 0.084472 +M 0488 748.833313 +S 0489 0.084665 +M 0489 754.000000 +S 0490 0.084859 +M 0490 828.000000 +S 0491 0.085052 +M 0491 799.666687 +S 0492 0.085246 +M 0492 867.833313 +S 0493 0.085439 +M 0493 897.833313 +S 0494 0.085633 +M 0494 958.500000 +S 0495 0.085826 +M 0495 1021.333313 +S 0496 0.086020 +M 0496 1059.166626 +S 0497 0.086213 +M 0497 1110.500000 +S 0498 0.086407 +M 0498 1141.000000 +S 0499 0.086600 +M 0499 1192.000000 +S 0500 0.086794 +M 0500 1273.666626 +S 0501 0.086987 +M 0501 1316.833374 +S 0502 0.087181 +M 0502 1393.666626 +S 0503 0.087374 +M 0503 1424.333374 +S 0504 0.087568 +M 0504 1489.166626 +S 0505 0.087761 +M 0505 1597.500000 +S 0506 0.087955 +M 0506 1619.833374 +S 0507 0.088148 +M 0507 1683.666626 +S 0508 0.088342 +M 0508 1761.333374 +S 0509 0.088535 +M 0509 1851.333374 +S 0510 0.088729 +M 0510 1928.833374 +S 0511 0.088922 +M 0511 2014.833374 +S 0512 0.089116 +M 0512 2072.666748 +S 0513 0.089309 +M 0513 2179.666748 +S 0514 0.089502 +M 0514 2306.166748 +S 0515 0.089696 +M 0515 2409.333252 +S 0516 0.089889 +M 0516 2452.666748 +S 0517 0.090083 +M 0517 2593.333252 +S 0518 0.090276 +M 0518 2683.666748 +S 0519 0.090470 +M 0519 2722.000000 +S 0520 0.090663 +M 0520 2802.333252 +S 0521 0.090857 +M 0521 2941.500000 +S 0522 0.091050 +M 0522 2961.166748 +S 0523 0.091244 +M 0523 3055.500000 +S 0524 0.091437 +M 0524 3110.333252 +S 0525 0.091631 +M 0525 3132.666748 +S 0526 0.091824 +M 0526 3145.000000 +S 0527 0.092018 +M 0527 3175.166748 +S 0528 0.092211 +M 0528 3191.166748 +S 0529 0.092405 +M 0529 3211.500000 +S 0530 0.092598 +M 0530 3169.166748 +S 0531 0.092792 +M 0531 3116.000000 +S 0532 0.092985 +M 0532 3095.000000 +S 0533 0.093179 +M 0533 3031.000000 +S 0534 0.093372 +M 0534 2981.666748 +S 0535 0.093566 +M 0535 2908.666748 +S 0536 0.093759 +M 0536 2781.333252 +S 0537 0.093953 +M 0537 2759.833252 +S 0538 0.094146 +M 0538 2612.500000 +S 0539 0.094340 +M 0539 2566.500000 +S 0540 0.094533 +M 0540 2412.333252 +S 0541 0.094727 +M 0541 2341.500000 +S 0542 0.094920 +M 0542 2236.333252 +S 0543 0.095114 +M 0543 2129.500000 +S 0544 0.095307 +M 0544 2029.000000 +S 0545 0.095501 +M 0545 2000.000000 +S 0546 0.095694 +M 0546 1866.333374 +S 0547 0.095888 +M 0547 1784.166626 +S 0548 0.096081 +M 0548 1672.500000 +S 0549 0.096275 +M 0549 1594.500000 +S 0550 0.096468 +M 0550 1554.000000 +S 0551 0.096662 +M 0551 1477.666626 +S 0552 0.096855 +M 0552 1426.666626 +S 0553 0.097049 +M 0553 1329.666626 +S 0554 0.097242 +M 0554 1235.333374 +S 0555 0.097436 +M 0555 1197.500000 +S 0556 0.097629 +M 0556 1142.166626 +S 0557 0.097823 +M 0557 1090.333374 +S 0558 0.098016 +M 0558 1087.000000 +S 0559 0.098210 +M 0559 1009.333313 +S 0560 0.098403 +M 0560 986.166687 +S 0561 0.098597 +M 0561 925.500000 +S 0562 0.098790 +M 0562 913.666687 +S 0563 0.098984 +M 0563 852.166687 +S 0564 0.099177 +M 0564 842.833313 +S 0565 0.099371 +M 0565 810.833313 +S 0566 0.099564 +M 0566 782.166687 +S 0567 0.099758 +M 0567 758.500000 +S 0568 0.099951 +M 0568 731.333313 +S 0569 0.100145 +M 0569 685.833313 +S 0570 0.100338 +M 0570 685.166687 +S 0571 0.100532 +M 0571 665.500000 +S 0572 0.100725 +M 0572 650.166687 +S 0573 0.100918 +M 0573 630.666687 +S 0574 0.101112 +M 0574 602.500000 +S 0575 0.101305 +M 0575 585.000000 +S 0576 0.101499 +M 0576 562.166687 +S 0577 0.101692 +M 0577 536.333313 +S 0578 0.101886 +M 0578 513.000000 +S 0579 0.102079 +M 0579 516.833313 +S 0580 0.102273 +M 0580 496.666656 +S 0581 0.102466 +M 0581 483.166656 +S 0582 0.102660 +M 0582 484.500000 +S 0583 0.102853 +M 0583 462.166656 +S 0584 0.103047 +M 0584 459.333344 +S 0585 0.103240 +M 0585 444.833344 +S 0586 0.103434 +M 0586 417.500000 +S 0587 0.103627 +M 0587 414.166656 +S 0588 0.103821 +M 0588 419.166656 +S 0589 0.104014 +M 0589 426.166656 +S 0590 0.104208 +M 0590 373.666656 +S 0591 0.104401 +M 0591 380.166656 +S 0592 0.104595 +M 0592 378.333344 +S 0593 0.104788 +M 0593 371.333344 +S 0594 0.104982 +M 0594 363.833344 +S 0595 0.105175 +M 0595 360.000000 +S 0596 0.105369 +M 0596 349.666656 +S 0597 0.105562 +M 0597 349.833344 +S 0598 0.105756 +M 0598 331.833344 +S 0599 0.105949 +M 0599 321.833344 +S 0600 0.106143 +M 0600 324.666656 +S 0601 0.106336 +M 0601 304.333344 +S 0602 0.106530 +M 0602 316.166656 +S 0603 0.106723 +M 0603 306.833344 +S 0604 0.106917 +M 0604 301.166656 +S 0605 0.107110 +M 0605 284.500000 +S 0606 0.107304 +M 0606 290.333344 +S 0607 0.107497 +M 0607 283.666656 +S 0608 0.107691 +M 0608 283.500000 +S 0609 0.107884 +M 0609 278.333344 +S 0610 0.108078 +M 0610 271.000000 +S 0611 0.108271 +M 0611 261.166656 +S 0612 0.108465 +M 0612 260.666656 +S 0613 0.108658 +M 0613 277.500000 +S 0614 0.108852 +M 0614 260.166656 +S 0615 0.109045 +M 0615 249.333328 +S 0616 0.109239 +M 0616 251.333328 +S 0617 0.109432 +M 0617 266.833344 +S 0618 0.109626 +M 0618 262.666656 +S 0619 0.109819 +M 0619 250.500000 +S 0620 0.110013 +M 0620 254.833328 +S 0621 0.110206 +M 0621 253.000000 +S 0622 0.110400 +M 0622 250.833328 +S 0623 0.110593 +M 0623 248.500000 +S 0624 0.110787 +M 0624 243.333328 +S 0625 0.110980 +M 0625 251.500000 +S 0626 0.111174 +M 0626 249.500000 +S 0627 0.111367 +M 0627 250.166672 +S 0628 0.111561 +M 0628 252.666672 +S 0629 0.111754 +M 0629 243.500000 +S 0630 0.111948 +M 0630 248.166672 +S 0631 0.112141 +M 0631 256.166656 +S 0632 0.112334 +M 0632 247.000000 +S 0633 0.112528 +M 0633 252.000000 +S 0634 0.112721 +M 0634 246.666672 +S 0635 0.112915 +M 0635 247.333328 +S 0636 0.113108 +M 0636 266.666656 +S 0637 0.113302 +M 0637 262.000000 +S 0638 0.113495 +M 0638 263.500000 +S 0639 0.113689 +M 0639 254.500000 +S 0640 0.113882 +M 0640 268.666656 +S 0641 0.114076 +M 0641 261.833344 +S 0642 0.114269 +M 0642 273.166656 +S 0643 0.114463 +M 0643 277.500000 +S 0644 0.114656 +M 0644 273.666656 +S 0645 0.114850 +M 0645 268.666656 +S 0646 0.115043 +M 0646 272.500000 +S 0647 0.115237 +M 0647 280.000000 +S 0648 0.115430 +M 0648 285.500000 +S 0649 0.115624 +M 0649 277.333344 +S 0650 0.115817 +M 0650 278.166656 +S 0651 0.116011 +M 0651 257.833344 +S 0652 0.116204 +M 0652 268.666656 +S 0653 0.116398 +M 0653 276.333344 +S 0654 0.116591 +M 0654 265.666656 +S 0655 0.116785 +M 0655 272.000000 +S 0656 0.116978 +M 0656 240.333328 +S 0657 0.117172 +M 0657 235.666672 +S 0658 0.117365 +M 0658 235.833328 +S 0659 0.117559 +M 0659 231.500000 +S 0660 0.117752 +M 0660 222.166672 +S 0661 0.117946 +M 0661 233.500000 +S 0662 0.118139 +M 0662 223.500000 +S 0663 0.118333 +M 0663 211.500000 +S 0664 0.118526 +M 0664 198.333328 +S 0665 0.118720 +M 0665 196.500000 +S 0666 0.118913 +M 0666 204.833328 +S 0667 0.119107 +M 0667 192.166672 +S 0668 0.119300 +M 0668 175.166672 +S 0669 0.119494 +M 0669 171.166672 +S 0670 0.119687 +M 0670 176.833328 +S 0671 0.119881 +M 0671 164.833328 +S 0672 0.120074 +M 0672 159.333328 +S 0673 0.120268 +M 0673 164.000000 +S 0674 0.120461 +M 0674 171.333328 +S 0675 0.120655 +M 0675 155.166672 +S 0676 0.120848 +M 0676 160.333328 +S 0677 0.121042 +M 0677 153.166672 +S 0678 0.121235 +M 0678 146.500000 +S 0679 0.121429 +M 0679 146.666672 +S 0680 0.121622 +M 0680 150.000000 +S 0681 0.121816 +M 0681 137.666672 +S 0682 0.122009 +M 0682 148.666672 +S 0683 0.122203 +M 0683 146.166672 +S 0684 0.122396 +M 0684 138.833328 +S 0685 0.122590 +M 0685 138.000000 +S 0686 0.122783 +M 0686 134.833328 +S 0687 0.122977 +M 0687 147.500000 +S 0688 0.123170 +M 0688 131.500000 +S 0689 0.123364 +M 0689 138.666672 +S 0690 0.123557 +M 0690 145.500000 +S 0691 0.123751 +M 0691 130.000000 +S 0692 0.123944 +M 0692 145.500000 +S 0693 0.124137 +M 0693 131.500000 +S 0694 0.124331 +M 0694 130.000000 +S 0695 0.124524 +M 0695 134.333328 +S 0696 0.124718 +M 0696 138.666672 +S 0697 0.124911 +M 0697 128.333328 +S 0698 0.125105 +M 0698 131.833328 +S 0699 0.125298 +M 0699 127.833336 +S 0700 0.125492 +M 0700 131.666672 +S 0701 0.125685 +M 0701 122.833336 +S 0702 0.125879 +M 0702 125.166664 +S 0703 0.126072 +M 0703 122.333336 +S 0704 0.126266 +M 0704 132.666672 +S 0705 0.126459 +M 0705 128.500000 +S 0706 0.126653 +M 0706 139.833328 +S 0707 0.126846 +M 0707 123.833336 +S 0708 0.127040 +M 0708 130.333328 +S 0709 0.127233 +M 0709 131.000000 +S 0710 0.127427 +M 0710 126.166664 +S 0711 0.127620 +M 0711 127.166664 +S 0712 0.127814 +M 0712 129.166672 +S 0713 0.128007 +M 0713 127.666664 +S 0714 0.128201 +M 0714 127.833336 +S 0715 0.128394 +M 0715 129.666672 +S 0716 0.128588 +M 0716 136.333328 +S 0717 0.128781 +M 0717 119.000000 +S 0718 0.128975 +M 0718 123.500000 +S 0719 0.129168 +M 0719 122.000000 +S 0720 0.129362 +M 0720 119.166664 +S 0721 0.129555 +M 0721 133.833328 +S 0722 0.129749 +M 0722 130.833328 +S 0723 0.129942 +M 0723 123.333336 +S 0724 0.130136 +M 0724 136.666672 +S 0725 0.130329 +M 0725 123.833336 +S 0726 0.130523 +M 0726 133.833328 +S 0727 0.130716 +M 0727 120.500000 +S 0728 0.130910 +M 0728 131.833328 +S 0729 0.131103 +M 0729 131.333328 +S 0730 0.131297 +M 0730 125.000000 +S 0731 0.131490 +M 0731 123.833336 +S 0732 0.131684 +M 0732 133.000000 +S 0733 0.131877 +M 0733 128.333328 +S 0734 0.132071 +M 0734 123.166664 +S 0735 0.132264 +M 0735 127.500000 +S 0736 0.132458 +M 0736 125.500000 +S 0737 0.132651 +M 0737 141.166672 +S 0738 0.132845 +M 0738 128.833328 +S 0739 0.133038 +M 0739 140.000000 +S 0740 0.133232 +M 0740 121.000000 +S 0741 0.133425 +M 0741 135.166672 +S 0742 0.133619 +M 0742 131.833328 +S 0743 0.133812 +M 0743 134.166672 +S 0744 0.134006 +M 0744 134.333328 +S 0745 0.134199 +M 0745 124.500000 +S 0746 0.134393 +M 0746 131.666672 +S 0747 0.134586 +M 0747 127.500000 +S 0748 0.134780 +M 0748 144.000000 +S 0749 0.134973 +M 0749 131.833328 +S 0750 0.135166 +M 0750 129.500000 +S 0751 0.135360 +M 0751 136.166672 +S 0752 0.135553 +M 0752 141.666672 +S 0753 0.135747 +M 0753 140.000000 +S 0754 0.135940 +M 0754 134.000000 +S 0755 0.136134 +M 0755 126.500000 +S 0756 0.136327 +M 0756 153.000000 +S 0757 0.136521 +M 0757 137.000000 +S 0758 0.136714 +M 0758 146.166672 +S 0759 0.136908 +M 0759 138.666672 +S 0760 0.137101 +M 0760 143.666672 +S 0761 0.137295 +M 0761 151.666672 +S 0762 0.137488 +M 0762 146.500000 +S 0763 0.137682 +M 0763 147.333328 +S 0764 0.137875 +M 0764 150.500000 +S 0765 0.138069 +M 0765 159.500000 +S 0766 0.138262 +M 0766 144.000000 +S 0767 0.138456 +M 0767 150.166672 +S 0768 0.138649 +M 0768 168.500000 +S 0769 0.138843 +M 0769 171.833328 +S 0770 0.139036 +M 0770 152.333328 +S 0771 0.139230 +M 0771 169.333328 +S 0772 0.139423 +M 0772 162.333328 +S 0773 0.139617 +M 0773 164.833328 +S 0774 0.139810 +M 0774 178.500000 +S 0775 0.140004 +M 0775 169.666672 +S 0776 0.140197 +M 0776 166.333328 +S 0777 0.140391 +M 0777 172.000000 +S 0778 0.140584 +M 0778 187.166672 +S 0779 0.140778 +M 0779 184.500000 +S 0780 0.140971 +M 0780 184.666672 +S 0781 0.141165 +M 0781 184.166672 +S 0782 0.141358 +M 0782 196.166672 +S 0783 0.141552 +M 0783 182.833328 +S 0784 0.141745 +M 0784 207.166672 +S 0785 0.141939 +M 0785 204.000000 +S 0786 0.142132 +M 0786 200.833328 +S 0787 0.142326 +M 0787 211.333328 +S 0788 0.142519 +M 0788 218.500000 +S 0789 0.142713 +M 0789 225.666672 +S 0790 0.142906 +M 0790 225.500000 +S 0791 0.143100 +M 0791 218.000000 +S 0792 0.143293 +M 0792 240.333328 +S 0793 0.143487 +M 0793 247.166672 +S 0794 0.143680 +M 0794 256.166656 +S 0795 0.143874 +M 0795 254.166672 +S 0796 0.144067 +M 0796 254.333328 +S 0797 0.144261 +M 0797 266.500000 +S 0798 0.144454 +M 0798 275.333344 +S 0799 0.144648 +M 0799 285.000000 +S 0800 0.144841 +M 0800 287.000000 +S 0801 0.145035 +M 0801 290.666656 +S 0802 0.145228 +M 0802 296.500000 +S 0803 0.145422 +M 0803 313.500000 +S 0804 0.145615 +M 0804 312.833344 +S 0805 0.145809 +M 0805 321.333344 +S 0806 0.146002 +M 0806 345.666656 +S 0807 0.146196 +M 0807 342.000000 +S 0808 0.146389 +M 0808 325.000000 +S 0809 0.146582 +M 0809 376.333344 +S 0810 0.146776 +M 0810 379.333344 +S 0811 0.146969 +M 0811 393.166656 +S 0812 0.147163 +M 0812 377.500000 +S 0813 0.147356 +M 0813 405.000000 +S 0814 0.147550 +M 0814 405.333344 +S 0815 0.147743 +M 0815 413.166656 +S 0816 0.147937 +M 0816 406.833344 +S 0817 0.148130 +M 0817 410.666656 +S 0818 0.148324 +M 0818 416.000000 +S 0819 0.148517 +M 0819 426.000000 +S 0820 0.148711 +M 0820 448.166656 +S 0821 0.148904 +M 0821 453.500000 +S 0822 0.149098 +M 0822 452.666656 +S 0823 0.149291 +M 0823 469.333344 +S 0824 0.149485 +M 0824 455.000000 +S 0825 0.149678 +M 0825 452.833344 +S 0826 0.149872 +M 0826 461.166656 +S 0827 0.150065 +M 0827 471.666656 +S 0828 0.150259 +M 0828 469.000000 +S 0829 0.150452 +M 0829 451.666656 +S 0830 0.150646 +M 0830 470.500000 +S 0831 0.150839 +M 0831 438.333344 +S 0832 0.151033 +M 0832 434.166656 +S 0833 0.151226 +M 0833 448.833344 +S 0834 0.151420 +M 0834 423.000000 +S 0835 0.151613 +M 0835 405.500000 +S 0836 0.151807 +M 0836 402.333344 +S 0837 0.152000 +M 0837 409.333344 +S 0838 0.152194 +M 0838 394.666656 +S 0839 0.152387 +M 0839 380.833344 +S 0840 0.152581 +M 0840 373.500000 +S 0841 0.152774 +M 0841 365.833344 +S 0842 0.152968 +M 0842 341.000000 +S 0843 0.153161 +M 0843 346.166656 +S 0844 0.153355 +M 0844 341.000000 +S 0845 0.153548 +M 0845 346.833344 +S 0846 0.153742 +M 0846 339.166656 +S 0847 0.153935 +M 0847 306.666656 +S 0848 0.154129 +M 0848 297.166656 +S 0849 0.154322 +M 0849 298.000000 +S 0850 0.154516 +M 0850 269.833344 +S 0851 0.154709 +M 0851 279.500000 +S 0852 0.154903 +M 0852 257.000000 +S 0853 0.155096 +M 0853 255.666672 +S 0854 0.155290 +M 0854 256.500000 +S 0855 0.155483 +M 0855 237.000000 +S 0856 0.155677 +M 0856 238.166672 +S 0857 0.155870 +M 0857 235.333328 +S 0858 0.156064 +M 0858 231.166672 +S 0859 0.156257 +M 0859 213.666672 +S 0860 0.156451 +M 0860 204.833328 +S 0861 0.156644 +M 0861 212.000000 +S 0862 0.156838 +M 0862 200.333328 +S 0863 0.157031 +M 0863 191.833328 +S 0864 0.157225 +M 0864 185.333328 +S 0865 0.157418 +M 0865 195.333328 +S 0866 0.157612 +M 0866 184.000000 +S 0867 0.157805 +M 0867 194.833328 +S 0868 0.157999 +M 0868 181.833328 +S 0869 0.158192 +M 0869 174.500000 +S 0870 0.158385 +M 0870 180.833328 +S 0871 0.158579 +M 0871 175.333328 +S 0872 0.158772 +M 0872 158.500000 +S 0873 0.158966 +M 0873 158.333328 +S 0874 0.159159 +M 0874 157.333328 +S 0875 0.159353 +M 0875 162.000000 +S 0876 0.159546 +M 0876 147.833328 +S 0877 0.159740 +M 0877 164.166672 +S 0878 0.159933 +M 0878 157.166672 +S 0879 0.160127 +M 0879 158.000000 +S 0880 0.160320 +M 0880 166.166672 +S 0881 0.160514 +M 0881 153.000000 +S 0882 0.160707 +M 0882 132.833328 +S 0883 0.160901 +M 0883 150.666672 +S 0884 0.161094 +M 0884 149.666672 +S 0885 0.161288 +M 0885 143.666672 +S 0886 0.161481 +M 0886 142.000000 +S 0887 0.161675 +M 0887 134.500000 +S 0888 0.161868 +M 0888 134.833328 +S 0889 0.162062 +M 0889 138.166672 +S 0890 0.162255 +M 0890 145.000000 +S 0891 0.162449 +M 0891 134.500000 +S 0892 0.162642 +M 0892 147.500000 +S 0893 0.162836 +M 0893 135.166672 +S 0894 0.163029 +M 0894 122.000000 +S 0895 0.163223 +M 0895 135.666672 +S 0896 0.163416 +M 0896 136.333328 +S 0897 0.163610 +M 0897 132.500000 +S 0898 0.163803 +M 0898 135.000000 +S 0899 0.163997 +M 0899 133.000000 +S 0900 0.164190 +M 0900 135.500000 +S 0901 0.164384 +M 0901 133.166672 +S 0902 0.164577 +M 0902 127.000000 +S 0903 0.164771 +M 0903 127.333336 +S 0904 0.164964 +M 0904 129.333328 +S 0905 0.165158 +M 0905 136.500000 +S 0906 0.165351 +M 0906 138.166672 +S 0907 0.165545 +M 0907 144.166672 +S 0908 0.165738 +M 0908 135.833328 +S 0909 0.165932 +M 0909 137.166672 +S 0910 0.166125 +M 0910 138.833328 +S 0911 0.166319 +M 0911 132.166672 +S 0912 0.166512 +M 0912 134.333328 +S 0913 0.166706 +M 0913 131.333328 +S 0914 0.166899 +M 0914 129.333328 +S 0915 0.167093 +M 0915 131.166672 +S 0916 0.167286 +M 0916 136.000000 +S 0917 0.167480 +M 0917 133.500000 +S 0918 0.167673 +M 0918 135.000000 +S 0919 0.167867 +M 0919 139.666672 +S 0920 0.168060 +M 0920 151.500000 +S 0921 0.168254 +M 0921 145.666672 +S 0922 0.168447 +M 0922 144.000000 +S 0923 0.168641 +M 0923 142.666672 +S 0924 0.168834 +M 0924 147.833328 +S 0925 0.169028 +M 0925 138.166672 +S 0926 0.169221 +M 0926 137.166672 +S 0927 0.169415 +M 0927 150.666672 +S 0928 0.169608 +M 0928 149.000000 +S 0929 0.169801 +M 0929 137.833328 +S 0930 0.169995 +M 0930 163.500000 +S 0931 0.170188 +M 0931 147.666672 +S 0932 0.170382 +M 0932 153.500000 +S 0933 0.170575 +M 0933 146.000000 +S 0934 0.170769 +M 0934 161.333328 +S 0935 0.170962 +M 0935 153.166672 +S 0936 0.171156 +M 0936 168.333328 +S 0937 0.171349 +M 0937 164.666672 +S 0938 0.171543 +M 0938 165.500000 +S 0939 0.171736 +M 0939 152.166672 +S 0940 0.171930 +M 0940 166.166672 +S 0941 0.172123 +M 0941 163.000000 +S 0942 0.172317 +M 0942 168.166672 +S 0943 0.172510 +M 0943 156.333328 +S 0944 0.172704 +M 0944 171.833328 +S 0945 0.172897 +M 0945 184.000000 +S 0946 0.173091 +M 0946 169.000000 +S 0947 0.173284 +M 0947 177.500000 +S 0948 0.173478 +M 0948 178.500000 +S 0949 0.173671 +M 0949 187.166672 +S 0950 0.173865 +M 0950 184.166672 +S 0951 0.174058 +M 0951 180.833328 +S 0952 0.174252 +M 0952 194.500000 +S 0953 0.174445 +M 0953 188.833328 +S 0954 0.174639 +M 0954 195.666672 +S 0955 0.174832 +M 0955 188.500000 +S 0956 0.175026 +M 0956 200.833328 +S 0957 0.175219 +M 0957 222.500000 +S 0958 0.175413 +M 0958 201.666672 +S 0959 0.175606 +M 0959 224.833328 +S 0960 0.175800 +M 0960 214.500000 +S 0961 0.175993 +M 0961 234.166672 +S 0962 0.176187 +M 0962 226.333328 +S 0963 0.176380 +M 0963 229.500000 +S 0964 0.176574 +M 0964 224.333328 +S 0965 0.176767 +M 0965 254.166672 +S 0966 0.176961 +M 0966 228.500000 +S 0967 0.177154 +M 0967 243.333328 +S 0968 0.177348 +M 0968 259.833344 +S 0969 0.177541 +M 0969 261.666656 +S 0970 0.177735 +M 0970 277.000000 +S 0971 0.177928 +M 0971 273.000000 +S 0972 0.178122 +M 0972 276.666656 +S 0973 0.178315 +M 0973 274.833344 +S 0974 0.178509 +M 0974 276.000000 +S 0975 0.178702 +M 0975 304.833344 +S 0976 0.178896 +M 0976 297.833344 +S 0977 0.179089 +M 0977 309.666656 +S 0978 0.179283 +M 0978 302.833344 +S 0979 0.179476 +M 0979 306.500000 +S 0980 0.179670 +M 0980 334.166656 +S 0981 0.179863 +M 0981 318.000000 +S 0982 0.180057 +M 0982 323.833344 +S 0983 0.180250 +M 0983 329.000000 +S 0984 0.180444 +M 0984 348.833344 +S 0985 0.180637 +M 0985 347.000000 +S 0986 0.180831 +M 0986 352.500000 +S 0987 0.181024 +M 0987 363.666656 +S 0988 0.181217 +M 0988 358.166656 +S 0989 0.181411 +M 0989 364.833344 +S 0990 0.181604 +M 0990 384.500000 +S 0991 0.181798 +M 0991 371.333344 +S 0992 0.181991 +M 0992 393.500000 +S 0993 0.182185 +M 0993 393.166656 +S 0994 0.182378 +M 0994 390.333344 +S 0995 0.182572 +M 0995 377.666656 +S 0996 0.182765 +M 0996 379.000000 +S 0997 0.182959 +M 0997 402.166656 +S 0998 0.183152 +M 0998 415.833344 +S 0999 0.183346 +M 0999 417.666656 +S 1000 0.183539 +M 1000 416.500000 +S 1001 0.183733 +M 1001 409.666656 +S 1002 0.183926 +M 1002 398.166656 +S 1003 0.184120 +M 1003 397.000000 +S 1004 0.184313 +M 1004 409.666656 +S 1005 0.184507 +M 1005 399.333344 +S 1006 0.184700 +M 1006 410.500000 +S 1007 0.184894 +M 1007 391.166656 +S 1008 0.185087 +M 1008 401.833344 +S 1009 0.185281 +M 1009 413.500000 +S 1010 0.185474 +M 1010 413.500000 +S 1011 0.185668 +M 1011 391.833344 +S 1012 0.185861 +M 1012 394.666656 +S 1013 0.186055 +M 1013 380.000000 +S 1014 0.186248 +M 1014 379.166656 +S 1015 0.186442 +M 1015 368.666656 +S 1016 0.186635 +M 1016 388.166656 +S 1017 0.186829 +M 1017 392.666656 +S 1018 0.187022 +M 1018 374.000000 +S 1019 0.187216 +M 1019 364.500000 +S 1020 0.187409 +M 1020 345.833344 +S 1021 0.187603 +M 1021 353.000000 +S 1022 0.187796 +M 1022 308.000000 +S 1023 0.187990 +M 1023 325.666656 +S 1024 0.188183 +M 1024 321.333344 +S 1025 0.188377 +M 1025 305.000000 +S 1026 0.188570 +M 1026 319.666656 +S 1027 0.188764 +M 1027 310.166656 +S 1028 0.188957 +M 1028 297.666656 +S 1029 0.189151 +M 1029 290.333344 +S 1030 0.189344 +M 1030 282.000000 +S 1031 0.189538 +M 1031 273.666656 +S 1032 0.189731 +M 1032 277.000000 +S 1033 0.189925 +M 1033 267.833344 +S 1034 0.190118 +M 1034 272.000000 +S 1035 0.190312 +M 1035 253.500000 +S 1036 0.190505 +M 1036 231.166672 +S 1037 0.190699 +M 1037 255.666672 +S 1038 0.190892 +M 1038 240.000000 +S 1039 0.191086 +M 1039 239.500000 +S 1040 0.191279 +M 1040 228.833328 +S 1041 0.191473 +M 1041 220.166672 +S 1042 0.191666 +M 1042 210.833328 +S 1043 0.191860 +M 1043 217.666672 +S 1044 0.192053 +M 1044 214.166672 +S 1045 0.192247 +M 1045 196.833328 +S 1046 0.192440 +M 1046 206.500000 +S 1047 0.192633 +M 1047 203.833328 +S 1048 0.192827 +M 1048 177.833328 +S 1049 0.193020 +M 1049 193.833328 +S 1050 0.193214 +M 1050 187.166672 +S 1051 0.193407 +M 1051 187.833328 +S 1052 0.193601 +M 1052 176.000000 +S 1053 0.193794 +M 1053 178.500000 +S 1054 0.193988 +M 1054 182.833328 +S 1055 0.194181 +M 1055 175.500000 +S 1056 0.194375 +M 1056 179.500000 +S 1057 0.194568 +M 1057 166.666672 +S 1058 0.194762 +M 1058 167.166672 +S 1059 0.194955 +M 1059 163.166672 +S 1060 0.195149 +M 1060 148.000000 +S 1061 0.195342 +M 1061 143.000000 +S 1062 0.195536 +M 1062 159.333328 +S 1063 0.195729 +M 1063 144.500000 +S 1064 0.195923 +M 1064 146.500000 +S 1065 0.196116 +M 1065 140.166672 +S 1066 0.196310 +M 1066 147.000000 +S 1067 0.196503 +M 1067 134.000000 +S 1068 0.196697 +M 1068 145.500000 +S 1069 0.196890 +M 1069 142.166672 +S 1070 0.197084 +M 1070 137.000000 +S 1071 0.197277 +M 1071 123.833336 +S 1072 0.197471 +M 1072 120.500000 +S 1073 0.197664 +M 1073 134.000000 +S 1074 0.197858 +M 1074 123.000000 +S 1075 0.198051 +M 1075 119.666664 +S 1076 0.198245 +M 1076 120.333336 +S 1077 0.198438 +M 1077 128.500000 +S 1078 0.198632 +M 1078 121.333336 +S 1079 0.198825 +M 1079 117.833336 +S 1080 0.199019 +M 1080 107.666664 +S 1081 0.199212 +M 1081 113.833336 +S 1082 0.199406 +M 1082 118.166664 +S 1083 0.199599 +M 1083 116.500000 +S 1084 0.199793 +M 1084 106.333336 +S 1085 0.199986 +M 1085 102.000000 +S 1086 0.200180 +M 1086 97.500000 +S 1087 0.200373 +M 1087 100.833336 +S 1088 0.200567 +M 1088 101.500000 +S 1089 0.200760 +M 1089 98.666664 +S 1090 0.200954 +M 1090 107.333336 +S 1091 0.201147 +M 1091 95.000000 +S 1092 0.201341 +M 1092 93.166664 +S 1093 0.201534 +M 1093 95.000000 +S 1094 0.201728 +M 1094 99.000000 +S 1095 0.201921 +M 1095 94.000000 +S 1096 0.202115 +M 1096 102.833336 +S 1097 0.202308 +M 1097 101.333336 +S 1098 0.202502 +M 1098 97.333336 +S 1099 0.202695 +M 1099 92.000000 +S 1100 0.202889 +M 1100 97.000000 +S 1101 0.203082 +M 1101 95.000000 +S 1102 0.203276 +M 1102 97.333336 +S 1103 0.203469 +M 1103 98.833336 +S 1104 0.203663 +M 1104 103.000000 +S 1105 0.203856 +M 1105 87.000000 +S 1106 0.204049 +M 1106 94.500000 +S 1107 0.204243 +M 1107 93.333336 +S 1108 0.204436 +M 1108 89.666664 +S 1109 0.204630 +M 1109 99.166664 +S 1110 0.204823 +M 1110 97.833336 +S 1111 0.205017 +M 1111 94.833336 +S 1112 0.205210 +M 1112 91.833336 +S 1113 0.205404 +M 1113 90.000000 +S 1114 0.205597 +M 1114 100.166664 +S 1115 0.205791 +M 1115 92.333336 +S 1116 0.205984 +M 1116 95.333336 +S 1117 0.206178 +M 1117 89.333336 +S 1118 0.206371 +M 1118 89.500000 +S 1119 0.206565 +M 1119 84.166664 +S 1120 0.206758 +M 1120 84.166664 +S 1121 0.206952 +M 1121 94.166664 +S 1122 0.207145 +M 1122 80.333336 +S 1123 0.207339 +M 1123 82.500000 +S 1124 0.207532 +M 1124 83.166664 +S 1125 0.207726 +M 1125 74.000000 +S 1126 0.207919 +M 1126 96.500000 +S 1127 0.208113 +M 1127 93.666664 +S 1128 0.208306 +M 1128 79.000000 +S 1129 0.208500 +M 1129 84.666664 +S 1130 0.208693 +M 1130 76.500000 +S 1131 0.208887 +M 1131 90.500000 +S 1132 0.209080 +M 1132 79.000000 +S 1133 0.209274 +M 1133 84.000000 +S 1134 0.209467 +M 1134 80.500000 +S 1135 0.209661 +M 1135 83.000000 +S 1136 0.209854 +M 1136 81.333336 +S 1137 0.210048 +M 1137 74.833336 +S 1138 0.210241 +M 1138 69.666664 +S 1139 0.210435 +M 1139 71.666664 +S 1140 0.210628 +M 1140 80.333336 +S 1141 0.210822 +M 1141 70.500000 +S 1142 0.211015 +M 1142 69.833336 +S 1143 0.211209 +M 1143 66.000000 +S 1144 0.211402 +M 1144 69.500000 +S 1145 0.211596 +M 1145 69.833336 +S 1146 0.211789 +M 1146 75.333336 +S 1147 0.211983 +M 1147 65.666664 +S 1148 0.212176 +M 1148 62.166668 +S 1149 0.212370 +M 1149 64.000000 +S 1150 0.212563 +M 1150 69.000000 +S 1151 0.212757 +M 1151 72.500000 +S 1152 0.212950 +M 1152 59.833332 +S 1153 0.213144 +M 1153 61.166668 +S 1154 0.213337 +M 1154 62.500000 +S 1155 0.213531 +M 1155 68.000000 +S 1156 0.213724 +M 1156 61.166668 +S 1157 0.213918 +M 1157 50.333332 +S 1158 0.214111 +M 1158 63.500000 +S 1159 0.214305 +M 1159 59.833332 +S 1160 0.214498 +M 1160 60.166668 +S 1161 0.214692 +M 1161 51.000000 +S 1162 0.214885 +M 1162 54.666668 +S 1163 0.215079 +M 1163 59.333332 +S 1164 0.215272 +M 1164 47.166668 +S 1165 0.215466 +M 1165 55.000000 +S 1166 0.215659 +M 1166 51.166668 +S 1167 0.215852 +M 1167 55.333332 +S 1168 0.216046 +M 1168 56.833332 +S 1169 0.216239 +M 1169 60.333332 +S 1170 0.216433 +M 1170 52.000000 +S 1171 0.216626 +M 1171 55.666668 +S 1172 0.216820 +M 1172 49.500000 +S 1173 0.217013 +M 1173 50.333332 +S 1174 0.217207 +M 1174 55.333332 +S 1175 0.217400 +M 1175 44.666668 +S 1176 0.217594 +M 1176 52.000000 +S 1177 0.217787 +M 1177 47.833332 +S 1178 0.217981 +M 1178 46.666668 +S 1179 0.218174 +M 1179 42.833332 +S 1180 0.218368 +M 1180 54.500000 +S 1181 0.218561 +M 1181 52.666668 +S 1182 0.218755 +M 1182 49.000000 +S 1183 0.218948 +M 1183 47.166668 +S 1184 0.219142 +M 1184 54.166668 +S 1185 0.219335 +M 1185 50.000000 +S 1186 0.219529 +M 1186 51.333332 +S 1187 0.219722 +M 1187 50.333332 +S 1188 0.219916 +M 1188 49.500000 +S 1189 0.220109 +M 1189 49.166668 +S 1190 0.220303 +M 1190 47.000000 +S 1191 0.220496 +M 1191 38.833332 +S 1192 0.220690 +M 1192 51.333332 +S 1193 0.220883 +M 1193 52.666668 +S 1194 0.221077 +M 1194 48.000000 +S 1195 0.221270 +M 1195 45.833332 +S 1196 0.221464 +M 1196 47.666668 +S 1197 0.221657 +M 1197 45.000000 +S 1198 0.221851 +M 1198 45.666668 +S 1199 0.222044 +M 1199 42.833332 +S 1200 0.222238 +M 1200 49.166668 +S 1201 0.222431 +M 1201 45.000000 +S 1202 0.222625 +M 1202 45.166668 +S 1203 0.222818 +M 1203 42.500000 +S 1204 0.223012 +M 1204 45.666668 +S 1205 0.223205 +M 1205 44.000000 +%%%%% ENDE # end of file diff --git a/Python code/materials.json b/Python code/data/materials.json similarity index 100% rename from Python code/materials.json rename to Python code/data/materials.json diff --git a/Python code/dielectrics.py b/Python code/dielectrics.py index 6438a9f122ec87422341e1cfd4acb566b8850313..959e99a9482349742718bcdbe4fa2fc08ba685b8 100755 --- a/Python code/dielectrics.py +++ b/Python code/dielectrics.py @@ -51,12 +51,12 @@ def plotDielectrics(x,eps, title=" ", plot_show=True): ''' fig, axs = plt.subplots(3, 1, sharex=True) fig.suptitle(str(title), fontsize=16) - axs[0].plot(x, np.imag(eps), label='Im( $\epsilon (\omega )$ )') - axs[0].plot(x, np.real(eps), label='Re( $\epsilon (\omega )$ )') + axs[0].plot(x, np.imag(eps), label='Im( $\\epsilon (\\omega )$ )') + axs[0].plot(x, np.real(eps), label='Re( $\\epsilon (\\omega )$ )') axs[0].legend() axs[0].set_ylabel('Dielectric Function') - axs[1].plot(x, surfaceLoss(eps), lineStyle='-', label='Loss function') + axs[1].plot(x, surfaceLoss(eps), linestyle='-', label='Loss function') axs[1].set_ylabel('Surface Loss Function') axs[1].set_ylim([-0.10,None]) diff --git a/Python code/myEels3cole.f90 b/Python code/myEels20.f90 similarity index 98% rename from Python code/myEels3cole.f90 rename to Python code/myEels20.f90 index c250a8e1403756f0be2175e81560656d56414201..28d953996d2b29f4bbc22a1ef84c8b9f790ae06c 100755 --- a/Python code/myEels3cole.f90 +++ b/Python code/myEels20.f90 @@ -12,7 +12,7 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, & ! * ! * It is based on the work of Lambin'90 and modified for the use * ! * within python * -! * (KMS and WFW, Martin-Luther-Universität Halle-Wittenberg) * * +! * (KMS and WFW, Martin-Luther-Universitaet Halle-Wittenberg) * * ! ****************************************************************** implicit none @@ -20,10 +20,10 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, & integer, parameter :: nt = 5 double precision, intent(in) :: e0, theta, phia, phib, wmin, wmax, dw - double precision, intent(in) :: thick(name_size), epsinf(name_size), osc(3, osc_size) - character*10, intent(in) :: name(name_size) - character*10, intent(in) :: contrl, mode integer, intent(in) :: name_size, osc_size, neps + double precision, intent(in) :: thick(name_size), epsinf(name_size), osc(3, osc_size) + character, intent(in) :: name(10, name_size) + character, intent(in) :: contrl(10), mode(10) integer, intent(in) :: wn_array_size integer, intent(in) :: layers, nper, nos(name_size) double precision, intent(out) :: wn_array(wn_array_size), f_array(wn_array_size) @@ -44,11 +44,11 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, & ru, um, dlimf, wn, user, ration common / mulayr / argmin, argmax, epsmac, flag + data aerr / 0.0d0 /, rerr / 1.0d-06 /, f / 0.0d0 /, f1 / 0.0d0 / + !if (present(debug)) debug = .True. debugFirstRun = .True. - data aerr / 0.0d0 /, rerr / 1.0d-06 /, f / 0.0d0 /, f1 / 0.0d0 / - ! *** machine-dependent constants ! *** epsmac + 1.0 = epsmac , cosh(argmin) = 1.0 , tanh(argmax) = 1.0 flag = 1 @@ -234,7 +234,9 @@ subroutine quanc8(fun, a, b, abserr, relerr, result, errest, nofun, flag, eps, d ! left to do when the limit on nofun was approached. implicit none + external fun + double precision :: fun double precision, intent(in) :: a double precision, intent(in) :: b @@ -552,9 +554,9 @@ double precision function fint1(u, eps, d, layers, nper, eps_size) double precision, intent(in) :: u + integer, intent(in) :: layers, nper, eps_size double precision, intent(in) :: d(eps_size) double complex, intent(in) :: eps(eps_size) - integer, intent(in) :: layers, nper, eps_size logical :: ration, user double precision :: acoef, bcoef, ccoef, cospsi, den, dif, dlimf, e, elleps @@ -611,9 +613,9 @@ double precision function fint2(u, eps, d, layers, nper, eps_size) double precision, intent(in) :: u + integer, intent(in) :: layers, nper, eps_size double precision, intent(in) :: d(eps_size) double complex, intent(in) :: eps(eps_size) - integer, intent(in) :: layers, nper, eps_size logical :: ration, user double precision :: a, arg, b, b2, c, ccoef, cospsi, dlimf, elleps, phi @@ -672,9 +674,9 @@ double precision function fint3(u, eps, d, layers, nper, eps_size) double precision, intent(in) :: u + integer, intent(in) :: layers, nper, eps_size double precision, intent(in) :: d(eps_size) double complex, intent(in) :: eps(eps_size) - integer, intent(in) :: layers, nper, eps_size logical :: ration, user double precision :: a, arg, b, ccoef, cospsi, dlimf, elleps, phi1, phi2 @@ -949,7 +951,8 @@ subroutine seteps(nos, osc_size, osc, epsinf, wn, nLayer, eps) ! * the present frequency wn (in cm**-1) * ! ****************************************************************** - ! implicit none + implicit none + integer, intent(in) :: nLayer integer, dimension(nLayer),intent(in) :: nos integer, intent(in) :: osc_size @@ -962,8 +965,9 @@ subroutine seteps(nos, osc_size, osc, epsinf, wn, nLayer, eps) double complex :: nomi, deno, addeps double precision :: wn2, b - integer :: flag + integer j, k, l, m logical debugFirstRun + common /control/ debugFirstRun j = 0 @@ -982,7 +986,7 @@ subroutine seteps(nos, osc_size, osc, epsinf, wn, nLayer, eps) deno =deno * (osc(1,j)**2 - wn * dcmplx( wn, osc(3,j) ) ) else if (osc(1,j) < 0.) then! Negative TO mode means: _Additive_ Lorentz oscillator with Q - addeps = addeps + osc(1,j)**2 * osc(2,j) /dcmplx(osc(1,j)**2 + wn2, wn*osc(3,j)) + addeps = addeps + osc(1,j)**2 * osc(2,j) /dcmplx(osc(1,j)**2 - wn2, -1*wn*osc(3,j)) ! Sign of imaginary part changed (WFW) else ! osc(1,j) = 0 -> it is a Drude term addeps = addeps - dcmplx(osc(1,j+m)**2, wn*(osc(3,j)-osc(3,j+m))) /dcmplx(wn2, wn*osc(3,j)) @@ -990,10 +994,10 @@ subroutine seteps(nos, osc_size, osc, epsinf, wn, nLayer, eps) enddo j = j+m ! we have already looped over the LO modes, therefore increase the index - eps(l) = epsinf(l) * nomi / deno + addeps + eps(l) = epsinf(l) * (nomi / deno + addeps) ! brackets changed by HHe 230915 enddo debugFirstRun = .false. return end subroutine seteps -end module mod_doeels \ No newline at end of file +end module mod_doeels diff --git a/Python code/output/I2L02 .png b/Python code/output/I2L02 .png new file mode 100644 index 0000000000000000000000000000000000000000..d1f637fcd4de3ff16debb3c9d12319385fbe70ce Binary files /dev/null and b/Python code/output/I2L02 .png differ diff --git a/Python code/output/dielectrics.png b/Python code/output/dielectrics.png new file mode 100644 index 0000000000000000000000000000000000000000..a7375762b536e69f0e0a639057d52cc5561710c7 Binary files /dev/null and b/Python code/output/dielectrics.png differ