From 2d8560f4eb78e2c169d8c94631e1a34b3309b69c Mon Sep 17 00:00:00 2001 From: RJbalikian <46536937+RJbalikian@users.noreply.github.com> Date: Wed, 18 Oct 2023 00:39:01 -0400 Subject: [PATCH] updates --- .../file_test.cpython-310-pytest-7.4.2.pyc | Bin 0 -> 665 bytes .../run_test.cpython-310-pytest-7.4.2.pyc | Bin 0 -> 638 bytes conda/meta.yaml | 4 +- docs/generate_docs.py | 20 +++++++- docs/main.html | 4 +- docs/sprit_hvsr.html | 18 +++---- pyproject.toml | 2 +- run_test.py | 10 ++++ setup.py | 2 +- sprit/__pycache__/sprit_hvsr.cpython-310.pyc | Bin 166583 -> 166588 bytes sprit/sprit_hvsr.py | 44 +++++++++--------- 11 files changed, 66 insertions(+), 38 deletions(-) create mode 100644 __pycache__/file_test.cpython-310-pytest-7.4.2.pyc create mode 100644 __pycache__/run_test.cpython-310-pytest-7.4.2.pyc create mode 100644 run_test.py diff --git a/__pycache__/file_test.cpython-310-pytest-7.4.2.pyc b/__pycache__/file_test.cpython-310-pytest-7.4.2.pyc new file mode 100644 index 0000000000000000000000000000000000000000..9ae29f4581998829665c71da6f540066fb8098e5 GIT binary patch literal 665 zcmYjO&5qMB5VoDfNxN-1AR#z`5J>bu7cL+ngqDAh&=zffdSe zYe^Lmgi#ge3EoC%PGu^{83sQR@pqJG0*}!!`3<>fKTpkE>rx+2gL^Ka0fqA%j-7xy z?Krf?0&hYYvO-5lQw)tl(G~iHmdNgkSR{fvvW_< ziL>(wO?VuCjcvypZI=byMn-;toxn{0A)o%kOa?h74v7gYPP?{h%lm>Fgn0J2s>NK% OCm1{%dfh}B>5;#1o3k(g literal 0 HcmV?d00001 diff --git a/__pycache__/run_test.cpython-310-pytest-7.4.2.pyc b/__pycache__/run_test.cpython-310-pytest-7.4.2.pyc new file mode 100644 index 0000000000000000000000000000000000000000..9a5fba29dacdef0c876e7b756f6ad908a5a3b5b2 GIT binary patch literal 638 zcmYjO%We}f6utILoJpF51rmZ42!TWvq_AOu5GrqwP=%TT>S!d(iCsF9=ONnzq)HY% zJ|x*7@k`#a;ulzPJ)x)``CQ*~$H&Ll)o|Dc9AA!h)enrnLUMjkO!jcv34#PkG^C`B zw1_AWDB`q3wiT8#NoDjM(yk;^7{|YH5A4qao{gZ#a6Z9lX9%Y~2Wuk9hM)>v!Fwnf zL0!zE6&%45?5>PuM>0p((FQ((z3Sq1X9*i(Zs3ZARY@ev$t743*F`Q_qTV<16`i7w z?!ISM6tz}|&re6fT4mg|o44x4ov}?L41UKc>v$bXa@_MBTbtbZ$TX!N?9~fCYjW-K z(t4(=xzN`4k!|+~<1cf2QZ-s~qiVxH=xQc3FGZnzk6Ur90_{hS+GF#od}d5#{6&6T znL;?OKG#~5!eM1hO@tjucfyg|^4*_chRFT7H@}~~MdLOzxmJtpb@1_2IFY@`-AOac zrf&{j9o~NVZhDZRZ63@`>V;2Qm&gQVJ3yeMPZGlLdrZv*B*uwP&Ms)e+xRPdd%MZD vS7>cw)F*5SUgParameters dataIN.merge() params['batch'] = False #Set False by default, will get corrected later in batch mode - params['input_stream'] + params['input_stream'] = dataIN.copy() params['stream'] = dataIN.copy() params['ProcessingStatus']['FetchDataStatus'] = True if verbose and not isinstance(params, HVSRBatch): @@ -2365,7 +2365,7 @@

Returns

-def input_params(datapath, site='HVSR Site', network='AM', station='RAC84', loc='00', channels=['EHZ', 'EHN', 'EHE'], acq_date='2023-10-17', starttime='00:00:00.00', endtime='23:59:59.999999', tzone='UTC', xcoord=-88.2290526, ycoord=40.1012122, elevation=755, input_crs='EPSG:4326', output_crs='EPSG:4326', elev_unit='feet', depth=0, instrument='Raspberry Shake', metapath='', hvsr_band=[0.4, 40], peak_freq_range=[0.4, 40], verbose=False) +def input_params(datapath, site='HVSR Site', network='AM', station='RAC84', loc='00', channels=['EHZ', 'EHN', 'EHE'], acq_date='2023-10-18', starttime='00:00:00.00', endtime='23:59:59.999999', tzone='UTC', xcoord=-88.2290526, ycoord=40.1012122, elevation=755, input_crs='EPSG:4326', output_crs='EPSG:4326', elev_unit='feet', depth=0, instrument='Raspberry Shake', metapath='', hvsr_band=[0.4, 40], peak_freq_range=[0.4, 40], verbose=False)

Function for designating input parameters for reading in and processing data

diff --git a/docs/sprit_hvsr.html b/docs/sprit_hvsr.html index 5122e3c..649902b 100644 --- a/docs/sprit_hvsr.html +++ b/docs/sprit_hvsr.html @@ -1269,7 +1269,7 @@

Module sprit.sprit_hvsr

dataIN.merge() params['batch'] = False #Set False by default, will get corrected later in batch mode - params['input_stream'] + params['input_stream'] = dataIN.copy() params['stream'] = dataIN.copy() params['ProcessingStatus']['FetchDataStatus'] = True if verbose and not isinstance(params, HVSRBatch): @@ -6047,6 +6047,7 @@

Module sprit.sprit_hvsr

return _peak + #Check stability def __check_stability(_stdf, _peak, _hvsr_log_std, rank): """Test peaks for satisfying stability conditions as outlined by SESAME 2004 @@ -6189,6 +6190,7 @@

Module sprit.sprit_hvsr

return _peak + #Get frequency standard deviation def __get_stdf(x_values, indexList, hvsrPeaks): """Private function to get frequency standard deviation, from multiple time-step HVSR curves""" @@ -6201,17 +6203,17 @@

Module sprit.sprit_hvsr

if p is None: p = hvsrPeaks.iloc[j][k] else: - #Find closest peak in current time to (current) main hvsr peak + # Find closest peak in current time to (current) hvsr peak if abs(index - hvsrPeaks.iloc[j][k]) < abs(index - p): p = hvsrPeaks.iloc[j][k] - #p = hvsrPeaks[j][k] - #print(p=p1, p, p1) + # p = hvsrPeaks[j][k] + # print(p=p1, p, p1) if p is not None: point.append(p) point.append(index) v = list() - for l in range(len(point)): - v.append(x_values[point[l]]) + for pl in range(len(point)): + v.append(x_values[point[pl]]) stdf.append(np.std(v)) return stdf @@ -7096,7 +7098,7 @@

Parameters

dataIN.merge() params['batch'] = False #Set False by default, will get corrected later in batch mode - params['input_stream'] + params['input_stream'] = dataIN.copy() params['stream'] = dataIN.copy() params['ProcessingStatus']['FetchDataStatus'] = True if verbose and not isinstance(params, HVSRBatch): @@ -7963,7 +7965,7 @@

Returns

-def input_params(datapath, site='HVSR Site', network='AM', station='RAC84', loc='00', channels=['EHZ', 'EHN', 'EHE'], acq_date='2023-10-17', starttime='00:00:00.00', endtime='23:59:59.999999', tzone='UTC', xcoord=-88.2290526, ycoord=40.1012122, elevation=755, input_crs='EPSG:4326', output_crs='EPSG:4326', elev_unit='feet', depth=0, instrument='Raspberry Shake', metapath='', hvsr_band=[0.4, 40], peak_freq_range=[0.4, 40], verbose=False) +def input_params(datapath, site='HVSR Site', network='AM', station='RAC84', loc='00', channels=['EHZ', 'EHN', 'EHE'], acq_date='2023-10-18', starttime='00:00:00.00', endtime='23:59:59.999999', tzone='UTC', xcoord=-88.2290526, ycoord=40.1012122, elevation=755, input_crs='EPSG:4326', output_crs='EPSG:4326', elev_unit='feet', depth=0, instrument='Raspberry Shake', metapath='', hvsr_band=[0.4, 40], peak_freq_range=[0.4, 40], verbose=False)

Function for designating input parameters for reading in and processing data

diff --git a/pyproject.toml b/pyproject.toml index 1ceac90..fee0a10 100644 --- a/pyproject.toml +++ b/pyproject.toml @@ -7,7 +7,7 @@ name = "sprit" authors = [{name="Riley Balikian"}, {name="Hongyu Xaio"}] dynamic = ["readme"] license = {file = "LICENSE"} -version="0.1.48" +version="0.1.49" description = "A package for processing and analyzing HVSR (Horizontal to Vertical Spectral Ratio) data" keywords = ["HVSR", "seismic", "horizontal to vertical spectral ratio", "obspy", 'geology', 'geophysics', 'geotechnical'] requires-python = ">=3.9" diff --git a/run_test.py b/run_test.py new file mode 100644 index 0000000..ed7c958 --- /dev/null +++ b/run_test.py @@ -0,0 +1,10 @@ +import sprit + +def test_run(): + try: + sprit.run("sample") + test_passed = True + except: + test_passed = False + + assert test_passed \ No newline at end of file diff --git a/setup.py b/setup.py index bfdebfd..90b48be 100644 --- a/setup.py +++ b/setup.py @@ -7,7 +7,7 @@ name="sprit", author= "Riley Balikian", author_email = "balikian@illinois.edu", - version="0.1.48", + version="0.1.49", package_data={'sprit': ['resources/*', 'resources/icon/*', 'resources/themes/*', 'resources/themes/forest-dark/*', 'resources/themes/forest-light/*', 'resources/sample_data/*',]}, long_description_content_type="text/markdown", long_description=long_description, diff --git a/sprit/__pycache__/sprit_hvsr.cpython-310.pyc b/sprit/__pycache__/sprit_hvsr.cpython-310.pyc index 111651e5a9a7a67eb993922fa868170852336271..d484c293d2945c77931d2ae9fc7578ef51c6de9c 100644 GIT binary patch delta 31816 zcma)F31F1P*=D}%ZtjEoP#}OPQMm#L3W^9xMIc5J6+Bor*-f&r*HJIKmX!uY&W3Qnf>=1w*NY#$8o14NbL`B-M%&(SD7i8}Oq~B4CDIN6-jhDq1@t zty_c9cq$!=w?_Cj@`KRWZbs(qjvc-)-J6UI zuVrU{;{i{oWu~+-mPiNJN0JGh#%!>MrZ~wwqN}FBc6Ke>JA)vvGyUr)&vGhuG511% zh8n-#=w3Ya;99}O&?lx)C$_DPrqZT}^{$Gxrt@pz=U2FIO}%Ues|I$izz=$ot0b>U z_sZ&%4+*hAf$?FAwG32x{RLdPWz$YSlx0!psr<2HF$G&Ay~f`iigoRS6!)rWr%Tnh zD^u7WNr%#DoyDNO9-rs;*&RK7aZ$oneU7_+`hpVDquvBy+$*M!{ZbDfiA$L+y^|?H z6bGljO7Jy;8wtKmkQ=0=xB_&O`@!_Xe4JmBdjH%RPZTVz;|rk;j;SVK@wA({(#so= zGHR<)-Fa0?ZJXSAhpd_Xb(SPn^fw5uBgnNa1=%U3kbjeV+aU{P-^?QN#Vt&Ii(uFn z9D$Ete5b!?)+klS(&$<}evA`On!lm~Do1l|mSf}FdA8QU(33uwDlU1_+=tFlqHD{q% zgYIJCZTtCxvo@Z?=M{iT+2q!c5j|HX6iYptG#ho!g12N8G7PEkqQ}$D}=CZxa z)|ou&x$_fiB0E-dHuf-u_9;7OGLKN)ao;Gb7Dq?(Rbw$&rFdPBVBzomj>>fq!UD2RUB(mJB(KV}{ zFE)WrkXA$+_nx|i>KE=ib*CxguBab>AjfV>qOFN|OEldXN}6J(lO2h4!no=BWk;H# z&QLTSOC*xusc*B7zu+fMdaHhm2`L}6l1X+;`z1>?y6@Mo5B#0Uz5Oe`^rP`b91+&) zU2bu3+SomOaW_G3$c^su;6k<5O$X0M@BSKGbml(>3*25{O4=gn){daaP_jWdocjWT zyejyCA9=OlTsnXR=#E?oXis$M)`%`bUf)JeqI+X-s(Z)g#(-_*>d=@a|GwU%hU#^+ zW@D_+xNAd)mXsmYEEDR6E}fLE)x9}1uIArYyWUpoz8R`p%h6zy@v-(z1pof?CVc(} z{2Or)g%{z%?3do_CR^689Ks^xmNF1yYb2Et0qZ}gxii}_f9r{>^6M$J@^LZ}+M4Gf*Qw^Q8)yLPj$u8DVYF$3GwtPK{dxsX@@Ok@>}9| z^B=Ty%V0~U@c2`L6j1Xd>jdh++{n!5g;J?rUBRwb5eP4ny;b<e4Qj`s~ND32` zTNRnA-gDH|k6MPRDhH0b)k+KXh;&) zi}`vDtMi58y=#o>F{&5bHXe@j1Vt2G&k{{i0g1X%mJ>HW&2*q8|yS8Bq`&SfdJ&VjZqCX(gKIABiDW?Wy z+5^#45E8a2+LDOIxjX3D9BP!s0=7oIl^=fUzScV(6!!PtL&j5(h6t5r{?MAjPO_|T z_iC9;AzhhUSsa{J>`MDvjn}cYGguL6Lw|#hvYU3h3${;F-R{ZTXRBEMdD~0M7tOqR z@hv0H9ueNWY{c1v>Cmzf;pNC{nK&Zciu|H01N!cC(b*%?2cPF$8A#pj{_x@%<8MuG zID5o->KuICG@_=xL8BQ`>|ZWE;!G|mQxJ`JYl>b*og?easO^S*u3*6`f`tT&2#z3F z3{X@2uivAEOLk3C-|GLbU5}~=>&Cz?=$wKf9 zx8Tx+g*T!gYI2Xc^ytA|Z|9}s|ASp`4E8u|((aEhow1bCJ9s-52Q)Fe%L@tAB41Pe zz()QfG-OI-sS5Y8M6$r0df7(stKQ3Im70RO<4+@KFz%g~&75-^yLLOl9sIVmGuo;l zJ3+(=gN^U|R=?!_VEgHmrpzTLRMJN7Hv09sK*pj>j{z>n&Q^5)22$}GaI;$eGCE^8|5DQ=N6nqZCSN*Pt6*YPXiwHqy~;p*)O7h&d+>l9qlwFLH-DR8=oNT9#^?dYL|PmdwNDBH<`r& z_K%d7M`B!OW~6MM`%IdU9@meZl$-dK11ZLo+p*A4(2&4Bkb5)=6B;!tZ2aLyS`~e? zbYl60p-GT*jC<5Ihqs#&i*Fz!1uIDl5pODh-(-**d{D4Xt}Q!|ZVc&AXUY^~H_~>9 zB|5aNNp#v2?1P{Wenw7}%@u+(?sM1F9I}wDLaYqMcN}+=wXi%pqLm7AS zwc|6U43wYCE4PK};$bOL(ix2prBN`O)_-6Nd4qe2=7n zc2>cN5&iwwj&~-ODgJMis!&Dz=lC2e!L|4Q+IzOL%9N$HMw7urJlGyerkaiW@U6Fq z4g|&%aXW6CI&S~Ei+jUurySYH=_?kM8WxZUws4|0g5B0imf)1Wwyt=qrsSXBP51EI zH^B~l$?a?AuR(9EF0huVtn38Nq$dy&1<6pKUGF0?f!GV)bPMk|DDVR5=w)~Q9rddQ zkVHd8kY{6|mRLmJCmEW?HskAzbf)rK#yeB8j8aW&qBD8S>3-A$-u}QHv#`6reaD*d zL3VXBL5P5l4SJI5umA3es@OQe)9|m(xHsOp`KUjl1)%94>^HfHp2;r`27rcJB{Z48 z2xW^G(C99?Ys>7-=!+N)VOFoqu~RFqY-LsXU;NKq>l;i3S97Y-Do2qpw5UIGl!bP$ zNSRUlz7K+AzL2-Baqas%PEuGqiQVJJoy>l3LSyIz#5OC6H@nyGp0jElyC9elW;#1h zaK(-Ysg@uo@C?a1R;du4xY4cHGd%;7riO%{(xXwMDF}tbdIP`ZKCQn=a1Vg-i~Unr zvk7#?qi4bRd^z71@-zRojtEGMQUnD00xX`K9IUVAC_X2+k~I~B@rCM1il{U>JsIn@ zlXM6xa@w(GVa+Rsg#(@w2)SCjFYc+%9K`QyP^9K4tBhCbNN1uu5`^6~79~yh;-%d} zPi3kHADBtB@C2=z5KJXC6s2&aEhI|cDJ-==L3Qgkoxldd(sWZgl1v4`@*^2D3X+DVePt8e8EJX=|ec zr!#3rkr%E+Dz8)po3>v_$qHDqS|jmvs9gr7=de%0@SEA7DQrzx<}G|y^487=TJ-;s zd;(N*R@G~j28 zR1p#D6AZetaemXCd*7xRW;At`C(1D-r?ST-`vaxtm;F8-xB#sZ3Lkrzc>j_6!ZY)+*ErAa zQa^MreRkpjqH_B9PBxL;WRG&(A3ZyH(-o|MQlGw@KuAE?T#-y1d$@!BxC$Vv+=_5n z#ws~3x|D(}mLr@_)+AC+F+q5;x~=G?U+0@ZG)_s=6oL~Xy7Rd-Qy~{7hBn9R9|^DIYsq0P)}l8O8*sknEaQBDdIrCjFc&B)R6GA z4qqXPnyS9yCamFV4%E0WzB=YiQ?xb&3sy7?L9s3t>Bd7MZpvvJsn-$gUGdi5;Y9hb z=`qR<2uJh;${Z9C8WXdF?7d>V zE3ww2&XC?Jo1_WAD#Hib0AQj|z~*QCdP^60pCy%K2#!dI(UG2XBp%l1vH2*$Hg<2e zd(vwMmWo-9)T1ACdtO^$$)_Y3PjIWGi;}ZU7eW#|8S*_a*Az)M^>(%-VpQeQc~593 z_f(7>7k2zp;GwDDSEqIWq>>OG47N^aZl415gp%PzJC@* zy{aB+d#DGcPo=JF3b4bN001g9sugWOadcmN{oq2w+N+iO_t&Qc>_X_D{)Tcg1rXV_ zkHAkLTPq?5Xs4HqLr=};J9z#H zyZH>kvjkVTKYnv;MtIo^eEu`P5h>2k#~pkW0W8O3NHFkQ2~WR-@BhM*v_~84PIQ?L zS;t?+8&kY5{fKzAuk=)A?Ub5!dgF%P8l8cnl2!VlRk?;`;yCs7(D1zS!p8{qguC z9_0)b<28f>it6M|(B%?Itn028F%EPlQX$y6gkCC$o=V64+sAW``#o!^;`f2pP`n$w z$QOntUNhTTP)S6PKU0xlwxX0?LgeHSaG~4z$pM)onJGN-2W-1o1^_v$T|DrXxE8bP zYuKK6ekiF!kA!F`QfnFKiIXHx5*j6`LW+kFZ7D@|CMCg=O!&77$${E3BW%z^ExSPy zFiy%E#-)x@Y(Pw`WJppO@zD%Fg}f>sd^%Bw;R}_QsM|gj!UvNkdqkv_PNk7Aj0pRJ z7z@~k!X&i5)k~48U`>>i6LCm^DKIEOc_VUDm<;$%eIabx!v))m-5>wq$Hl_quL6SH zX|pEvf9H>zoRjzjqistANQx$MEg8O8A_Oy~sf?l_SRg46P^>kjv{Q&Y5lUO93ymlF zt+~yg9zLm!7bgDIkV;eMm0kv_NV)RyxR-J+LSY+={?@Ok!nVB41Gtb;!tc`s5cet=@oOdew5Wd zPVgkb(*(B?(8g{>#h)`J_JR(M{|zpHUlP3O-8xcDZxK#gf)mRB9oL-#0;fHUCd}0U|_~=(}~U?7PW@T zLIDEY75$LhR_Lt7Fd)QJBAOn=RzypmX+>)h#K*Aqu>{087O^$ZrG{ht3a-RJPirEf z!zR$%nqfk}$q%pQJh*&c9DzI+go>Ni*DzD2ZrCKqZa@P-ZifqVY9TTqG_yOfy=DiC`k99|Xy(Bcvm}cEdQ2L)eJw-(SYhxImykd3R1x)2H{b7w4h? z{8=FxI!tJhSSs2V(Xa8{9`BcGq=ObZ1Hwl71`&q`q1s-_NJT~f2a?hGs0)K0cS zOOGXV5l$xakz0nI!qij((hQ?9TcaR`h@Ql}?F4V~1;Y!;Kd=I$os(fZ;trzYd@p;A z$hlwU;|A}xscMa*3cNo}RbNuOyrZjC{YD{q!K`5Vr>wJt7M~s&q5cxPdOQL52)LJ@ zVf}(#*-4J#t9PV}y)l@=;lLDO+jy^6t3wuD&6Y1ExB>tMxUBWqL4w6fp;gGu~B4 zsp%QEY2$~6YS*Hk^{}0QeiOaxE8m)k#}l&?9i5@pj;JY<$HAvkdCI}UBjKI|%`6lf1M~FIBOKRA>x;blr7AV<_r!~g@Im(RY4*o4-v2CB z;{%TnU7va%E>#oMUp@aawNvf&E?TAzDMwIIyfxO9ibBiX>^-nd9o}{d5mdR7Ztm2e zCJ28yD?LeGNlXZx)RBzL6}G{sMSv~pR8P;&)!@BU1`1xlz}i9G^L^?@-;>!y_0&vwOBwB=4Z`FI zJ1}{`{`kI#hE#nALn$)o*UKDF37jDv36b=KzD)7p;(;`V)d5K^)w}|01@pKn(#=3x z{4lpMFTD=22&wS;&^GH9BF{~ooC53ELj>$$Ly#x2Z3ZRU5OOz7q>YVDtKdz;rW(%X z7xtj$bEJ%TZIWI^tQZz`dy9mY@n<}I+87jFcjlqS@kMIWdM~BjxvaTeynF!D+ z<$&iDqb7AUtbaxtyq7=-xI}u$RZV$>oo1Y6)`Z!se2DYwLKOYU=zw%VIz#2&xPyj_SR1s#@*LcDyMY)WVTu zo{1Q&VR7SwX}8mOtsB(b`NrAS3S%NNfYU^*Pm;5l%-?)Qkxqn%J>xyPK^@^-THyWP z2KABZ@jl$BE-4flIPUY#Jxxt2N}y-Z9G)4t`80SUtAzE-X0(N5e6sh$)7AEj@WMI& zIWOesz*@4Nns_1##n9R}9>puMh3lU%CHwocEbbR?l<_yyDW&PKw<0d3+oP;Hn=@?W z7@%;eZ|kwTc#XG#r-*Xx2+4LfDCG+|1m9{ZrTu&ILx7IQ3Xh|up7mqpP*pHwmUUsE~* zlXNHJc;i$>!AU!VTd{;qsijgIP9|&%I`RnJ#kL{bLt`aZv&wf6c1J8L>;9JP5-QqQ zv+%C#yZK%I62fB@I`s@Rywr!@pEs$GoGLZ&UQpc`$dCzJKW-^gA9D7@qnTzclsl}9 zEVTG8r?o67%w}CA91X=02a$b(or6YXcw|J1l`V1<@$odP9Lv<_#0+;aE2VPZ%GiBO zJbSV+Pe#F-33fy3W$WT&m93%(f?{JM*pCFu!N9;eg|}18;+8-cUMNm37=YZ!-@S(jW8r`3d;K}JJeLBu@5wZ>fIkHd!AtjRk}0p`eVDML(e zIHE81e%-Fx)x`sAI@Ax!OTUF7H|y)s(|X3SE+0b_$O`s}1fZW63^UguC>VTi>lNk2 zkG#=kcsW*i(QZ|5ZHy96CAx-K21LCR^ML4hVj~dg;3?L4Bj?~JOy%_JcII*Q^_A5w zIBL<6`2UVnz4Cq`Tvk1{o#FT_G3geq=ACe{I>2+ktVULozxPcT z>REdPv22(TtF6p}q4u#atLinZ6j8hqW|8ZK_@9P#r9dxXS9aqAEsmgJ{ps`Oo(J)| z&%5p-HQDPvPwgpXpX}H!@#dYcCS}B1mBZUO;wH~^W+$*FFUHhwk1I&ezR;gKR?jLS zrekrgOSKQ5-07g)SIOxPMz0;md-{AeQ%(1-x)41qyFlG6SkrfTbr-AZ*;z49B(qsH zHuQPi5dQXgt6%m0ae*5dDR!G1L~wHpFOiNylzuvEgx__3@-EiHMFe2+i{VaAix(SIKUqJ5odL@Q)2t6t}-_H#5HjJ zNy;8dDa$@UC*9|mb)odC#_>MA1URp824-HWP7ECXD7#LL06lqwjM2Ian4(Bx)V%#= zN8g3bT_lkAE7|ia3%m!fR`W6lf1^0P3hPZcrA;)UdkkYhq!r>C3L9W>N$9QFP=_wK zqkf0xOyNqf>@>v9kND(Y1^-Hi6@f%f8(=#kr^x~@rxT|ha$(q$ zB;WI{xmKO0zBe%8tLn7C5rRIEL|-OKe@5^Mf#{k-V^0Ga_4*gcM6{Z)e1$;UG2VmM ztLd|3#VL_{!1{fhj4(f!26L^R;RS9`M-+;b_oUBTd4u}OL`vvb0g*V?3W6Wgy~fv$ zU=Katd*=pqa2ZY)bcWia5RQ6{H{nJ#B{PY&Zy?(8=i5#&U*%>>45pNiXp_x`yz#Wc zVD#Ht?BZU6Hwk`C@G5~JcufE?HZ;KuZGd=3u?t5z$pXWF7AkyhkAJ;Xtxa)bZII@B7F$naVs|(H)0V zv%r_63Uz9w6eoHR)^6Qhjfj;NqA5xyAPGUcH%0K|h@KB)IP&;0WG_VF3_cO*FxBJT zYDhLVB_(`om}$ybLGH7`C;J(n>}O2b&zKrc4&V3weiL}?_XjS%8H!C6<#tolw7R*j zVa2*S{fu|dEh>~5kD)Z{=ZWNUVibuHIG!TBR^(x!MVadT$-@(fIFXjc^1?~M38hzA zolJym_2fu;5lb#)N{kRMW0GOb7y~$rif1_9cfJX$)(mIhjc=;&6pehDABjze6*Zve z?swINnO`tlEDUs-^#$g1rXpaEq7@1^6hAX|`4L()w+@v0Or4>{2QBe}Bz2q}hS+Xh zH&}VcvQk-YGb|2dhN?i*zeoafsW&)i+wP*tbZ$gUb=CMO?w>vv1F{ZvhSmXyB3q z!<3nLhHnjP7R^^CW-SZTkpb(@pI9$(OO_^wP+;nyg(%vIy-7n-XA}&R?+4!IZH0xScON0)dld32G95&krI^Mf?sY{&$oq?|1>QyD`r7rd3Q+B() z*=qCefm~q5pv>R8!*NiJN^h-zINE@hMGh40Gw*?Owb~iD;~tgt zA3*x5mt%Q0T7%@=Y@E8|v28CJP`9efyf+5a8nw==xlc`3&wFd{Q->{n5T*J?ZGzMg zUe{Dtx9YgAwzi0#RRa_}#7>VV7+PCEu1$Np?^DaiG7XzPPR)?~^mNY0_dMr*wZvO~ zzuK75^U&33f>rEAYVhkBD2L4xQA_Xv)Et^UBQ&Qqp(9#W$spmvh07KSA;4+M0yD@C zWQVwq2$HoRlzurLr6cyDlWlDAb44ZmAFKdUo=&ohl4_IA=4LA=opRo^Fsllr;7k$k z8`T5s;3>PQm9vbj%xBq$kt`TU0c%m5$AbZ(3y6SipBqv}{}t=Y+9ZayL2K>a$!~?m zXw;~gV&f%*IAn3lsu+@|Ea?z63@*xU9~Euz4Md;?r#oajDd(G1G;XuuZtd`@9#FnB zWD6_Au`DT7jCh>5J!7D@fDcF=S*4kMx`_BD1QZ_h9*f2@e^5+t!(kN}bUpQw893`fb%-PExVhnWOGa>S<6vDmr~!**Uw5P3Eso=S z3#x=4%H|>TCnowzL8Td)wO{j)6b_Q(be}0@baN~uv0`>?bQqcPZR;T}z}+N8aRUGc zd>y6;k#7h#fQ~R{%WqLK`*tB zo^OO1z_ZgJdr-|-<6Zf<8kbqY`bB2pm~ezHf--qrgyg3WWg#I+*+1wP)|=QTi8~_) zZpU#!Q=U!WpcsNEa%3V3^8~@!!8!MZPI6F;rqJvZc)MJnww@BOx8?~o zaSTFa#5n~lz=2?+C%m2~)a(rtwjuZ;-GuWH zg3)2T{%oK3w7HE3t#aL65Nx+3yDcYR2VR|E)cr zg?6%#kWC$tY5ne@mQ*6vMF~Hih*KoOsS*qkzIyQn2=f>+3hI~Te1h!Ef=e+%Ud?vT z8N%u~yQsxGBwil|OM?biV>1*#N+CHs-2#JVI2u#vi|JW`V>8}G&#ILpMc_G}-E1J3 z;P+m8R!y8H^hpnnCDGw4>>Qsr_BplQImb5;eGX#1U`sc#2J{iPnmM|e4_sw%Um!*v zL2#UZE+GPv8p$+M&eVS#*`coJOQa87Z!q}yc>R|ilIW|x|*%l6BEGFvYr`OCYx92 zP5d@Uu$dr4&_eKq2|!*eUxW!b4?4osFnZmN7avz*UtSM<7frxXPRE=1idyGv_YI_9 zQ9G0`~Yu4%2q(_4^tBxjyfuU#PXtxxRsE zud3&jGq=F|*r;il+ei>1Rf^6a64x&=_wRn=>-_KjA_Ph^70)^yfMJTzWJ{wAy& z8a)(Gwp3yXC+E91rP)&e;vpzSph_?eXEDSu38BO3Y4PS~l#Zo#zEI6x9uud<~N zy~p;dg&7H?mxU;j;Lzlt{)<5Vur0)JpGiYPU&5!NWLEytU*XzdH}(7MtJ5k1dsr@<5*~66K>Dq(Woq9p=ejH1L~Fn{!jx=h zV6n)`MKWV~6Ux4ilsFG>^l#Kc=i7d7^>5TCvp6VyEBbxAV~X1#P$B{dM}z#KSh6G3 z+<;Xf-F?q%eMc>B&*8bkcJpHPaM!-wB+7?((}533=PCR$Rt37<=YBi`4WQ7v8RG>-DZ z0*TwVBDA03wgX<|W*_V#>LK44S{*gh&4j_@>t53bYJNt@omPW9llh}WfUF5QSupgq z6`%HJ&=Ec*J8?HjjN`YM_-DwY0i_++R82&8!$Oic@{nXG7%pqTy13Wx9r&R-V6HGz znGjiqx!%gwkcV4#x!32leW;Fi@O$@%2#k7w-MEKfSiy&o*EJGbG&(!N=lGOnj(%^T z>bGiXAj7=}y6LyLYzqfBPa6BMsqo4GQwjqMB)Etuj>Lyim_Pd!yNL4*2;;bcSjAJ8 zoslQejG9cdD_XAbCV?#BU2<>t0=3Yw&uj)d}9=JWmCFtP-wzkr#?zh z`cyVQOf@@(jf;>XdWA5V9c(-&8VmcC)h%R|)3q-ZyaO{?mdnqsU|!}*!m9|bCb*p7 zHiFv;?jRUq5q4u-JbddV-?be01ExMC_$|Tj2)Hp@MTZC(l*Ch_6&H5FVJ<0LsWy0d z4xg>KDU#*}QuO zfx)x}8=c5okK-9uie?xqphfW^AGiOCe50- zh?@10Xu4)>9y<~%k(@}#*#tzR5Icbo-sL1M*}+8u$dMrB1F+agG^S%eYJVPE+&xm%<|Y(dB9y?pg6f?pH7P4F9n zcL2aANBF!Ki=0i+}|*Djj70i@Bf|ke5E`NyN`mk#oX=cQj&?95COTm+GPfec$pi~~8RAc-RD#WI^B|br z{I&>#lktgKAE9r7v;haw8Shs8z{a8ctbS6rSF<)tlGvU?Nx;!!S?u{H#C z8Z-WZPGP^7i^Yx)#oF|>D1i(_rXvD~fkAut@e-7!6?j{7hT-dCq4TAgl$@&2*j^Nn@ZW^(8qqDHL(C?HhJ zx!{fpN1-ou7Ua?#1q*45&O+T7VXbl;oFN&Sp{z9{-i(pp4%laLNFps79CHgIKwm;Sg*7=3q425je1Wz5y^MUE&zn2W zIo=uP8|WS9O!dpDC3~-V2YVtWIG_H6oM$vgJ%)f5p_+2id`q7|G=X$2u4(Ggqfo*W za5(|YC#VJJ;A00LJ6Np4`0YKjGDEp(a>fR{i3d8f@H?>bKxaiD^d#|PS$h~4(rC;a z!Z3(Y@PusmO}GUIj|$KOU2Pi=_n?h0+#bVGd~GzOMHCt_YW6rpY+_`LjU4ZNKHZ5p z=M{LZGn{Eh$kn^O(pK8T`8xbf@|gY(a9*F{dfhb_f)3W6Ep$XTX!%&a#o_V7K)_s z6`@)=4phztgr)Se6=BjDiNvOn6!J`03tb9nnHErAPg0l)h}%Z4bF!9jPDaTAtZm$e z5ewZN!5vJt3txQEX?sa3me6R4;Aw@?Zxw8bZi;f3hf<1ikR&a{VYA)VsOXnJ|X8@5X~XOCFO$L_uBT^Pa7BCN2>C z@cf|Nzba6o8M{$-VZmma)B*J5{VEcMh~!w59{m*F`TCM?f7rxr5aR zm)wPvDT9YH792b(Sd@x%Msf17NB@MCKSS^w!4m-3qY>CF*HBw-LVXZN%5VUck->KV zJosw%P*loCS%09Lp8k}fCv>Qt9^)Oh$T?`KSfhnme_2Q-7C}&J*o^bpa@32_O;jaSVNqNXU6!-(#^UZ=r z%0zZ&-;P1RjEL~dCP@3ZlvO&mHqmM`xo#ObgcZ4m)nCXl3pb;^OJ7VN<}W!DcQ^6i z7tsrbHt3G{LqD}g)kD0Z=^`lv(cJ#3pDP~;-cU!bU2d>8y&snh^~@U za%qgm^>F?8CXTAO1+)Y9%CHZcg!0JtkUgA@))frAafDN^oVo&U<`U=0jKtPR@Xg~y zHI9*kq3H(Hi=(}ztPb)wqMp%XIEb@p5TOW7bm-e#3>>2HrvyWW0UYHRPN^*3!P$rs zaE9jVSNPGeeY6{|Uqzr%stiSnifUCA_Wcs&E70A6f$NWSuJdKc{t>B9rX?4F)N7$~ z@^@X2&yByAE=u9(dNrFnfk2!WOGzK!2YPisfryg$&IpNIpR^@k* zg3qp~Wo@s$OH$v(cYQ3y`vLUb0DT43v#ag)a~Ers3jw&t^cBx&lmF(;v&-=$GRA}a zp0nIpn%T-q$aSr<6TtOS9WYJN`~DSx%;rvuF#i663aIV4Q2q!0ZMr3~i9z^;&k9BiYRK z;SQc0Lm|cKOQj#Zwyqp82?_m(BNY9r2qOIP!e^K|O^i!q7;-=v{!jfXiyC$uH1;+|>T(QiXWipk6bMW&GN!6TW#!bhSY)Psw| zaYb=4@7!sYKOx|6X#p1~{PY%b#mssB*Hh7(4 zC24{#fM&gj9XJ92d$H^limaTzCkoPLl`awS3bz6Pg0CCeR$LUHFU zj=T2C$$CV*@llIK_dIFC___tZ-PE$VCUFWFk_#^QlP0gmoNJfhDgtnW)AHDX!N5quez*!TK z<|i!pO+G%v)DY{^vUGDiEh8LDc$aW7;awD)L>{ui-r!cafQdlPw(j8DXZiR7`|u*c zPYEbS;e>CioRGqOQthp6xjoVv&zI$S+-0PhoTFWQ?BZ?Dxq^C#<6pCKi50q#O$&GW z6<;ipw*+q!$UQ1@1sB}YHex6MJ*ytRtw7diER-O2b)okL3D&!VX!@@Vcb1Y;Iti6GIPjPCMGtC)z=Tv8pFT)6BxJZ9Rvimi` z+XTNMc!%I!g7*pDBlur}4*=@%k?iO!@R`##8~+Z=y|#+hUY2g-Y(@#@pz9EwXcf^5 z2^Iy8AwKPeb1hi<@`p-HId=+#TtEQAI)iP7Nmw|d!-c81xBg6Lriu^roawACv}e81 z5dN|S$wSt#7&)b*WQh9vY_<$Vnng8kuyIM^2ct2f-pV$NY3-}peopl<)$`|7&jt6Y zCI^$@5L16bRU>|N{nz@+q)hF+d3^8aUwg+zoI^8$LK*fM9QK(6n*d~G;Tony7tT($&9g5S(3KvAaRax5TGF45G1~(VLk7(p zg(cAE^TvU~NN?VjHs>Uj7R5ds>xg)dp^D7&Obd~2M(Ms$Gw08nQ`?3)n>laJ{I=?u^HXx^ zZ%rxxG*}NdlT4_nnAV1BG?ITfx#~{ zZ#aD>y-N~KyOUJjuA9|?164`qgg}*`LB=a9>Md_ew{uKJmY7&4gm7Be3H;-Fd6VN5 z!pRYd>|VmBbEHwM**_DI<$^01>+uPq31D$*kb#<9c~Oq5;T z(c7Gq>KeFbo6}J|c4+b2l}}Adf9%cK=~R!C@ksY+n`*X`*y@FMg44Vd@GipdF=33S zZ!~9VVO#Yw&Jt!zrb*Tx1)uChf3#eSVS;j!o$_wK_YT_0v~UghElP6O$io4+1+>XF zi;8l~S)lzhQ=f(w5Y>vsz&|%6qN}VA{H==){28Dumx+MWd?BKv8=qb|eBO2QL3%{5*~C#!MLDExy2+o8bdQqP;*nF8XW`dtE)Kh8YdRt{`Z}wXp@w zvf_9tD0x2uWD2d{>mhWcpx&;CWko%~AF_eP9JZ{g!;lByi-`hbsDvH7=|dW+2+ z%ty*^#=GZ2XL81r;yP#A<(hFK&z=%{kA+19yN_E71_w_!D;2F$WUTx7;s*p28R0hX z>cpQI79x;viScX-|3>J)vtL)TP$1chYm&BTe1$6*2w$6$<90Eqc>T?Os=L?CAl5i@QkM|Pv@~s5k zOlJjsd=wYiB&Nm@jOEOU7m2Ir|IB?i5=OX1laYUx9>i+?1Dm|6} delta 31579 zcmbV#31C~rwYKL7AxcO{oP~t!c488e7?YSKmqo~yY>UW}GLlyou@)#h z7kb&t(iXx_P3zLqpVHSyN(&U)0ctzrM8p`qGMIMMY@={=4DXh0(Jfskp01*?o@`Eww^+Q99h7ibod@1gZxDO9lc< z2Lj6m0?X}UQ!2E0FjPGlS~3_~YC9WCXP~K;(v;dzgZ6`#-hdyyk$@3;JwYviH=?;M z+PpgyOQbX5M0505#HIZJ<;q{IHzDAJH;!|Tm{r)fxcn-gTH!o7aysDtk<&-oh3Qx( znhyH)8sxK`k(J}d+r@32X&nh?!lof4s70qj-Msi8R}$=N?) zg*wBTJn@yel6MPZcM}Y$i=Rxmj)8dJA(K8=N7fFheTx1C@}XZm&=GS?c> z&Dc(W`e32n=3Fv;PPO1-u!voxlU>_l>5N^>dN;wm(PQ$FM5dCZGCkKvCQi|5-N?ePjb!|{#-km6i3 z<7BCNukwmoqnU6fqjMP4m!kO4K09M(9$6eRRku4kW-ciuJ?eS@+qr7yxTD+nNnFZo z={*byqBuBx0l`HCml0e;kRPOkxB_&Y^WMz)e$FpRz3+%wPx)4@V1m#F#}o`0Jl(`t z>E#y@^3-lmbq=XgYO8k^9kzAB7nqY+(H9e3LXdA;GO|-jA^vjbp2LNtQfoKSGrYOeZ5y$Db@0@Ul10CVeFpTg)n2~H&#HjLApPlNMR$eBLx5OsRrnt2ae z!7#Gb>nJm~9;5_kZ32f8>l?e^X0_^ee0!x}KRT1{ie=jL*@(z33A+enM5@*|BXW4# z+nje7Zcu0R9e%_%HCHA_rtE)P&ed!*?%cj;Y@i*rHCyf@x2Q1znT&lX*-)+X6> zp-!=}8wmCgXacD}&CoCmFw1II-O+bz^+>f=2(k;2o-Y;8c>Z)-vTNuHx65SV>u9D! zCx%zn?Yz6hUe?2A&Iib?Ce3ATS7_lMp@Ww^$LU$RZq7Z(7*zVafeQ3qhKBX{0_T@Y zSE&YP^0M!ZI~cNVaz0#khg#`evwYf|49njOFtAEV`ok-{$$4@4;i%9$^3SsmQlVGM z4FQuF4Y!|3Ox@>9UvU{4xNpUb=?}2JhY21!NX?Hp@2x=1eG^w6rKTOEYR#2;g{d=j z)Xno`G5DL3$+(t^zxlsayuw+wYB{-)n{9zcvA9EgBv%IWk<$ z-ET&v*L@`@UCB569d_Yh*m~Sqw|;zt%Xmc>qz-Oy+Rx7mu(VPK|zy;BkUHP9!AD?uKhO z9H(yR+p*yUU-j3J$t!G*C(}{Szo#RX@roi*5KXd2itbf*acf5`q?5@U1FLnsjn9?X zfp(BqRNKxYo0h64oZoJ;mF;Y(n=qAQw@YKq$wX5u(;iOQB@Cz9l9{CKWa?J0w2Rxr zu|zzXOo5I5kbQiDeXugS^;0Y>`9nq+$?anQ&0Mw4J9Rq(uQIr=Z^O~wn^4RVVXfZi zl!RuC6IOK#Q}aWvb=HNJs&!5#bRl~8m(a3PUmM7<*XNbCL^I88ArYHiXZ#xk`v-EE zdIVv=LR__l#9t!-x+7Nt+LIl+IjVn)kKSFJMCZ29bm!jF8v>@4Yr|t#{QG*38m!mS zO{hxhvz=|>dByMGvk?>OzR;`2xDCCkBv-R@dw6{C-&eiER4v7wMfVgtpM~dbt3jpM zWGHRQrmp{eJv&j)C-_%zD1{ec$sDxa=cJl;)DEH)a!WZVvpJehi-7eX)ZA`*=oB^| zTQjttG9w?SqT$`?|Dc-7Of`8tH4mvwgs%Ufk}FLmz>ejWQ$PgEPj|#KX_*n-%Mrg! z@Hs*EkXgwNZH;@-f6&?;1Fd;Qn@9UX~GtUbC(>!jAJS*(chBih5eqoG(P)Ep0|(>?kzT;pD0CYFguy&~+7@i5kb;O4I^Ac9~W zLxPMJgo4xcZy51Dfed7w^W*3wb(8a(=n|~WNi9nO>suzMSDcoX8LNn6uQ=M>9HM}4 z7crDOD ztJ*oN)l-M`S#38L1&BSp!nrGXl4^23PtGlRixKY>IP+600@Ju3Oz(@OK32_w`$_um zOcoJE(3Z!dU_gYS2bqN!@J5G*kaDOg4349J%G94RbwqB?P7LcrEYYgBI^Wf6)veCw z`jo89{-JDsApw!*l_YhnH6$aF9Z&YY(M(>;WPdo(v(;1Gp6US$PDG;JArVREGDi@= zDfDZ7xLDfym;!tExa9pY3b(qfl3cj(S&BesXM?ffA%m!|UJVrZzDdQ8kvoO9HFEpU#_d^;dy66zb-G2K7rF2whigi^8Q-SKE-VIv5| z+0e7Cem(nJ9B$r24jt7$Ve#iV$`Z<}Ax^s&h^0f2w7udplCcE03|+#Z?qarpsZpQI zf-gI-^vnc>{k`X~38x|ga@B!r3fs%--UqkHWOCw+h9irE(~7N1=404ujIy=u%r3MQ z;b$#BTbw0l&QLAR@n9JcQn6Z*z z3Bgi=WdzFsf+hbt9W7jS>15U3ckiYDrPji<&<0*%m*{Xh)r1w*<5@aw`*nIrq_>Qb z3#*3Ykf8oOnt>Mj+XCml%OR z9Mc;qt3vi+u-tIuKj5BYv3+4PcK-YFS*zBu5z7O;L&-9`xiFZ;h}TCl(Wk#a#Fd(M z#9|;hs-0O^>;elq`--_`wr|sBo1ory9=>AsLQa5A6J%JntUcDOAtOPQiGVX7n5w_z z_^u2FT;{&XslKv4{3x3BN~7JhJdi*Q-#Z5xpy%Aed|85X38X8WV<_i=_GqRp8POFi zABcmYg|+8zkA@R^4f56n%QdS;EBmjUqO#84uAHGdo#|Iim~EGA;rO|cc%|*(44Gd% z))Xu=&E3tG?r~1ODx>x~&tDar_GLtPC2XN97Rj^$N7Ahve|Ap1dT!QXnd=CyCm23v z#!*`K$_?g?reXR=4(g4pqkKidXeTTV9cfRM>RBu$d~}y#gIru*IcyRk3_UEem@MHJ z3nn9H1?w+Nr{b}UD5z%3+UiWaX0FOO8?KqD94B&3RhF88=Sw8JxK1*Lw-#oSk#LU| zR-9t`xR8)nXvfYF&0QgMxC8c@e@G9KP3crmT~Y&e=^?zv-x^Og$(BhT7;P7(ZuAqs z(CRl>+VNexdb(^2!u>DeAVsqMIZ}-_aI$mSwc`V_VeD~wuU)=^)Z&#L-&hL-MH{hi z=TYqy;;R+wkM#`Fp#Bm+Z*tzgb}Dv=LT9Jyc6K=@XGI>9c{`gKNMu8(i;K=4C0pSE z(L)gA%dr3~5yqvXM|pXi3L4!5W666`HnR%!fdAIzx%bc4k~ZA#0a|*7I={Cex-w zM6#5&#}b2KWXy&2Hv}i_T)t=-KL;|Gu!mdf(P7We z|3<43s+j*Qzh&gHHGOZn&x{;pm#25fQlVrb)EZ8u8|^X93tt=KynOc+qU(V19Ph;M znLd6$=S_AT0kY%#s=&G9o)cF78U6K2M0JLxBkGl=h~FqSP$MjYQ~Fyv63zOr$T+kP zXVKo%U}C;%@7BdLm?UI1lcAhA1ooppMZ++#gnP|0AB9iw{yUwrdk+mf#!fx%EW5XE z<6Q`%p<;-r@o-Z-s=p!;`c3@UU}x>=p)DiQD_u=Nr8(K2T6=OK>H%+m;@-K~u;0FS z>x4oJF#&=i0)E!(KNj>IbKkM5#D-xX|2Emqo%f%9)Hl&#py^5Wo6JMMMesHNG}%fa z%%MxJ25pFOhd@C86t}%$+dJnTrwZ(G!!jcQ`SpB`#*coWe-PS?;60Dc zIA3{i-o`&MSrEl;1{VwwRk2|f7^NUD;tWbRMvV|oxy>2<(9Cu_kkOFo)A~i$?F&aD zdMhiqo!}~hTLHX6v3!cm#S*)r zF%tL$&X91ne=u)7fbGnGcxv|V{5+V%JPHlOI%AQJa2zHFaw1*8N@N2siW)X#2O!O2 zUw||W6VWD^dnU-Bl7KAj2 z)gL}QE^C+Y+c3ji5C;pMYPZV`l0!^S+vQj0*$dd!@BIPsp8*HF||HY@$3CI*YjTQ=_WaH*(z6-0L8nA_rIF1~w+P_iBcg z5G*AqXXa%LEhjjVUP{4|=LI`G>z7{UmZ?swn*FGYPu;CO;jvb}`)RU)Y@P zlq03j{0mdiisLLUJ)ushxx0ueq$xu|Uk2zTxQSp5kW3vCL19HJi1d}t&wsoayGh`g zOVx$W)z3`ImN4&A1hP%!XZq)e^8#{&miDY<+6!Qso4X>Ozby%|%x00333d?(U1ZUi zS1c3Ho>SP==Q7S{svaq3_0vamPrdwHx12$A5}`hOOxwNB+P_GGzLnN6$%q27ur@{g)6&{eW4 zt-pgfP{lEvG2%dPMaV91(2(hKWGd7Wtn4k>iIrW&f!f6zU?QLp?UigxCR$^ejtC;R zg<+nGMG(0u9__?OGGSNL(ehEJqu8es%{_V;`Ti@)Ezs-OApvDhoi`t9m-d#;K9YUy z272lmfo5~KEDGk`L~19&-wK>(yhAq#F>Ygo?6XA(d8NiS)E?HmWkd7=u(|MKGa7?c zJ_(DSSE$eE0E>LtS^Ub_Bdge6Il*@Jhg9qp#+d>6!imuM+ILFyWto!6UR0vi|rIF(5_@zWXR^D2HG0Tk%%xyKEvqTp&P~xvzQpMJRi)l z&;sN%4*RVBBN4!*X9rS|+9Zsjze5amKr($D3(LlGJTuKF85KGYygFf&9l#ER(w=HY zg~=}4dGXab^XrLFV{)VcjOHwFBtwm~b@4MA$%%=5b6!){t~3i|2s#LaSUVXSTu{*S zmm{LCzOF7?k-PHjhopp4E7=@+3ktH>q*$^UC2j0yJTV=#g>OfSU2-4|iJ)~b z@k-BT?$Hn*nxGj)Lf}F>Wd5-Kr#m-^c|`WUSJ=sM=mz$=_%0C<^ttSV2oPYch7;Pb z2&ZH}m)QswYj#JUMQ}F3;6&+Ldb?l=%@g>(q+6hcj&)z;~$Hr;@1pQ-uKMTBuN$rQbq6QA5zs>FoLal6#BV|4^Vf+V7D^Hqd{mkhZj+^lpugLqB8WS6AGrhsp| zC%t;GvEQa~j zuz0O=(w`2Wd=D{n2~)mEa5=%11XuUn_NSQ>17bLcyHAZ&XDy)R5+<#%Zf=7rnbP9& z0^KpCV{xd3u+dS3^L%G?B=mW1vQjM@ZTr*FROX#hH26*#TkUqgS*g~I>0lFQ5tyZS zxbAd|N2#fm!)%p3?$S|eR^x-{Wl-5A;docLCmrhU(IP-z$_g$cxPss+0xAzfgh~-5 z#ddHnM|vEA9H;Jd?;E9NHi>-sHhKnQ!y+c;4{AibCa}wlNmkU#;QS$L_3IoG7w8bh zJJ01srmS_37_AP;4znV1_wb4vVAv71`68>x5Io85i%d%+o;kdbbzQ+vS?gTkqL}L+ zF(ep!har)IUtvguk}}Gi;f!dN0bPjNZTE%I3TkO_i%zzOnswSMfq5X27U49nmCrW2 zj-dD@R`)UiQ4C9vG0~-mqxgB?4pL8ZGN~h8phtc+&ztp634X@pk|-p&P)9ZaBk<1| zD|0mr>%!%!@aOj-IOk!r*$C^q7&c3wDrefT_d{{U(vRb?X**-|;jHsI7X3a$f*4W8 zMoD))Q0RvWhE9>^FAw-#j=TO0Cr7(9@QDd(v4K6Vo!A*!a=a((q44LF7AoP z+F{a6DV2U$VWpMS3A(Kxz1_gI9SbQ@gZ0MWa)ic!WD(<`3 z4>KFFk+nnmyPo5&Ve0~F-ScLsud7l0OJ}Of#;X1QnYxVw6# zx^j%z>&jqB3K=akr02V@tyC+sWRx(ca^iI<0S{>)NRldGbxlVk!tAI*?~k$8amX1Q zY4l@C7kUnWUECdlEhd(Xctvapd$AY>A=c*G&4fE!;ke2eAkK5IKMEwsHjV4hVC`De zw()Eyphx?txb03R63GS0w)SvyTg)z(&w;Np`7%kKwL2+93+)J-wD>ziRgsYFutC2# zjc^l!xDk$z#W^y(o6M)@XILNa9s7?x9QUJDYNDF!maJCisj=>5tJPr@a40V%mZ}L^(Vo1(PS6RKu1I=gLo&ipx*ZN*Z6TwX z>pg5t_77^YUQr9Ahs=^{{Wr#w!Rr|e%>?j@nxdJmXf)vkBH-KT{WK{+oN)p-Vm&Zu={egAB69FBK1+p z)p=R*ady1a9IuMFM?$ho2385=EGz@aXz9icSQMC_jnPgVeq}ahqn(UykPj2y4vjM% z*%9u_nK8-vj;D&akV7SlwjnC=G+JB=jRCZwVdq91r$Gje*e?_khgDc6A?n__UQNu}g&8s6mdFoxKD$&FMF!Tb6nDkQ<`c{U5-T&NrV)1| zGQwF90;l10Aw)0H5E_Io2lgkSA{w7-5mgq4`ols;e*|H)YPY-B2fTc_3rbi%s4-cr%FbA7JMj>O>A7%GNP=@OSF}?gsS*W&NSR z-BhcBYHfeAR=qq1T>j`2)sj<%%QrD6xqL8yQ~PiN2MCm96CzD`0kd$tq5;#^vB(Zl zu>m81INFg(Vl#zD!IOygR5(}&n$LvenXv7?e4^T7y|3JvC#j{QsMsXqFpec`KPioU~uSR-S5~%;Z%>h`%Y5Jt!FLwA1A4g)Oh!!UFx!+ zNWd-V5nR{Xu*TYOWLE~iFK!$(;_FtVL2G!l|0^~QD=TFekbS9wb1>ii{>kdhMKv7R z!v6;@*nWzdoV}lY6EV1tA=%dDD7{b)p1i_Fo9Oe52VnYMG#i6HWB58dr>lD6EVT;>^jZ?mGf~*iOJrVvEWs;fg zh{Mq;0}&%_QunY&U^~;tnPa>VZm4N8I(ww-193TAW2{}6fF-q2la+Wr%#qn{2KnJ| zL?7zD)v8+6q5a3Usc%;lx1;xs`YfM2s#6`Z?jw%91^`S+%#C#%_v?&(gJ3@a(+#L) z=vM$^aK@uQE^ymBRb5t;;KfW5eM8IvqWFE0m=f(z3l z_Ik^!mmIZhCH}u~<;s=%dB%(O{VGHE5xholKLK<^SZf|&P&7rV`l2rCUo!L%fhdcp zIHxf_%s4qqn#RVb6YOJZ3q!307+nFZii?<26sg4wz0XYNFg;6fE&w3#= z>Q;AyL-&{UD8E|x8tbSb_$9&q!B2@BoDP@M@#`V3$p2}<8~#PXcmGw3x8@7(Sod9w zOaII4w;_8PxeDgx6eAi&^tXwDcM9C5Gu4*t3dHmlo_TiY%(G?0!D^$+J=QA)8;y1w_iNdy{=sxn@`c|VWV^!_Mo-3{xeuJBCU2qY z5Dbd5qP`dn!aztngUCqb#(7a9avtetlS_=r37TSL%DFJuZspl(Y}Ou?zu+*Cg8iLW zV2y78q6PO#=TWKB8=?tlj%?@qRZyY5D-T;6mAhciDm?6AnKCA~0q#y7^` z6+VJU9(w==%?PUv+aJmF(L=B}kXpl;6(CFTqW&w+QKu?(k^5d&ZCWO@lg~1k&#R!I z8>0IaH}CHZtc`M|FZ8+5bJexhZwuVNoU0BuKI!Y-`RA!As>wb2JaveGRYB0m>Y`^WFMO)D+h}UyV}h-M#0l zY3e5Thv%!RttkdT5ad zJTit}c4MrpnW|W#3hE(`y`lZjHOEmmvc~o9t_#)Ftz@=_Lll}T22#U6b!RHGeU5T= z5hUcq_H?V6sXQjmv)u1psH!deetMyLyhyrpV}X0jC2G3+(nV@ymit-n_`z=FL^hGK z=8B_i@5IBa=FJ-#pIaM2)59Wj!$LU3znH91>}%6MVa<TJJof=5^nh#BMbfMHl>67qFfs3`{@~^6*udsX*>y_)=*e+FUhbTLAqpU# zT6D1J*jq)&01i;_&PV|L0%-DDwFtgnXkBOM&mjEK1a_i(;28A8f^$PrJxnY~y*uaS z0DlL4I$J8L1*6VDv<&0RF5HMClr|&qh*WkL+YzZu$h(4Px;SM@M;CUt#~VE#V6U5W zRPh*lB*G!seBs+94eKwvUv$(Q^`QHRquP!X(b&?L`%t0S3uD9qSwQd*vn%Ge)E_Yk zW5qQ;x8XW9anpqe+W}rLiZ$W>K4e-g+wS0Y4L_qaedf7C!--Ii^S{mtmJz3htlbB$ zQ`On;GggkmljsVLuuoxxu%gAxJ9r9V7s2I)4t64?YY`|#mV_8-rYVPa?0g3A44`*9 z3dqZkI9b=xmRPsG04>=?a5Zm?9fbpdW_B{?XvRY{@gnnlnG~+eV0(r4!j8eKlQq&^ zf@=4M>(vFSx_|1IlpRO=@D=A1~^FztrVw{ka=)yY!~AcNA!4{416$j^2C?i6zSXL-A z<0{fEjFpLy9aogb-!ra~-~)mp&cPi8V*zW}u7FyGexuyIcfyGE8`Z!6PW9E|QKCGG z&4v}#qv!5@>I>PMSxp&B+=ftZU}1YY3YI8(-(8G8iy6pxJ{r z(BRpcxL6z4MH$g+OnsGLkbCL>;U|gO98QXY!c$wXydwiq5~6Jfh+qE~qO$sYqA1A?u~AMOYW~s{-5V(KM)IQ!_t}o*#mLqyC85h*>*#Q>u>qbJ(>1 zQn~Lvpf0z5srq{!RKHMdq*r|p$H;EiH5$cz3Wwv9ae9R~OU6lUQG_Z%0+$5SaB_%8 zF^z(T+u7_KXLoALzMJ@2_ z75p7w70BUFoXv%VNPNl`5P76dQ6>~jU`#h zKJ;#&=r@oF4)O3Zf@%|-m%IG4usNxtTIfDci8iY*f;^_VyPr>@66Vj#A%?nzbjd8F z7v^0O(Y%F~JgK1hdl1@miv5#Ni}cjU1gvI=u*seZA9|kc(lHh-QzO8hoHqN_Pnbv6 z8g0+|XUrH#ceG%i_e!ARcK{2p(x1j`U&BDN#%eIQ!x_7^P=-(q%?`n7#tlMNksNIc zoy*oEJCX1{qK583xt(D|3#|N6+=a@t$xczhlyod%qH$i@7K=x8Fa6m-*GYMt1=ayN z#xlRB&aeX}uVJH{AvxO+PmN~TevCJdFj$hL7S8(&LlNYb`&GQE70!4iY(0g(;En=K zlaN8p!%`JP`0Ob=Z_D~zHlSs z1ardhWSbCLJM`RjAy{H7@9_PQH46oCj9wv5MS;ycHz9TE?=TZtilNJIv-c+> z&Nz~h<+iq-0-jcps}(oI;3bh;ye!5=V-QSNnEQ zI`zYK+V;o$;>Qh8gI}fCpoRe3*DPR zWw7SSnFk(30sX^-&%-e>o=!{<(VR_(AleEsduLRQZbJNgCIK#bI5Tly(KxS(+60JnEs0Q0v) z8En))V}gk1g0ji~0m?i-oD0)N5Fu%bYer`>j-~Z68-7s?d_QN!FA@BIVs?tbtQ?QE zupXjUzf8cxB>g`GzaaQAz~H$Pyq&VxwRgjhJQ~L_csRY=9H9iTa?FuRxO-dU2JNxv zD+I3+yhiX#f`6I>#J$dhHwZWvdOtx8$w5FZLutf)3~fBU!>xKronWPX{oOCAbCoDR z!4dgYtcC7`423F%=!7(dt>uaCV9Roh!O9uzfE)`7wQ$YYBY33Sfzvagv$TGOc|KLA zlM&$52wwFj!LA%Q79Px+qc-G!j)i_gkYB81dbk?fuAgW8uwi%l+?W4HZL>Ok{j-0e zo>SJvmiq@!&B(rCA+zuh(E&tox&!3$)=f9^zsJ_S_3ShJ@BMYSW-j?;gS(qc;2r`w zOudUC5s|*i&{qgLfo>2To(JM#v)+YX!Oy9QlhDN;itv~fn(5~Esmzdj5JfT!u2vly zUUnF{H4QJCNBq6$3ThP@F5l@#5Nns87)xY%DMg-U!5(%UaAEGLo4pt|?60xA;!JlN z5hZSFzvAe4#HoMH(A@}u56i&vE!@fML8ju_^2Ty0()o|c2n#ay58I; z-c%P^CI0@pU#T4`D?|JDraoan$V=esN2L~Fn{T#QU-AQt9? z9htEjVpKpa`#y?qA+);_exsIJrG9t&Z`AMSa!~qp^cx7Y5f6V2B0p8p0bD%GBATn>sQXu3IQW=}?pw zFT75`x7EytSvTU!Ft{LsGrfrrx@pFbyC=S<7H5URsooDUY|GM+Wza=b43;*6E->rIIO*ul?rA&ydynJtk z`47P@JG#f`p7Flg?B9a|&T*?Ixj%Vdja@vA(=vfzSoZ0N>llR{7{3L6%da$J^u7IK zKTxXzCv%U1;@J+|47|^q*6-sBo5$OXEQ39SmRr#Zi4{+E$RnI6#g5uOjL0~A<15F< z$c$|F-cP{&?s1>`L^Wn*c?WvnBcANi!mvah`-q6-0?}pY)qp$o&5RJ0P$o%?LLwUE zXGLZ}=0bK!Q8q>G%nPw1JqRm&L)6@|?E_*Do)+~}2GTFGc@b@eM~ZB5JsTH!qmm(E zAmaC*mvn`-%E}e`$?M90Vtqr{DG~A++S4Th04Y;O7-<+-u-C9`fHJph&ck&ZmB6E@l=|TdFz()`uC?Y5( zkZt2gUvnJxE5pYZ4FneJG0Y)PAB|(Eh7%^BmZ4v({%8KE0?I<*=YLX{7L}1H8I9^@ zx9uNl#S+Q93b6yn1vF6O^1%q;2-~z6^Fp}jmao+3{_G#>>FgNR#C6iBzmFopaYOiy zSa~KfLtcp&v#cogLU=-|vO5dq=1GtQ4_+FL#f1=|4Px=a>Hy12HG^hbe04Nt6x+Eh zo?cGU89sM`vgT&xcATiGqtiGvhHI@EiSZ)HIkL&H%uap+3frS{Ss{4g=2IMVYV+_KbwSSF1l7*8;PU?RaJ05G)$ zKKI39>$L3)@Cog`qdk;Jhf>Kd-a+G6v@7>j#;;!UfpQ|ygVAD113% z&?E+7F?S!r94ThEa!5sy7ITqT7Ih&X1i?AVKyQuDoz>85TMclEZK$?ky=xDpdUL;v zbHCLxVYZd(h4ir)|EjV6nURPZwvIf12yK3Ow0mc%^#yCWzkgnt^`bT z`zB%I*bO;sG2IPG?0}Q_f%#!UiQ%aW93&V|b2Kn)Ds-*Sy?&%+TeZIak4IYH^an`h zdaV1vSZmYx{Q40}Jrfwwi~Vl#IBQ!rkKIAS)62OOg+MtK+*A34ReBE&IP%0QR7|d< zspxv8hFIx1A%3L6QF0kaSFs;lvAT(QNx#NDzn%#n|7FJZloxmO|E*P+> z;`TyD)G@~TqC_9cz&Sk~StM|=z-nW(tzmrvY8hJZbC-;_He2g`{b!H2rWeY^A@h*< zTJ}loKz?1#0k0-lL+~-CJy>CbyqkIjq6vhnJuB3$S0RV(wpHVe zr8c|JJnTBGKh|2md&pF4E`Ixuo@#9fjQ<>Mi5&{Jr}z%cN_OB5_qLhV#3kzy_I#m8 zH){N+$6eUx8^%BY3}B$67mSB5B*wvqT_=2LmKC*rV7aZctQn0mkHTI~!-#`Ey97fq z{KGEF@evFxcWXb7x*F8~Bi7HQ9I%RUJE$BGt@EM1P z0v6v#qL`!!WP~op_PeqWXgI#l@E1(jT8~dCF)M^EfLB2!0 z+D$J>yyA{5o;1y{*KLH&)UJrvCR?%nGJ4=?-oVZA-3N0cjLaF_~?3mUi~wel`)jDSQVS1I_LRYbu6O$ z5adnCuC^FnV6DI~ahFXllWiy3E+U4kwqCQ|>izC}i>&F{8(GtO))XfoQ|~Pkw}L7N zer=gLS+$L_4SiO!`G^y$U`&rC^%PE^v4`Mw4&LC@A0x^x5yK+`h1@IVzo7P ziD2$9j#DXk`pb`U$HxNi8XbTn6d_u4{Mbbo{6q5sw(#ZeH=TP zR}{NZmtBryzj!EjeGht1v^^G2wsz}(umIj}g4#mxSA2T%IGa2GZTJjTL_F4uOQH1G zH5=-jQ;B%dGVf>ofli1vknVZ7IoM9O?H;krI&_uTnuVcHXZ->52!=h19-)*QaI^(7Y`H36-_@M=!_ypHFC_t=Z;epT?-GM89;?AaSX?W7HV*e>VUoiX zW#9`Xq1384Z_wcqc$lLq;o1 zxLOM=;N~A$o!*6daox6z)jPqOfj2*gWe3-LH;E5Ox?;grgp zt(=WE0?yE4y^V!x*Z~2x3?-5Li?O+lySTWTj;ZuX#Unoc{@YhtH~X{P7`&1c862mp zvCdFAIltb4;&x#V&l-i->*Lwn2?TO9aV?vi%+E;#B3q)I5%fLPYj$rDG`S1t{md_` z@=b*B#zb+`wjT4+qy81s8rbUh82T+hudixBm9hTGo^mKtx%>2%%xaPU7A-Q%v4OPm zO`jWBXRX4YS-}T6v@tjbBH?r!Y;&|ofo(Lgr_BVysMfp5b>REtwfcROVW(iEOvwYB zCn6gjW65&?QAfxHyvn?lWud6w`f%1M5|(fXS(QIwFSs<#jVWOwf8^&Lf=>uO1+cRf zI042Z$@GE3>reVGO!+gxp9nr95Xum0q40vW9FL`pjss)4h>p9_12`LD?e~D|Lo=7% zFV5};V(ok-#}2d!H_tst8alXaGh04fE(#e8je7N^%ln%!7<*Z_y?F&r*Wr89LC z=QwX<&9e~Fn5L9x2YO6ZdGZmW$^L;OB>Mn=NJlmcYhlP16BBua+#bd`rY`4HRRG{l zSR;oCJzRc6Sp679QD&RG{veMJK+HibVYr*`;hUEuOUS+iMRpO>V{4&MA%wM#fG^yQ zC%j|r7&_j%Y22+$#erW6p|Hrxoel6Ak%%0wx6gQ5Kd-z5+eomib^K_ARmCx+?mv_Bg?LfFMk z?96Nq`~hGl_dYaVWuqH}5iUgC7DtvB*itRQB`hZPXeuIJ1)h+~-7AK3k^X?GayyBp{2)8s_s6zlG6=k(%tzF2OmTdE!B5*(DK6L_j&0pTB33L1v_7E#|qDa3JAD z!j*&+$SsQ)ZbZ`FhE7Z6j!e6II`;rC<@`uY;yx|AK_u9*5;P{K4)y= ztB0J44?79}3agOg@H~_GfUSIu3E91bo$Q(Ln*|)&_n5qpp(6-r`801?;Q-lq>I9%g z)x}z2SCJQ^3?+!UT&P&~oby=re1eMzWbugtBr|XU;^3;7h02bXd+~<}@uw#6r_GE! z^E-43*OK9s29`s+4WEb|2#=JExcdd4%Wx0&s4qzMtZ!8A?HJr&AdBXv;Nw{Q1ce0R zD@g%emk`{%og*)0u#8}kP3lWn=}5*^5=cJIhFy|u#mMEEoEjDuP|MJ-RR8@aSgVzV zz`G|}5BalH<#3jH0;l&xf|Ce#5!eJL6P!wL3c+avA%Hp*lD(K2KC@aD;NSdwvQejw zXYSRUzI6nf&|wHYG_B|EME)fm&5&WO#XAwJa z2J!_*n@#m{c2T5jj>0a@$@X!)Kv|m2eZcy!uRuQRpb+NjKyDC1p}WFzn`kycSbE!F z%dA_37vhs0xK!A($RsSWadaJn-2qQ$8oH%wZ&u1dER37b$qU9qtZYC>-9(FZoZ8g? zy%y{GqES@hjm7s$_W6r`w_vX?O|w6+|%dq>Y+yl7!{OV!${*^3q~ZmF8RI4y7I2FpAMaNXduu!;=Bx*%Lo zFOK+v@Y4=N_)xMFtRI1NJ8(vN7SBBd!LX2>0pnhHky$V`=&F`gQ|Dflv|6DQxL4e+ zruL6bS;qt_1r0J@Ss0COqSIQNm6al<1R)%G?-vb`?62lHg>dpjB72W;**s~KX!g$z zz!xaI8!r;bZ^|ylGj#Bn4jE$58(G^m0A2t(2X5L*hVjhtPiDle?uITat+w`mqswY5 znLN0L%>t$>BL<<~$?osYvnEf}GAbDaMVCf1l`EC|?s?XfBAJjY{qAMwSu@m8?gQsp zlMa9^9uEpS>!f*!Y-Hg8Vq>RKZ^CCpH4s;KHt|goj+~5JCrQIE#G{t=VneF~xiF$$x3IES0%2r$KGoNl`{O4*F~h+V zrO5bsFGok}>%{jJ*~;~W_-+*SuD%|UB^iBPBOc~e^5$W%5V1JEhji138rC4SQC_nW zk*ZOT!MN+tFDB?6M1mjBC7^}+I@G|*%wGoZ*WLhv5JhwPnxoGIrpB?y01=;eE{} 0: + if _peak[_i]['Report']['P+'] == sprit_utils.x_mark() and len(_peakp) > 0: _peak[_i]['Report']['P+'] = '%0.3f within ±5%s of %0.3f %1s' % ( - _peakp[_j]['f0'], '%', _peak[_i]['f0'], '✘')###changed i to j + _peakp[_j]['f0'], '%', _peak[_i]['f0'], sprit_utils.x_mark()) #changed i to j return _peak -#Check stability + +# Check stability def __check_stability(_stdf, _peak, _hvsr_log_std, rank): """Test peaks for satisfying stability conditions as outlined by SESAME 2004 This includes: @@ -6033,12 +6032,12 @@ def __check_stability(_stdf, _peak, _hvsr_log_std, rank): _hvsr_log_std : list List of dictionaries containing log standard deviation along curve rank : int - Integer value, higher value is "higher-ranked" peak, helps determine which peak is actual hvsr peak + Integer value, higher value is "higher-ranked" peak, helps determine which peak is actual hvsr peak Returns ------- _peak : list - List of dictionaries containing output information about peak test + List of dictionaries containing output information about peak test """ global max_rank @@ -6158,7 +6157,8 @@ def __check_stability(_stdf, _peak, _hvsr_log_std, rank): return _peak -#Get frequency standard deviation + +# Get frequency standard deviation def __get_stdf(x_values, indexList, hvsrPeaks): """Private function to get frequency standard deviation, from multiple time-step HVSR curves""" stdf = list() @@ -6170,16 +6170,16 @@ def __get_stdf(x_values, indexList, hvsrPeaks): if p is None: p = hvsrPeaks.iloc[j][k] else: - #Find closest peak in current time to (current) main hvsr peak + # Find closest peak in current time to (current) hvsr peak if abs(index - hvsrPeaks.iloc[j][k]) < abs(index - p): p = hvsrPeaks.iloc[j][k] - #p = hvsrPeaks[j][k] - #print(p=p1, p, p1) + # p = hvsrPeaks[j][k] + # print(p=p1, p, p1) if p is not None: point.append(p) point.append(index) v = list() - for l in range(len(point)): - v.append(x_values[point[l]]) + for pl in range(len(point)): + v.append(x_values[point[pl]]) stdf.append(np.std(v)) - return stdf \ No newline at end of file + return stdf