From 2d4ea278629ce7009d8d00d7c7ec37b75212e4dd Mon Sep 17 00:00:00 2001 From: wpbonelli Date: Mon, 25 Nov 2024 20:07:11 -0500 Subject: [PATCH] fixes --- .../test_mf6model[1-prtdry_drape].npy | Bin 31584 -> 31584 bytes .../test_mf6model[2-prtdry_drop].npy | Bin 35184 -> 35424 bytes .../test_mf6model[4-prtdry_stay].npy | Bin 17184 -> 16944 bytes .../test_mf6model[0-prtvor1l2r].npy | Bin 24384 -> 24384 bytes .../test_mf6model[1-prtvor1welp].npy | Bin 24384 -> 24384 bytes .../test_mf6model[2-prtvor1weli].npy | Bin 24384 -> 24384 bytes src/Solution/ParticleTracker/Method.f90 | 1 + 7 files changed, 1 insertion(+) diff --git a/autotest/__snapshots__/test_prt_dry/test_mf6model[1-prtdry_drape].npy b/autotest/__snapshots__/test_prt_dry/test_mf6model[1-prtdry_drape].npy index c7af1441f4e4c30ea3706aa2ffba4148c3cf5a78..bae56af52fcce27f891e3b32e8dee981cdd9e7f3 100644 GIT binary patch delta 238 zcmaFxjq$-Z#tjZklOM1POuoR#E8Jl1zyJaYK@LY1HWXzXo#8OKfze~K1JehX8WEV9 zgkXnt?hmV`0M#%~W>iq$EWj+l2va2hQ}rOgf$RF>0^7;4%-%p_Bqm$1+<++(ZZLye z2Qvz2oxH ZgnWz54h0e%n=e#6UGHY;ysJMBL;B_8gM~*f=K7Nb{J) zEU>vEZ36q`3;YlAzR2oc{sP+H=JGMIf diff --git a/autotest/__snapshots__/test_prt_dry/test_mf6model[2-prtdry_drop].npy b/autotest/__snapshots__/test_prt_dry/test_mf6model[2-prtdry_drop].npy index 933a6a989c75e998abf53e0ca44422ee360795ac..34e1279aa13e07f81c5ad20fa1f6a63823a0dc35 100644 GIT binary patch delta 152 zcmew`iRr-RW3mG~$K;4CkI4^M1QuIl3BbiXsyrqKunBC|sJg%m7v!ijn4D3qvH3-9 G1rGo|bSg#w diff --git a/autotest/__snapshots__/test_prt_dry/test_mf6model[4-prtdry_stay].npy b/autotest/__snapshots__/test_prt_dry/test_mf6model[4-prtdry_stay].npy index e875a76d491b4faf9847e0c74a6239422b37c787..71917592aa61f1c4d041018cfad9e782d98b4a01 100644 GIT binary patch delta 39 vcmZ3`#<-z{af1aTyRn6irjA1GWJgB3&9#iXZ8mSPJHR}7gPp?U2~H9K0`Coz delta 67 zcmdnc!nmM~af1aTyNQ90rjA1GWJgB3&9#iXZ6;r^;h21YQ($s{ox|h}HUgV9>@G0F R1qGZGCObG9On%|$0RREt7g+!R diff --git a/autotest/__snapshots__/test_prt_voronoi1/test_mf6model[0-prtvor1l2r].npy b/autotest/__snapshots__/test_prt_voronoi1/test_mf6model[0-prtvor1l2r].npy index e52948316f1b0df5420c5130a5102a1955525df6..ff1e38b14d337b44cb7610524080d0a767ab7cee 100644 GIT binary patch literal 24384 zcmb`Pe|(i?6~OP968=DEz(7E~in5_ouU{6b>0(onTFm^)K)h3E~u)m^(tAmV19M)XX+a&8y3~~E~~4qtY0v{ zS6M?{buE-S&}T` zJmJW5a_#hEG2w}n@?&#{Fh86t_2q{0!@sHL zb~(<0W+!7;vJ7wMzFd*mtL0+(tznja;;M_$cW;ORTBNyYh=CZB; z<%Sc5n^JR6H*&F^_g~I_)sMMtx2U;AMlPyVT*2IvrGIZ#bI%9r3%x^^U<2271dDML4PBr(;2;3cKaUAXEI127kbB7waI1U|v`3B|Zze)VgQgeqH zxyU;S#$QrD%-mFknmgRcMSEvLp2YDXE#pG0Qq4Wf$i-oI{W!LFAdd?zdhXdCxAZv1 zg+w{qTfEWDJIZIP?LEiHrE&jCo^PD8(fQXKdTu-dcPq@V$e;h$?Cf1p6}alc*D2>3 zxu}-CklR`2vGF;9t1g^-o_jD>z*5`4=*E3veuXl9!bUg0ds)vNVdSD(%T?UJr?Wp- zR0poQaC^^p4<5I@7)P&9;Q8}V=El-$?gbuqw>JRDrAU69#8LLP=(!_}T#BOv=MCsz zswTMG$yj(hml(O|PNQL6hB7Vw+^Of58oAiNh39eq%DOu{S8eYoBNx@WVZBf5GU3K+ z)Z7b=Tx{NA$X950Qu5U70Kdx zYVO5GE{*#TN7TQ4*xs#r?j;ep+oXRBBz|kv_Fih_()f4r5x=5pI`!Oyk&Aw{8s>Gx z|BuWk?^WA-na8bH&fAgqMk#au#oYLOHTQBO7u8b!dZA6vO&Yn_&Ie(=K=C8%g^C4% zt1f(8xWdS#ejU!^(OcZkv5cB~rICwj!7%3jj=3p4ceIg<`C`|Z%}(B28@GLO;*nx zZ{(66!?;i5FU&Wnnsz;Rg2$cZIj&raaq`9)WBwh=Z$NPT(q~= zf1gsJ=T0(m(cUENZ&1F1{SB(7RnNW7z zk7;j8&z)@K9T(tKnf9}?Erx>}%jNf2!;HnFcqZ^D|w0FQz_A8m^H0ik? zHFD8^$HKas{P<5Y{_43O^SCp;fRRg){5XvZBmDR5i6w!%AKuP4M&Kq#aDM&^+nd#M zZ!&Umyx$A^D%j7-KQ}x1Gd=faBNzL1y5GO?rGcw1yq!O8M`7|rKkTp;cs<`y*u?tXYXf5OPc_N_W-(IFfZ)i=O*QBNutM0+;;xV6$7VysPI<^SJZ8fRRg){5Xlg&(Fc~z}*jzqv>2x z0ZT0x?cH=X+xvU&-zGix4kH))_f&t}@6>Z=7`d3Q77yb1-N*JORs^oP@OGYQmzWt!}}xTC*?hRho0MD zCLx3gZaR$vAb+&!wyQI-I-MDIPbM)<=*h$&cmy zMoQ0JV&r1G7Q%TI*S>7;QayL6k&AX7l5;Y0-mFE>ZS=Tr_Td5+g`VtrH_;=0seJdU zz}*k;-(~K>T>l;qw*}@o)W3h_d0n@jyF3E72~fF+y{(YZ0ACMAF1DifvYaupI3R@-$V(FxX3R)WP;9 z^xV}ExCaKXKTH2k({mq;!0mwZNR+RB%l2mV+=q-@RBM;{q?||EqUSywf!lTg_wV1? z-gZ5AjggD~ydCCs7*A>WKIJ_bZ}5 z+&4W3l1q{Ncnmq^kT-}2d6M$S+j{O3MlQB<7o10;IFj}6z=s1@UAR9#Y2>0>6YL`q z7tUXioS^46Mc`)P+%@rFeMEA$p8F{y7u&hY-$z=j=YHDAMf)bhcWvZP19+X%s^_lv zxc&KBg!r*YDMxPK@f^QHn7dEUeagtib{;XB=U4ri8(R~&>cai`Ge$0|1^#|rsh<0^ zk&F3aCd3it4LR4B(sMr>f!opJeHpfUa*3Y%jFF3K?ZSoglf->a&)wj0mlSXT3-M!- zau@yiX~|d6&J<6t=(*2E;EpWe_=WfER86;@`<#)BYE3ZTp!S9Dj>+CZYXetZcpN?N z9?Y%xwe5_&D%b}me}?l{#2v5ael7yHUd9EOPZGCU&uunx(az`Ld=ByCJKV?h+>H^q zGp}LpF#ay}MLoC0EU>bajcaP^XRU?=B6V`*2CuJY$89jHKk&AwHXc%+5c)ro5 z=YGxO{*tdn$h=CVl%s6*+xZ7dgsqlK?F{?RBxV1(=&`_67w*sBh`^nC8n^Rn zyw5O3&uxvs-3<0p-V%E&^xSVo;O_AIcdefLEh87>XAH$0Sk#Ak#ZONcQ1T@ME)#)f3j20{h^VIYOQ|! z_Io05)rF4>KXMP|#s#gt*uRg+c_jIJstG;!*jZ=fwX5BBELw literal 24384 zcmb`Pf0Wc!6~Jdj5fx=+bzOdlBd80CFgk)jn>@`>{6S+1fuioZ?8uK5hMh$LQAbu0 zb=5TsSwq>P2#abBBTTC`i|Cx~98Q*DQSE6vL}j&QdBEt$z4zTa?7m~?d%0iu=keZr z@B8fi@$UQaeKXI@n0b9>PD$eKl6x!X)#Ym&DyJlG8C|(3SJzNEHd$G0#fsTiquvXKf4vw&Ms+aY^Yg$`;vxvV23HmX*KovI%v<& zugTRxee}ht)L5{~eaV>rPAp0If7+1-x* z^Ue#tOibW1+Rpt+Fp&^_kA&MXlDU5o?i4lmcq2Cr|MHdH4&$r!pTeE0<_<7&f85#G zxoS+eLw@?-!o5PxJs}2nZ>rm6@!zpLehXKsxhERADAqE8{X0&I? z-tAz#rr`VZUL z-Ko~M(#S>L)^R)^m2kf5(sM`0;0}a*g?>)=V{T!-THi57E{f$Y3i5Muftq{1k&9}+ z8SX(aa!lBgSevc193&=tc7ZQ#~Qg5$LDgslKkAG=U!;!qMf(F z{)qGhE|u#Rsr4NfgF6bi$WQ-`xzp;^-0?;(=Bpoxo&V0A9B$ z+?$p0jNGy_dH($k=c{3P=evuxGrzzr#+ zcj~#9d&3+t7wwz}|I&OU{;gzWo+l4vVZ$8bN=&aJAafrD&bPgrMz(luiGK6NIawGPBU`R zzx9wOF}~80C)@PgtBhO}%k)}5Q#U<&?$t&v=?UiyRLc8bcenH1Mf-QUk&E$Mo@W0J z;C*+Co_mdvi*a>uGV5C==L~x8wE;J!oEIoBPvX2$%G}gF&UY7W=j)7I${TRbK;r$G zThMc_H*#?t&4qOx#TBgUX#Vfib7#ch9)fip>YaU$^VNvuUjIb2oj+#eqJQ6p_hqD~ zyf4e?xfvrDd3!HqI}h6E=l|FB+#8HsH1AHwb{@Xb+q-;)+Rh)3!7aav*Ey9NeSFpG zxi`k(PMg8pG>@ZpJ@=-7dz&(Epq=;4u~^GqWb<7(CA99L&CH>2myGICL@ zwAVQybdB#Hl;*Bd>0?rb9$#X8_zh4SRv%ysuV-(7URn&TfVuAjAW zh52fQ>=)pD3hJMIQqR5F$i=ui1p8m&75Be-?kz?x&A;+~@BrI6waWSKqU}7_$VELT z%KOmb`vpDsRwEbV>CK58&!_PC-L2=|7H~_0je%69Bd22C%X;p-fGhWOx%vOz_r`eI3Hx8N@6UMO zUGjkQa?y6KGjfqP8T2J_Ip4_WxmhC@#d_piMdraaJ@-x{7wx=R&iUouW{;jb-^j%{ znhX1ViYJNZWWo9FqW!zT$i;C~2KN}TO#g|;(banHLL-;rIlM2L({mRE+_GQ-kgAmY zcmm^jNO(VXy`EcdMMlOma!*$LiJvSF{2L}NpRVn##Dq7&aDC(U(aigP#YCU%;cT~cqmW#Ol zu-_;D%6*l!dTv7u?wW8N+pg#4V{j+-T5q8K*|+rEMk5!;Q5EbLa2%)q%;PAr+WGFH z^Ybzz7jb!6*U7o)Xg&9CBNz4E8m`+jdhR_&F1mA{oac<-akN~|U2f!}Sn^sPzvGzO zrsu8*xc&NZ0~b=2Qa&YpflKuib9b%wbw!?BNrH#JU#8`vSm83(cPMjvRy%A@D7g3f z2aEgtBjF|{@Hmot<`qrOcNaa5R>k19!#O{Wr}PP|uhVnyi@~ji{UFWza_&7}&%Hke zw-nCXi7Vf0YteHbh{5d$&!2bbxrG?q9NgO^eG}~8E%rbRZ4z5fqHI$^)tmYeBX-7@p^8Pk&9woa36`p;XV@Cd$ykY5D6a6S6Y2BPp*Oc zgtX3im$?N!_hBO!#a#Gaa*Lk3#>hpx4)1l2MBVJxb00BsF^FHK{onyTcWuBO z5(KVPrR2w{$UuHZKWF7VZ1Q2}KaU4NE z>ztnZxRFcxh2Ljiuje)!xj1g;K)#}UA?NKKdhQcOE{YYxd11ev`=pVJ`o0M7_etNu zJRg;=alX6gcz(*rMX}vi@;Ew&xs&wVmVkRQ*QJngE2Xq!-e`h#ESBl<%&pdQpEhzS zu3*1F+%uWGR?q#UkxM*S$I^Tx>)3WZ_ZcG>#X7=!6npgCPZ_zGH@1h|#3Rmk7ah-^ zHgahkh1}74?m8nE$L$8-l7HoUiy1xlStA$Y>Y0$cT+dx^98ymWi0# zC0|I%%?Y&cE8+L8cImm##o$iwwVtGIdi31S7`Yf%TjBhQ@`s$qR;+ctyXbMW!N^6i zQkZ|y&smv&ou2!=k&C!#;re;Lp8JB4OL5gJKhrpVLeJf3OlE{gSp=Y#_ub-ugkc>bJ!uy~H7)fdOnK*(39cky?7 zQhM&D7~HWJ@Vf4OJ_nwy=YHPE#rT~bzQ-x(xnD4HF@Cp&_kp+Qxtjy-v&wizu@rD= z{79bMt>=Ex$VEF(hWB)oFXZ0!0X=t14DOgAtS@|Tg5rAQW6pOMJ&s;9a*?wS?(@?= z1-?&5dQaDLzhva1SQ75llYf87+@*T%myKN1Z-v+yz86D!uh(<82HZEfE`_Xbq?C5# z9)x=Zq_6yLK!=|D6(bkL+RNF$a(`pLp8M4p+|*#^{)FTCl*gU#E;^pSX5^B82Qjyk zxfAu=b|V-4nuq&26i=WpjiVYp_a!42#hlz1mG3=0rssY=26s2)XUbFZedl&PcU!>S zbPP9eA>&6%`4sK@j?BOCdlb~ZN6&qkJCbQg%SEwFc+Q__cD}plc>aceu(;0Aa#7#m z@IHm|2E5;=xEiD9ZZ~q#ujS$IQDpSoZyLF%@2YSgx?IoQVdSD%7o4|~edTu%+w|O5 z0&a^07`ICtOUca%^y?xxRV~C7iOc*uThINzk&Aw<8qNNN@6%B?1wHo% zMlOmK;GB@e<##Z)=(#&d@NhhbaeEwp9=4;t8x}BE{JUGv{o&E%9?)}l9ZhcXDd)TE Gn)*M7HDdz+ diff --git a/autotest/__snapshots__/test_prt_voronoi1/test_mf6model[1-prtvor1welp].npy b/autotest/__snapshots__/test_prt_voronoi1/test_mf6model[1-prtvor1welp].npy index 3a4f3890ddd06160dc10cbd4d9ab0cf8ababef66..c3a6a65497d530ff1045552acb9c263b398c8b72 100644 GIT binary patch literal 24384 zcmb`Odt6mj7RRreju}m)Wztw(M-xpllfp+jp0Y4ai^fOMLq*_*iikWk4JRh?QIat_ zmZ*WRAt@|P%j~IRW%jY62n_2}<^ywT%BjYw(Q?k(>$f4!r{*r!iGOxybJzZ|*7~ir z&N1)jnKd;dwee@NbCPm$vm5Var6y%h&1kfjla-q8 z?y_+Re>L;%)!oLxUh$t|uV&T%?oQY3(A2D)tfY)_Q?pXsIP?$3C#7eny8GE_Ntvnc zb+_KVdd0YLnH}``)23N~f28ZIRxrTQ!m~F1t-ZR!;r}`|h}Zpt!tt8J)PBw%w*&a$ zxz?Xwhuc-mJ?&TZx0jgPQpsIjOh4Zbilv{;S2%X%7;`M(HneZgkl1LCme|E#|gXa*e++w>j+hHQ-$)=AN$Pntbk@fV^IZICm3s+jzKX zHsRQ}y0i6$?SSFVNdxbC2_{@0m)j z(Nn}U=4Ic+x{i|T`L}+!>vB2HuO}eSc0Ppf8)y09Zn=m3+;%hSSl%}xm0V+XYd+#T z65J~tJ}$h!XM4E0o;tGiV^@4Uz`B})IFEv!FBfy$D!Hb9hs}S73_!`-)p<6Cwl?!uQi;`>1eyrd+>X?T7%I{+@R&q_B3_i*JRe~EOuivgpt}%<< z$o0E+BJ%_C`i=H*r+9Q_`;=Yr@c=(&Z{qw;pM*He{q3gY8nc>B%o{qH^E_NWkGd`rk1kPiO&t|%VQ$wc?01yh-yR7yR&ZKi0$T;mPm+t>3qHaGY<>W`3RgIpqo^m(23x zcQUtGPLuO{rIJf=t}A_&xtnr$ZNhzBoQIp}iL>pGcE!g7R7csBud%L@bHSD0=k!r> zsg7#f>~f5GHv zK74=17nt?_F0%sy0L!y0wgV==efBZu+I;pmB%cdx7^vq)cywjm_i|0&c=rJNIchHJ zdS?j7h4=SrIAGnBa?Sn5y`L~QW-hN4a_%)sF8Q0d?L+3B1Ap7yCH8l)l1u)6bFh|m z^+N&k^W@xXm0U9iqd!3X6kwHeuTyfVek*$JWj`wm*zcif;yA}Exn#B~p@#i@9&xUa zbB8Fol;3uHtGRw2gTGPs`8#|aU9aSlo})%obDY-~IOdp*a_&$im-3t5zMA9wcmex8 zcoOpkbAOl*Zegox*3)aq?;bh#20z^Mz&(sOJCnu!4p(x?EI)+&2GOrt$hjkwTF zDscZJ=iZ{^l38Z|I*U-!lw2~~e-q-o z61o~C=iaI0Qaa-=MVwu(IVMHU9joM0oCilB&MUypk#ol>xnwr}^Fut3mZOg5$hqT{ zT#9qt8;2ZY?sDlmv{=qfQgX>Ga>F6kb1`(aO3t0&gL`Zx{9T7QKPBfTE4ier7K@PI zr=hDC<=hk{m-=MCnTL2T>_FexE$5~xxzr~+rW|5DRYF%aa&Ab;HToS5?#IaQK{@v> zAKX?qAkP1Uo{!17X+F52YoO~ zVc$L{*cICW%5U82I>%gpUBdn2h6Q{rcwU%{gDg(UC9|T1xWD?og#8^W=cX&U)Hl}7 zuH(57DCIhyCg)C3aw)%wxwwyQUCR6dIX6SerTjjejyRuH>X>7e$+=UNT=KU(6>)A~ z%DQ?&&dv0}JvtiP3&DLw&YkArPO!qpzO5g-V%sbmhakV_gIgu%X8GU->~qSwrCdjK za&ESgOLa83C!SNzMjd@G=jJH6GzUNGg8W9Hjw0@N__*+WGS>&UFbZ6^eo5cO7s$EO zm0YUdu4h75r@`M|a_$Tt-1r|4v!ACh7Y56@cPqJM7CeUMg%g+yW8~bKN-pVn=fT6= zmyeh5nkMJYQgTVpCu{LM=+0krOrD&3kCJQV(L3<>OT>ARoSUcQ8h>{ozyHF&gI+D? z&Q@|Mzu)YDp1*;fH_N$y@^I6HbAkLV-0I=(l5_7>aw)&NH#Fqe^s8DqH($x6e)ak~ z==mS0qc7y#Kl|XuFGrjYBfqEQ-20SVGK+i&dj1r3)b;^~j|)GK=I{sW-xEr8MRAUw zhd6(PICqhA3w&_P-1lb2`~$?fkDNPK$tC?Boq;&-MVyDpx$~43GWGkdocoZHOJR+_8uO?MeWUe4hmQ+C7mECF`@!Ez=qgIi z{fm-IX5)J!zk4tjddRtpm0Y8vp78f=_&ZR}U83ZYp1b#eziyn3t`g+jzxv=NMq>`T zx-$MJ%ehOHT=FlbGxS`IxiDSMeOSpQ|JDZK?}v!<0y*~)CD)j_T(j08&MV~HWj?sa zE&}%${H>64i#^=-&9Q;S&I7yJ-%@^ucS4-q=We6t?Q-sNC6~;KTz{!whE8JszAxvl zP;$*Y>Igmm06o{qxhs`ia&vkV{QVJe{!Y$)RLLc?gbwi6zRx+M(Bb34&!ZClVEuDS zQ*@Qr9{#%T+sXgXg>r7Gl1qJKPCNM94DahMmvdJsx#Zrnk+{EVj`wxf%DJnRT=I8` z`#n49y0``2Ki?wfu2FJ{cY161>%N~P|3j1H+_g$BVai&fPqspw^W@xhKDfu5!(W$c zjwzCJ|EA+&CLuor(L&x8&Ssm0aR=w%jvtU-yZeyH&|0 zf8(v5?fu5ra_)0}xVCS&_gAd1MGhYqelBcNa>=adBh-<5A8U^3Am=`B3>x$-)fJgJ z`%p*jbBfVd4>|V*C70q_`vLUqJ`b8W50rCXRB}mIaeEu?ugv^SkaJ&Ba>*>R8uRzR zn7_$#?#oIp<@4Zs(A9T{bFQ5GijqtD9QrPF^>5tgERb`zE4h^4>F+>S_Ivgfa_$Zv z+%c66^VjHkgPi+!57+yCtR;T$-KXt@D&*IFzhU&eUC!O9O|-nvd44D2{e1e|hv)77b3EtXbH4ZaJ{-qgn-~wYs;*xL|f_R*!)(V|w#VlF@7|Co7ol zj`HL*d~4?0-gGm;_V}yV-ppU$yTdg*ZAw;7R!YX5Q?kanIt+{%nv$L!bjP#Dr(_1* ze$Re=`oy}r+!y2X)v;B8U()uqc9?*sji+t?uI*jo@V|~OY33E>JdsbYOC9dViMj2S_8MOFAScee&=&U-b5Brm z4X*flws({Ro}Al3$u+vxCxSNwI`>6CbUcSCxn`P^1pi*+94MuiJgz!=xX;)?!?^9* ze8X_S`1krC^TMx#pXJ;W{cx-90`~@PhsonPT*)=lGABXb8y#LI9&KpB*{N)AwEP@< zmXd39eqvzQp=2 zn~40CpC@}Nxh9YHzRdRJPJ*5CdeuwGHTISlBc9V+%)LO#CI1$;dzE?Uo)-;zll^=4G`711&tv;5xiqb0>syYQUSqlP^VI+|vFZF=;?eo{;+b@A z*uS6Ss^A`+bL9RVsN_|RsK_7+V?{>sm>gFM_CkG(dYY>SU`v*R7a^9*q1=c~a= zF2&WxZDp+Qu9>JC^8M&?4|li6-rw^t)s63y4l?(YS+G~WAN|?G{m{d;`q~yB519Nd zI>fvW?{j!u&b`9J-6G_Y&egkWSl`wUvfa`0`iw97>pgUphx>av)7M>k>@dgG(z(o! zACKii-@C7*1>Bc?zR?nHAoml;OddCv`RQ`*Rc4}&>-4T>{|=gqx*_KdQF2Wlb^Czz z?F{aIIrnNG++o4@nOi@Hb#9b%uTgSoT2aP(%qyA0I`^{o;Lve3RLQ0MEx)yrxzB)` zFXvwCgIoOeA?B{0gL*a5;pIZxdz}w%!+=Aq@1t{=S0?9PujJA`*uJcSdF#PVw)gMQ z_THf6n)Br93f6ZgxX0z(8+dE9jHTLE~=Wmd|;r8Ainmb&{ zHF>li_MSAC^=rzay;aFIan%4BIW0rYE z&P`Nusa{Q6@;Tf2B)EI!+`lNfr1R0Z&)H7Pt&np^D7mC_!~S~Kc`fX%m2*cbxil?e zUOnqvw37J^a_((DxVbmfJ7&7&wwcF$K|J5C%z>KrrOa{I}-W6Z>+bNkik+s4)@+uCu8=Qj7(aXhVoy;sS(Nj|vc zW1#cum0T}xk#p}*a%o!3#n9R1QXMHxmUI8=hZ|PMcCJF*$dq&MRC39`k;iJ;&L`mC z`{mqZC6{zfDzA0SbeBu{oWEGkP4U64ezBH$`LMS@&K;}dlFkjwp|e{j%`zM1+*Bo( z>ea}HYPn8rg1y`1+;K`S`7>`S?A?L-{HB~6RC3AQeW|rvuMR-x3OP4T$u;}YXvEcr zs2eqM?p;bQ*&9C`{{0epbWF}2ujEoZmtPM5eh2^lAm>g{a>?Em7eVI}@b^$19(H)S z(D$;5N-ou_#Q7`Jh3TMP%NPWHMM z_MTqIadkk>%~Enn=g2{@Hwtl8Bj;u-x#nEc7u?S9Z-bnhqvVo%yPw|_&t!kx33(1L z7rIVP^}*d74gZEAu1=G4r}^O4hvPoj4*rdmbMH}d$-iGT;+_&f-EidG={~s8U*n$A zh<)K!Id_JVYtHST;hyq8>4AQ zxfEBCnc%wiQeLHvmUACga>>7gCxZI{bWWFZ=PSA7-%&yMw;FkrE9d4ZxyIfU*!wZK z%jMh!KDZm*=ds4GHKHcAPj#!=bZl zFU5D<+j8z=C70rAz)*1AxT1dkM{@2GC6~-=b0xU8ULBWnA5n72zfpt1t$@yL7CO9K z=zU?S5ANo8aLZ7yqU791{czpi7fc+NBA$E7xsNHire0l$bIt+Ojlpv6G9{O$mGngZ zTAgo_bC>(zMw}0w401;-P~@-BYWRJDCe$Ha>>qyu88L!p>wUAyIRSmxC-k6dw)V6eJAIxQF2WlodSE^ z=ianGrFC57@N%K|qXH$DFx?|yuYHcxP0lU!!Hqf@_O`O#Ic>X*>&V5qJCI61>0DIfw`Q&&xcb$?;`Makr>~)`uQa|oqId{F1 zOVf&e#_xmf?>MA;{$e@zDL>pF@jG1y=v*Y{KCR@M{B4B2F4yQ>Eaz_UaF+_7n-Mqa zd;I?E>P-E#eRA$IKDd?N;CFU+Uoifyl5?L`a;a`qe+8ZGbHbx??sG~m`7`wh;_4{& zg`ecy=apR2b<|<}-u)T&--yKyFBf`0+Nk83`uquWcJJM0ne*h_7kqHHe}s5;xh9_D z<=nq3xh8F@q4TG(_c}TEMJ1PX4*wW=^f}IxW8~aTez-O8ul-$Pik$nBl1ncMx&)f7Iu(a_$Z#mvqi6 zLtMFk=Vrb~j-0zw$)&j3_YUIfztDNUoLi#g(zJvFh%5K+87QyfR>`@$d~j=AdyT&C z{!6;&zaZ!C_Hez=soLRv@3^tgjVpc|b)!Vi-J|5vw6b@bxaND5%ei}%T#Boazu_94 z>*U-wlw7iJb|jvQy3dh}pBv@eeM&BIhIMY@lKp9&9&vcN(EH$Jan9-6#Qg=2A;Y}@ diff --git a/autotest/__snapshots__/test_prt_voronoi1/test_mf6model[2-prtvor1weli].npy b/autotest/__snapshots__/test_prt_voronoi1/test_mf6model[2-prtvor1weli].npy index a0e91bdb042c158ab85ae7f65c29a314106d736a..9ec69de0ee02a177b877f280766e8e09e2c7210f 100644 GIT binary patch literal 24384 zcmb{2dt6mj76NL zlu{OHnPnD5Sy`0Xgp+2PmNjOEpk)u!$;oQY+UvJ9&gWBi-F4%i?LFRozO2XI=bU?4 zJmBVQ2M%fED`+$$E-gJUHZQJ6^so+b6LQn@;^L#@dS!Nx>lh#H9th;;HaL(uF)gd% zf!sXLJMwb#8@?keby~v%W7D%dCz+m;o!;<0f&A3`$$^H)^3qcSIoS=4<>#emxkq`j zAHOyG_SIA~!M@^O#lD(;;{QgtVf%CP^7B%&N9E+Dxq0Xj-7hsOknSE2WTfV%yRSQR z?b0RQ&C9IlFn^ji^6`&!oF0J*7#e%VhW~3{b#wS%-}K>ePoH4kDVS>C7oOV|esI^u z^T)WcV)xWPG~O;^w~5kST1G#|`XtcLx^9jg`E)k=@Hutz0w%|B! z`_8U69QK)b0~zcMxf))3jKTZrbezccmDDwDv~dEvyRJdJ9mQ@7rE8{b$mTr9UJFlN z-o5?-@$h^zQ{T(9hTN=$G#st;2Dz zEaZ5P-GF`aa^yRByk~mc$)0s&*N+|Xb|0_TdF+(m#PQZ%CU(zKx~9IZZ%17RVc(Kp zkIq)QyieZA-s-`e?=ty1YNK?`IvTJD^}Pl8y%b~c`u<7jnmn&8VaFMY>tVduja0fE z@4cMotYOIW#bWmykK5bR&%bqDF#Y`SVs?v%!|fz?+bUfXZ_54HH%4In%J<1ArEBWD zyOiVGFcR^~_p3iET{A6v3CFo@6#Ioec=QGD=W~^=S-(w}qQ0Zi&ppNNc}mwzD}9LF z;i){PB#GU29(RJLpY6VEN4(u<`f6x7_f_nekn8&wrE8{bdzjrFaHFnp_;kVRd%n^& z>vzoxcFM={SeU@0ps%8R{A2Iyq8RWSEA(Cg;=F)^1N#e$GbWMeI>u2j8nQat#I@+oZrk$_EY3;2ah|^ z(^qzXv?Jc`qj=ZltYt5Le8};3RJx|V=RWJ0>EDcJzot9C7VK9SD_zsiBi5n56WA-1 zUk5Kyy0kw>t$mKYk}U4)b@JNv-alD`BaK5AD&)r>=u35jgz0CPd%f{ma$%6U$n^Fp{$^8Go{<6h#a@9*s!u}N>RdwKx< zEPpP$Lg|`y)VG@RTNFT^<@e9sm9CjK@=abxE%M>Y*Kd;2HSy-W#cpvv=Q}~Zj(T|9 zA)fl$>wz8dcAvQ(-M5kJ+G;Y~IEVi~`1>I|2DgKP3E{$%H?(}0-NY#&*B4)4`s{P4 z)TMX>y=&O3o)Yppau8RjF5Tx;wBO9r2TtX_s*&%{S1Mg|pVQ+#j`#3X#4G=t@G6fx zO1OT>9nx(J*SB;UT>0m*$ziy=KH0+g?KYkB?UO$b{!Qu9v<;7LW%q;WoNxF0SUz3w z{i-j&@kBp=WET#O+n3&TM0u{td7twsZN9B2M) z^p*U&e~rgoVmB@vx8Zx^rG0Xnvy;8#TOEG8{B!<)C|xtH^DfSFvqIEYzK&9quGuF? zf5hu358o4pa36cE#~mBWKAE`NG1Kb`dGyKqYJk$EX@!-aa-K8iaNK_R{(POsP4UG0 zTYW|SJm{=Fj+wr14v#+hefvPAOZ6=)|BU0EI2ZBC?_;m`xT8Js+Vx{cyxnJ>2j}i( z_sCqVUwM6R@VLqL50r5H_wroPW}jo;Hh&(EKKbVxH{$$uy9mWC7*NacUNoP_GWql1 zO&)iO{RME`=GTrmZ1P;({9CSTi{*|vUUCPIg3k$W#)&K^b*Y~t4t~vbJqPY;xjRVd z(*At*%0ryrcFTF6TrYPAd)zTLa2&V!wIdFjJg?_Ai zZ@Zl1{zUE$3ByhP^h}0bak~!JR30 z$0}W_Z?oLB`g zzuXM3{UC|&C3_y^$H zJolHobCoW|n^**QEAo7c+?^MOTR9W%HuSS!?#@@bW?DX6yMFWK?rln!`gzs({f?RL z?pL&rB+Qe$3zRPP)$VkxBYXY5NABJph8sU>KR++H&t>K`<#P89rAz%BIT-uZx9F=i za<@q7QhlRR&{xNBU05%7@ASBlr{Dkvd!J!P`@iHCC8Mv7qp#kTy9<>rwQo#M)Yo;X z&uVtc-D0In&UuN*v-`Oy#b5S?++Cz}Dc%+FSVup@JtBAS3d2o^ML+)nw@I1 z%sTpB?v^WEGwoUA`2g-G8!vMBbiwa89#*;(UvUNc`E%UooGEu7QMxA2<;b)9IU)6T z!ufJ{h0>*EQ~D70XV<0QqpXwMeN^e1X=SLd{amz{+xQhdkm$2xL97d7$rm%9~8 zm&&Vof z+fiS8pR-!-KBaUi-s20fPr83^X5xKG?ym8;*~0q`vp>&AeOuu^XOrB0TItfX=s9q0 zygTIXGhw*SEUcsExPRU!ch@Ris%Pyq)YrX_HTC^Y?mnw@sh?LAz_rf{zPlVgUGV#j zbxN1wTM)oHvd?9$YL!YJRR$| zyWD+2=~BFzRGrEBs$9DQZ~esZPUeJu>P+Krd|gdY)arQF@%alPN0G{OISk5gVFMxNZSOugTf zyRV1gmb!W7>A3#BFL&Pv!%h5^YvSD_cdL~y?axKSJn!S8c6DqV`Vy5i*bIF$FYc5?S^rAtnd*X<~GHz{3;H{n+< K<=vln()~YhF0gd~ literal 24384 zcmb`NeRver6~-5VD2fu)qI}h~v;-t#bOHHNrs5YE6bhmd5MoI9%Ey`wfgnf_Oo;>% z1kor0f>2P@XcR;dBT$7Jd{7h#U>mgt0#vbJ1eLaP?s+emc^=)9Gwwe(GrMwY3FDYx#48!8bRP$=#JW)cgt_~ju#oQ=gdV4uR^TJqmwU~Rhl56y> z>gVjJa9;@ycyjJJO0J1-aO@6k^9y8ez>`r z9N+K&)Nyw)_W~u?_$n=%InF@Dxto~VOv&Xs%4MCEgX*2%XeHO^9X*|SX@kL)`?oj@g2;qNBv%i3)yu+$~DuX^O>7B6!j~wql=VWli%&LnYZ9h z=(YQf$aQqFl55t5N9J(duDT1l<$nH+l56ztn9F|Nkcc?T>$im;?x`Z?mEDc82kQ<* zj`Jl-uJLouBj61~etU|!mwLFhcEg49wyeB$Wb)g6A;-Dp2yo?kTPnF`+6#*uGcLOa zy5+uVrR178w<^IpG?KYNc^$R(!;M?a@!dEI@s{s%+9T-uK7unBk zJ=`p7B%HT?vLj9hjIY8^v(DUm;j2y#zaN?VTU^NYXDQd@cj|J)HyL{6>(OOOuF*Si z1@qQ{TP@!|$0)fbzlqN=e5Y|7{s}=!&?Bw!+lfOuMD@-diGVj3~=Rrva^y)eKO&j^_<_@4AilF-+rZ%YwGBM4b0o0$^2^h zdK9nZ8b9Cu2J#z1-;l4rS1GwPE$@?!%-tSB-;n$HY7h5okFV^0&W`wa!02tUiFNk6 z&oSroZCR@8dB+9_KDfu37h!%UN&hERJ_exu3f#xkhidw>hu*S?ZYh`aRu;d<+%=QXujGD4VCu~!KZghHV4XuI^JVhn`TdiR zI47KVpYvKf8FB6a8IkjQy@y-n$*-*+JL2O3@^f0Noy;rDWnP~AKD4h7Zs9g?+e~3@ zK;9>BP;$+*1s|}V-kZYluE%qg8$I03o;d$p9nn4~D*PeGIqiPbliXJ~DY?eaIlGwK zW-5;f@;-U9l52d`_ao+2Om)m<^5pCOEgtSZPn@lf?TC*DNbiEok2${cr!lu4_o4li zT(j=4+0A|$`ar$8{XN_}ZNtR*pYxT`yXq76bL;8Mtw)@1RdQ)s`M>rs_t{Gj~#^8gR`HKDI)T1>P3?C1S?@U^^tZ}Y(&^~$Hr8#06SmdV$_ z+dbTso`U%~U(x#8@%$?0R?Xmemvnaca*@}gfl4m*$=$;~W1SNoM4V&A+(900KVg2& zek1*J<~E!OU&;IApFP~Z!Z@40@xm9(E1JnV>ruahm0XH*{_-!mZm*k#`jtNy9pd3G z_vp2UcEra66z2_@2iZ@jWE4+qc1CuAd*5u&-~Ib4U5$ z&iwcg^G>u_T7Y|W3G04J&P`Ersc-Z*Zov+4kHc3R z<=inI?qkCBmwZ+GTn+atSFgEDg`7K9$)#!WMbK+~wNK6+=Yv}{9^6AqSoa}0cf68o z_U+f!aKAbLU!9P1Cn&izt+Z_o*UuNo?>P@Ue7VT$ej0!9bib;PuezOG!<;JgjZ5U* zbS0Pi#^9p|nY$nPZ7=6$D7iE(e)mE4(;novi=3OOT$t8JP1|YxxM4T7Pxj9NM>79D*L5}ZdsH3%V?j$9b zrd4%Bem_T?Uzc+yd$_xBYzO;%!H)Jk^^L?9XY`tUZ_ zcOE;y`Q3*&SIfEgE4dWsf&&K})(Jl!mvg5oxs=bfA0Oa;^%43;!#Vst1kaJCDY>Ng z_-)0JGB7X2LbR=`)C<=i|am-1V<99&ngso&mm z?hGZ@#JL3AorrUSoco}XOY8okB9E^U<=mM{uBoG0;JS5W;yhN)ou%ZG-WIvwy8Aj4 z_epZ@LrN~yZ*Cg;q}^XVB+Z)4I4VQA?7$D~^RB~xr^r>pcOm~0BG2?JK_fa$9%yWKe zpQLG}$E(@LF4y$QadPe=C71kL_Ek0W8bj}7IrlLom!{>_RP#C*wUozvIrniTm-

HN7ubjKY2RC{X@>`2_|Dc?^ zRLQ0Ol=vFB?(cLa&d254WlAo^wS6h#Tmx>Sxei}0@_yrKC71N|UIVWCdy2{Lg>vq4 zC71l1yb^JCpA(um$H=)Wd~n-6js9$(Pj;4bpHXs2-~JNxNvpTFoV!xVrD?^FqHnnO zp(elm<=j~ki{xoebM zs-t$rsNe6}P`H1rc(7R2}eNoA!xE9U=w<-2lyXD+c zC71ks(+v3eT->L8Dd)bV*NBz3b>&)@Za_(zNE=>!K2iIMXOn)9C=WbAPN#BaGsNe6QH(Abo zUCA}|n}Rwz0lgtP_YEbN^d={R>)w~qad5hvyHUv{z3WDyj@vhhu$hVcZ-rs({{SP;^|m_zm{{i z^2F2sKk3cmFXi{hh%;Q{=O5+VzbUyC=N0aA2CJ9)Na_VeIF7u}c}K}L>;8%}T>2it g)^hHwxil^G6RznWUC(g;59Qq0VgLXD diff --git a/src/Solution/ParticleTracker/Method.f90 b/src/Solution/ParticleTracker/Method.f90 index 89d4f7e4b33..23528e68c75 100644 --- a/src/Solution/ParticleTracker/Method.f90 +++ b/src/Solution/ParticleTracker/Method.f90 @@ -229,6 +229,7 @@ subroutine prepare(this, particle, cell_defn) end if else if (locn_dry) then particle%z = this%cell%defn%top + call this%save(particle, reason=1) end if else if (particle%idry == 1) then ! stop