From c5a2dd1bfd04caa9bb14ed83521d22d9a487c0a5 Mon Sep 17 00:00:00 2001 From: Eva Rifa Date: Fri, 10 Jun 2022 16:22:33 +0200 Subject: [PATCH 1/2] Correct path of the experimental dataset. --- vignettes/AgriculturalIndicators.Rmd | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/vignettes/AgriculturalIndicators.Rmd b/vignettes/AgriculturalIndicators.Rmd index 3d9300d..b155d3d 100644 --- a/vignettes/AgriculturalIndicators.Rmd +++ b/vignettes/AgriculturalIndicators.Rmd @@ -66,9 +66,9 @@ With `grid` set to **r1440x721**, the SEAS5 forecast would be interpolated to th ``` -S5path_prlr <- list(path = '/esarchive/exp/ecmwf/system5c3s/original_files/chou/daily_mean/$VAR_NAME$_s0-24h/$VAR_NAME$_$YEAR$$MONTH$01.nc') +S5path_prlr <- list(path = '/esarchive/exp/ecmwf/system5c3s/$STORE_FREQ$_mean/$VAR_NAME$_s0-24h/$VAR_NAME$_$START_DATE$.nc') -path_ERA5prlr_CDS <- list(path = '/esarchive/recon/ecmwf/era5/daily_mean/$VAR_NAME$_f1h-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc') +path_ERA5prlr_CDS <- list(path = '/esarchive/recon/ecmwf/era5/$STORE_FREQ$_mean/$VAR_NAME$_f1h-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc') sdates <- paste0(2013:2016, '04', '01') @@ -182,8 +182,8 @@ For the function `PeriodMean`, we use Growing Season Temperature (**GST**) as an Firstly, we prepare a sample data of daily mean temperature of SEAS5 and ERA5 data sets with the same starting dates, spatial domain, interpolation grid and method by running ``` -S5path <- list(path = '/esarchive/exp/ecmwf/system5c3s/daily_mean/$VAR_NAME$_f6h/$VAR_NAME$_$YEAR$$MONTH$01.nc') -ERA5path <- list(path = '/esarchive/recon/ecmwf/era5/daily_mean/$VAR_NAME$_f1h-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc') +S5path <- list(path = '/esarchive/exp/ecmwf/system5c3s/$STORE_FREQ$_mean/$VAR_NAME$_f6h/$VAR_NAME$_$YEAR$$MONTH$01.nc') +ERA5path <- list(path = '/esarchive/recon/ecmwf/era5/$STORE_FREQ$_mean/$VAR_NAME$_f1h-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc') c(tas_exp, tas_obs) %<-% CST_Load(var = 'tas', exp = list(S5path), obs = list(ERA5path), sdates = sdates, lonmax = 353, lonmin = 352.25, @@ -292,8 +292,8 @@ Here, we take SU35 as example, therefore the daily temperature maximum of the en Load SEAS5 and ERA5 daily temperature maximum by running ``` -S5path <- list(path = '/esarchive/exp/ecmwf/system5c3s/daily/$VAR_NAME$/$VAR_NAME$_$YEAR$$MONTH$01.nc') -ERA5path <- list(path = '/esarchive/recon/ecmwf/era5/daily/$VAR_NAME$-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc') +S5path <- list(path = '/esarchive/exp/ecmwf/system5c3s/$STORE_FREQ$/$VAR_NAME$/$VAR_NAME$_$YEAR$$MONTH$01.nc') +ERA5path <- list(path = '/esarchive/recon/ecmwf/era5/$STORE_FREQ$/$VAR_NAME$-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc') c(tasmax_exp, tasmax_obs) %<-% CST_Load(var = 'tasmax', exp = list(S5path), obs = list(ERA5path), sdates = sdates, -- GitLab From 6c928c5cb9c334207ca8896a98108cd78156613a Mon Sep 17 00:00:00 2001 From: Eva Rifa Date: Thu, 16 Jun 2022 13:07:04 +0200 Subject: [PATCH 2/2] Update with correct results and images --- vignettes/AgriculturalIndicators.Rmd | 26 +++++++++--------- vignettes/figures/GDD_SEAS5_Corr_Y13-16-1.png | Bin 9961 -> 16049 bytes vignettes/figures/GST_ERA5_Climatology-1.png | Bin 8885 -> 14335 bytes 3 files changed, 13 insertions(+), 13 deletions(-) diff --git a/vignettes/AgriculturalIndicators.Rmd b/vignettes/AgriculturalIndicators.Rmd index b155d3d..50df881 100644 --- a/vignettes/AgriculturalIndicators.Rmd +++ b/vignettes/AgriculturalIndicators.Rmd @@ -206,8 +206,8 @@ dim(tas_exp$data) # 1 3 4 214 4 4 summary(tas_obs$data - 273.15) -# Min. 1st Qu. Median Mean 3rd Qu. Max. -# 3.63 14.38 17.89 17.65 21.24 30.21 +# Min. 1st Qu. Median Mean 3rd Qu. Max. +# 3.63 13.97 17.25 17.29 20.75 30.21 summary(tas_exp$data - 273.15) # Min. 1st Qu. Median Mean 3rd Qu. Max. @@ -235,12 +235,12 @@ The summaries and dimensions of the output are as follows: ``` summary(GST_exp$data) -# Min. 1st Qu. Median Mean 3rd Qu. Max. +# Min. 1st Qu. Median Mean 3rd Qu. Max. # 14.23 15.78 16.50 16.50 17.17 18.70 summary(GST_obs$data) -# Min. 1st Qu. Median Mean 3rd Qu. Max. -# 15.34 16.85 17.72 17.65 18.41 19.60 +# Min. 1st Qu. Median Mean 3rd Qu. Max. +# 15.34 16.77 17.22 17.29 18.00 18.75 dim(GST_exp$data) #dataset member sdate lat lon @@ -482,7 +482,7 @@ summary(GDD_exp$data) summary(GDD_obs$data) # Min. 1st Qu. Median Mean 3rd Qu. Max. -# 1195 1504 1687 1660 1804 2055 +# 1195 1483 1569 1583 1730 1874 ``` To compute the correlation coefficient for the period from 2013-2016, run the following lines @@ -534,8 +534,8 @@ dim(tx_p$data) # 1 214 4 4 summary(tx_p$data) -# Min. 1st Qu. Median Mean 3rd Qu. Max. -# 287.0 295.2 299.2 299.4 303.9 309.9 +# Min. 1st Qu. Median Mean 3rd Qu. Max. +# 13.83 22.08 26.08 26.22 30.72 36.72 ``` With the prepared threshold (90th percentile), the WSDI can be computed by running @@ -554,12 +554,12 @@ obs <- Reorder(drop(WSDI_obs$data), c(3, 2, 1)) # summaries of WSDI summary(fcst) -# Min. 1st Qu. Median Mean 3rd Qu. Max. -# 0.00 13.00 28.00 30.65 42.25 82.00 - +# Min. 1st Qu. Median Mean 3rd Qu. Max. +# 0.00 13.00 28.00 31.22 46.00 82.00 + summary(obs) -# Min. 1st Qu. Median Mean 3rd Qu. Max. -# 9.00 19.00 22.50 22.91 25.25 33.00 +# Min. 1st Qu. Median Mean 3rd Qu. Max. +# 9.00 19.00 22.50 23.16 26.00 33.00 # compute FRPSS f <- veriApply('FairRpss', fcst = fcst, obs = obs, ensdim = 4, tdim = 3, prob = 1:2/3)$skillscore diff --git a/vignettes/figures/GDD_SEAS5_Corr_Y13-16-1.png b/vignettes/figures/GDD_SEAS5_Corr_Y13-16-1.png index 43e291cfa9295e5011912f10b02326cdb4762795..5d05be945f1dc68da843d2b67436a05a23315cee 100644 GIT binary patch literal 16049 zcmeIZcTkjR+b7yaP!TX8f`I4%k`<5~6ck5LKynZe5y?u90*&vWO^Aat0G3 zAW3qToO90m>)BJ^cm6r2cB{@-ZSDEGYIx^mn(pVh@9VyaKsb2q z@5{;r!cI-nKjb^{oAbM=eDOc~tp3)rB@n0#NdN5kAeCTBAe!Vm^3#7aD~G1>buY_F?didM;msc}Ud-dG1BMl+<8tv(#XgkRbWo zEP)HzzZ;YCrt>^~PG?a3+Sn2fsF5f&xJNX2YOqlf<#DiQUcy*%SnONTaB>`fcr`9a zYqU!j*E4mQdMDwX?Ej)4(&Xue<&<6etk+-v>3lNXQ>w4K^7G}H$U>`trH%)s_EViA zTwK38I>dy9oob^b^(?lH~QP?yZ52><>}lD^1&QxcFn3O+S%@#YwYao z@}WH392_=%&v}fO1hllY0@)O&+Oy47RaNmu?Hr4?H2u;W;euL|5rKj8Q-#wx<}Ez9 zu$C4fQPBd$XvrK=^Tp{N{S(o@eqB#uNYyJAb(m=NXApbpO z+O2fg)_HUbf?mD4!%1f+sv9loh6_@+k@X$UNW>o`C`8&cKTjPpxRht!k~%YEmGrG8 z_d!p%&;v8EXi@vbYk{olVu}&M^);CjU4<>T&z+S>VJ%NoOOlq6k>^oQYYAYI>h4!B zw-Dsj%J^Pen>O&OzrQ~-GgH1)HmZvKQfs0br_sKqM74&To|jA|gF?1L*X=ba_R3gT zOkj6@-ODU%zJC4suLu#_Ur{HCIo-51cg+8q6n2`qt*RQ-RC^?a?t&TVo9}f09tfj8 z&B0OG*YK_Pgl@Icn+%R;U%!6c+FaK!u&WqrAB@)T1uulBua_t}a9CAzIp%Kb}+^0fvT#T{_nXt^2<3Tbvzh zN-DG)RnFLK($H6v8Z0g@ZZdct5a2j*ZTcsv8rxq*Mn*tc#hU%*)C`5^ioH^6~^~r%i)*H(7 zM-CqrH?jBe@oB6L`uzFxU{$Doo|VD;P&Kzsj_Rp#3X$wW(=cl*oE#Ri#l`*HP2viga(%xslx zdIu2^k#=mvwe^aj!9n}2jn%o@5wobH#EOcFO4>cUci$?oyJopPky+c=yFM(n{ZU|P ztV!PT*Uz6yVpQ~k*Dbdvb9;h*bQL-tpqFDU$7xEsUNBv<-Y6Sb(TGi7WcgKHy%^NA zTI;bDr?5%-oTY_@1xk-qEaCk^QS7w7+hS9HfRuZSuH)*;N@J2n8d`$wa1FXbM!U!Q z(qygpJgW8j#fy(hOFcH%7S|WYm%mX-X>Q28fB!ylusTvKH!shrIav!^Yu1@JYni}t zzccT={a7P<&|WgK^_7*&SFRkS6Zo5f(^6vWE!X3j)^tOU^^P9x-1`hPG(*F~$B!Hd z3l4Koe_W12aTs|JM=;6rXU-eac7?OL3+W3#l;kjG(8)ex#`mF z`{PaNf8HLKzxigzHg4$sjqv6IdvkkDJ&(;rQ;+&7CFj0pZ1`=C!=(27q=R3GP5eX= zP5i?o;TdUR8e;qj2c0uvQBn8J(@{n9v$HqXJ?5*0)C*H|aup`yHk0`~s`<;h9D7n6 zUW?kBts2{1P~89Z=|OsqT|q%X?<`9nu4afjp#7@RH5opzwzh7s>2Y6AXhiw(NNJBsNH`oQH*_%tCN@aPT9y zqSR(qPNqZd$|r+w$&xqGkQ*gg0uB51v$l0nf@73FMn+6@rc>ndM8`BfIZU>9<=jiU z&e{+cAZfN(dUMmbVXGsl-n8mVV04>8u6Ar+o(x7TXN^6zV1arpmMR{B)huXfgY2 zgHHT)UZUY8%{2qZQ{>r&p^|GOXT*$JZ92pYt$H7~$C@=Ij){FIwco(Ubs?xaJ%&Q- z{vQ}S^TV}C+`L6Ik8|d;&089;;c#RYTwJd(U=|W`>?v{EMMkN*=IZQ>a;ouEPtYrN zseCupRY>Y27>a)MLaEB}a)l0)vVL^&_h(zvQKj~)1_lO6>M51Kesz|5cwkgDE@9Q# zBR$q9g|crZG!!p3>)`GsaJ<&1OVj)!a5~gS(>CjFKF!6lA7Q<&~KXiH*El?(XjBQ0h@X$A;sUlq!FP@uMBS4hze7Ssc^OHe>f|O488S($GlJ zzF<@tR5(>2Z@Rp=7|DnK!>;1v>ubqf6e08S_3M>3<7msajKR(2o;2lnpY{@)K{>}_ zO(Qz1`E98T`?-N@*}AUJt z^2Oaa>3+EJl`B`;H_whOC4J*BcvDW+i?w-k!TkO6H+pGOE(;&pESmN9o2Bnkaorjq z)r5UmJeBmMR50~v0fB{3efQazrt@LfOG-*A--U4NpcScXcmR5|ryF{0zRA`Ph`@mV?nlNbrg6baX@NXPE< zH9q!R%Xfn9KUdcr_{x`}TG8@0Iaq8)I%t&EK&yuqH=nq+IGz+q#UOh3d^P~dTjJ-> z7cHYZhQit`)HC(TW2eMzop7!u#H#KZ8A;51-Zl2*(h)WVi@KP9s86%8U7a+vc8&-Q z4PC28L$)7ow#5u>%QQA$whXyemR3V5#_@|=8}o{TGw)pLW2yT#R_Ax^rTAG}o16Ip z)yP40`plW>!fDBD332g_>=r%OxeB)EhjZ8LuQvSfqYownK4?$tmG<7TYwrjB&$$+D z%5ygUl5T67w~Q@~$M986*5lp9g$fBO|E`~IIgmqJP2%S2a5TD)iNO7ioExU>TwEI% zPwpGDSF;1^O|mj~i&gR;MZ0X83Mw3qa&0VcH~httKUHQ^`AF>jbxt*9XAZ}eC!gZP zA3Xm(v%IXq1au&xK459;b#uMdoDle$N4#R&8kmgWfXVUk*0?6(rkHt)Jh{O++QUbV1eXgv=;m7GO}rniWpsWKV{WXUgl5}APn|ob z)b*-aPfz#Nm)qxPDS&4q@!m%}b{#%^I8pWE?-Dmxd+w9!HL*&|)f|ln>DF(p-4QYR z%%Jo>I3rO!SZqnU^=6HbtZCzin=)<>7e;Pg<;*U>%V(PQB4?iayB2>)Ae*jI8V3(g zl&Jlf&VynjzAz%Sq@<+n(69QblXaBQ4z_aS*7HL*7Ohu%vhBy4)%wkTJU@|Ex>Z1Z z2l%imgxhs0e^@S`S}4_jNjVm@flDXndRs8QFzy(qI!5^5(9nM7Mkk(YI*K@46H<{m z+LM&j!An`i6B8D`9aF9PC1s63p6$G{)`OMD#8e&t4J?;#>lRAM#kPzGAFnbPJh{KT zguQU#-z|?p16IS9ud&LD^w7q?eDUIw^ZdQRS&S9%gycw3!ygp(J->bM?J3%7&4UM+ z0qB}i0IZWj=AQBmhNVY{^OE!P^JBQ?x^KD+)}oKazkhE=wTYuB>O1I&O8iP;Ptm7FnD3p)-y1`I{RX zuL1)*vdu$-g0k*b93B*`qfGi4tCatr_Z(cvk(;=)n-<0P`6`p!J8``{t6sF~Y$gWX zLI*V_u?dXjGblQdAQBQ&>h8*IvJVy|7@%{Qc#btJhfr^|nhlC$~4JJ%Vdl!6GqOJ8qwA z{HQVvL&G}qZ6&^ z{pjV1>M4mEtDv@_34+X*uTVazjOSHfB`eygUO9HSQ!QWZ*vyX+5KqvNV42s>&d!)* zD+>*Z^0Err+HcA#9>pgVi@kAqQ5H`sb%rxVz|qT7T}`b{tAn8pzkkd9cq(1IT{YWZ zPIXc^CdXsB%WCc~kFY0LGQH<=UI5u-@uV}v zf4c)&kH>bg!*DHVC>4`rvMz&8W9`f5&nYI|%z$XEd)$^;{(S+K-Y`|IH7T|3f9k@f zRyG{%k;=%BzufIQssGtg=luEeW6rH(ZDjj1+ll~YvPdt1MRDI(Szd~Osgf7=r0Evk zOaJ^TtH;OUB)=1ik_Kp`bndNWS;qX&<6mI=C1dpKKo^-x;L zT4HjGOCkVDoW94VnzO|@+AOv&1tzj=A6^HlL`u8L5COp-0)l;H`5-jh*lokY%6hNN z=NORWNim1yc9UWMm?hrm5wvbmY)&ffs0;lbA7(!}$O|$xAJkOeJu1b+qmX}ktY>>G z*X6cyTFW-tx`X;<*9=z^lMx)3jLgh{lOk_rlj|aOr85nv?cUPUgz)8XG*KaMI`Z_S?tu`pJ9;^V*z46@yIUlS^of5VS zCE?ecn1i#j)M%k}87V{Pq6-mYw|Pg_d7+#C?L9}6n`c#9?qGm%Q0A+xPB(^yeV zJ}1v6{ja$C)>s)bG5L|YSkEQLgYr{6IWymnC~6i*%4vDd%=CV1_b9?V_syijX4s84 zbH`^kRaIG7#+aKrY;Wg8(SN|D{9U}n z6rrcwMsyBhXH3hnn)lZaFe*wcs>aFs5xY)h#X-PnNz-?%gp_flo~b#)pT7 zhRhTNfQcr7yRy72xJMRUlB9Mw#9a*l?Md=Tl+rvyR7!h$NcyCxoe3)OUSZ)^4vVv@ zJZ|ZgB7Z13$K`?t#ws~`>EqDjp}^#1ZrlF({114Ba$57*Bs3~ClWP(&N;~*N$h)iZf@xJ(e7)bi{mY!Huk^n*<$X4>o7VG zy!`a(+=BX*9H3tD{I9zQ~JAYDq$ zdh~ZZC<1#ZVQqa~$gJ@L5Icx9&tSk|i_W}suxp_Ft*tG5Qk2kxsEWI2e4%R%C(%Bcg&Ag}QUjeXnt5`Pj!cn8mbgkq6!tBiX9P+9xWd}W9W%a4g2Mosj4D?66 zy8hvoZo=?{#&d#+iRt9YiScpz6DQgV9Um+e3}A#Mso2FG7UG zF|?Yx!>LQnqhIo{2Ev|`4kw2k6Xx0+8 zwNW*dy7SJ}z+3;m=~+`pUg2)|9GHo;PR#XPEEmbGKX>i9NXwUTc|cb|pG${GAbdRZUqAlOUxO2!x{r=& zt!i=ld-P6moVw#lUr@N4ocXf*^@u|GM>8ERpSxH-tF@lqP59gI)PMc>|Hm3s|C!pE zJhMFKd+`3q6bIpt$K14^t>ZrGD^kr?S4`5c{BJui-o#oIA2b@*hb{Mhi$GBQl(9$0 z%4)`P{dM&afQ6_w6@hRp(c;qU;9#TDZ=($fUI*z(qBrw*+ggcXC)k;)1;Pu-K_KL* zwn0I+g9q~L2wTmY3sradx`3xivJntzgz^boLd>m8dBAK2-=0Q_*s}ih*Ih^qGEXV$ zN9}47DF~1D|G7;KvTFb{^5e&k(%uPCE@Kc8p&L7z!S+=p@{z?>+n-1hDRBZT!B$h-SAGKqW2C}ho==UX}h?HvXc{LOpH2C!>CZbPXjaI zE?*d{LWD`_HdIaq-h%`}i|tH^DL`5)G-%89BSgX@xm1fwc@}Ma$D4K#9_a*elLH@4 zzZE4AtPZ$8-A6ce;W739^&@+(+q`iH#GWu<7i;VErjMnYOX19q*0XjHIQ|avwE9Zi z7eBp0Vz(9x^YgwwJ{cFP0IO~MgetX<5eWY5w$wG@LcIL^nL&{$dksKx%b%YxG!FqI z+95eOVlJ&pZX+1SK#DG7wmuU}3*_w_Z0P-12E$u?;RA?C4&)9PiA zH=x5v9k}qv0uGD;;@R3z-^Y4@N7v}G5(ugN1E(I2#Cm}e4tE;uB>WMbpxZ*q|4ayVClHQWGk|#X;4fQx>Qc2+KrjwM5m~#l`H`E>L|uzbHXI9#$nIJn{{8 zZW$ptl-7Nhi}jco8B0KhleZPgXW=B8ZSN+;@NwEw(@{Ky;R=-h@zbYI&=sbO=KWKc zRCf|$=4@x~hN1eVCOIfd2n0qAQWl`AFc;bx#DTx%{C^){K^lGW<1J}S?H&td1HDr`yFYO4q)-EdR$`Bkw0+-S%n(`9DtDaN z|HF?Ht{~;iCYEQ9VTFZRcgFVL^cRf0~G}T4V#56rPI5;KE!SL@96Ab)+p&Q$Q zk&!8xNg(%Asr3~cUA&FE2y&HGRpF4M3+%@!nIvz;%LSuk5f$jeMzxIEnMDr_3|ug8 z4vKKZ8IW*UIJ;>>q=xQ=$bdTMz5@qryUo=0?`9^bI;^g);+I{p1vE7^(f#gs7v($8 z4?*zthiC$O+$`6y+<(-r0RD*9uc_f-jrW7(%sY4P#P0WwEhcjPB0Q>`gzA4t8amAx zAt4P@sN;&HkcDX1KfB0+Lf*bT&BPQ|o<4ZXpGZ!6_kxStojZ5z?6Pvs3kax^$69>e z$BYBQN;JVK770DW?0;zIMWbA!DjMGq5a0(TtFr$dHV$r-_A%+BG&Jwuz1x8Q4Sj8P zcDC4c*&RKmk=`=&p2*=lzSIQqAL#9C%hTIG&jd6ks*OYNLrc1Wrpn5$LxizB<*_v% z4LaQNQN znADx(<5P}&x4ZYx7f)zDNfE3rkWpZ!!dAZ?E;v~&)WbDsbJT$++S|v+W==bwUh3O! zLewyV2DK8)okBei=OIxnPqe9!-&9h%Ef@~+o~^Na3;T|3^n3Q~=A8_Is*0K#N@#o( z6j;xDWy!uBK>UAR`XjXwbIB7us4f3bt%aZ$kqGH&6O{wEBVNm&cswV(LE`)nF<>u*aB01bV6cTCfd)rPw_rSDZQ9z}uzj3o&pP!!IX?6&EN+3l z+9t@5C>l8x4iqj)#AH}Of;jmdGoagRBbHVZ0_$|+CKi#eC z>Th!I9zKL0ayZf&qjhz4sCNkL`Op}RwEYhgmwZ5QvF|-uKLeU zW{iR1R8WdSPPh}bkY3KH5}hR%%yDe~-NiB}9M>$n3e*y2>#s897oPqsfm8OV@!ua5 zQ#^!=EKN}Qp#Bf0uMzU*sfmPS+%@C;~!)d&b=!K?)u z615w7v3bJQ#^zpaR2PmJn&zIpd-v|%4I9NCX)G)jUh2lryuFna*+v^b22y{$%jX8= z7AK!-|Ngb1@CTKRHw*;Bv7{0vbukDX>_r9O1Fl|Op6qbMY=d#`?^gk-h*dFaW2HX; zc9OWTu<+TlXL)$0dcByU&KZ0~XpCDFZn1VkgZ)J7%?@GM!tq;ZN%8Rl35lZT-g%Ny zP0fbxJHmSQz$)@3ZBTcyd$7E_APvLX{t7I~K)*lc*RMyz=ji-!YjbmR8H?_WA~ts~ zUHYfV`8dxlS~@y2^B*W>SeeI;AJ_SHjFB;QDSvxw9Y6@ZA|dD9`}ZLq!`mz*Q95X{ z7Qrx)^E|exS3sRNp1wnWN)lz*{n?v+gkzK=nn5=CocNn>0J`=k>nOaX2I5lzhX0JX3* zAzKVPg`b`xmZL*IEhCE$4b?z82uIqqz680a$>+dtT&v9VhYlS&c1$f}6IiyTCTX+L z-_LJxdAY?QhmM$2aKi5gOiS~QGiS~KO-H_ZrJQ?T*U}YBX=7`fZPZ%%4zrqU?_RTE zB-|{1zq^)subPYJA0b)_c#SuQ3IO9}Wq-8jd!$UBbDkS`PI^&X$eq|7oB4Dh2WAt4 zM>!#0hWaaiXOT1iUw>U*KO|!K<2l6CH`Vt7E3$xai?=+;icq@A?r-S*NLh=?=gXb8 zwQzcZD6xT`S}Q3iDNEcp3)YHQq)5g9X66V9cq4EjX*Z%YTIdFsrJWKz_?EkS37@OG z`xb=mY!&0`H}Se|I$>G&yF@2%L4D7&9}nSU=jEM7M4(Aw6K+RVRcD5gufrI?qsjT9 zoT$cpnX&Qlbrd+Gn8Pm<0Y)kJhcM*ZKgPz!8o8$tj4tC&9s&fn*Iu3Umn;KZ&thoQ zH3Q9Q&M7Z1*ExAvAndwsfnEL5>Oi1k<42n4b}Fe2D+L9GiaEIkZg%!-lhh|pgoTDC zM4AcEg(F7b;O*zvq>vqa4WnGGV{>z}>8_gQ<;XH-AurFme3=4-&L)fC=EB!$a`^e{ z7njJNM@Oui&x(H|MbI|q^M|7{N6y{-`GWS~V#GJp-n9-(L?Ym>w$nr^ECcqN>p-tH zNqUYo%M)@lcRAR#5 z5-@B0`Bny=Z2&3Q_+xW@*?gOX+~t-`|3&>)4h}hU!wA_WH!J1d@4_MM$x%(}0#;C0 zQi|CUfl&iPpVG%~zXtTU&|5-AM_K>d${uuq&OQCjIq+Tz)=2 zX!q>7^jTt5ytneKp!lOHz3TF3lt`51X`nG?*Sh5->AtaoNM8LG7Z$2;@(yNa>Wwb| zW*|X@m~NPjKr>JU@CmG^1k zScv%L9OPYjwP3OZOB+SQvbD2S=_hqdt2o@Kpykq-{Xsj=HaDY`y&Ewi$hpDr*hB^Eo7VL2CJGn z?snvty^;3{I(2!^8S>uq`}gnXa(JN)tQzx5sEd*R-Y-kzqtfE@zxbCP6j&zZC_8i@ zR@5{v$5ph{ZcZcBew#CJs5bh6u}FJly0kaAtk&~{rXsGv8@akd?{+hD(3ayXRoo&=_&O2`q|dBf(=3&++v1w5-+vlfgjOj%;~`iUf`iL0X}8jQa3l zK8O#b6jFTA<-0)4IeP~lc&vw)ao?ul)k~K6r>wj@A*l3}ff_bAW<%KgXw>(f9uEZe zFlDqe?p7dF_z;)_BLD#8g6Z#fNaB9;qRO40cdC_Kw}3GOFKi8F_y+PC)s>aVk)*rH*n$R#x z!*i<@`CO0f^&UhyimZC2ot##7A@G2tE=PwgM4LAb$49WW*+Mq!KU=n5*#AH6>8wbB;B{J267GJxky{4*P0H|q77Zmt zyN?Ame$>iza@0^|r@irHiYGSXCua@Rz~D0Oes0H5C9Q#TN2Q^!HQVb14+~xA*yrpy zQyI*O>3HMDjZMf5NXLi0dNntelonNN#Tu{ZQz)NwG?nIVK_RiS?^m@7=_7?-xixmN zsM4(So!D={!FSA(%f|7T!YN&966VKjz*OtbGF5?jjzLbsQZzKIMAFi=9Wd_HJw=F- zt~g92;d^;``5QO7KrKPH6x|?80`G$JQ!7J7hHB z>JEZqju=0ALqmg#ib_kemhId?1=3J{C-{c|50@4f?HzTYf(=trd?m$zfHb?iyRkLx z-=7_+kRm){aQtIWRFqsDK-|1MfGIj8L>7f=@4ecMu5bo&Q`qOEhWE8b#zN zI4WwZoXI20tZ6KId&3xx?0>JTH{s#qQ)2`|@XyE-#Ke)2kwa%>ExHO8QIr#PvEM+_ z;KBw52BHMt1P1n#-a6&zMHZU~B|`2rqqx&2a8S&;B;BZ<`f)pU`$R2NfWHgVRmmr`5ot`0o+M(I>;vO0f=XSZ5%mrq)XR(-F~F5H?vkUNx!tD@PKvc zfA_OLVq#43e6NPaiim)!rsl}xHJX+5SppZBn0f1XJ@S)k4+`CjPKgDOJXzV<`=CvZLkagn+OITBI1I> z4f-!bLx-oQr>8hL_B@p?Q~EIVJoaCEcz)%EHi1y-`5)%VjgV6t0KnKt^TrRo)6>zf zUkke~Sp-(tb>xFnibP%rkXX{((>QV3|KJnqKJf1T@_>MUk@ADOm@hqBOSymwg$3H7 zHn-1iY-}W#u!wK^_0AODfgFml1GD@-W;UW#LcsY5ifd?$D^pUHfwC~VE!3l70z_lv z)kKJ>Jnh*v4IKsupiB2kei#cf#oD~Qyv&#J0ZsIm57M6de3!2+`aHlQc`w}m_3g!x z`nWqBj-p;(eV7Ud=>?BUBNK&X$H(tyK7IN$57H9ngT*2ZWp3i&>6-_+W+%LRr>eQB zFfdct=O@&yB-(=0jEPuTb___e*t%b~Om)f*H(PNR0nr!K(g_V0kOcrssywy_z3?LP zxKEDh&qCT7?C)>DQzMSYmGvCmsp}pdJ+&X(<@~pt}g?;52ar&jefvnz9%?;qwXsV6^s=)YMi; zkR3ldu8md0Xbk?`3S|DQj4Xlr%4(TA06M7Dj3t?y+_{i@*Jg zNQ>N3&TQ%SW*YTrLBVFV!Gc%cWDe2Oa}ry#OapUfAtbpJOENR}=Fm4Y8Uwnv`;P3B zri!ts{1qjk)8Pb38XFHB5fO1xhyQ&-!U5(2$LXmxI!^CzGO-T_ZbXT8|!<=e- zAD3X|bXJZOp1yeTA|6VSSeUQrjinW_c}Lv`ZK_cay!SHsAw0D3)=alDJ@#is1-uPB zx&ksz6Zl5J^fXb}X23g%!yynT?g1t>6m2&=gfR8|RMEU*hzP>W)NesU+JAs%qq^`g zPo{Bg`x~Uc4W+#??UOar@8Hj{zVskrjy-P1ORYN8!s184=V+wG=YF1{nS}g?Xk|KE zx95-E*iTTw^dHV%39MN7LdN7YBJu1bqYowkq7v9N_HA})ZyY8lXDZ~e6N1aasK<<} zm(yB8AMV<{8;+>lw!_>}(x5)hrTOXhm&X->Y=wn1!+Ga*@OeAm4e;v?Ap~2g_k~O#@sqxCu#l16x7Tyb$6+ArgM@e}MT&F~sN# z<{|Rt4u{I0z9kt<1_qEjI>F*|Ho;^dWvMXcy}Kn`%fLv)R-o5Qd)L(|`CgN&vh=L0 zsv2)ll#)Oc7D!fEMFmW}B}E6hoU3P*DU6tHd@4qcIVKy_ z?pvB!=0V11vv80Q-Focl3HPRd;WVeLh}Y?J=RVVVcp-_ z8(TMOSIHiSDv+@5BmQ>62qq=JLKSZnI#k;mePf@RIi>1o*uGYqWgA=TNhK)^Cj?UF z!4rJ%FsmVAPc$X9K>Y!Lgp_CV=ktzG)d@aIsw-FmkUZ)>P)?N06$o(H^o)!QWZ#0r z!^OZ>J3DV3Z<;tg=&WIaRO}SK7eJjwN>t*C^V2&>;9R?QjfDDUzERcUd5_8VY}6G9 zb^taTaxxElOEBWvq$kOm5}WYZLFQ~AIdp;+PC$~#hmkX5Z0p`U7coLfCGGhoW@6;* zYW(%K;*SyBmQ8_*k|L1Eu`;IgP*he%;sW+FLwX(Du+gpPjm5JW>Qg*YkTZgUjt&}r zZE10YQwTSb^ed<(x;h1!!Na<-ivTKPb_FD}X7rdkt(VtV3&GNjxhgC-z%A*a6f88L zvzX&l5#n%=KY>>+M|n^)6H80C(c@67-RR(uOmOX>2cS}rNWO19|L_-9{-|`PeT%NN zH)I^mEKY99qN~9w@rr22#6FG?Lz&oy(V%NIY*iQpkahLR(_xNisE#WSdT5*!w)z7} z#%<;c8LCP5roj_yi?bLJPFP|*F9u)Eq(eqVCb(nVsZVl#u86a!xEMu~iK7n?j*jYt z>+kLDt^8tNE8)%+fK(Gm4UxLCtTZtqA_75N7-|#>TvwpH5x4kUj9%;Kj!_ z;m2zGT9wF|gRWsx(@dh;rM+WfP7Sw=-Jzs^?^3E|Mrm$bCB zBX8t%YkD28KuUsnIyXPhO$T6zYnRwedPOm{aSC@7LwyJHvHi#P#;}xZTslOxO3O2H zAVeAzpZX)n-p|;)f3-0JJpH(5ikLu@w4zoAwdv zc8W9j_7t*nlOWHnJg`_&%czv-Q2LlZzVpZaxNbpHJcJ}y1-=wSGie$4T-MNx^co}e zpAg*Rer5o4l+(3%7Jhz}XT(dMd(}ong`a(su}1M@$)NGE$cLN_^&1`-IBQ%9fi!)< z{(EpV<^sC0*If69y0^ODdShP;Tfe1#3t~rVGup7Fx9}X@e&#oVmSfLywJydBSY=-8e+8LWcvkS;a0>PF=HV} zWVHIHuML#|g}6=^w;&}6lLF!q8w*P@dUj=(lSHVY40LPIkYQLGPa)Lo-i?dDht)MI zIHtuiFr##~=S2xbSD%!YCnxFmV{pMCLPXdT(*nh225MlQR*Ytmy`Q%y>?FO<+qZ9_ zR^>pZJ9R1s{Y_DP_B-i0tln&mEHPGc`+EV!2DAjipyO~f*Dmmu{~P-KPaORJ^l#f_ aI}Ze2v9<2)!!zH6YgZI5r%K;`{NDgqBXKAI literal 9961 zcmeHsbx@mKw=bp8($Mw=iWS#F@KT^yp`=jUg1fsFhhl{`6c6swP~6>9q{ZDG8r(I( z<%GU7_x^L{+;8qXbI#0}@0r=N_Fj8!S-&NFWfH8UAo&D92Ef3;cp@z&4#vQ^54meN z_tBiKr?t^&^5}z=a~l~M<>loaA0LN>g?W2>lai9Esi`?QI9y#_ZEbB$&rA;u4J|J( zCnO~N`0=B%va-9o``fo~()b!rF)&{6N{frBdd%#kZB~5e!-G7EEi`V~;k&$SHn+ta zBACbA*{+w{@{Xhaj403nNy^I7Mcm!{AAJY7x&X#}TOme1QqJZh zWP)ej2xoYxJkzZy2M?2Lh#;H*^C-xR_n<-Vbl@#Mc-{O_EJJ0|@UwOedEzKZf`dkf z6mq{rY7T|`HU_Wg3?;~DBTINFXF+6qJGp|jsy{GSQ=Kf|pvKFB#=^(Tl3gL^T4H&t zVhFdx$9p>UMcX#7a5muUN5WJ^2$VJdH*zVICCkilPWoD*L0!dId8SBt7XEBV&ZMz> z@S$3KMhOMS1)^JmxEh?X??iUt=&V9e34h8W%=MkJ)6vK8HOD{}s9)O7o?FrU%J62r zx|-R}Y07P%hn}$1DV0O>B7XMDM(8Z#>E~aWTQjg-CbmXwtC1rn(ao+NnDo4;!SYY&4E^_Dk?h<=#POEvmFrHE9NC{sXakL^YIKjkXZ6l`QQkOTdNC%rA}^Srr}*$>eXZj}N&dn2 z6YR-MiCSSK$aynsV~MNxC~+fq$3QqHZ;8viSFe}WV&ieU5J%yVi<4QGw_#jO;=@9Kz$7ElgF=tF@T#rAk(~6=*{0J=~CLN(R1V9Eq zD-$0z{9I{3Ocp$H-)A*yz%HzkdGfAyFU|Bzbn@rk@|P?d%-5?9BqH4)gBXc@Pyd7> z$y@{ONo7i07oKQhIYysJ<9hw9@ifydqkbFI=YnXcYq)K7@U3AumZM7L%}kRA8~#Jl z3oe7yWuCe@89m@Cb9Ez-Ih`DwfBoWlZAzu`OE2lLE}HZmFsU!I7_||pewCOS{A-S z=|;&^z84AG(6Q8jjfZ{;ZEjBD9B@>IRPQm^1c|f|12P__0JJ7e9Hlv@62nVl)j1i_ z7N5v(FlwYxQadF=oi;gowjXnUztM&d)%nG94jkqhn2`CeSIeX=P*x4b zYZi*hWp)10Qcxz+ zK#P!fPMTf_v}KUr{T;lXSR^S)b3UGWf*mfB0-mJApO}&$=kdEQ5K4Ov}GGrD#LS;JBb8U~H-;KpJC{eWN?w7@w9OUX-Fvnh8tl z+7>$Kmwx*`(;M3(hw3P2WIbJShBAF(a9 z-+y|{`+Z#Zyn_lXi0$dXyOfbUBsu*Ecf@#kY1KuJ02J$p^WL*gsi~>=Mcd}nuozM= zWt5G3DuiY2WM0W=cnr2~-72pSA=eilE68s%IVYf4{0>^c^L#6%)A#ussGJLnDE?zg z<2#7obkjlfwna$5@^BKZg{*gA#0*n%23Nvxe1pjZlXJtFZvQPDeRL|VorKtV!wT64_stGwu zNr;`-RSj0@nW>X-t69yh7)&L}E;90yV~e%}eDz+|AGi>`r%jR7hX}z9)jXIDqW|Eo zO%*vJQh>f~7Nsqa@<`C^&2MehDNYlNS$Is?d%a_Cd#Po~a8F+;by$y%i5kLU&1pC8 z!FKj!B#L8k5|>f~LC6khyFl2ZAmZ z>9imBHM}fL%6(1N$@?~lhd_-v`+A`%68KmeQjI^L-?XIW!Srw!ZNRrbcM`X(?lGFr z=JkR?6J1kSuEum5siG$O8j~+ZPI#@Ya%#u7dB07A5&SqZE52f`6sKw~2^t6mPE8`J zUbMecq5ZwjY%1?|sw z7gWt>HgqG7K)rWWJQh~7!~UV1@y0?0qV1#UR>#3R!o=JC>JqeQRGV%4v9;hp^9NSv zX`T|7X+)uh)vQ|A4mb+%P8B8?;nE)$buKrp4N(?iXCp~DRVaJ^GGjH_DxNz@R1Be> z$tHt(^s5U<=2w_;J2&T2cC2+^*3wb`9(GyE(81XlC2#S$3AdqlxTMWdvNJONLfGJU z(IteIu{t~!ZE_zPmbrE8LUcGZ97pWWw3@$C?fG+jTRh<1&4U)Zndi~9W|I~>GVMzF z@m5O9cQ3sC_NLCO$sZQVip9{%VSy+MyQIr|&epsuf*DlBXT+mVK+sfh$_-n8&gz@N zX?0iurIqI7mRCMINl)EJSc3IS=IX@x1tj3dAg(@Sf2G&OI{^*yYv=v?B=l3PP2ATo z(c{{iz<=cGGNDS_ENI?%ddbklR^_}m`u>#|@!X?E&g56Ss5YnIP3UG$&8{XGQT3+N z!&x?*t@LW|NR!{7t{1&?=TW=R4_QbP?MGy=yiVUcwlq#c-P*h|E;G`{htHKE!Ssa( z!gMYmF=Fk=(Ww~op3LD2%WJuP#l+EFM+<)Ay|Ka*gpn($l@6rC!g6_B4qC+uZ>D&k z$TGE6bWO#HMMTv@vVSH*R1vsw-*(2B%2I0sr)2hrjM{bvUtSQ`+r>~{ld_QUCrp2)FEnE51&oDZj!l%HXv zyLHG5R8gpoI`X2tz;sMgWiyLMTHBW(ELbnX{D89c3vPx`UQOcf@J`w=F=cZJ+;{ zkMxiBz^uqWSP+7x{wNB#V<@BhXa*Vt^d#^Ol0OVIKz|r$xc)HE;Qf2~AEf?2p!FUQf zkdXvJzD|_e-E*ENs?2}l5cdcIUfGxP@}%8ijtwAOIO;YWxulN^Ak104dQnvRf(RVh zCdv1_=)o}vC}rfB%-NR1Rh6o%VQBt^OmkmJHhaBE^U)P(wfXe1B9W9 zeX5EflwD@8sOT{_jkOaC62-@Ne6`1kbs}R#4kEG?2IDlnf}p1&LExZfeu$Z2M-ZXM zA_Y*<4^xU3hy`STV4xcfBn$7ZWknKVqFVr=NG?_j=giU7@UI;ifImrDjqvg20>C}~AQ zhd4ipenr%xV6F43euxB9=tvKrl?KxCI&H9&1(_Wpt||ub-z75(owBZz6>)rsss!P3 zGa5@|;Y~673D@Yao-%HpWqBQE4OJ8}(xJHw)e_>6U;=z5xTJzCB5Nq87S4Ysj8{Ip zR~#lD;HC^*iigXwR=u46u|L|-LL4Ray2JLucaLR>op5HZcciX`c=|(qHK%cPx3u^c z1jIvj0AAflI*F2UymVlTW+%)X{xq4x6#<8(#GEJ%B>5R2s_?=eT*VfMa#Z~iAkrt* zcQ&Eb2At7WvyYs3K!ThaV?IM*c=Fg+pE|I7PJbdEv-&qDGpADy-l2Uml5NDf|w*f3<;uj3v!?3I25N-$)QWEOv- z?uCtqzm_C{Ge>PV;jv?1CBCJ$ZOp+>^;tv3#WPRbJNKb~ zNw76u!^d!k&07jhrbPAo0T-jF>^FYDlVNzDCcEEh08a%Qnnl=SG~12fnmGUJPl-!-Y%9dy$~u^eP(#n)#YUXmhhdPT-abfO}})ZEGB|isPXx(*WGb zj2{&|N-p3kIo{s-i0ZLLKrL3Q^>t(87Sz1_Q~-n3lUNtIcCv_%@}U*=A5lq9BEHwx z4f`;W2&Fzeq{y>R<2N#Y`W_J?WvM?FNdMvDwE037h%EJK)ls|XS{1`{A$r`bpQDI* zRh9FuxT)mJ>PY^Q61PUjo|or7S@@(D%wekDMs^KILF{wbv5z0V;765hePl$3fLrqF zZKbH3HRbpeUmY4$8U0*-YohGc5EDw)VYBX0qmT-kt>=vIQ{mTk5rkCM28}WSn+8OS zSQ{82vb?)C*)O&Qy8VAL)A=m!PiD*LFe&pb);EfoTIy#(l*G?F2;(olN|(Kag!1lM z^Y-S)l^;4j&Y3?gcla9UnmMvLkQNl*c9==4GqD_JBV%MuF48CHAd4`Tmv;~rw8nOu zfMR|xV@bKG12#aVw>!n9mEtvcH#3lFLOdLT?xA#?LS|~sH^x#%)EZ2ATU#$>Rp!@h zw;J|tvQAMWTtbz)qcHU-nYqp@X#Yg#OBb)a82Mf6cv{%*LU)rFqcwlUjkV6cVm znB2nR{sQjR2ug*4J zr)UA5-MW0N$=mk4VC_GWw%`v2YgCj5Pr`@e2}{Pt(ywfvPv}qG@CTUIToQk8)#lG= z-+<07S*?_44+^qv9U+8e>Xtd|2^fVM9L>d?YO_9M=bU(tbSK|fVz{l!c?EX{xR95*Rh%6UPSia3iqQ<@b zY9GuuruPGpy)@}^RrZ^>lhy@#pJf1!Gzk!aOX*VHc^~*;gzxa_b+7kME~^}{$izn9 zUe=&aU4SS)Z_S}Jivp`lkx}Q6J#(?pq{`zI+^Lbj>R*(AZC0b&HWT^{zU3;RvSw;^%8MhU# zt^h(x+@2z|0GIal&BZKd9p`Ov$U4&6c^J(xnfwsdwA`BpicfXBr&iC4NKo=|`7ylZadgmo;FP{EJtZ20vBdPnrZ zLp{@wjO8EYi*adtwjN}?D;K@&z_9_EhGh;1ifq%5sm( zj5nTB4jJ8RmhgMgD&Rv8xyDXAJh^dmQ*dZ79pJoU@dYOjo5qn7)oiosL zEJb21V-WOHmT7P_JXg`An`v4JwRs5?9$!wKjpHMmC|z^w83m?`BrFjW3lp7rqHQk% z6GApft_Z<=i?yio-1dGtwOwe4ev{j=Uo$%dK64ckHEq|oE-J|Iffr6>zQR$8^Ws>)pR$gk&9XD!>MVVhN+$-XS5E zPbL4XLH^%Wh>};Y?zU)yLjg(jjF7%4ZaVa`$mKtKB1)Km=u+y+OC^yqpR))`*cOdP zB?_6Q9wEY11aAK6HlLm)aZS+j$ETDDe76Sztuk@y*cN172uSwJ9kLOI369G zlG;ztk|PlrT6FK==@&WTTMQ)+7MCt5m7Wv3j&CpxW;brnb~Jk~X)Mq^OI8n#5274! zFPKZSKAD%w4v7%H> z{nRQLAz3tA9J){$Jm`4qnXno1(apxaV*9<|xLY%a$KcbLQ|i1bg~9z)GCpr0)BR2-{lw&=`nt_kgoa$-?$ z8py4Q6tusv09~`fk5r^^kOrXyIatl#jZ~nO-2N|I6$^*Nz_opE%;@sy%;0r#-ka%L z0rc|Y~S3~j|~>TJc=c6E~t8A8!PL*$QmB3|@LN{{%qmxfb02jmZlOI=UL^($AEY>YBh z4WUBmJz#*i=f=xc5yLBrv39^3|G!i(VY}TuaINd>&~R)AI-}LX?T~{u zoM?l+oL59ihlYV?%*0ba0jE@j)@e8{!ZfKE(38c<^KC(zVmIH)gx6;FavhZQu~UW~ zcJN&x`eL3pj6vPUA}4xcCFR=(MrpYKGdbfYhae(40+=6@&B|r0omHmPH??n7-gvwI zamMb;FutFY&zt>?)!kxv_C6~_XNiwdu{@ML;dhBpT(xK*DZp!N9MN|^yEt~SN*X`?RpUH=o z#7czn1>c8KBt36sx$;R6rCCscXsXBsf5(m+mNr>o=6X^at8ym7sjz z2SRNu>UacugL8ph9%2(lKlnm+mlY^;NSa_QhKP$6 zA|MU#F%%dX$oT*J|h9t}-3rx%e>dg%K ztsSKdIB`=2?h)tqIRb|$i#)vjh|UOzg2G#_1*{sI*E(54N?(7^?mu{I)4Ibl%!L<^ z&o4!GkoX-Zj-URO-~|n>pan`YuWaF`mN`4?48H=0{BRnUy8&H*n?3*tR*typZH>pH zns}N{l>*NYl?hI|4{sU- zl!zFr>%YN2j?OXl@PT`5tMxM%azG*Mm@!c6+K@;&nRb3v;raB3rx!j)v#@4{aQBqWp*VUoUctrH1Ah)n+%iWY&r zzTS@+D_n(>XBiBZtTjPO--Q?%td3g62sXfBFb!U%67JzISlIji0d1R71asJ!31&|M z8qRK3S&-Uwnv9zEWbbM+vok*C#!(a|2Q*`GK+v8rP3 zXSIDgA;WwqJqp%q9E&1I+xQH)<;T0cbc@cil1)W>;jgl?0eRAZgv1li_b&Ei0oQaj zt(A{X8Vnm3LW6#L7E>v<#Q+YI&#j`jiX(VlG?2nrCIPoM;$X%}5mxHP7-TjJvViAN zz>!l$VPFMLNb12tsB|-z^v}{kfV6uQXrAff=j&V7=q5f+WM$OeUh9j8bR46uX=6Sw ugrKcwq5;h{&{jh((#H=z|I1`K<2NTkFJImJC>e{p$52{ALA(TH81Nt0F1R27 diff --git a/vignettes/figures/GST_ERA5_Climatology-1.png b/vignettes/figures/GST_ERA5_Climatology-1.png index 90128537ef623ee0e360f271168a38248999a2c1..e5c5890f6c04015a960e0ef5236046e2da096967 100644 GIT binary patch literal 14335 zcmeHuXHb+|yJZ_`B{fQrOzTmS7?6yBv;s#^$&z!<0+Izh(g=sB0xFV0a+aKxC^=^V z$(bfL(4=|#d_U&isrfPY$5%BqbE&c|rTcx~z4sH=v)1ar6`skEUu3w5Kp@B;|M^G> zfjF&B{DV9V|8w^&Wgz@@;l-aCjtB&$0rAhNG>KFT1cD9m_>rWFYxL@vo3_f%anlBl z49QW-cG{9_ErqRY+DXpF{_$D5dLI5NaW0YNidc9oEuLE@KjxZ=O-V`7 zu$>8EEWH*O&xUM^o?uPw5OYyKYfQUW9qUa#u_b0EJ~7!=IiF>uKrID#|8ZOmhHrkVX&wCnGn>BivYk)q7FxVWgOsJOU^ z=}?K|Uyhwed()w1Wo3leYce4>r%0{tT**+)G;N=hm6z0b7aqR7R)NDW_BM5PcGlG; zJ&h8$CLJ_X<>68tgB!1^(Esr!KAz3a&JOPPmYngepC4VHhTY81_e+}ZKY#vg=NBFx zZqb|Kv^vs1I9TE^+hQ|X%zV$Z?eEL%EekKGrGz~R+pq{mL0j#7!-g;XR?dq(Cf$ir zySRzlYyK^F_4FpLdVP)-_tcDP?d;6uC@3v0%~H=#j8E-oZr;G-2X^D`8c9pWuzpUD zjEt0`x5HpC>gwu}q<627O-@eEcgA&!8aX&Pn46oU&lTE^5!xc}X`WDBlct`%uB4i+ zS;WzhWQ?rk{k=X}ch2()IkO0#n3&Fl@U2@pZd-Hv&LRH(95Wh*hKA8GG5xXoH<}1L z%b9Aqu`c=2N%`+W+5MP95;Gmy z_t9b-+89`*wLsL?>R}zF1oiPFsSCBesq%Emn+*Jx_4V}&W^L{5M@@QW(Y&=qE(fj$X(M{*Eg@(pP3TBb*|R;R3S~>tTn5t^^YkA&h^F zs7na5=zuZrY|P84`oFuAWMrhI#9jz})2MW{M;|O$4(7zz9W=3eXW0yK1~{w?<v5hHllw5w?ot%<#o`j^SsmY4x zufP6+)3+JO(y+8#`s)md*ZwSBkr*eZ#rp3bO(D#gD(N3azL~avIYUB`7Rs%Xb_W7_ zyvoDF%}ob_e?HDzH{Z~gFkWUqE&b?G{fFyvSkz_qXDfrban-Li3QR3Viy*MfhrcBY zq+@(9QVNSuH;z9JB2mxN{|n!_G?YIL8<4A0frou^b91v;C81za5M7q25OrR0g;l{B zp*B)h8xzYIZ{k7$RY^(g|GhA32HwFquC+h3)3 z9p|28-@S7ukBppMP*9MK?Lo}a+FG?`e+HcP%5Xs{B+Q8R#>U3=gHI9n2s{`YtGU0n zM^!)F()srJ%ZCp=#Wq^mbr5`8x6s_&neaJFlBA6qvwKgX1b8OQn8aU|ju?zrx@B&L zMMY`3w%@MaOT^z?TU&!|_iY%Mul79H=3c7xB`tak;jVcA$6=5^Bj$lOfBd}2vOhXH z`W30`)ZgrdzSa}fRqK|D0Z5Dt9TQVl^<6I2jPQ5w;)}eu+66~E=1z`rCw+}Sd`PT< z)pE4EMdYJ}6s;s&*W*WAw-*!@UspDFC(F?}DV5lc6-pNjGf`82Efy6PHZe91kC?O` zD_4yovU>t{j%IdLQP$`|FX6jKp!t(ZQ}>^(!|x<^D`HNR+MJ-8ezJSb~Ae z;oe4G-W^ufxBwdNV(a0}2sPa@&jaW2 z`6E7-9v*%q&C|5Vyc;tiEFz-0eE5TeH7&S7mO@@$-l->?Q+a!pFwJz1dC~s+WaWA- z>CzuiML&4j-qNj9*E)tu1$Dq=<7+}ff>w!bzjtxQQh%n|h9v#v%UeIMh0bX0L%F{x zzAbUEa?i$R2-BPLR9Q)>F!ufX_ak44Np9XtG9-K;D8iZ)9IpxR;lgR907QqEukvtjpT3viNK7O$*$1miqNPcly=z zJRzXCvX2&1;|5G8imcSe>5YtxT#i=P#wzx2MY9#sz3J-aw;t*mB>wC+(CPJ+q<;Ns zj*EmdmQv)scaM}{+|7#rEM#w3E&p|*Bur4aaQB6_cJzQbZ)72J5ILj3&4q%xQ+`P? zgD7h1?eAZG+vmHX;)FCmKN$suF00(Cl!+WJe<@T5{XNqI{Azc1 zH{mi976s)KPPW1YH=xUneiS8OQ#f0BxY??C08y|uQDY|%1x+*8=304K*;u)g+-jDD zKR+Z(l%U;>8#io+@+3v&8d_WB4%Yig5>Sb%na?L)ZM8#vl{j2~NJ&YVedNdL#qsK^ z)P1X$8&eI0vTe!I7n8L=YqFuXLUsXq=%qN|H}ColpP`j6hgh(|3y!-!R1C*r&-(m` zFW%eV*Dww$ZsNt=C`B1Ux>MdKFV(4V);xi-Pctr@d(W(+w7gvNer|5A9aQwISFa`} zCT@;hxBcs%?R--1I3I%l%ZiVoGe$I(Bj8v`FOc~&VT&C$AqvME$buaiwkPJm$!6|%i+1n z$pistm@N_9nZZz6y1Yx5Z=c`m^gd1!u*_BRqhQ`{INXXm=@YT)Orxv?__Q4QyCY6Q zFKYYv07oa@<40MG>C4g4^6vIc@jf}$WL;IafO5Dgkg?5fFl$wdz2Hk(An9{n)9k3< zi_6>lq_XuT+Wmc&WtGS|rl+Fn zhG0O~Ls=>kr~T958QPLeUtQRUlpn@_ zc`z5MK!JvahLm=^^Z9qEM1#7wBhvuQc%+~Q6k7jmZdPQ~u5wppH2w2W|Kf&g*oEHI z(C0lJF``A5j*u`84tR*3Vq-Hi*88sOsc~LCrhT+If&eCu#=9YrEJAl=%xtnAhuye) z_n&B$^02hba28OG5pfzRvC9=_5qc3nFn@rn8SG4~ah}R9883(MX^#^8bWdZhH{(Ogqbj*bEYHhuOy>Pj& zPZ$axG|>c!b!C45X6B(g1}SzWpTCNDA6Me{CQgoBst0>6pbPCTSn@SCG}sJ(GdVdv zf;}+YJ~=#AWTUS6)KHT0BRuMqE@>%k%z=-w#%u z81LwR=gzH66b}gsiiU<*>|JCxQN1J`1@Y|6+_1MXT~JWa@12;)0ncaCpP?$*uyCSb zY0$tNf~p-HbO9h$VA@^|F$LY$9x%7W@t)B{wI@N+Qp~^#qr#f}T87E9cU&4G%O|=S@J4C0YYn9~Cs=>B+Oa}{M zQGS%Hx@8Wt@HBy?C{YR-Z%%b61#lYaEb5j6S;Mf+hlhu}yu9#M3TCdU(_1&RCs!F( zsY7e8v$JPXLLb>$%FG3@_lcG{NIHqYOIuPqb&;IBKW;?%5#XR6ik2_!a}dFuU-*}y zwbNnh8yluRwXXJkBm%@ zNn2cA&h$GcY=9K>uR4 zB?8KV)d8iCkI$#udO4jJR2u?lJUu=8D?TxczfzZYTw`x!uW6+C3))P@=hbri=|(`LY2F=Q_**UcHGKmE<)Dqp$B10V{wXIXCo9{t_Jfw5USn?{ z`vxzs!;d%Tq?7K7i+e$9vLr_Rj?4o9_bEK>hVn(e;+}-Py}dM?MK zBRwc0@m?)yUWXeptM-FAI)PdbVfrc*d#9+VNH}oJQP0PCFkQLAfi)j1-=TARwb{%b za*qj(b>83l$fH%r9E$OWiN^dOw{KEiHAy?dN z7AS(_)x0=f?Kx^V&MO(`yi{$jI(`sHJ2O&bweU$35H*7H~1)yL0FBZWkk&$2MnM zi#$^L(iD^6U;yS8EWqE)J5R+e`_hW61~Y;IxX0Zy3n?ktEw&jgiW%>8--%u{|AuqLxt-Irjb7nvKP|6@Tphe%-Dt>|l7GNk_s|277zypY$ zJqTvi68m83q+@LbZhI#u<~#bY-N~%s{&BE%zy$NXl zZeU2rfWV&k-ftgBG%^OhpPijJfFJi*> zdv4Ad@Cfl@-g3gLB(=XeR@_$42e3^ru4*AEDb++J)=TtZ$1un&fjcfJMBQ-U(>#BUpstM_ux+<~Cr-H5rqKK5$qTBoJV{bQ5!~D`mOu(uNe>hC_Hx8R73L^f*U)< z4{qWvVIDqAeG-{?lzyc0&Zn0(zn6lrv>UjDb4wXalp^K#!=eQ3#t(M}^?XtC5=YxI zte#&bsI!SB=~S1i25(0vyaDf$o9O z-i8Bto~=ow1O{`)0GJi&J>_z&yBqn0TP|g9qPX(huBwT3J|(RK0NF z0`NAJWMtdx>uKuwhR~~ErGnN&?{+`=K=ka)gvq;mc&Gv01t41Pe2dV{eJStgLBn>RjZiy%`ZQ2}b4yE_xw=*8hi}bGbhQAo0#Ix?QF8)d!Ks`= z?48HQ#}IMQ<~&yNn+Ur&Rw^ngMn24p>*xPU;f|U1IJ>QQ~PCcH>O-hDv(Hf-T^2INc2=k>QKem`1?kI zC~l`q0{xh0Mp1R%d?GbIoR zf#;R_na{IF&Kz(AeE9G|smC3G$V;C(X8>f(B=gktG{}Ud7QHvn^#Bj#bm`y|ik3}C zM@L7P5$kAqG_6Mp(*ZYn;8!%h89%n*VL~9dIb+!FzuJ58FXHlJRFtUWTw9V%C=vat zd?Z7>R;95UG}i_C6|e%U=l1(yYZ#=OBu);u0Ch&lnjjEE-yfcv1S*0kQtoe>$jf)a zJ%T}z(!Ue_ocQ|(7kt*?_b@O6Q9_L&V-#2c1#Nq0ho`?5R^V$$k}zzo3!uI)UAecX zXK66^WY{bYXzNynh*JptGYz5|bAW;YM0rm7()INA>UkgSz|}syLA>V8cGbcHcvZ5a z?M;%Q1YLuQO40RXe{23zB_9Hjm5VKTG75qgC}b={acpN1LZwsZCV@fBRv|$koE0nC zn1mg^Ml69Ad~1*X#(xrnS@60^*NGYd4C+8TzDi}e@0#}|ESUQ=;%DZNEdd_xQ<^n}+kRe7c&N%$EAn&$WMujy^o- z%x}oPhh(;1};({q$vbtKU0rlxa1VuAQ(6Z{Q!CfJa~!1p3etjAgzFk^@F5+`m; z>8SuRvlVF(h+@Us0%;iW{Kf#J`{(C|D36LG)s3*FN&_9konv;NHjX&?b=PL2^wJjpfAzs^r*E(TA+7o;6`#X%> z)xhNV=`p~kRBEGawE&Pu3fk!cFDM^-UjrCoC5;=9ObQ~BX-hE1!xPEb|Aj!PGnS~I zLHrpEQT2bpOZ)KnIdxe%IZ=0kaHaEvis6S~zy(tU5a1)g$V9BHtN_IXAN}JI>h$ZRdkn!)>R_r* zsLqG!9D{T^^9*T!eCURL6n5h&i-Yi`DL4_K@+yyA&BCuKDLGmt&yXASm*vxVVC1@t z0dc_3*Y^l!c<@@l89aPT7J0|ujdRj!;!PQhHPz_{k@g@Rfi8guiUeAKy}f-P-EG_c zj0lm(V^_8F^nVOwYXLo?koqhJII$;R_|j!US%!y)rSw9$u3bxtZ)MkYo23~g8!#QZ zdiFK%G`t;D92=k?p$idBBDuO^HXj(UMvXz>PRMuvvH%6Mwzd|(HYN%Le944Yg@3Ce zjl(`RE>4d#v}K%FlQ}$~!8L}m*4U1fL-tWH2`M1$PmVpM6t#79cGS)uZ#1zIiH*Lr z_Y7fD8qYx2NKH+xsPJ%iFSltlbb(nm{4w1vZFA7-6ciNDbcdV)Q;Um{A3r+7wm>hS zGZDUY9!)%vgdpz<3l1&Cl*mc++9-6Bd#%{y{-x5)w5)AkNLqEF6FV2Gk>AHn#7>!(|+;?=M5} zQDiG6%YuET42}{yUlCaa=q+C#hhb55Z^k#Y|w?VdsuF)qB>v0e2(!cW(9wpv=t}|{T-$*8kd}0I^$#Le=rAy;p z2df}S=#;-~Zf#xemSK&#V}JpM;h3o#bX`aiyP-U^xOg=z287PrRohCS;RO^i@gqrZ zo3?)Bftb9ta~G6$VwS^$|Jw&JBJ@YF+~nrA0?!JFsN$~c`e46NQu^TM=NB3(4;$-_ zU!iAVS%YlFG5qzdC7vUo$_g=K$Qub*sSp-2-4r?q+$V%MwBPN`&6bh^04gA+XxUXc zL3+NSRZo_UNLNlJ_755^l|0S<2W3*I(*hAOF(3ls0I^&^-xm{8R#v8w0<+V2k=3V~ zni{G3tgNj0SPxt1%~Bc;0KFhg;B9g>=JDFoSJ3*(aJye1&_kAW#EL}$*8B1WSg%~r zfkEK7Dk!Kn(Q7S0$t0BL&&h!Wg+d>!T%a{&EAT7C=1U=m`mTlEhfg^G90QX@k(6RM zsN_JG9YFyEEw|sj1JbTBh@RjU6CHh??8fTHb_9aoLi^D#C1eUA_p7t07J^~&tQoYq73ZexoQ5<&wT?y*#vLFsf!E)O%k-)K3 zP&6Ue_aUo_K0>DH)p*4?n!Qdm8syLyRRw>QYUTjcad_53X(>feNN!`zyAyo}_8_H~ zmAhJ7xhH)5L12Mf*nu%dE~Oohi3#}64<3Q8$f@6RzF1^uhAJ6E^9baIN}v-nUw3)?Zk$;uHo#D8%}1u_Oy zz9Z0Fkmy&hUVR7)dwN3|EC)~rX;7$me3|QpA(8&q_%$|m`O+nLCDk;A+rg5HB$oTq zd%;vkd@vVyuz@ExdO?FNhadurQy&BoR0ufLzK$3wOcvB|VBSaKR!2)tj&>(j6*<|? zt=&8`^d&5;7n&G2aI_WtK$JgDDEs;U-c~d_luhp>Hxn+ zGA}(i|Ia_yL>X9F?Q&%i^eF7>y6!lkYgqrlzyf9TYic2l^Dl&mdO5(yzVu(ew9V#E zPh^9u0f1u9@isKR3m5VYJs<4X&R%|9r%0U)D<(wwgXJIr&g~TS#@-(WN7lg> zv}rUN&B@70J4P)PDdAmxoHB~>UtfQr6^;z2BVJ7fOkp5JrP^(5ZZbft0Ne#a#{2N_ zm&b?ufGB&W{B8pF-6_#Jv;F|q9Pi-hNW`s>X>^{nz>R|c@ji&i&$%yh_`!(^ahRFS z2QvyTb-v-$hW-trjnu`{4UnH9tz%=JAWngLu)VRt`Be)*HuM+osrm~?w;4`zTm=AZ z^1$7l08v3Rz1uo75BGQ+%Kqw z7l-k;exQxVTk~Bm;6KG!5u7_W1SDfXf#3s73o$ubk?40F_p~G|EDY8_O;5kS)j2v^ zXhB6gzFP3HN&1Z zb{rRZR1uPtl$4SZuu56gG8cOB%{o8?D68YJ3(&kPT&0bX{)}Kz~-?Uba?L zYcoqro8Rd#)YhQjLhE*C69U*&Xf9#h1q;F@dc%_}0@9R{TyW!XF=+7i&d#M_e60p3 zfR)amKG1Lp)f^KXsi*Ke3fs}L^G>fXh z78&7RSNnQHsA5H&4x8Cex^pL{ea?%`uIc8-A6ru$`|@CDw6NYrfeH<@2}^%9A?3$Z zLsB_uJeSF?++h_3Z!wLGbNO2 zyibm!`Iky_w9CE$V1o=g#*ZGajg9*G?p(aEs8O_}H>ilhdZYhmZ1lWQXs!9oSKw0h z_K2_paW*Yd%X>+4HS&-+y>kQSeliJc+GJO%r&ib@|$ZY6dTwtNgRF-*ILbfb5o+D!`p&^3?y z4(Y)JD|b$p^;)yTGhtJLgCE4e<2{?XdXe(~{O4T3XB18WXGXTg99 zd(sb&22%A|5|YshlxQFX3fuuS9b+73>*5Ak3KT`emwClCl9ee47&n@2}S>*{<4(Za&R`ghhv zkmz5(eu2@nL2hz()@i2Os{kni0rM;s3yBuGauY9{crcMaTCDgl4$x zNX5kiM|=BpmGtkdT<g&Wc0bue0k=5df2i@^2y3<$P)qx+6WBVc<@3+{tCRGfXtZWVEDQ+NFe*zb<9?Oh z0}ud{%yZDfKsf9K`Uq;MwTABEw?=*N*4h1(2S1r?kuxLb!gVO?tuCL*1~9@}to1(Q z{9x||uai8~Yj7A}qNe5phYyrXt&eKuS$BuN70Plb1q24dk}Lq|RV&+#mG^&rd>M`H zHYYFzyp;8sv$0XOv`mA;*`JrkX8X9)Iec!&-G(Zig{*+TmmZcdyOsnB# zprtiEuw?xZab7ESfMiw@_oZaqUGC^G3tAxnqkiSAqoX5OkoG-lT?ehzZZ|(rf zgjj`H0OSAWc*BQ;k`}X#n~>n=S5*w(-rfcSCA$6{`k}SfO5L|5ebCCa5-`4Sj6J!d zZv!?j&6)^Af#YKl#=sb1$Ht7=h+CNy#jmm`WvO(NVIPr1EBZGVoNm!9d{kb1|0`+= zrh;Ode9G1gAD)MqSZzrm%RUZp9((v~9X07yx%W3WrvvE$vx%k=3`Fpa0hseaq9tGU zVxvFeQL}KNmea+Gr2`ED_V?5`{eW~G;H*4gn5P;`OuMxj%Gs5RN z9p`<_i8TiTF`J9$99&>n^ZE>7=WoD+>2!4|U%$fsh);;a4&94<0{r2Co}@NdwZM=t zR^zQ#GcF{QNr9hg2nITMkB6NloG6~v7^Jj9Kz?JvaSs!)lCWGgJ z6n%3|D52V?@3D!mDv4I(^9!OR`7ulqUfR^_1kdh%Ga&)UKeQL9QF!=D8E&C`f zWEPrrvVg1zmmi6>)KswGc%G!{bHm^t z9}9LQ-MTUh0g;1z`#ep7XcLARYiDPN22%-uA0|eDoB|bfqdXI<2#$NA1Z(M|;aOm` zp2dkDK_`G>1bA!&a5cmvv>O}TlB~wxBpyUY01H_1n@Q^}z3S+9?@-{ofqDZ&Yz;U~ ziaE#eZKhk$T{2ZO1)P>tJO&M&8W^&ymwhTKDqu8v_6)`j4*|Pz_~>t_2Y`JX+)}4L zkQWDKUjljkfMbEN3(O_`pORsO1^a?x$rNDlz7uUIkigHbD&ijIH|l%uP5KEt&K(kG zg~p*Fb?6}IdVsgtcTvlnQ_?*EAbe4{%B>KajBl4-auy6Ixx ztUho*z9{bJS*~dTHaY&oU0_T1`GEqU9F>t(|IvphiFouC(qgBK0U#3aGcz+_FS&?@ zZ(A7G@yhw}kSu@BJX;}s;`c*8O{bz5l?a9_f^$;$9(?D96b*f4MLCiN&Sz=#27Ke> zrR>W;U}htLl_3sH`vCrjD`+#hvwwW7*_nj6$M>hj>e5QP_{rwEn>Q2L$v8nwpv6hVmaG zyJ6Nl2sdHA&S(u|?1cMrHhc*M(D1@2ZmeP#Jba{RA_<)xuv%}->j|b@5b1@4hQ^E3 zb?iytUhoAahhLX&QRS1oXpuEmvJ*;(orOIwV{G#xa4=^#pynK;{ z>!x0{wwU(}_!po$&|+RfBQ)k6t{n2)S@&sgZvGokR!~r-*HQThFVvfV{Lz35zZQe3 zFC$;XecPOoEaRovdfqV1q$weTA5}G)>nMwVlD@5jf39ZAk@)fTYm zQf4xey}j}DOQB`Itk`1P=&0z88;$s6@k&2dYYn~YL;PPtA$@`IE&zHF(i>{X27H@B zz-9!#Pz6JeHLgciulfA)(o#8S7L?+9@$di;b5(9zreL)c2OS6U4!kVGS*}sTz-yoq zp;`fvN-KLuJ$--#CLSKG;#Y8D;P7hu`4b%AZ)$2zz%I|i$jEr*iuu1PgjVtM3&82( mb50=;6$bxh6}0Gt#Ff3DukMuX7JQ=*@mT8Fqa2L>oBsjjkGin{ literal 8885 zcmeHsXH-*Lw{}!)pkf39f}jM5p@k|)5fl@e^xgtSdhfk?)KEnX9cj`_XrU9DCrIyI zx&i4Bx*+w!d(Qjgjyvuh_nhzh@qJ^QJ+juEnDfJ^A}aqZF>vV5yD;!Iq3kkxet0BC{d<3c-RTk;IlQ&G@BoSvSZsm{*MKp?PF z;o^I4vXfcNNj&?i1NF&QWfTg<%*?#Cwbjtj5EvM!q@)Cc!P3*ydwO~#B_(mV;+UA2 zgM)*`#l_Rpp{c0}HTsJK#=sv%(pX7AIhutY003W*gFn;o7+XymkHO}Ep3*cwrmdZ| z-9MG|N=t*?j?6>KRSubV`{gRV&*9fCe(HX{-M@N$ITdL%sUdvmTX)hIuTQ9wZ2*jO z8f(x_wAQLhMU8DiRe>f zJ25eP$ISBGT!QEp;%@88^7gX4VJ zv{vv&RwFAJ*YW}kb0zA2+Wl2@^dr;N8s`mn8I{-Lswo<-@yH`)v9^jx*5@7#+q(O(E^?cIJaUBtcA812LFwb7HOg-Zy{@KC_Ps_UoO!@QxlkP_iN`0I5 zt;Ubi6-CJ)t+K&o(eRDJnThP5iwVo)I&|1}*IwRHb-;bNz0Xb!NtuDp-LvzyPq9G_ z>eqHF{^tzgMIZ%rm11Do$4ju9w*F8r4Wt!rVp>59w-wn`P~LL$%C;SZE8VlVOW=MM zza2lz`a;874)cy*jW69yZK87L!}rK-7LL1uPdvSc##3|`MW+=lA~#DlS3F~S?2rtU zDOjt4t!(lSyp06?{?_FI+e@KJH%5m6**xV?1WQ6g$;Luc0VdtsN#3l$K8b?QIXgex zWGmFzZO*|&)QMhijz1+Td|VrHia|pXqjY{^TX=O$d%S51pAF2)DfpHJ5~R}2&%(P? zqFm_?8-YG*?e6rcEYNKg=l=Bbw+OLC3EjlH{PlXY3@EL;=|?Md#}KKM6CI1oaR!!T zvz%h)BrAKbB)QDye-plUR8Sbvbz*Z!U$M|gH+NV3ZmKbnaWU^|^9lnE>6W?SP#N(! ze62uSjWdO#(2g_w}><7#w;_SL)ZPwn)rJziEu~DIA z*01_(JWusPP2@bM+UP2a?a7c61jtU6t4Q9&dW2gCzI5`CIeDR@<>4{JssOh@r_ zRe1Hblk=gM$f!ykRwv#_Z;#FlpUk5^-Gy=GukywmE*J2v^W$c{!e3X*hH{aHoHM%+ z6(eZHN##ZjLj_goc?tnb)TnTWzbVa!`K;qKO~Y6bhlVHI8Bf@r7&l70FWSNpT(u8} zn9x@7z6pd?q#)=gtu^a3a)P}-z6u?^?_IFd1VocF|M+N21&Pdw1@~ykJ8F2i7?o~# zEh$vfj(Lm}c{zSUr!bceCt)=kzQ%;dvK&JoIpq{tHf-~gMJ3-R3&MO zpZK*%alehAImDVgTgmM*CuYg!W#f2Xq4))(>kTglqW6rE{A>9hz{Ip!DSf+wiew30 z12(md=L{XubeM__t?Tzdg@?Sil=D*;6|ts!dY-?OkH1Q>;L&-tBuq%|Ew)=KZgShE z_gCCj!08VzSs#O-W<(#M3gMY+CI|Es_9G0MCjy9`0k&9?{>|1N9nItu&M8#`HuW5Y zO;ci7aj0ILOs7uM^3fvmE21`6ZE~)IpPjVE3FYpC=Cn~w{!v=99p?;M5I!K&rpK8} z$D?w}f>C^Js!C??uVK8+P_>_nRi4%N|?9SRGS%&2j;j9ifwdKL)) z&po`S8Idw{r4CYelCD!F;h+*wKQ5O4DE5~Fa;O4@t$Ua zS{d}fBHN(E>>GH5Rg`X*ou%Ly0jQu6>E~N@qI1gFGi8WheEX`%bLfPeT$zi?b5M5g z*Ga44VcSfFjP9C5d0~Wo3*ktI?^Av_*HCaG*{C6Sze^{_x0911Vc^c|;U2r8rlaDL z4@IAxC6(0+MD+!uCDn;1l}PDo&yi%%_m{-5t@KS)Foo6ON}PrhW8%-s{^gyymAEOD z2_ee*I!XEP>2S1N%Ggew$W;7~L2Jd#Q=59;pVFCC=yrG5{_s6*6aQ>2kmu~INaGj@ z#8k?6XkBtN+e|oAX)eF!vYBFivn@28{1*ftMKypmgMcnLHW_)|XA zwjA+;FL)ro&O^T=RC~A~v&zK(vI4!QuPBq54?yF&MGLaSe_I!kU`l@b&arN^#}Bq zx7U@@PtrrrkC=bydhHMH^EaXZyvpg(X+>Y&7fmVIFsgj;i%0h4QiT8XC9bbh|9t+x z^uW=5Z}-X>2=(Qcg8KS;f>uX2Jx^)5jrChyF1iK;>&NX*dPHL1QoL^|pjMM*M*oc| zgiwUtpfjIb3M3t7U6N)-r+t(nXV3U3rTT}xROp$02B`dL8~U%Ci<0c_g(jgkJJ7)qTAz-n4E{>jQ@Wk&cpjpG||+{apd8S!%2S266QzP`ui+@Ci)-g z|L2$DAMF1@NBqS~b(!h!m9D>6U;k!V{^zm3)6D*7ul?0lwedGw-~Xi|L7DHs!6dAH zmUd$7Tr-X=e>%hui0(SddUIB^JL|**VtKCFpLJ$b1c7)5IJe5gj;4mN8w_LnNLJZr zoJeXRSNLG$@flk|&7%+u&MQTFpX-mKfsU5qV$Ye;a^yuwIK0N=P1rg8WGb0XTt@yJ z4%}=Qr#iDuq@sbX;h#7Ps6RW06I@g3g-7A%YEM(e+V!)y-N+T7K^|r$pRFQze&=6D zg~LwTitvPU7ux&(`%v)Tdxr9}!wtFnHXX0m1B#kdpik6MrzZ=tY__MT6-1fJybKk2 z&HaT(8OEo|VY`uT>K>x5K=fP$Kg}&_`8IN8Hnds^;e!;rF3M;(=7qc_G*~KFN-j|K zn?huCq#BKzs-DQdkd%QcQ`jqfZQ*;J@aQu;?LMhdlLc8S?Y1J?z&W4eak@D#aeaEr zmeBq!l5st$Tb5lr(-=Of0yQ0$hT~MAY(j&;*x#hsal4Ev0_cOLsp>)*D^Arxwj|Ei zUl|43JWof&xh`>I%I}oNPvi zN>3k5gYLq}_95l6!^-x#Aim=wH4K!|QaG%BiNu;}b>+Lh1!ksti(#SK>;3aRhSa-X z^eq>hR!n?^BR`2(j$i8oOofC>w4@kgBv#A7ZM0#xZ+@?Xl?}{42ZQSlOI!_2N6$DP z(&$If2b(=G;XINXLn`%`XVf@ZiNcHNs02cVSlCDM?|`k=X7ILj)4M`L!Hci=?s%H= zA73(=sePfAj&P}&z6O{Ir8kt5OpMe9&wXTB;|JDQzB2%{^bQ|XbAHQQs#3*-$X43D z(*nr!Sz>mvZ0#pG-%Eh<7w!dY&h}Vx=M0D4CjZXczxLbpQb$w>#f7OCT>PfbC@>zo zMVbu)IGcAC0Q;tp_t{HGq{5QY6|bV0;$t<(^Sw`nY)N&afN>t4eOh~&nH|zWS2nal zMx9A5E^(*r#9m+Hz_s5^S>TsEcrt^!?{6cm%YqmubwHGqS))2)r55kZ0PC4UV6fw< zxzOtH0vId`0`@Ic06$%*9S?BJpAeeemY!sNEbdU?F?sfyRYHQQmg!?+O@kTjFv{N4LycC6do;ER$1f&u~S%|SuF`FkP# zl{$~c3pHE|@{>=x$8LY`+RO=p5A0;8@L7t*cT1EC1=IU$g6Y5Lx5W-)=dU|4`e!Kw z00KY5^DCw9Meu}8bu_&n*2uZCVMZ%G+vDb28 zO~n+nWBPNxK7OP6xyNjJKr+yg z8CIJ@QPd!J8vC9O*@E^rY<}l&cqR?bhqWCr=rJLYg-uDUy_jQdUnlxX{z?%;j< zGHpNc&q>vRneo;i(|?6em|%mSl%0sT0-+l$P(nyhwT zv~7r$mWJ(~t#@p4L?To^@`#m|Lp9XinjeU~;)VBsgi_FOfP!r+En}qxIP#%=5wk^8 zA@3RU^^mnyA<*v2u6IMb(|Df;Jd~eSP3v@x@)JBIpSr1}-GV~srmRUEjqOelJC|pMG~u2uQLCP&wo*I&WA3Q*<|3rsmDPSKHhkFi zu%M!XnpBLF5!=^AD_X{mAlcjr=);xV$K>NgA62Ou-r$B60u|zyMR`yKM+k!TOVn)W z48#<}r8az=K6-;@&l#oOp*QYlfrYa^5pP!s{SYl-R5~jynq(tg6`@AZz+pVD4#9(wM2(R6dC%1A45{Hyr?;QUHa?3RyaZKJ2H4EYB*>QS$c*ZZL{P zXfziF>^UvZxDAZDe*9={6h-O(u_W-rel`8YJQH_KNRug?9quAFDMY13E#`Y$*nvkT zYeva&5u6rZW?LpJRuB-*9Thh6OiNV8dq}Sx{$cIJD7YlfXD78f$n8d^n?anNM1OqR zkL~i6x020TQ65a;GrPSSjN*`cw-?PRkK`uTkCCePKBeGDGh@!FcXC_$K11=l*Rs3s zm|~b83%O3i{&mgu=lq<$TT<+r)u4^puwQKiQh|NHFZ=w1<3{z63c1iBHE#mQ^V+E?SeFZ zm?m%xokUgEW*ZZoHEC%cfX?XdaYE5sHQ^SdYTD>ixa3KN2iGv6T%E%uJn}Pzzv3=Z zm3IATNlQ_V?+Ho3tHAZd>=K1hNS&alCIQoQtfHPtGi0Hvc!Ey~N=y0xSO|%*1_n>w zSG}VbRY4706Ma#eymd6G>8wv11CH1=1q$U;AVV}l7xUd;)m2};lWD6x$~#9U2D=gK zGwz705WbkmDK9XuI18(Vs-s5^Vq6+Kj3dC)*yAcRMP%xOCBEwr4h6Y)3xvFebdtTD zsG%L1aA}$}JVd}CPHrM`?(Dx)O6Ni-k_#Nq?fTingE<5UMC$DU+f5h+#cPyaA!dFQ z!W6Q3=pO*|s~tN1z6~i%c$nTc4T_o=MtfGh`_YV6ter}bKt>5M#hBd>!w8fJyxpEm zD>Hx&WLd7TeOzY1Z!~h&S-;7sj{zoe={d;e&5h%G!DLO`+QQ!;7Ewyf_G&Z|#r_^7^ z2aqyoXkPh)CrwnQeIoUAek4&^dLBNj=lQn2v@uigS`LMf5>E=D)@zk+D5gRa=M{oS zLwXXs{K8e3>i241uP3xZra})DS3U}_g)w&5#21{Ty~qP?OS5%QPMrzs#I4Whv4;SV+^x(JKZyd!%c+8cssStg;#Q?FU!j8!xYL|d8?G; ziwJa);sU};iC>qbYKpJpR=ws3u3y}*t47JjC`exE#w@p44TR_N?G8jciYiFJ2*ENZJR2271me|0VWY%6iu9GP0UX72 z<#JX@cp8`SFCAE0{c}|(57>h2T8}x>VuhY3k>zH?N0P3>a4JQI1$RQerqG12+Jf*V zmP6jj*IgF@S!q*BqBbqPx3o<(iW-GHx}pBq4fSnJgm4^E%ZlUfn-=eQg~30^4m H^38t$It^N& -- GitLab