From af8231ace69119133c4c8a906e98946ec5aa79c8 Mon Sep 17 00:00:00 2001 From: "Marc T. Henry de Frahan" Date: Fri, 6 Sep 2024 13:43:49 -0600 Subject: [PATCH] Fix input name for offset vector specification in the inputs (#1137) --- amr-wind/utilities/sampling/PlaneSampler.H | 2 +- amr-wind/utilities/sampling/PlaneSampler.cpp | 14 +++++++---- docs/sphinx/user/inputs_Sampling.rst | 24 ++++++++++++------- docs/sphinx/user/planesampler.png | Bin 0 -> 36119 bytes unit_tests/utilities/test_sampling.cpp | 2 +- 5 files changed, 28 insertions(+), 14 deletions(-) create mode 100644 docs/sphinx/user/planesampler.png diff --git a/amr-wind/utilities/sampling/PlaneSampler.H b/amr-wind/utilities/sampling/PlaneSampler.H index 468f284789..e3c3542345 100644 --- a/amr-wind/utilities/sampling/PlaneSampler.H +++ b/amr-wind/utilities/sampling/PlaneSampler.H @@ -72,7 +72,7 @@ private: amrex::Vector m_axis1; amrex::Vector m_axis2; amrex::Vector m_origin; - amrex::Vector m_normal{0.0, 0.0, 0.0}; + amrex::Vector m_offset_vector{0.0, 0.0, 0.0}; amrex::Vector m_poffsets; amrex::Vector m_npts_dir; diff --git a/amr-wind/utilities/sampling/PlaneSampler.cpp b/amr-wind/utilities/sampling/PlaneSampler.cpp index d097f2cd7b..67322b734c 100644 --- a/amr-wind/utilities/sampling/PlaneSampler.cpp +++ b/amr-wind/utilities/sampling/PlaneSampler.cpp @@ -26,9 +26,15 @@ void PlaneSampler::initialize(const std::string& key) pp.queryarr("offsets", m_poffsets); int noffsets = static_cast(m_poffsets.size()); if (noffsets > 0) { - pp.getarr("normal", m_normal); + if (pp.contains("normal")) { + amrex::Abort( + "PlaneSampler: option normal is deprecated and renamed to " + "offset_vector"); + } + + pp.getarr("offset_vector", m_offset_vector); AMREX_ALWAYS_ASSERT( - static_cast(m_normal.size()) == AMREX_SPACEDIM); + static_cast(m_offset_vector.size()) == AMREX_SPACEDIM); } else { m_poffsets.push_back(0.0); } @@ -62,7 +68,7 @@ void PlaneSampler::sampling_locations(SampleLocType& locs) const for (int i = 0; i < m_npts_dir[0]; ++i) { for (int d = 0; d < AMREX_SPACEDIM; ++d) { locs[idx][d] = m_origin[d] + dx[d] * i + dy[d] * j + - m_poffsets[k] * m_normal[d]; + m_poffsets[k] * m_offset_vector[d]; } ++idx; } @@ -80,7 +86,7 @@ void PlaneSampler::define_netcdf_metadata(const ncutils::NCGroup& grp) const grp.put_attr("origin", m_origin); grp.put_attr("axis1", m_axis1); grp.put_attr("axis2", m_axis2); - grp.put_attr("axis3", m_normal); + grp.put_attr("offset_vector", m_offset_vector); grp.put_attr("offsets", m_poffsets); } diff --git a/docs/sphinx/user/inputs_Sampling.rst b/docs/sphinx/user/inputs_Sampling.rst index cde86490e9..a564e738be 100644 --- a/docs/sphinx/user/inputs_Sampling.rst +++ b/docs/sphinx/user/inputs_Sampling.rst @@ -124,18 +124,26 @@ The ``PlaneSampler`` samples the flow-field on two-dimensional planes defined by two axes: ``axis1`` and ``axis2`` with the bottom corner located at ``origin`` and is divided into equally spaced nodes defined by the two entries in ``num_points`` vector. Multiple planes parallel to the reference planes can be -sampled by specifying the ``normal`` vector along which the the planes are +sampled by specifying the ``offset_vector`` vector along which the the planes are offset for as many planes as there are entries in the ``offset`` array. Example:: - sampling.plane1.type = PlaneSampler - sampling.plane1.axis1 = 0.0 1.0 0.0 - sampling.plane1.axis2 = 0.0 0.0 1.0 - sampling.plane1.origin = 0.0 0.0 0.0 - sampling.plane1.num_points = 10 10 - sampling.plane1.normal = 1.0 0.0 0.0 - sampling.plane1.offsets = -10.0 0.0 10.0 + sampling.plane1.type = PlaneSampler + sampling.plane1.axis1 = 1.0 0.0 0.0 + sampling.plane1.axis2 = 0.0 0.0 1.0 + sampling.plane1.origin = 0.0 0.0 0.0 + sampling.plane1.num_points = 10 10 + sampling.plane1.offset_vector = 1.0 0.0 0.0 + sampling.plane1.offsets = 0.0 2.0 3.0 + +Illustration of this example: + +.. figure:: planesampler.png + :alt: PlaneSampler + :width: 800 + + Example of sampling on planes. Sampling at arbitrary locations ```````````````````````````````` diff --git a/docs/sphinx/user/planesampler.png b/docs/sphinx/user/planesampler.png new file mode 100644 index 0000000000000000000000000000000000000000..d5d7e1d2c3ae7a23b67f570d37236302cae1d0ba GIT binary patch literal 36119 zcmZ_02{@GN8$Uj+TIsYYYp0STRQA12N)l2bS*DYaFhh33jOvt-gf!VJg=|yyeRO0u z#KbfhV;N+dVJu_J_J3aKob&zte*dnkt1ISxp69;r&;7YS_x--_gk3f@+`dI@3k(L^ zZgj!m3JkWD34^W2ZdwO^(mrT$6#R3*;ldRY7%cDz3>NYi47Ln@6!HTG^Em;7P27aR zv|hnrq8>@*<~lIgny|~4uALJI1g)*D6%`c}3WZLm8yFaXOAa197#|=1+i$=5`T1E` zSU5R3J$dpZD=SM`S(!IHy!ib2oRiaXLnA9Wg|4Wy@ZbTvtbFF8F+)dp-rIZjx;3|} zYbiU2KRat>WwmM3CTC}74Gj%A9M0{hE@fmbJPu!sj$6tv;*C%6=jYi~c&4N4?2Vgq z9zLw3*UOE~9AeizDrB6=PO2S=t{j?Q;$}4exRLm+ibVIw>><+T^S@3yz9l^@q@r+? zrv6!ECi!tOwXuJ;e{ylmEqJFdKNm)kspd{jg5`BZQFL@#EF+LUm}r6du&tM`%Ebn!#fY+ckKA>Gx3LG zbiKMq;l3Lm(yAyk^K9>=mh-4`%f~gcj=22Bp}#UZ2FB+*`suOxq_=P1?%%(E-@bh! z^@r1e>(|^eHZ*`e*&ulV{AaVr1#2%DZ0+@xzlrtL1>lE5-bN<>+i>@}u}ahy{zqy5#ox+Uct)8lc!?f>vx?vE?RVhTjD zg-g_zD#RYYV+Hy5vUB8&FlBwob%QtSM{7Qbri>X|Xl<1EdvHladn|Z*0>azRKLa`{ z2_Fn~`xcsO?Z(iLU@(ue8`Gv?2nTAz5P`u0t{L#@DI%mD>@~30mU0I%NNR(4I2Q)HEF0TMReU)w)!~qF z@&-C5uTcbc+P&ympucV*>jJInRw2)oGPW588~A|;AAg4R{A5xPFDZhH*4@dri%p$9 znF@A$mZM6EJmq?wto37XAs=-R-=8~~4lOj@cU)#EGNaH>r_#ySj263=l8gH)&s+z) ztS%#(F(FBxDJN@2n<*8`fS>KKAdlP;64I{p^);orM@z#(z9L3N)5Lby9d9oY@mcb; zW1ZI)6=5Va8cMGGHeGzKflw!q6ZTm$Usw#FOYRCtQ-T({fh?PBqa8*ImQRODlMY8R zzAN1+G0-#P9@K}y&X}ydi%CgjEPnD957pw9`nC>5;BK}Upp{pDpF7p+qM8`0&E4&w zZ%@NZ9e;5&0(k1LZ{HQ#>HbU6lFCf#?mGVu^3jlq7YIfvx}AfDzLkyC8_zh<4u6Jm z-v@5B#zW5Dbcv?4;8FDn;YG^$4g5cIS2X@xq0<|pk}pOKvO?=OcFP`~X@>UH^qdNG zbwk&_A_gw3S87pC=}Om~wx{UOM_|>>)Szg z^DlY~S$#b>U!xDMYhHcwREf0LlJ{{}vwCb4UWXt{4wy6I`($sjf^6}Ov(?Me38MSI z3~Hz)J##dnA-5L?F8%4O>B*gYM#<}GQ#W$6nSsrF%vTJFa@VG8~maXW4G`w7g5J^ zZdvPcIprDVGG7O>RWz*&*Z96RW_l*#*-kb}rq}bx`4Lv?c1+>D&~bDslemID(X`C) z)`CBJvWz>WWSC~>i9LLN7yOdgWOY8Ym2LRAHPd&F+owUccax4+ChP|$S#wQ0$e z@!z!G)A;!S17Lw-owt0Di&yjb+S7FiewqQ_X+-pS)A?i-VNusVe7sE7!w1eMfaDB(_g(QR-h= zd%6Y8=$X>+cUzSORkI)O*jq{ITrb)7iGFNlqc0xrS5u0j7cN8#vO1=7?ud>!yt++8 zRP>bwF0Hum*4f6SJ{`HL>FeMV7V}nTF&Xd|NmC5BvjEQn58U^#TabD;X!z~CG*4?^ z>Ab(t0JX60Y-bwYha#Sbm>JKQ+-`G@pMM1g!FlVa-*5XCoa1Zfj%MCbLNblAgo};1 z3EpdA9!))6AC}q<=f8wUOfE%ZhSZjnl4js-T%x?^9@nQAm86aXBJueBh>2Fko)We0 zbKO_Z@tZD;%2XL+DG9e}pXpJu%-Wdq$_tsH74`w8lZLudu$fJ16Vr5&l9ws6`d4Vk zZN1Z-7LU5;|Gr@iH#?#~i46Q|aZq39PEpI^7aw)IEx0SVobi`_YA{;o^%5WF!1vkV zsv#hjljsGHbkXq}+rFyHIP*)cB-EiOrq189s$(15zsZ!TWfukf=ow<=&+ZGsJW z9yu9k8`$eZUSyYntNRlMFClF&e&Hsgf~(tbAw_aG~ZfO!mt}C zO-B-&^*IesxqCWZSDpGcp5|-Q6=fONQ++%_^q{_I)k1E2M z@w~)=KeI&c&k-mGw>DChqGl*vm-Gm7@5`!H*=<*f?_Z~d3K5;ET6pfLs2`t*`xNvo zNjBmk_wek2w$L-ude2=pxM5`WMRlh$hR^ev<)6{S~6 zLT;Lry+}pXoS1s7vHjFnKO(e-bnuDU#OET?7b|q;SZcY`(yC0Xo_odAOR}z|&9q9l z+apm3m65q2i zRK&5|QNM{x^y0cMmPx2HaZdd1k3n63>QtRoH|i#SflDb{io$LTPpIBM7A<&N7#*Gr z2CJ!K(@NvHjUv*mI?k+6^%6{_c-3dF(6}ZeG4w&D%hhL(c2;%mqG?N|*zUY@;%>il z*!hB2$RkI`30}PCXUppdlr%NOV_ufk3001%a%CCMsiQpJZ=ZF^?+z{BJA?egh-3ao zpV(r|$};sXDw=j__YL|KP_%JF59~_dNHw|x6LMAt@dDqvF*GWD!i;^~_7W`fcrR5KkYy(+$d%s`d z@rK(LhEqd^E8_8nLUVNqs#KD%4!1N9m6jY4S`siu*o7%6g$Gw4o^2FCx^2AB!Jj7(q`1ph(y3VPNN%<@#8`;Wj9BC+#U*}a&(K~DPo~qcLGE7w~p4XmK zkru3%pw3ULAQbslo|&~#StadW$l12UVEgijLG+uXlG*RFK%rF$=!Womuyml&9R&6f zalzBqQNH$BDBQ#hyNJUj$Y9q=;&h|G+id0yCuSFQ;xk#%FO1!RD%qOT ztMLsiACTcrk-0m&XSaCN(cWcWBCgT}D+_~}rjsis<-KCZxwPc&%vTw4mR>3@-o!pN zJrl;s>G=YlWHGAFFiPQ`xruit+#s8*(>Iv0kMrh`VE(YcpWyB~d9AW&O~?TQ{??A~ z3Tg~UPsuWI_7#ZFcjjEsNO#zQCo#(p+8h)M087ekUAGM0tjvF@MH>H$)?y11Z z9{{>4zjqur>@+`WydX?8XbjWKx0RiO~Ij!%{s8!0+vgoO+Dsw)IN~3kkO6j5~J?>a`Xysg81w%o1D>nV@9y7o8(Cl)z z&xDkyhBr^Rdz3ygnYU_-N)WT_Vc$AMQ-ePa5c3Ph;X(QSjkV`+)wY!kSq8d{ZrW>i ze3B<3EsT#^f>&3}oJb8WYpBKWJp&uDcW9;Wp8ViN5oz3UZF!;Gcau#~R1-9xd}UN{vrATd<6K)tO9O^mXV_vk zJ1!+6a0$GiH6`LD5N>>H{=(lN?WKKFZRBD8p*v4|SSYieTKf?^GEGx?`ZUqjgZI5w zx)2e526dyML#hvt3Fj{xs}%;;(WIt7XL8@*#>{jB2F$(4j};5Wcd@-IEl3&Ft7tHTtUG@R>@m3Cvc(M3`vm~3ZxKa>&+x8}Pg z(nZ(t^0G{B$>5ca;cf?EcJaS1I>vw$&R%k#^!!X+O`t#08Oj&o@9A zJQ;0izv#6f#5mHBYcSi}2w!v@{5_hgx7Zdfh={z8VcIz~gkq56&y||1Eh|wQ#s&h3 zn+IloN?rSVn5#}wdZRP(_9e?wdus9ak#m%JW8+UJXZ4qdN(pB&Y_6f*Uu}CWcqI!u zr@cD`&eGK!dWJ~GbZnKP{GVK2wuq=;t}1FwA2X9Ec}P$j8^rBLwNTILp;+pv@6A3g z)LP|A1{Kd7X3+Q+HEUFm1+J_l*3P2lArH)g0lU7j58XOQ{FI}+AFjt|s-)f?BgnEz z-Q9Y6(V~ede*pU&QgE{aaryjJwq@P{5@c5W7wgLJ`s8Zne@c0zl!sjw)JxZM~91z2)!gdSKD;K8CUCzZmW5B+$*R$9&cfP zs7g%FPsg}=c+j;W{|K{Wqm3^2V*&Y4I}=ChpWl5@|Mc;KT?-M(zWwp;580VRrbK00 z0g*#|Di}G8v%-#5<4-cpQkMz>4iS``hfx^49y2{EXAr9<(61aKlTws$2g>^o3#Yil zzdm}QE48=owM~2oN_?jme)d%@Lm?>wS;BBF>iz%G`_svd8CDZw7Sg{n*Zp-$S6N`R zlb&&g6C1vx(mRhkm=_y6(`BNe?!+8lV{C#+KO1jJLfdyQV|qj3Pr1f8|6(`K1dNx) z+wtOPugF_-jeW_gf0U;nUl_sX&bBA!@6{o64-4Z1mtLNbQh3McuBG!;YZ=>|{mK-b z#?b-DD_&*W9!O^0h(?Lo%No}Otbu8{otPN zkEfXJ2p1)%aqBau@dt03;W1N-A`6WM{gfreo$NIS;}N|A!i0kB!^J;_sE*35>6iGX z0h>Py*^JKG9mFuh_3dnpLTd6d-sI0Ng=W2RDOpl+G>i3KP;CdT2rfINuq(Gz;AVYb zdoyuXq&ZdgER7MDt1ZFY4Ab~9Zp3G0A?5T5f{(hR$#yA-t49QHxvura0YRaSse0cl z4(gA+waHiJ+7>SFne{6jJSQ!~WXUOw!R#9JR^$xuo*M9PInC~qlnIB&#PcuPSTIhV#<4o?#GJ)S?UmcK4qS9d^#UbHR1;!5@3a*xbpCO2I#hij;NuWN9JS6}nr&6B3zT`Zm%@O1uI4ViUa_kk^P?KknC^RCrCpEO#au)2qm!*J zhXluGMbso@nd`pH4ko1QSWPYxgFZT*alci;q{Kxo;y>6UzK;(|3kPIy@j9CkS}wFy8UAcK6A6;wSyt7Yu{w5s z3cq|4|5xR)Sj;uP%y5B%-*3<<(bY&AzJEyWmB+1uwf*$?D&+_taR(<{xsJcIEVF%X z_M888m28}q_H`7q(mJY<+N(AThAj8pK!=Cm-{eVVyEW6`FVS(Y`X?i@O-#kuIKy)& z9l=#F8vturtyg>HuF<6XyjCm^V?G%(8+}kO` z8d2MRWB$)Z)Ohp(+WeHLcX6~jvpsxQ@deP!pFSXdP|yJW*6-E&5E;)+?kXQp((es4 z8Z?L&IDOCMqI4D8aK>BK5PauSNeiXL;oK2*X$2tk|U;G|Jz9lO33z%80 zXZA0@wh2}gT#v^1Hc>-E1bmfQq-TskB^s~UKO>78xTeo0{N0;gEdS;km_L=%ZGz2rv)2q|{&=puz)pXrYxsm)K;%jXaT9wDF=G)dCqWg#EBOFZ zy~l>zA+?1^^VRJ|+%qK&Pxq=L@(Wa+^0klVamSd$}Fey^Sc#-eXVhScV3e#^=~2hcj|G&%C&e`kvxk zQPU-*(6RBj(s0%D@0A;jz}!*ohdcf-m^^w0RjZj{Zd7gbC~Ooo9oi}AIE48mMwg|i zTTZf;$Gj17&qudjsV!iy)nc#rl}EVe7}sZnb7P}hjRUdG8M63j-QpX`2bt}r&hsT` zFiJb{+~{w9&V2E>RAT2=A4z#8=9<`&KA99Q^nhD@lE3yK~9xB2UNwaGRcO zEQM<4k}@{{fWg<)NBgmD1)({u}!g9h|A{(So4z`*rUZB;Fb^7x^y z>lJgiKO(jd9F4{>bF+dph!Nfjqw|R|$^x5Xk=!c!@`fb4o$MfASH=lfFi^RFX_trf zx8Si@)uO6~aBiw8X;1Z%_e+*Oij`ce%RO`u4GKiJ%R6FQF+_KZT4lPBlDO9 zwqLDk@0oMc_3Ie>H@++v$>Smh^&L3*1R?6hxeH%FT5wD8j9 zhaxf*lhY{Zuf19==%`@uM|tjMGf379$B^Iu#|wC3N}lMA<@t76pPacurF;9Kc@dJE z=dtdTvw&xnx)CZ9I2y;N7LV$DJCtAL=qK*Jy!)rGHenHueWx>O>rL8ZM}}i zY3b4UgA2u0kF|ooA+nb{F18n7&+&(DMq|+AVo5V|Fr#`?X_dmpU#H5AWun;nz5F$P zK}:JUN?iYXTtv+x%@(bgTszQu(ZvUxV{oISMp&Z z-C#+A2|gH?@;n9IPPcN0YdO)B7!{giH+<5QjD?Nw~hy5w)IR8onH71qp54B3O-lN>FO|A^8z?VU6#uZ=lDwU#rq|H46-7)xd9wPZsLr} z%Rts?lP~%NeO=;rUPYDeO}aK^Er9Tx9%ce1(O^%I`4 z7NND@gw|aan%mlsW)2#b{3e(SMQFOr^M6p{U0AJt(%0$)?V`^D{CakI56Z$1!{yQq zD_Q+RlxlxLa<-v9)aLZzEB)vK2FAHx&?Aql}pb1@|ElpB?(Cggmv6No;4- zcou7Q9x{%OkFy9RH-#Yb&w0G8icoDYn2B+g9U~p@b&G98ItdbFsMbdW>nRi3fc86D zl`#iRxE2k$uG4dMLgI(TS=0}uT)Ge=4*lUff6zFJdk5WvVL!qbEJ@5&C%T$u#yw|$ zNR`BmUNf9qJSBoM9d*xN?rQNce8x@w%uK=ed%JkZ>}Xf-^zPY;e@c$4kP}0j;tXeI zRGqsGF?RTNg`MoF#49d?^m@!%z#U1<0j_m{oWT%v6M3m%#~nT18C6% zwOA>3e~|T`Yf31=B^$JJSn z5a+W``Zd&I4+`9WNUJjNo+=4-H;e7k->J7R)H;?|{J-ZcsgJK)C4SP3ZkKkt-aYLS zIO&iA3n^z;9^Qc)o!mZHJ|C^mz8U*I<@BWJPQ4`>rAoSbzL7O58Ig5JFcmfJPbwZ2 z5Uep@1eSEuNZTBN#~UMO%)V-Wzgx%dC|z%AscpMa)gW2zbU5x%IWbFdejJ*^SfjZ; zj~XJwf@&wHK{k#3h=__kVVfU0IxNzCOevT*|Ino0T4T93(h##ak|F{Rq=kx9g~)~2R*MIM7qT-0lio-ce#l3xk376aVwi!C6X01(&QDg9m9Ck20qDHBl$Fn#f>Cd@_dXhY|F~WB5vZ$EGUh#0PD+&u;j1-1ZE8 zs{ek~VU?1WEwSphQVuhvwa$XSd6Vzu5$#V&Y6jD}fAMDtT{p%_e=Vjfc_!`*aLvvR zz0 zq)vGpQ-ldmrx^_i@|T}9)u%?8X{3{dd%E;pF`V8#7=fA9vjkqGUPI$K0%6@yTilqr z3d&EWS$#k2c+vT~pnx}B5+--k6MNBp8LU^n%VfTAZ`Dwm>^(4Ssk)@ri>sQUAD+V$ zP9M_ipI7&CoWZsoVlu>0~o*n}_Y3yY59J#n?uX)4GJfA31ZgdmALa9{-e4 zaZfOgl;`sBm!3&j204ceC)mBWJX~>2W?_T4??jH(Cy_?Z))302a6c*>^OfhqZ!bpG ztososw1`=9e~&j)n{TgH`A)8NA*@Z+iZQm3Vz*>HFS%kd4em64Of;hDe>oGXe3X~A z*v_xkVhR>h4~v}CRaWa+)*OU6SsdjR+-++T2KPjr!Te@o{1|kR*m-{UeAJCj0fT|A zFuhW%9ce#&HFAI6-r|Km9sZ{$AKnVSRwr)DRbe*@1V9V+td4UvGMra_v$}4r0YKQsSAX~MMq@u3i zr(<1_n{nIKjMOG&nky6K*>&`y3JX&S4QDEAVc)wN5)-}Ar3BK0S{H1ttssJmU`i`V z!oSZpe<1iWd%f4dzO~*aw#m9G=4l&o6V(I~rY>c@`qX`FL0#vpN#s?HqC9t9*l80J z(>ct{xb}!TUO*Pg&~2j1^)?Hn-7j(RyfI_HQcHYG7$G5B=hiNv=bugR%5+`+pRkGP zIgxa|<7CQQ2Z|_eMvBAFbGiBy2{wkKZp`wxb=CP!9`LNP4$(QHN;;+Qu~_}ekuRc$ zW~KQ_q4{huHak@@QlMm}mXw%3LBaPPWGb=7r|yz#jR<*yVuz-s*bK9KQ4)&1A9IXs zAIY|-j7cfh^y?+aayV@|yHXYSQ-6iwD`$MK;#0ISCEiiRCVP{veKM;^8N+iXvYpEe zg|Y5|5AH-{7g2+6Vz8)E7>5&djEW0nf$MZY%d!-!mV8h@qf5airDve71&0VuM*Z*r zsipX}YI&5t&ecx%$&`BPH6cE$Uw{0v0s3IV^dx!c(qJn|FpZ;Ewdwkn-l_B5-p|(- zP@04mS!E4dXZFF-vz`G1g89PxaGGp-f1d8!6dm5((R1|AwKA@VxnN#X0b8iMTw+Jx z7$cp7(qqp@Tt$q>?_h`BAw}SM^u&OWr3vv3_n(_%hVv>bOt@Z3UQwb0v;6`>MXG5t z7)7$u+NA~lsC{wgG6c6H8U-cBtUHWK!dj_63yE=l)JpE%=m)18?>J2JB9K8h*9hkp zDha5F-`Ijb$EB8;URNTUs#6p=3*c1o;B(HFVOPls0#-|^)#N>x(B(i%ALcuV1?A1W zSha_aUm}~+n@$;0;&c!_8LTdw4U2y!!{;u8F-S1(!TbL zs~EOp%i7;l$M+n_j=*<{O!f91+L&2VZ{q3v&DPJnue0T?TQ!%?tn)KOsThS4t@~SX zVBrab1_5a&tZ|lfu0rEtYGVJ*p1p(U2ExaCq}vraXckgUhL0CG8#hGxO{K2QSn#>m zJ!3&%^Laz!-Wf|llYey3O)5?5WQ{}QGcVsti{mt&%Y;&M`$E)3vl2JC0Uox+tc2qe zJ@9V#`aj%TxAex&R*-+Vl$gJZ zCrA^_Ql>cNG=5z0={Lf_y5Ls6rG@=DZ96O~8ubMh)iGr`^^s~-!F1tZCZcmAM z|2Pg@IiYVCAM7m4UYL9{x^f-o(ly?OZ`Pt3N`#^Gg(bsymmKm@73~lC`A@ZS^h^XI z#-u@>l^|9yJ~5^39h16GW#TWf{q3%;&Ar!buzl!sDZ|oy-+cTS zW-*oD-r__hi-s-k`IhwHbY@h2?VdwI@5Pr>k4Rq!*!GGI&pcC?OSIFEt6&RFN|gtJ zTW>cd^CZ?hGdOeLTjN4)06&gdpzY#JaV`_uc3C`kyt^@c?a( zYwq5+l@dt@{N7qUO|qO4BwuCQ8S8tBM^qli!9J$g$&%{A^F*l@Qx1?^YL1^7$O&G8of58 zj|Qv+q6JWmPz1%XR9=E6wlrUgY5fyq>vJmmRjO zeItd~YxyPALdHwXXKll&-hob}D63~{nI3^F-^)Z9i<~Uz+vxnwx(IPP!YJre`l5eJ z#gTNXT2?qpx_pYNYr*DnmVKLfq^i13+I>OHD+3Q9RpoTGlk2jNlq);G8q& zIvYun=(uiQe0@2=Q@!U_fIJb;qs9ks zHav5H#`+lQTUl-H(~~*pLXTcIcFWSfA}Ad7@==S~`Aa3e>ajxm;qHYenQ*h6g1ftt zE&WJf&UZ_`ZQK21wxhe#-+(o* zi#^RC2VF*|HmN@dIDs_b&n@3gi;Tuw?I)@;Hs4U&e~nxcIudHy;JU^|X4tZL_~0 znHO*bJyXWlKt02;i&G*o(TJ)GdLNbUMM-HOiTktU8br7DvPA`gjpAOM%zOIg4keBH zcDujdD(L81@VGi4@DnT{TdTjkH;3 z$4=}LQR_$louqK~+MYx2@P*1+ zM9-r}a&U2bcdvJk+Y7fWq-ZbweiCa!({>0x?uuBHkpPY~bS1g#`l(E;e|F2nb-5_^ z+kAK_Z@S⁢-f4q`|vY5zDKsmE~uv1Cvqdhd11Nl^e=flr;jibyK%ph%YZ7tlm_d zjVrEBv2;~Y6mRwL`hhHpdcH)c{*rNENB$3bslaC_*GF^U>*((R#i8kTyQ@Y=pXvIY zuV}m73;%;vu5L>iAcDsO-k#)SozS>5q|Vuy8ne6cW=X^soeMYWf*tLI{oi%}@sX}r z+~|UQ#w~Dtg2W_XIE@(n;?c|t%T9TMz>X=T&v&WTY6I-v^>@E*mOYUh&^vl@#u#V} z;cQjv8KR^DGLbC7mabNC3tqkjFd3XMdS~!2LHP-4qU1O2C44LyJF|QJ-%Uk}{OsF6 zb!D=BYOF@n2w;Tygmah7FLqa9g>hM8b20rB#@>W2ew+%$*NEXbGWN%dcQqeOKu^Q^ z?wXL(Qq{`467n@KyFD2MpGDofF;+k>zV7Sw^I<3A@44e8v4}cJo4J>7ICgl~T?I-V zpbTVKPT|QayqfA}-i?sP@QG@j#DfCDRBN4oDBI0}hNoL8-8WV0v_K0B0*lXn2>$*z zMYBL@@R-s4UE;_aDfYGH5*`|{RIWBv;#=6vu^*&XgW0&Sh0 zGFu$^aM58=rc(Xlm4SCR>OPjRQqv19`V#wa8#iU!m5qMC;a47<@5LW1t}pKWQ4+~b zJ~K{yDeQuFxj)#<`m6pmMN!zQjq~@EvG>I*2cDCFa;!Zn=T{>_GT<9EO5t?cMkD8K z4Edc6cd?XHN_}>cce61A59~o4j&1tuG=!MjJy?9I+Kxm<9GE-DR&k7t%|pE zakJwwTIsjj+X*VOEaU1!DkL=~6*99kW_P64oB9-M>^flwscF%N-mc%0H!$+A($Zqq zhBR$nh*Z6(ZNP848NhMC*(}Q3F3j`rW4Q7xFFM0L!|uFzclyAMx-~+b3GNZ&A}H2} zvX0|D$s2g!IaVKOyiHn#K-*}9_x1;29Xjn*ze^CNrZeG8+VvB?iu12axC(Es2hl%g zk`tb}oDg^R8($);zij(>^rVhI4I~pFh)o6{T~8+@PpdeS`fhNwPDE0x<_8`r&sLd; zRQCPI`0g>!By$JTh4j6lD^kxV?+67Mdsz+KGj^E|HBk{sPek>|tLx}+{f7@J+3mS` zmIgQFMovFQH)3+flnB-9GY)T9Yh+Vb)u?;mK+gFRm7}69n=zO7He@)@SxFE4$uhwg zG+m)g(=@}+u-{CNO!InkFH3w3Vl_{-><2pk{^4XF_;Sr!&fPuoWtEtc0D;rV=BVbK zi}Ohw50mD+O^mHCZ*d8{0ytdN%1T&kxbT>n>l`br?M$L~e+MZEB;n};P^euN_U>Br zx5uEa7Rr9YgIn$R>w3x5WYuXGnKJapkGynM+tikZ7Ty_#5)B3BK`rfSiWdW!eQy0~ z-}Xz;uh_C$tN+5jd{;@(J@C5cD&XO>cNdRPd=sKDnGJU$!`Yi|5R4*Tzb*STDA+7h zfU6$;EA7CJwV(sMnF{i%=A1tuRxO!hJDym2!Zcq6Z)Ow${~vPS!vZ1;rUNlwDI$lu|N3K z0UWf-%|#gyd$^G^0zWBO~I zaiK+b`kFPi>wi*;23VE`ExY0Vjtf23ER+BI?zBBn42=hpgX?788W~^lPI(|k{CRYU zj83P>Oy1L)T$}V@)lL9v+6+w`oW1Hw%v#SMIa{AY=>nP_iT2pG^6n9s>F>$G@G+?! zQx}c4tdA)U2d&d-m?(5I7o+gDtCt->se$bH1MUUY?G-x)b61x)ixtTMe$D`g)r!>m z{rUVl&T+`iV0*wk5F%jrPOPRtwnm6V4A}~}!yGJAX{7fNa{jv|iI9UT3lE=5GtK;D1@<|0!I`RFCbYf+ z*m?iX9wEriAHdE7zdrlUJR@dzrZVVC7QpQ1 zyqK65q_T!}*Zx5vo#gHj$s*iS3kg!T^XMgxti0Y$?{GeiOk zB*YIu8I0#__1Mz%PwY|-OMeGl_0K~>jMD@dCt;N}pzWeTG9-g!!2U{Gb;zR-?3{u8 z76u6nVcYIT%v*EWZ{KZz`}c!DLa~1fP$K+Oqu&3Oa-HXsFJ5_F4+UbTdDY(Ac63Mj zEX;Nz6kd?Mz|;Ta1hmdAVDJhl$sB73xsti~rTl;3w=&E67b~8tu%GEaYII+t)6x^U zl3}bH2gLyR8c6`T(1`)`Ix=_Nor?eW*T1l=!>*q)_5jxIShW`N zGbl|UG4U$@=Fq|a=FtCxFMhH~b1#T(mad83N8ra5hFVRzpHWCj{F&;j&ifmcdtk>9 z5TIbHDHYNJF$_3Wqd4T|cLj{>zI6}mWMqy}fvLY3a3w|hL6Sj$g0Q*+5ZAq>WHNRj zMILwrXp0r<2KGE(p$rH(5Di#kQSm69e*;}ukO#5a|B*Z^`^VC>m;b>iF>}{iBS?;| z@)(HO3P^u}+a5sdXC45rL&OL2xw05=${8(AxUoU(%XdF?HmIimM9WIHuMXxk9%X=J zfzn?G@Z_qyX6!*W{?~|Spy9BCvK}gA;0(U+|IOS%P{_2*&O`#xuQ1_%;n%$izjqdW zR{s-2Jv0(3xtLp;`>(i>0U8KfVVagvwVURK8Y=uIFme1*NPR0A_7g)mKm} zAN~^55Z`M^Uk|yuR0SaN9}89rZ>Kd&e6>ig7Hw!qUULk9uK)rcM2(O~AaDajLBRms z-~UN^zxuXcr1yVe^f?p-s18Fc14{1s^(Pt0EB+SwD`Jp#Tj>B*(^q2##T$ykAlvxgpeg+- z{AOvUn?PNxp1;2WhZPe1RjLXPk3gcRE9EXI=4VcECR|N5pv;v8H#V&QtB8ID0T=h~ z^gd9*SAwt~G^0=<3vd7&x^A`DtP;R4H~+VSuFAy~$NgV;lV9cM7@smMsC6E{7B|wK ztZvpl5E8M%6u`x+X9p#C1Uh0@ibltJDA|9V9q1#q|745=$i5XFqYeRRrOSfYT2oeyY37KRc zfH{rL!GskLtQMZ1k`m0SApr>C6Ve^P!~aXz|9lv18)NB-^NJE0IJK(VnDswF0ujYc zDZAEBe(D6ME%}`Hi#2Be0sWssp%D_> zKwx5UzyawY@Okhws}$3jucl=AXNt8~e84tzYBt{h)s(V40Y_(MlgH>b?PBTqr$34lESrpUmiT zlVUCnNAaHSnkC~;n+}5!^p6x!^?(A%O{;(Q;T~x?D$CSicI}6gtK9*RKG|!-U6D90dsuNkWnw$nlT%ah0~w$&|EEvrifV*AL#FTfsk8z5^3Vo+ zEp2fArb^&9NItHrq$k`WU$8I{!ppy{LvmVuW-!L6Z}Fp5!1j#uUQQ||aa;4GG2Ba& z>H7+0NwE`!LUne&3d_l}C`v>@`4kj=7PAHVwX^x`yfqJ9M5)&`o2l@z&-)h_No zZfu$tZ4>yQ^yk~e^d&ZL<9IX=+`4=0J2I14J7P%BmFT$Zr?#K*(n_3i`}v%#LX40U zU3F-|uN2c)*T*XrEv;kdvwhm&2-n$a-XA~op#l!|y*i&1{jL!Z3Bp1$8GorCUW&U@ zF*&=d3b7z5_ZU-^;9KExM?oubJI?ph&d(g8IKonD$tm?WCBft7b-cetA2n3g~?k=GnszWW<*Uf!ssWM~SMn`~1HOw&}UyT1L z@1?^ok=pMqcmAk1&kHbN54AlE4*WYf0M}A{NHO^B;s75$GN>2X zJm|(=Tf*nzxLS!VaOEBfyIoW1pAwjB(P!Se3+|_~J)JanwPZY?|BVMh9{HEyh;1#6b4uQN;QE_`F+MKbxp3JwToPu|cdqg*flQ}Yhv z#L)OdLIKy~&5@gnN4K&AaiXQQ4E^a%-pUvm3Mc)ddiL2NomB<(wt(*0YWhvl_eQ)Z zw*OlHv_>#nmRy~x=kBW3I5Aiv>KEI(W$p4Fn@hs+muxHqmsy9zb1cr~zMyZmBPcrqUTG3nz|fDV2dQ zxH*_t)f~6wf$@JvnqcKGQ+KtJoT{52bc6Q;CK@&5@Eo_Lv^ecc?WeJ+qxyHiy9%k~ zl-uKRjg@mxoZ~Ca1B39~b8ft~EG|#!F#r8znri6ueRtil#-BnLw1OhI%q_1Ac7DsE z%Y*kK&IMe`t48^paLrBB59(}=W6*z*b zW8ay16{=D`We*cmXzo_jYXP(A4hDR4$1CnH`gwunW$%ysn&$hktCv7L zQl7XV@vj?ut80QW0+EV96fXtK`HW%G~u^&WNBebCzwj<94mhYzK1 zEee)DHE+EvBF;uR3QVh*n95{B9WtFrowZ#qR`Q@y1v*Tst>^E*P_8uzBx&%wGd}ex z$)mYP>vc>UQku?8$6RdyH?^*l(Vce6-0`@yi)K;MYRr;me)2}>h?27l_Ok2}^nDJW z?;M8lsH+3ULe%-CyzkGV?vyOaJYyf|_4Q}bwYhCl>fOQ^lO$a?no zA?MCs%i{Ce58~Z*MZc`TWVLbblA}NG9-B#2tr(2j&9<4#HeCQWwup5ehLamLdS6}f zCEh=TX*!Y;9;~Y!Ra!-Dnpc1i^y4wIY$fymOFq|G=JYz9>b-VFsxoR9TQpAvh15Bo z(B&K>~aPts9Ec)U$G~6Ju*Ks;}*QV?NI*$JfK0e6k5bL zAjHe@-lY_wadWrmst_vrH^^>ozY-}dz@^x>)+t8@PEF~s9 z3%aHKvm2&=zLcP{^1i`V-+$gW==pvo7HVmRY(?tpR=lI6JPmt&Z zyzbqr(RL`i~SvJBvFeN>__^o5|`=!l$QZR<6h0yABY9M^6$Ctsp2` zaI5NnRluLB-R_ild-bd#i5c^Q0v~+f%IgC4tv_EEFv+PxWGMFqY$db~H40^;o-azy z_E0~kt*&-2epM6tc6Z4k;IaF36|1Jp-8kI7e_jblo}uHY5%wRBe^TqN>LKj)T`cVW zvTSFl)1zdIua{fB~yqr zUA*2Ug{i3z{iFiU6Ir-)*5j_D;Y@h@O7i0yF`uD(moMMxv_N+qz>KH7tsHgNQnnfT zHUV83=s|oNQm_2kWW75j7Q8z)sqYk}MA@Rgi%->+p1J7{4JHf(3vLd(vwE;q5I-Fh zYycx>NYk&@_RMsil3{kqw(fdH+Ic=Y-~c@@ow%4lkc~AS@)eU&mpsH2akjagM?#RO zK3|BK5;VEFp_Cb4jJ}WJLosyvDGq@a7dn+IJB?dE+|e^1+}8;2ugf>zc_xEHFy1vk z{~0`CeXnbxUs4~-m%(KS?!}n|_)5*H1Y(m~k(NrZpR5i%*ZPyopBEO;K|>XuWtYMZ zOwRvSI*=KSF)`zSuj5$#?tx=Qcv4(OMM+m`KZa7BqI+J4a+A0}V~lM>RGxQy6SPu$ z;57W5UlP4o(QikAL0J+zeaxRU<)w+B!6OrL|72T_&CUx`4r5H`yx8?Gx`^Yso(l;z zFVnSA-Hb`H!CT!(lo=6TA_>SLqnn)8D~=bAP%?V{jMxFv5e-zv?LBArhD zXb0~KJAIv`bEPOlqx3W4>b2-D`Bw>!J>VhxTL?<)*@Z)<%y!>Jx9=-kZ*5R#97KDm zpw?JuNBjEU?-lT+(aA(1_>z0LK(uakbdx(I>00l;YVC2NQ z9Qb}eqB4IPb@l*5kGTXMN*PKDcAaD#7dh|PkijnNuyydMdNVuoMV$ZHPii)EXlbJB zuUz{91YWQmpEONHesKoC15IKBxG|CEQ&y#V(NtSBThdX@V5#j;RO&hYt;=YCPWfFf zfsSZ^$Du2(Yi{VL=a}cmI_iZNaRL{`JjuvRBD_TI6JcxL|5x63Mm4!@?Z%4Rh8-&? zh^Q#2h|~}Z7B)h#frt=0(mMeHi5*b^5fP;;YE-0XC{m-M(gakLj-dpQ8VHaO(!Tj7 zLD>7-=iKw%JH{P%jJyBrVR-9$*P7*do;gO4Lq6i8QLN*Ad* zXV6F;8F&2&_WcoBa!BASX2lG(0S;gT3+s{mKiU#-++|hD4Vg_ z>u{@D|1UdZ1iC=%>5uM)A*;-qw&g=1se7aoamPx(rY`xMZsY0`e1MWxx8j9We-MG= zWJzmyUDDOy@}a7fOA7g!UjC|S${MCn5YY0km^AQ8DfA#LV6**8-gwXMvK6( zP0uA5E%wfD98;8pgO=*zun6;QvGf_(1L9SR*!b8Ge5nu(p+K#}5(KLd0aVkoS1dm| zL$j0=1VLo@-N!kH?@UCX{qJVbw#CUc`10v}=Wvi-gDOpgQ^PR4t#_Ws@y{%Hx7pX> zH`;1qk!wI_)YZ7%1DC`G8R&-)PXskeb3q62qOE!_I?t5Yt;YyFVY2M90bQrU*mA(=$0jB6c4Ly5)v4Lk@Dz;U(VR|Ag$kT`wMualtS`q<9>W#};owc=Y zm6Cx0uIO1pjSbO)P=+XckngDqk8u+@cRKtS0Y3uGhamZSDJJtb!cTP!iJSpun1?#8 zwxR?*F9?}O2wEr>iKM>A@NNPMJbKwC9kB!Ze%hb+{`wY?#sA@T7=*6>p#K? zn=ktos*vqh8!5|G(o7D;rrl zhqo6U{SZ%^I`jfW1K$KY+vEca^8d%s`Y*`$%kFnPJo4yNRu-xp54^yPUJshlyJw{t!SgUC!Ba1F!m>F*YR3rEg|dU6r7j zbUrB}=ET!I)KC%eb5gF+iXjS)9e*}9T!MXtyQ6>gDC zI1niW2|y9i_&RlI=HBgVz|#C*;5ZT5I$!kB*?)q)dpM3d(*0LvT!Ke?Pq*FB@iO7a z^0x+KUKyp8qb|WL5!e%ME{3UX9zFo9+Ep4tqK;`25CRsV4G{wN zFC;TUL4p(r;T5{!D~z#CGbRF|^Na=X^{DQ{<1fHJKv|+_NTLvsEU)deog{&afmlP@ zxB!`nh!i6U(a?-%X(aw^wb+SYk?3}SXEGDfsd>OYicyp?RRy52WBRO)-^45^3`D8` zgo8p%zKAwP-ZL@!1`W~@hY942%ilK zJcLpkzf?ptLC!x42hnDK72k=KgAo5Bb5bPRFsN*`f*Po|2tF!eTY5xRZ5>!OB=8aO z`lGTnx&dz%{Tg6B?>=eHvbDo zQ$)y$RY1vbvGI(oI3pAb3J}VDyzoxsUdGKxkreA^gaq|Ty}^Yn1ChOn$mTnqm53s4 z%+8J`3+oKLpRb-`yC2D1#$`I=m4oBeR)j}L&Io8<_h?M?@a7LKG(@(; z3R#cf@R+q?2mlfGj81HX7(vuD64wX?38R{MiM;rKqB#+P?`2sq0iv$twTXfP#FLFC z{C^h}|J#`D?@$P$utP5P@1PJQWxFtQR9gJMiG{5GHwOAY))oJgeeODADUz{D0(XfL zS0u%OEE$Yy)cZpcDV12EAc}59p6Vav#7Pp-*I#V(KPrkvQ|_`y(5GTxF4I&`R3uHVb$csp7pNt3SAZ;6pS^!9 z;rZwI0UUEy>Pizt=KkLl-Tx>g{-5arfd>=m0y`O5HHGNOUa8;+MZzB>$1ru7?1t{f zJI;z|)zI)=wsM+8M+fr|30x|dG)xsZa#e<86f?GEMH^>_5pJb~jhgVI;R?4;lGSqZ zvTQvJ)r3R*iEr1kSDw@v`*Pu$H;l;|TugMCM^WiQvhd+G3I7UaCTOMOxd36W&_g?a zkmXrmPLwwoM{Cd3PrzUx?lDN(qbg0&kDtF%aBzI(_f>tomAgZPzCoWy}P?mk^{A#<ZgqjENE&i4x z6@F*~9MOUkuy#yB3wkkR)~tN>CfWLMR&s$UQQl{)e&zddx0AT?VZD0x9R>SVrPVO{ zxu$cborD<*%FQ_&vSmM)(=^sZvMfrZs*ABEV`;(O1))aV{BN@0&8uhH?!>IimdRa*ZWd{DS{())KWHm^Jm;bGTlg^wU z4%i)Hz$#m>$DP6jC-Nq<{wi$W9`MK=`TYvu{Sx5HEJ4%uS&rz(r@u=R4iKtBD>xm& zC4`|}@EJw9jjqe1fUK?tYp&-=OOdI-s`u-=7>S4)p5n1M-VJEVO5{~#?v>(-!H?@u zhUcM?SI#P>lNUPpuAC{EucJEbk~QQ3iy|z%qORl>p_(SXvRG#ThR?U7d(7k6CYoF2 z65zUXmZ@AJBVjJ zpirf()JzT#zOvwuvoH77(}^L+H~dz21})9TEWvmUGt8Sw-KGML0nw~qF*ZXo_7Jk< zeOn*2bct#kw`&La(DGBtq)tEQ7SWJmlVUQar9F@t}D9K@!m1;@5`s;(M);3;qRFXK`B6Qyo89 z{3$&nUNljGZAd>juDt8CWzgEEGNlga=2h3C7d8~=%66+tGO+by#Y6g=ecaXfcRtWM zZk;7EkxWV)zdZ8l&Wri<=Fp2|G{$^X=?DNo=d(=L-rW8|B0AptlMsUyvLEyY4(e3lb)0JrX2*6HFX`dl?(wIL2&;;O0@^3z!B8kb6)gfN1v{Wx75bBN zxm8ZDbG7jN!xcg!H=h+rl{=+PltzmF=D>!w=&Yg+<~Io(1}$SEHlRE0KQtmMt#Lf3(UHwE%-3b&(x%?@$eC0}(A@mnBL``B zl+GJ&mM$Q$D!S8TcKMOt5U%(ak_g;p_C@2W=S1T!!0JUxcZmAPnNlBiD?jCAKM>eg zDYvHIXw4y!ZU}9)n5vaS#jNuLc!Wu*}vPqyF^73|0-?*tDCB-$rg zV~LUlvI8hUKjV%+MYX);c;g*kMQ|$FuqC}$-j|e;Wk@#nP@x{URbjfr*&QpPCXl2c zu_WZsz)3g0SE^Yf)vaN}M%4MFn!F0z>Aa|;+FU(AjjO|8`mMGMuP_gPnMC*)XStv> zcYv7obOrbO4bggiJH}MC91q}qHypRvVS1`uDa4kJhNFPNLc*{GNi`9+K5h0`c|W%e>_X2cK5)ua+7+Z zRY4|hO&wBB1J8`4;mhw9Vl8(^`R!%Pk(=hmHjN0DEX~BNt%DLtaTQq*WexL?MgK?y zZNHvz^w<8FQ(nT!1(y37osm@zEQy^wCT)m!Hnk=L9ms2lN9WU5Qw!y;AtByw0D_?^ z(s{c)W@w=!tLO9bo+Y*Nq)lQqJXDaBT`$B2^7X2_>oA^lIS+lo3Hf`E-+p|pQm|w` z=lQGi6vl~&!r&hV$X{hR`AGGnZeaEL+3E7kL)Jyz7VNw>%0A*_(dIv-m7)x~2thR4 zchSjZP-X{Vx5Jn2pTvlf+!o#8QeEgm8}e8}Yb{DXvfLiNdFA`kywc`2&8knivhLdQ zqKsrE5Jrm-ka7FgFiZ zQ)C`;J4gGuvr<)GZ~<}fP$%IM1E098+v4*is0~(R?p=_ydWwp9im6MAVn}RO_Ha$c z=2<*`cJT8;IcW7XjqZ*7X$AbLqAKSXUb~^`@esLC?hgD3CN45bSyn4fupDl|Zb24Y zM}lTnl3q$?+-6~aGd8j#H?m}@S>5@D1p}+nKs9`w8pvuf!Y|^SMoK#+Sj23Y`~q1t zmz8Bx5;@uq#vr|SD3qW&bztyQV$#PS%U;&oXl7xB%X-%`3m!xIiqJ#=mTRIyVn?_m zE)=JPhINPT9^Yao?)|cpf1*iGQM!S)>CunmmD1blkIsjsaVD^vDSqeRxH$7_AwrRzLvBV@+H`?c?-D_wABH9SxDj@PzIwtxFv zvxMqeaFxPKv~(jBg{8sH0&%b|B&B1rMs%f11)uXM=`6mdOLFe*u8Pmp*Q(Tv!8vOD zkGZ_~6>i-3S_0OYeLem6@VWCilMvkJ{K5jxTuwmaW(s3NGA5@ext-9#wYeNs=~o_I zWSLdkW}hweplETnY|5hgWlnGOd1z67{ zR35V`r<%ly-Ys~dHbSZ3ES%*!Rax$4QT$8|w7DCceTFMAVzyM8&@p&qqE^AeQIp4P zhPn)vi;XERyFQHa6(nHj4kVFsg?gv2hE-OQbB zn9-wQHr++WC~xlPP;WbCY?WCgbcoOW_O8ArN9k8kq1@3i)-XyYBzEPu8od9T@5@sC z6IOX=g+4t(j1@IZfk;0gP}DHR02-!bDKMpAHc#|EoS^wk4$g}4b!)7_*XB@={wT4x zPIY06KGo~JxYa!?MA$8cE?RT0^YmT%YBH#TT+Ds{u1{Hi(#|eZAc1o9Q=RXE-GVP* zNC512vC=0#@)9uLQk)2U0_ouUI%>?KN>7VA-L|#{tqS8Eif1Q)NOZA-&%_YSRtdNT zDwlVsuk)_BRM%wGEO+u^VE^MhG14=gd+Ne@P*&>i>62C$&tgA1|C z{eRj0>|WTzYF78G%3&)h=7c&;U52G%<%p-D{dH^(O~t3qFXZS-W&!fPuT%=OU{17z zP@KBXY4*N6%KoX^w$~wNYSs5e4)~?~C+%|5H{06UV%_9P&XA}AT>@ZH0~Kx-ScqL& z)SYdJ8!&bHvM!C4i^(qcN(khc)p^v~yPNI(s*=9HAI?I)x)_kT#1G(2X8O`K+OWCZatk_h=fX6gAE zX_v<$mL?D!)YU>P!mQd3SRR{El$BkrK+fglVshnsS!7N|Fps^eH3((|WdXt;}L| zp~oSEx^d0$g69HlwHDmvUAHE8(GNzS5p+Lk z&gR!W!~D#LFd6);d2#z)EycF5KTi}Km5%05>L6nuu`K68&#u+72cK1_gB{Hwf zsZHpVfHV`K;z?YuhCn+^puEsr)agMC=Cb6O^UI7W4Eii9A5*_wbU93YZJo9?G)<{E z;(z1Y1!5~Zw6zKr65_#uJ$6Li2Mk+O^hiW#MKI<}RQ(X(CC>?}l=5_MI ztH>RWhPN!g&igQ6@6F;I8p3LWp^A;V6?r`V;6$@Tn`TuK1Iy}F^&_2E?HPRj+K-Mi zVnGWT)g^N{)eAW{|NsCSkHsn)$zTfYungYA7TPJ@fbN)Z~0e^h~0#aVOC=j zkF7__*n=11=aq1Ge~pyl?b)1LF)X<7k$ir^jC0vnWe4w3S1-7UMK4yjs8}`YTz>aG z-{s;`&v-v8}5~|iW#mbn1wyyazxb@DHMhB_^)u^?<1sS!*>Q!pZ8=8K09^y zr;U0r7PBNldXv!fE5rVED*GIB!N)k>7f%}^PPl+IveF@1F> zf{Tx2J1vUmJ+qhoGWeLgiSLcC3vs_w(Wk-268(N0x77S~i$KV!Cidc{a1mFjOR}*H z`U2fDvubX?He1@2>T`20K{uJ2yRJ60nwBe^xuAqRkWG1tPF(FFEz6>%UhBfX4-U0$ z&bEFN%nEE(S|$@VYrBQr7NOUnIul!f>?YsY*69l&VaoiF*bB#&;VbT@{j^8glrVd* z*g%sK6f`OEk9!2uehuOQpp^JYS=JR;abFfq$;|YhW+7fC3R~gyYmo}nN2c<~FS<}0 z4`>Qljkw;nLdHg;$;i+{vnQne-((%ISM5G-%~+KSQ`E!6WwtCGGAF>0t`zfiQ3$od zhd1!Hf8)Z*d8x0YY#nSYd)ay^`fyk5D}uk$8K)v4^hosO{PZ+K*%@ z(LWx~N(A+c6zcwb&L1}XN?C?A>Vg!69f@_dr!8O?!HFTUah3KxuM**K$Z!rHc|awF z63(DdY4GxHhrkF2-vL<1*wP>>H=gElTdHX=fn1pjW3a%1)VFA z4G)nF;{WZ#l`c3@x%q6;aBnbX1HUgABN*hfF;6g8k1=JhZTCH)wF`0}c78bx%VU!L zWBDT`Y&B{G>A4o;9PthThKEv4!?DUl%cb?ib2YZwG?{cU0aj!LWaMbQ~OU0)=ALcyZ=^4xb}ef##nib zI=hEWBk6ER?2?zHn5c{1M;&rU$_5Q{y`^`p)Jut47ZMx6Y7QiNVTIU4 zeg~J)9@>+tz|2sVY7h=*N6X3>e-xi@59ikgOH(>=Pph_@qtCnU>V3F#on-tVv4X9J zkLvJGbEbFwu&1_WWb3L#-C%E;+&rnI29hz#)yh9YW&ap3T%{6md@4UVQ)-5y3SL5Djo4<2jH^SG<4L72dVqKG==)92Kmn%p7 zCNtUCnjzhdBSxnW@;GKB_h(8WvD)0GrmcQMNyV6?7-tSE?@ecWEGe(fg(Eef%IW6H z#tgGqkzM6xZv(nH7o&Nl4>*b`1-u)j_{?`t!Q7T zkss>pH%;jL1 zlccfbtB7!nJm$=6z7(r~K{D5UF3c|<;g}eLNjc1T5oNW~#ma4n9Cmv!p$48&fcBN+ zkVb~Ep+W`buBjHJ6Hb_?C8bToVH?iY1>XUbnfW7oy-+*ak1XWmtJ26ayEky=M%7m` zni4Q1ZqFHwi-ppe_S}p70afKHdW8ER=VSAhAuvtvBP_#c8=buJR2z(Db6}6S>>gUa zv_7%$OqX^o0cqTWZsYcGyGJv^Q!F8%htiiH>9Qr&sd8-!aI0Y%tl6+`HM?YTyyT^Ep;R4oPBhE z-E_g*d$_%f_G0ykg*6X926-^K2e50PV5k4)mWe+xv@U?o8YG#Q4LbS&x4h__L4J|lT$ zBs|_tc25xGYtQFoLnZaq%=z14j&y9Kl?u3Ij}V&2fx>fK(e++ex5D*sBugzv zW?)mYnTditW&(U^*hln$SK|7dmU`!8ql z-j?}xt;VQ->zt*n_$hu*Z^*kF7*qAlJau#^z(v?Ps=t}jLQc~MDSM*k)ynhr z&%bYcOtd}sI>9yVPG?=1fw|#>mw#Sj+v33`m@99j`g3i~fG-il9phc4=m>`9eGs zBc%TLLi>z9c7M2+$l~{U%0r9GAvYzm#kV=Cf4zeaDQx)p^1a@{7~b0Hd<;XZOtuk=s(uojvoHSTr6*y*~G>|Syx|@O|Zb0PUY9+vto!J zK9kgzWGBps3cz`0z9=ehoonN28CIWs#V~F385mM5m44-V#ccl`?Vi?B?Xw{ym-!|i za?urEbxKLHvYtmATlDILubC}er-_(W;=1F{8h#>+Of%CD*6cm_0#;`Y9KGqdwL=c= zvfSS3KMLBx~^jy72-saL59uJP2ndj+^PE2>udkKPZ~7e0XGt@$b?V{s_m=#xBN>4OU+Nu6YGTrTGN& zjgDvOxI$hN;4CNEt6!a}0_%0J^al-iPd%ucLLk*XCkIc7tcTyrHFwpbPWLAT(`1$A&LY#-6OF^rT)J?PVG_WnsNNts$Gt(3u)VkNn6ScWtRSMN}tNS$BnqDu_5F+#$ zZu#*unvGlc^3r_o52k{SMU$MYG#@lzQZQz#1!pps1UXdRq2>x$FDiy}Ae|UQmo{# z09?MTv7P6U7bkyAfOf4$Uo876>#sPbGS+W+uqW0s+1fRB07~zxaHL$!$$|KD=rfdZ z!Z1sZG_0k_%PVUAd;s__HtHx2Oi!{ut z!wClybFeDGB?Et3qLD>n9Mgh>Io87cSZseMmCkIahlX=S3QQ8B%g{7Yl6AaKSM z5B0iQwR;?e^>s@!yY3MS#euD2def)C;YbCN>aVgG|3ox@Pw&H&Z0cv9c>+$|5(oS= zarCK}-uuZ@2XLvE9oU@u4(?2pf{Cxl=O9f<+;un|q%peDaFna3@e-qfl?BT13UtY8xO6%3Qb2>5P-7%5^n9y;e2MD-7R?M(Tn>7gN!2{=sPB|vXNb*ZI7kwH$@=++0^fVcQ7KeTmUXhaA*-)j65?TK1oZ zd-ndsgd=;MppIXwF&{xM>m8KEv+W%(1ucWKw&@sAapdBTA|^34lvnY(Rdi%jZjFtA zrH7YAr!QC;+7-8O8hl9qCIR92xw_Nuf4p6^qhJ;aMYtc`|C4_gf4Z`kc$cf*{lZBJ zhFRt@YmOkBIl-Cv=8+)6$*m<&|N2)$MP-rqfjW>Cuz6oIz7c|jfszh&efXD zT7qnu&A&1fA@f!#JHKiG|0({?cUtGHk(;=g6MI`APoLpQx!^H%InYt*tAjN3DRJ>J#O1?iUnl>lKuOhuvOZ5C2K<*~KVSpI73$m$3rp{sN8YAkRp+DDHIA z`eQR5P=7EMWctHn>DGFX`m-{n*HIff(p_N#1C-ax$~Tn98hlPl#ev)zC{+K#=nsBH zy$vj49W0QQd2J4b3~8s{^uFWlHD_~r9<3cQScg3Bq3nt!x`VAYLv5vmPul><*zOw` zM&STCxK=wl4|zT7xV5FMd}Amb1OVmjgj>Zo+>t@nT&7r66*jD4v&~FUlN2c^MWLqc zY{n&Q?*4g1gDs^Yjj@}}P_^EslEfw*hr{OvkfOU-e#bli;WBq6jPhah zB$SEUJfU_jk&LV7bL$IuMbyJ~zFZ3TG0BbV7Z**q4zs$<(l=~axuK2su6rxdDT3o9 zEcz76*G(jPs_XehRgdX_r+*p$^-+>sSU* z4?1xc2!5$FvK_UfJ22)1pPJ_S{K=>b>jR0_FGo@j##b!aD8LV>(C}PJf{ez~0BS5D z3cljQs*Y4efo)6Ki4JOH#Vw4(XRQc9`-x!q78mBP)us&mxVui9-P?B!MRf4up28HI{lb1lm-$;^EuaZstmg1xNA*xHxUl0OkP zLZ7MQ7eGo33+3TyceL;Afic~oBM5YCJDcSex{a5&4p1o4JmI;I4M#HWZ5--?B{h;| z#5>zF?3BfNf;pdBK?*b&a@1hI1$Ge%@X{$#uJ*f;6c-3J`*9>w#&&R<6?Ctfq z_5B%Htd}n}Z1R(*PDEYs$u?p=2@obtdIHOA>zyRPu}t+``tqUUb#PyR}Fqixb@ zr`w%eT6lv;nz0rJ&VzG&@$34_h};B>ibX237x6>_$GJZ-L}xKh3PO8 zP4J11WDFw*b9WeDksbzyKWo{B9(}n#c{;B||H?MF#2;SO#|G_T63ekTgBfg@HxqTC z_Sc$aUo9@1H_O=;QV(!U1$C+u1wAO#vztVrHm9i6@|RO^65Hp^3^kR8V$Ania14*K z998V=os24XMU?obejk{K+HMkPntbihE`l{P*AjnEfplOICo(#y?J1mh&@Ez^yo`{WVYVGekxZ60Kb#p)I;QSBgXoznz*Y!p?oHHLe p$N9`TZPT-^ZZ=MiC=>lX8#k+H{0.0, 0.0, 0.0}); pp.addarr("num_points", amrex::Vector{3, 3}); pp.addarr("offsets", amrex::Vector{-1.0, 1.0}); - pp.addarr("normal", amrex::Vector{1.0, 0.0, 0.0}); + pp.addarr("offset_vector", amrex::Vector{1.0, 0.0, 0.0}); } amr_wind::sampling::PlaneSampler plane(sim());