From 4f94243eed254a6c093443ff6b8130833af624d1 Mon Sep 17 00:00:00 2001 From: Phil Elson Date: Fri, 21 Jun 2013 14:24:31 +0100 Subject: [PATCH] Updated sample data NetCDF files. --- generate_from_data_zoo.py | 14 +++++++------- sample_data/A1B_north_america.nc | Bin 1761859 -> 1824028 bytes sample_data/E1_north_america.nc | Bin 1761858 -> 1824028 bytes sample_data/hybrid_height.nc | Bin 659594 -> 677731 bytes sample_data/ostia_monthly.nc | Bin 1700161 -> 1713386 bytes sample_data/rotated_pole.nc | Bin 12402 -> 18809 bytes 6 files changed, 7 insertions(+), 7 deletions(-) diff --git a/generate_from_data_zoo.py b/generate_from_data_zoo.py index 2bf7cbe..88615db 100644 --- a/generate_from_data_zoo.py +++ b/generate_from_data_zoo.py @@ -23,7 +23,7 @@ def custom_file_loading(target_dir): def hovmoller(target_dir): fname = os.path.join(DATA_ZOO, 'PP', 'ostia', 'ostia_sst_200604_201009_N216.pp') - cube = iris.load_strict(fname, iris.Constraint('surface_temperature', latitude=lambda v: -5 < v < 5)) + cube = iris.load_cube(fname, iris.Constraint('surface_temperature', latitude=lambda v: -5 < v < 5)) iris_cat.add_month_number(cube, cube.coord('time'), 'month') iris_cat.add_year(cube, cube.coord('time'), 'year') @@ -42,7 +42,7 @@ def hovmoller(target_dir): def rotated_pole(target_dir): fname = os.path.join(DATA_ZOO, 'PP', 'aPProt1', 'rotated.pp') - cube = iris.load_strict(fname) + cube = iris.load_cube(fname) # XXX consider taking a 20x20 window for meaning cube = cube[::20, ::20] iris.save(cube, os.path.join(target_dir, 'rotated_pole.nc')) @@ -67,7 +67,7 @@ def deriving_phenomena(target_dir): def cross_section(target_dir): fname = os.path.join(DATA_ZOO, 'PP', 'COLPEX', 'theta_and_orog_subset.pp') - cube = iris.load_strict(fname, 'air_potential_temperature') + cube = iris.load_cube(fname, 'air_potential_temperature') cube = cube[0, :15, ...] iris.save(cube, os.path.join(target_dir, 'hybrid_height.nc')) @@ -98,7 +98,7 @@ def COP_1d(target_dir): fname = os.path.join(DATA_ZOO, 'PP', 'A1B-Image_E1', scenario, '*.pp' ) - cube = iris.load_strict(fname) + cube = iris.load_cube(fname) cube = cube.extract(iris.Constraint(longitude=lambda v: 225 <= v <= 315, latitude=lambda v: 15 <= v <= 60, ) @@ -136,7 +136,7 @@ def ukV2_in_userguide(target_dir): fname = os.path.join(DATA_ZOO, 'PP', 'ukV2', 'THOxayrk.pp') sa = 'm01s00i033' ap = 'm01s00i004' - pt, sa = iris.load_strict(fname, ['air_potential_temperature', 'surface_altitude']) + pt, sa = iris.load_cube(fname, ['air_potential_temperature', 'surface_altitude']) # extract, via indices, an area over the north of England pt, sa = [cube[..., 290:494, 190:377] for cube in [pt, sa]] # remove the temporal dimension of the surface altitude @@ -151,7 +151,7 @@ def ukV2_in_userguide(target_dir): target_dir = os.path.join(os.path.dirname(__file__), 'sample_data') open(os.path.join(target_dir, 'version.txt'), 'w').write('1.0') - + global_map(target_dir) custom_file_loading(target_dir) hovmoller(target_dir) @@ -163,4 +163,4 @@ def ukV2_in_userguide(target_dir): COP_1d(target_dir) lagged_ensemble(target_dir) custom_file_loading(target_dir) - ukV2_in_userguide(target_dir) \ No newline at end of file + ukV2_in_userguide(target_dir) diff --git a/sample_data/A1B_north_america.nc b/sample_data/A1B_north_america.nc index 28e80e8e6a0e9658ada21da3fb258f35c78f593c..a77fc154a6873f7524d759dd9ba79aa69abfba65 100644 GIT binary patch delta 36213 zcmeI5e_Ryh-uJJC1r|kJ6k#1Pw?xI?Xj<6uBq}N;DJrKZEy++ZNy#YDOt37mJVmXw zWM-IX_@fglDWw?}DW(;b7O5F^$jCG~9nXC@9iONBy1p~l?tVYy`SZSC&+~F#uVr`N z@6TM}|2fy9zZV&|e|Woj$m1QcCiap#J-&)JzV1R0dr` znc#L-wEsGD)U(#OJHt%TCP>)4T@D((pRE{>Hppb5@+n9@#*7;X6>Zy?Q#5?2scWn$vqh=* zK{lJM*%ab8YjAWDHaTigP_8gXmG1-KN<7VlZ(bJlo!Pav~8=5)&)Z7)Kd8wnTIkW5uE;FHKmKsToph9{{ zv=Kvl$+vPiXg#0n<;e+lM~kTnn&r7gj6=i}$oHPfLD8Nv$AC6Z%qyKYX-vw5&}+Q8 zhm5O6%r%qZcv;hyaZ{&`6&-u4@#T{FyG^uH?MZAc82*tO#nF+b&DNZ6F2^;Vy6Cy# zx&WJZ@iwEk?RM|=6~>G%->=v`CUBe;i_LLrR4U%PWQQ@M!1ei)J`s{7glN>tC zUa)Z1J$FBlHzRk}!nq5wF0?nt-7;+aFQd~_l4YBg)73#Hv*P!!H(eW{$6k5p@uz9B zsjRwRK3Jvd(6i0vYnH!x6Al_3;#9j5;MbW()CuFsrSl$I8{ESEPNRIhEUoJUY9L=(*ms zZVOEt746gTs9_dp{k>p4V+*(M+Te(3+aa(qlyC0pP1eGo-O4g6e}CQa zK3r}L4bJK(ygSQ{Bk3!<9gU?HLzURnYqAK>-XjyHoAHLS#9+5RvtDw&8+FKN7~pE(H z4%Gl7CC-Y^>TOzo1NIGIPB>u1d0*dRw29pp zZnkIhudq5n*|9!->qzhM!^XhQW6VDNECAmTaKJ(^X@ok zM0tPQV*I6pIdi{L+ZZ{Ns{X*&q4?Lqqu+m7EQFgjPKw4noc7DOl;qLA9*oBBcyZx$ z@9n3Jq<#t3suC@^`b7GD47`^;q^Eswy!YH`qryA>O=FaI?Kj41+d%K9Ta9R&+dHtq zxJpiOnQbO6jtddy{0$f3RJRblZtirEQXy(G)jBKt7Hi|{U$m(~dpljGp>BKENDl5B z-qlkW)Y%(WZA=cHKSX$^RU1=M|8>5giCV4t0K#UsW|y$foS*f8?f}eLsTY4=YqG3% zTJho!OH~OGcX;;OqP^9djnwGNO~(gDY^L7Y=w)CQUMe>(YM>2-d=6=^iVtQRdPfYg zpy1+BW>_r%Rg79Dx)TlZ7FQaEH|d};wm>f&)y@MKR(ttggs49>ta|*5|4Z|LdLf(L z?0sa*%$s}PtQl%Xd@;mKn$!BI6-cngw`Rmu--euGYhYiwoGw_qB>tt@*kyH%9?|#w zSYyn}#u%G?1j$FReAwk9L_R|IG{%Ib55AguQZZcUE)%;$v+Im|%v^xFJG_B9$^2J1 z{R7dXdpR9J)LqZ%ABpB&ND!S5NPFq9futO(9yG%;`-;8#6h*nP}o8oK7IR`bADt zi8j2+=`BQue#YrUqD4P*I*I71HfG15KSdhRgx;J^Cc0t_r&Ea5&){?_(f&(0O($Ad z!s(xg9;@W^R-*Bra(WxlrRO=#AX*o0?iA@yaXZntUYy=RG=DUw(}*6ui__^u-4Anm zC(%4Fr*{#p-Noq)qOoQVmiMo1ndJoZ@QJiKGt-h1fSwy3j za5|f4_6wZOA-ZoTr`bfEpKv;tXx2}h-b=J9EP{F>{VC=VjfmqkhiGOJr}q)9oX+Wd zqK?I!E+9Jnc}{bQR_x&Pexmk|IbA5XnG0;@Zl&xD^|r*n$0*0@)h5?E{!5SbKD)(; z_N0Y4BWtYPXZ;z2sC?IbW*4lfI*RKanX`&hJJIS*oQ4pMdXLjkqS>c84I{emS56&7 zoiS?LBz(KD;fl7&x>i%~->QvXkyVrQd$WwTG;Tx=rx8Rmi#TmZw9?v-)StCI)#Esz zv!4HBGvus_h@AelMrKt`73@&)FGjuYH;ot>Zxd;cHaF$E+>TUR#t=?B5#2tK(IZY&5H=EPp)cLJUyeUjR9P$w^A8q8Lt$al6NxZ3DL{QZ#)erwi;77?wR#px=dagTEPB+>l!oIXYL z=-ZsGChGo@(?1i<`#Yyk6RqvwWS`<0qOtusEhd_q!s)X_YqB_fj;L!nryim?FLSzv zX!Tp1K2J323r=4kn*9r>UZVTjtFKf1%c6v+vmd7~63rUV=~|*y_i*|W(TGA$*AdNJ z$LY&NE2}tNPtJL@U~LWS?RKQF~ubUnQC`j?;}qx6kCXOj5VAfWu7` zO@E2g%|y%haQYh2OP_OEPBg8V(=9~HBGh>U|DxDRv?ZR?ZDvP#&f9L~Id6OU=pY|X z`G}N{j(g5|I|X!<`aoKLXQ;P2+oaAK~;bL{rvs`UcT;yE&~O z+WZ-(JBTJ<;B+U^lD7Pt&s{|O59PFyXki+sZxTJ0!|86K@vAu9Lv-mTPOFI49pLmW zqH(7=eVb_hubl2BdNjH-TcG=hy07K*9in-YINeXQb{?ni5{)h5w3=vc8K(z`)>L!) z9#PlVoW4&q=U<%G5UuX4z7O;_%MXY~4dL`4(d>zw9wNH$UQTO?I#+Ufm}u5UPCq1C z^)9DJh(?^^^eEBH7EV7RS{X%O0{XU`9}{&9=JXSy({JIlj%dYPPCq4Te}dD`h-SRX z>E}eZ@8|TGq;BU)4v$kb{hyqELA1P+d2&jx#$OV>G>FrBqG_p|enqq_o6{3STUKy- zl4$A%PEQdneTUPpiC+AQ)6+y#{)^LZh_36{g*B)KqRj(2{g!C*1WvyrS~7>zzY;z7 z7^jV9Pg#1V)5_8_QS#AQKBDC#Mm}77O3!rh^^_4p?@H>UZas?|spzCqUi9}w*X-l; z4AG`~PJbXe;vbwg5iO21&*kXaek9sBfYY-?6L049C!(uob9#=b`x#Ep6V2Pk>CZ%K z4{>@y(NVFDCZ#1M|6H(M^t7G_ThHb(ju2P+&?|kY*(9#?p)?=MWSL!o~=R=qFE^sWX_gyEA9n)(doojK;Rv9p3$i3|Gf{>Sq0sSev~ zHgP+RGj5xw&A2HkX>{H%NSKFglrgtV?SA`BU*I5%QK!tOMK?UHKlp)V$@|Td+dYE3 zC#vP?^vzq0$bvK0$!**s&BO_pc%Qpd^4%GKwQT%FE%yX;rDcuG|&%T*ox%6CRz zFKDky&N(&I{P4mEv0RPP+dM;An?GoZ*sPQ+BSbz8IkR!G7^OgWkwhmEXFRq1Ns+0D zEJS_QSrIpS#CTQ64f0I75_8w@5FaV9TF_-B(HGC}J^Q!!R>%c-r>%gdG@MO`uE9J@ znW@UK3TM;Lri@n^_#xO#U=uQK*x3_uB9_~R{QlYB-nPqMhUFP->(ZUBGm@XI(0B1v z!*k|BQyf{7h=!zdyK*`@0d@F!UBEbRj<3a4&pkT z-O;R5uj6vTu^|`yAD9&)pZfwa0s@6RApQux2>gxUi^0olG}utWyu3gwiD`K)1q%qX zfdypB!H*7b>SYD^G2p8Ly}Cdw83NUT0qMqquLVC2d>#1l;Ohgts$Zl)pdlcj)J@=< zz~2nM8TRZ`SX-HdE)^zJN(4YhgXdgP#IE z0sK_(iP~HBi*yJi>43%m1U?1)t>DwZ-v&Oz@&$T5&VWEB6mUEEZ18u0&jmjXe7@-m zXe+1ObO;ne0e6Bg0)H3yV(>G-mynn9U(bO|2$VtrcY`knKNEZf_%)K10C~0l7&tU(?`^I&VhWn|TPxExFe_iTfH%+N1%x8`A~R2t`Mm?YdHTuc z3xpE_9emDw;Q}A915N=q_(&K(9=!QQzYfU;p8x^#>u((pA`!g##j@Wgfj7Uh_4^d? z=9i(CS1Vi^1k5iQ{RL!zH^0L3`%Lg^M>TNB2Cw%^RfbmhTnLy;)L%qCc)D`YuSFsF zZUIicECPR(-Ke_^ZK}gTKb-&H6 z6p#Tv5qu{2;cDyWJEGKEJ^a=Nzx7FgZRi6~OCtC@@JZYkhzB8%!UN(V@M+)|fzJS+ z4?dH5xquf#Ae#k*SOPv5{KMe$!7l}080b|4djtYSfdT23fiDIh&`2p+0=^K&FAeZC z?s5o}2LzOQ6nq8v$G}&CUje?_@6Bm{90Ile0GTJi*MVONz8-uL_y(W1=G`g?H2DIi zcoKXw_@}_PfL{$>bhdVj^s@Ui1ROeG@u$H%!9N4u1-=-(+w#66%4Z=E4+T62J^{Q3 zd?NTY;FC=6Kcajd0x3|y3*ghhd%_Q*7lVHVd);*S7l^+=z{vyR4e&1T72w_AcYu#)Ubd;75J+GFA$EaJ1YZe03H+PjQv$te)4L&% z78sCj5BLo5Rp2whzXd)!z|%VQHUx450!r-#pAUW?_(JgSfG_fUs$cAfK(RkS=3Ve5 z;H$xxf{u#d{E_@C8iqKKLr|HQ=kke*nH#d#iqN5CU~NVDUrX>%rH8ZvcN7 ze3RvU+td#s&K1F%8{PhHU0f96q;7jls;OoI>g8vG9w)E8g^aPxMKrR$; z5_~@RQ{W51e+|AU2DYiEAy5njd;`7&d;|DW@ZW+jkAZFKcMzz60{#lV3VdV0HZ?-% z<5PXzm!8j!P)aqF@;&%k@Mpl+abF;QfIvMDh$iq2;C}?)1pX}eX69v^`UwIpEFkZe z0$#xN8|T40!2b;18R%7;z5oH270_6Fy+|^kJG-(?t=6;3dcLwvy|PVp!)knGn|fuN z`aiQxz3laYPI>0j?Yzs}*>qZc)o^t>G_rOor|LiWVj>f>IaU9z5)+xbgj02Cc1&dX z3X|%t@PcQOyg^U7V?_9#L!TIz!hHWUD69#F0Zk!V z3-l(U^+0bX+5~h0(H5YoL>*~-3Ee`}1#}|Oc%YMrCIU?(ngVn((F~wdh-L$wN;Ds6 zI?*Daej1vGd2v&{;%FfX*gb4s;IDDxld!Yk|%sS`YMIqD?^O z5p4mQL)0;a&)@rqx`56n8V_^<(L|uRL{os?Pc#Fly1B$4TdRvxeEqK+se8EF+5>me z#VPUz2#)XsA5s^ncqWB}M2SZ7a) zyw64L1U?JHRU*7zQ+0_+_sAyw)|$MwRCo)}2%?TuUSK<-E}-qHns`k;muy`<2_g2# zL`|%k)OIeWLkg#KFNsbi9p>%Q-EGgGy~{0M6-eRAes*}j%X3k zoJj3Zg|oA17J@^a-NnKvxp20$M~gcOmbH zRuL@(`Xtd}pidDk1-hDO1<*eetp@ru(K?{d5N!ZjOtcy3vy!@;5AYB197U0Oh`J%W zhG+uN=ZPi(eSv5iP%qI;pe00efxbwz5a?Q>#Xw&oS_*U>(F&k16Rifio@gD=SBN$M zEhX9vbc3WWXC9xuuTm80Mxt)WmJv+=x`}8K(9J~CfWAgF6KFZnT%cQs76RQ$v>2$m zS;V!rJ%_$DmWFy)zG1BKuB|X?yocnifPMX=J|5Nw()##y8qus5N!s!LsFOXLB42qQWWVfqHf4m5={X5Ceb9IyNRX&-9t1JXcf_1pl=Z^ z1o}47VxW78mIB>Jv;ycmM5}@BCt3&eU7`&@tBE!PJs_#e`4FGI?@<)#`$XN4ts$BK z^aG+vKo1g41A2&PCeT`E8PlDeGveD;1vQKWw*>V~YkL4cWW#6nuWMj9y=g&2N`NU?$(!ED_8$G?YJw|GOHJTdqv;Of4 zlM12ls-WI`=j=3Y^nO}tboHDJiSV2cp)#iT_)p`g@#N1{$PT^I`jYc^ewe|oy0wC& zUFH@;PBIO{t!Mq@+{Qir$3pI2?(yo&()YQCD@hDdcf<_n3iWjk>IWOiby;4lSa(<2U33$L zuA=e5hI_;VG+N}8;1cnheA!q0yJqzPp{{C|zgP3+kPi@vv7Sah&3_c+M1KuM!z*bVO#yyoW_Wb`n{Yw zIdMPs`+29FxJ+9wYWMmle^lt9ojI@B?0U-Vs=@yzyz>tkZrfb%Glz`08}-jzT;^8= z{?m>2ty_#3ws7yXLq^vgda5t?^D8!cTkpRQ8Qo)@FPRI}W{9CMyQ|x{7CSP%onu2(?b&rw7%2BIq;*Qf@bhzz+VjXvi??gdT0p@ zNIw<|xCDM2cu~d|(RlFo0I%wAb*BeMKtTDM)SrMM0{qPZPQ7%3p8!6}@2UP)cY1L7 z12pO_P(Unr)h`7OZt#-=rFIQt%7F zmz7ypl;{b_g+MtJa6kC%;D7&4557G|wVxTDVqLCeH}CY2>U*nBq{oRg{obl77^+NP z($#z4VIww6WI$5ymgo_tgZJ%0S%FOOzFjKIXD#XKSsrjphiu6EcCxH|j`rp)9h|V4 z^6hn5@hI?m-^J3x1-_jY)^jX)xqA~h6bLs2Is|I@Iu5-4UPs5`!AELgJ@*Ii+l?eI z*Q^8x$g}-y7z#d$S$WF`@M=dma2NqTI#8-LI|%|YfdTnC8N3U83ivMIQ^6Yno_dxv z2$(x2{!}{nSihoI8Q{BtpAP;izo+^|CIq_s1N15jyc>Kr_^ZL^fWOA)&H6-`b8lGdh3AoyaIeY_#*Iqz^?}1*Ydi4Q4E28P{11S>R>Q% zC;>kpP|DZqzz;ONuYOSqfk99}8TbV7<=_W{-wu8Vd4K()0s_}U0hQo~g0BKU4E#Ru z*C{_eAXE*3>jRv6Spz;1d@Xo>O9#Dq!!I}%LMaunu~Y~0fkP$uyg(^mS8?C^mJa)P zpw%rMs=+UU3TnXTgRf;?F5txwILZR8Z|P76{$UvZ82F{&>jS-NV2?oHRA4~5W#AjY z2Q*SjHi9pN@tXoXjk_EI=K=ysJqo@V{A1uRf?ol?#qZ5&e;fjr`~fmgfESz0ZIIsM zTM6C{z6iX-=dF3S3IY+nfGM5??*#u8_$ctJ!Mn8A%kIw*h}8j$KMmdu{u%Ib;ETb> zTi$o*_$&naLjli$PXO-$KNS2L@QJ4PA38n{fe}!^3*eK$d%-7zF9DxIUe15rh+c$1 zDip95d>Z(dz^8*>2R=i2wfyx2ybOWqP{4Zdnc!amp9Q`Ye75w|{`3TFfItou@GAIR z@EgJBfiDA}zX=W|_co+DOz{i^2zfJuZ0&Xba z6Yz21>%hl@{}gxjd$#uwum?X8|Dsh*#H0l delta 3281 zcmZ`*4N%m_6@UNxJ~+5T;EuaJ;iPbUADkYR#xV#7NQxf_sA(~LBzFi?6;w2uGJ=yf zWTF!^Tl0^gPBYq~G8v`OUNRFw3z{ZQ+i3?` zV>op9sycQ2ZMDD8`zQ@aK({PQx7%zCy5h2`@8I)}9`|i^mwzJMepFqB4;qz4?X^z?RGj+35OS9*m?c!V&bj7UrSuJCisyE_{@!n@kqL!j=U!Hz0yjUexABENU-4X{HR~E zMsvkZ&uE{ew>9;yG25X{3m`XX*J#!|5IN$Gubo+Pre^_Mw=PY3Kueqsx=Wi~Xj+8) zK`Y&LsFbe0s%}|q)eON8)Da&iJ@YG%LJbO1d)nwCMD zt2VP>F&jP>X>5ro50>jkXo0C(N~H%)`S?3*ayZMiiXOO443FAU>1~rkDyK81d}k_) z5VmNwo5hsV5_7&;Scs)1IKxVpL^^{q^s6`vJF`Kn zqihN`iT#V`M|Ku_Y}EBOjXD_gUcy_da93 zr1`o$MQib(^W;e@S>mWE+k1iO)cB+isa`ITi_N4htCe349xb$1{?g>$Ie5#r;7$SO z9z%Fx(G&Z!Uitfr9f97CfP|t!(W2;3^eBc%N5EL;?=3rEL@}Y5QKC?yQ7nPV&0SrJJ ze8mP!0Px_@lfHm!W&4fO3tXgIX#wrzITy92-;7?lX!Oa7T8tY^QD;Veau?5CL=6*& zep5zpb@Vh1<@{!`IR5W*PhaI{h4I#N-X=N4Myzoxf?9a^Yc+!X0P@K5MlbvK)~)Hx zxuB7YDx#ZiHe}-Y(#BjlKBsLM7TSk~OfNu*`GxN@dYz8EYbwoYbb&|7TP)6wA6HMw zl`oLYi0w4#*UD4-$ZKw9PB8Dipc?#{IvZ%Wp;7m{oNs*cPxh1XZ$bMxvvC1UF*dp$ z?S1Fq<9#c82goY3pF6QiZk{07^6>%ESM|;F_bw|J#LRfI(t#J`rpjaPg1md`oz!q+ zXWRDhXMEEi9NP2heC~roF4~<|)P1md!mFgIgXdHu(v$?T$}!)#vg@azFyTjDy+!|D zdoV0&>;8JQUg4!@yjVsFWgg?1FAAbKYWqZ2m4Y%blojAf1{+gXo`Ze3*}83B)v(`j z)GgpA`vfV~UKn!8^%LYuY>)e?O>C2vAj^`J_<}W{q?xvsws6R^wmlTy=2`A5D4hNH zDr2|&c+&+7RYo~~fPB}PII5Cv;yV-d!l0B!x!X$0qePOngZ!R}EHki7t^xU_K9Y`G zmo6^3tdKp7myP~YzgC_(LTp7M({Q^m#H+pdf9$%tu=xLQo3y|-*D0G{ub?-U%Kc= z{u3o?FpZxMIG%W>4E)My@BZIcLM?u#!!N%4^l7Eux#2`Zbpx<{ z5xBeO%jdv~S=<3}HQG1eGLB}$zsC<>18)Kjp64{a zS^p+Uiw>Vvc}Y6M>t>Rqm2`UAXLLmWml?q(=}hv~H%Us1$X1e0_|MQ9Zrj@0nQ}>3 z5T4Yr>?SUUq!Vu#-LpGFyYVT(?Fz@58M;x?U3al(+ZYXalcO=$DGie!ea5?3ez}n( z(I&l3-XBHmvBA%|JE84XRuTu*XwMiQYHaz-Wh)nK4b$}|4j6+8CkmC z&;N6uMerHrq*-W3*1B{ym#sE`>fAI`oQ*FUVw3)*k<*e$*Whofh$+&si}V=mitEKo zlo^4<x@o(5*(#XzE`C0$d3j{vc;o%*~7fV I7)(k314|^ZQ~&?~ diff --git a/sample_data/E1_north_america.nc b/sample_data/E1_north_america.nc index 09bf79c597afdb086c1b310aeb7167f69e36f9c7..a3bafdec97d986153bd59b0daac6cfaa89fc3ad3 100644 GIT binary patch delta 36143 zcmeI5eOOf0-uKtY+n}g}@-lAeI3_A8D!x=yo+KqCr6lDPr4=bEmL(MxmIan6l@&GS zP?=GpjW3-@NvTY+NHMLbtflUSt*ZLj1WsW%Y zxtqh`=+Hv^WF-BR6TdG!d7haZgY5>5v0o#))8rYUrmqkkOP&4ox7zIY9KX$ z3YmG)U5pzdKg!3u_Uq3>yo`Pq9#w*+sEb7c32?jCUA^vT)N zBG2|M-f#V3!a|)CpVqW*%9=Uz5)nI84KL@-Nr#A0OO$9|vmiqa;_O3hb99vIIl0=~ zoUkz36yWfcZn1_&T^*AmM1=4S{Ko3*n_g>W_4|JLj>(~eFJGbuZ&U-MM8?SaO#aDN zt*kQlCl3$SEu?0IcULd+(R=>l+vd%^;r5d2i|5_8aB<;@Ub>u3=S=zAq`d4*xvn_aDEB8P67>x~;bu6p;0#%2ZnvTHo`kQa`KC5@KTHu5@YR|4*8c6<8 zb3-;a`j-Rp#XBl8*D=Uj7CBn}eI+unUsf1tVXvloF_&nN_nXK>N0#?rk)!(!vKwI* zP_v$8c_&7V_Ff+8^xhfeitcSMAHBU3BNM$Fqg-<6ol#S|WraQIOi(FDy!Y{_MDMxI z6z{qyr+2>7tzBOD=IZB7!4kFM@eVay=Rn8-AN^WnPkQ^ZAkv zh;5w5>1{+i59uDsWVaAKK9p8^VQ8D8!&>jRbb~U z+-{{0oHl3ugL1Hm@A%_?f0gpx?YG8<%^e|pullWn-G|kD-%Lx)VBZt_t%1YdPS%y` znfTXn`qLD%hzFjEDDoXXV2z5Otb6pYj%{Bio>r}%C`6Pm@;z&~)uKMiURf!FI{!4zntVb4&UK})-_}8l7gkEHiniF99{}X*@7G7z$Fe} z?LI5c*S^X6b8kIyzf)Tn`KVTXL7+1Uu4sq8|GZKN4{e$hO}r=fZ&}%ylLEaL4eeU9 zeu*#tsFg7+&0Zf8Bv)f(J}1Iww{P!J>kD7jYu3fSs&A}Sj#0km%~pcL<9mObl^EqI znzwL4QHgBx2FDznAtz$=6pj|rRJ#yEruIEasR-3jwQ9;P#NPM@XKqfnH~XGudEC6m z$_yJG-Cu1I=7-%pR`{-Iux8{OTN3p)EtsnL!VzvyED=8EmcrZ3I<61OBUWd z?|L=rpN`d8-MahjK$pW*TuA)k>u+Y+Pl#tfwXQJgPKM~I!A?hun7(DTf6au-n7SNc4s7ZN@6CZ`jK zdRjRBBhiwRoL)q<$)zWSscj-Fv5G}rx(~F5_-pA=AqE*juI+^J49h_zo&G>}V zOMH$kmdAC1(=5ui_t2}WY3~%G6H++MCR&-v>7_*53OJogH2q#qrx9KCG^aU4TVLb! zGNR)@=Jax+6~A&io#>J7x*srOkxMjfD5qBtT|Sx98AMyI=X55~)TNx}5iPIc^iM<& z*Kv9!(Ugxky^844UpUPt+8nKSb<9{?O*DB3r`HfIoy6%ZqKB^IbT(1XJ)B-kw8Y2h zbwr!Cb9z0|q-IVFh!!8`^ai4hPQ7Df#$pap_h3%v5?ypLr-eituH|$d(YU)golmrA z4W|o;?%Kv_5mDENoGv6<_%o+B60MJlp{~e`#Z5$Gk~v*Ov>=1in~B!V=JXb#&Lx~K zCOZ2mPK$}wzRKyXz9$aI_4or$Z))I@iMkT?=EJlhTG8%>;|%q6G+0Ao>!+K~x{Mw)Y|J7~V~7@1aN3h-T@9zbs2=BD zL%sjUVaW*<6Fd8BgLG9c8oc&DIE~e_CqZwxP2>AeHh(OqeTi&)qCe4`F`TNykO8sPmvMR)(US`~ z{R7eL6`T$ry7om*Jw!Wpb2^Y{=2x7aO|2RXyUvfHv=&FBmno6{_kKQjcV=B3dzx z)6qnaEZ{VaXxeg4#}Hj!&FNU8ExR}!M>O>dPRA21{|Be%5Ir2r_m<8jnlg&h^N22; z%4s^$=J}kSPyOD?^h=}U(c z+#1P`nM)+~xCU`}2SxKHa(X9G|13`LB6{jBPD_dAuI6+J(VALL?PsOJ7>gA|Q$}*SlIYS)IjtbtJde|dh$i39>BB@z*Kzs?(L--=x{9dh zb50*6TJm>JA0yh-+r`G>aiU47oK_Mo&gS$9qK$=|K1tNQj8iYsMbB}%nrOr8oIXV~ z?lVr;5H0!}r#_;)da18Yg7cz^sA~kLPZKSi!s#*sR%EYX;9PS+AGSj*{iMC4$v?fNKAqdWj%|ts>INhQ<%HzJCE06nn$!Bl*bjfF|eD>LS z+}Agxqcq3R=1V}c)!9Nt=Ul>z-b%Fk22Nildh%XQ|4cOd8BSjzx^@SrwM08U=JZvf znI|~iMzpF2|0;7k(bVyr))6hw<@7b8hZk|WgJ{Y_obDvLbOWdLM4R_=`a03%qny4$ zwDjMczDe{@LOh$HyNG(marzd~lIfi8Cfalpr+bJdRdCusw77=Ty+j)uIDMO_`)f|$ zA-d>4IBg`_5U;-b3oe%L5{({;2h3h$ek7)fKP7e}| zIl}27q6M9tzE89+j=u5>Y;!&!>Kw!AheT&z#%VLr+J&5cr2Ci^oPJE%{1-X>gy`1Y zoF10chG&hLYdY0bE2n4bJ{{QH;2}#l*Y@eh7N}OD9iuq?mT2ZQPQN2swSd#V5Iz0?r)|2Y ztp1_z%IY8DQ-g9(|fhf1@wtmzVre?dJYAu7w!J{TA?PW4LeQ-Vz1p>ww{=%W}%?z=Ei7bp}jb zOJhaEjF6|xyuHFwXm@A9tj#Cr>=+I02zUD4+Hbl0%x#|cvnff4#WUCcs?N3ES7uI? z>2C|=)i`rLzLj66Ax~e=qOai+9xeIj9fk8U^tV|ezuVsqQt~(JT**pI9_sy zt7(|5E#B@~GqQ8(ym2?757{WAw@f{M`^{XjAd69_$7dz9K4wOI?|qrO^@;33-F#m( z$W!MVHd(P{KiDU-P3LW1|FcKtkx6+uS+X*}8Xw8knT$7T{_PQus7htIs)Jwo)}-?U z?N!Msryh+zu9+Z~sX>P7Gu%z3qi2YXO35-JS z^?Xue819?6%^KpHyu-Srk2#I2_7Gw`>QpXWt>fLA(9`F6$r|SSe6!U%e9Ye;yV^Ht zuXWyeE?tb-K%3Y17mvG(manmPz_DA1Nvg6*Uw?9|5Ia}zw92D=lRmV{%go7Ev(n0p{1JQw_=~_-f|nO!utyd1@(QdBrsV|_ zEFkof6_W9TpA_QM+gk9G!Pkd+bp==^1R6pE(p>_+349j#X7E$Mw}g0AzsQC_Ye+z; zOTo8;p9;PM{50^LK~MFI90-WMd?C9Gyc7K8;9cOSgLenKu0L@3rO9~$0Zpzj!hTHw zKLdOk_?h6-jkoI;c@W4j0h|8`d^Y$i!RLa%3Vgop%glP54}k(G;A-$i;I9E+41N~) zQtiuV^QPQv2$Vws*MhGAe;xQr@YjQ{A}^=EnF0k6sD=V=0PhDs2YfB~x!~)SSM%RA zpb!ELP{2I!P2lH)Zw9{re2et7{Fw$6L7){1SO~ry{Egr{z~2PEvu}T@KQILrK|sXu ze&A;CPVl#YcY$9F-W@lcS5OQA4-_Cbp6rnV{x6c)YPT%((EwiUgM~ht!0WSbCNWruW(d#@nepYZDf3`L9Hh5TXsN>PkZ3wi3TgNtiyb?<2}q^_G>42ea0>%6wBw3J{cDDPVo9P zQ^1!A7X*3-oPOg5AF`uO0T1}Vc-!TqfY)E>n@Z&JodyB@HMa=}kq%yefgJQ1;PqFt zL7xp?f9Yv^wSMJ7K!2eaEFd4e{;DqM3&5+r)6ho|c(WU-GPHgdLqN~bU=gL@>FUIw z7Ukg23UTUf1^7P%z3yKsAuu2qpm$Z^J>aXs4+QTAe|Es@`b8}S1_c88O+ENz@D1SQ zYRVo>;D;D%*DsnOFw_L>*B0<8;9J2D1K$pQxb02-q5}f*(2qSj!KX4SKLqy~$&3W= z1V2jqK>fl6fzeQa8+;mg5BM?QQ^1cUAFN-bL0}vdkPd!4_zdvpfX@bhuJWOKRk;v2 z4+_W!pANnN{P}8Q7dU`48#(-T2EW}&gH7h`P)j=a67U(^mx((dkj(?)PVl+l?*g9> zz7%``^Ku3+fj|)p2yr*~V(|BXF9p98e0iu>3)sC7s0a;6cOUpl@F5E+C9A-f!|>H1 zo`zip0e?t9sr$j#f`0&fJ^1C|8-iZ9{Xqyc1p{POfNuuB5_}8z3h=D~Z%@02AkZEN zXz?)k4)BkF?*zXJyok4Vi_E-x6ar2Yu=!))UEm)F?*?B9-edc~0pt@9NPz;L1fK@p z3qBqEYVaA_2M-{hfTK-G}UVwlf3V0EGE%^1|>%rH6Z-|Ej$PEx^ zf&w;zZwCJo_!jVf@U8K10J#YQ?NGpG@Ezc{ux+aOGS-~XFw=Ppl+p>MYy~e8?B2rW zFN1e-Unc$x0T&O5SHQc$*Mj$ee-(TR^KzNm27xpd5Mn#{bntcHGr+$FK0DN_WqJn$ zazg{s?F63>z8-u5_}9S~g?L)0-he=HNI2L*PXs?+4xo@6EI88&j9}^_-y6X{5K8w3<9}Oz~|ue!MA`f0RIK}BI#-QGY$9>0>x0k zSKv#*9|2zu{%i0RiLgyQ3V})};2ZE&;9J30gZ~!1KM}U6-$9@j3iu27dhl%_+te6g zj!(@gS~H#7pp*tE<$Lf=;C}$$%zc^o5dtkdAlkvVg8vD8JNRSZJD8Wt)Xxy;WC3~q z6Yv7A)A$9v6a25>U7=nr(P87Oor8ayc*j0xG;6=!HZ(fleUmwD(iY z2>y|%8|XzuQ-Dq+nhrFBXg1J`iRJ^HM6?L#WTK@&Gl^CJy@Y5L&@7^Upi_v}1I;Gd z1oTp(EkLIdZ3jAyXeZDdqRw1Chb|-P26{Qs6rj_IrUT6-nhjLlKO!+UALtCCML=f~ zEd`oKv;yd#h*kl;lBgf(RYdE7<`ZoKdNt7&pw|#>2Re&rC(zkMomcSbdo58n(Cdh% z0KJ}QI?w{5*+6d~nh$gi(ITL8iIxH_Bw7J<9?>eG^NIR_E+ASDw1{XE(1k=>fZj;7 z9q3I&JAp1D>YTx+@6ANrKyM+M0(3FabfCpVvw_}9G#{wCf5d1)W*B65xnWdcF3@nI1wbQ+76XkWS`IWy(XdK)fA0%X z@!sqRS9H@fdlIT^O9n)@1U|}ZO9sTYn~!$8dQhpIKx2qHbNI~YNz@Irm!i=rhI&sq z`g_wN!ei47v1?KrxWw3Ops_^rf%YL<1hg;FQlN1}D}crmtpb`r)DJX~XgyFj(I%k% zh_(Q=h_(amPqY(g5>e-6eA=Ex)D83xL{oqcAes);Lo^%cK%)6T&n8*~bP&-}pvgol zfDR^F1#}2eKhU8>>w%^aZ2~%sXbaHcMB9OmAleBum8cVTKSmOD106*)1?Xs^=|Izn zW&<5VG#}_#qD4T*5iJEeo@fQobBI;}J(s8-=y^ozfu<8}0;+BmG1?K>-!oks?l$UO z%p&L?Z>J(lftC=h0D1?}Dxh}~^#i?&Xg$zUqD?@T5N!c^H_>*W_Ymy_x|FCBcB=0s z>IQlr(G;L%MALzm6U_#?jA%a4`-v6-eSl~w(B(ucfIdjH3g`-=exNIf)&s2|T6`Pt zh#n$Z4)kH7l|UaMS`Bm+(ORI75^Vtb7|~{+j}vVLT1m76=o6By+l2b zT}?C%=u<>9fUY5$3)Dxn0B9A_VxUhGEeHAx(Mq7t60HWhmS`=|=ZH1{T}QMT=<`Hd zfmRdk0Q!QYZdVDPyf0D|>3X6b$kq@|1G<4|2GET}bAi4@v;e4|Xfe=DM9YD0CRzzp z-3(%!W9Polt)-E^6|Y#UeNWd~jlO;IKE8qe(HsYx188&Hx`j%th7z|Dtp)lr(FUM@ zCfW@26{4*`Yl(IMeN|Go>kd9^wow%6cA_50))7qu`Wn#;pgV}>0^Lcp0BAkYVxX@R zEeHAr(Mq6i60HWhi)by-w}>_X-A%L^=pLf2KpTj50NpF8+jS?Oyl+z!={rO{kZmNI z2J~H`89?_D%>}xjXaUeBqQyWD5G@D#9??pm2Z>e#Jw&t?==($)fPO%<8R&;ZTY)wc z?Ew0bq;A(;eDZ!wQKX*`^+5J8(KMi+63qbm8PQyzpA#(r+CsD#=oduGfqqG}66jY% ztAQRNS_|}Rq76Wg5^VsvAHAb`)cT*+Dd? z08K~OMgyn8z~2+i0Qv*bT%bP^Edbh1v>515M9YC5BU%acXQI_Wj}xs0`U}wppuZAr z26}>Mt4-NAs<8#ohn(pNY5JV$3H4%frYEH7@V~MrOrRURo1YA!@z{LNyz{+JI7WEi zcJzsfkf@t{s@sCA+anB6kudr1{7!WjfE!8|>fbPN*s(m{!5!A%H9M^w-;O<&%VU0c zLJg%J=1vLvldyk+;@)l6P;W!TQ134hYS`I>{%?cHpOujPc%A*7=I{J~Lb&S9$}+a= z4TWqeJVT=&T-Km|SYbib>24A~zRtS9 zzIA_bvL1ArudL1*?whyGdOZA8RYi-hR^A{tO#V~9x}(Favfb6Do&?83%ju>F5m{ee zJxZucdd?Q>I(ttN@$zk=csjl&Q=}4Ix#WXvak+duSZrwhDqpCp@RG#Z>vFH7dn5E0 zZFjWJ73!Zd5=HEcuqyGJd^=qHXXC29LS4ize{kmWgCC1w4tb-7eqvYo%X8#dRoM4r zhzyendg#i`EU+t7ixoX|t0K)uU0am-sH-s-f!NK8HdZ&dyZNZcwZ}m7SvNPve3b2! z&HQP^uiIoZ3mije!8gC&qe6qKULE5IH!TZSOZ;!bck_PB<5=i>Y`>L!k@?Mv)BH%l ze|Hys`yL{eBieV(eyjf=)9TaxV2UHWhwtC}tpQ1{XZ5Uggqtym)W4e`4y$z_-m&x7 zE!Hz(z2qvII8sd!`RJS|O@4UbA9BokfXg$k4~N$C<@o|#T;Ojr_zMF@zd6i(*E>42 z@IaS4IvfFi5me9$ej@lb=F7U?(V?9My57;@IQU6WK?nHB;7^8nS$~&1I&_8xq`w3T zI0Zfnyr|)`XbSl75U=X*az_VeNI>~Z)t`_@4EU)bPQ7)3p9VfI=&AlLcXV(E12pJm zP(TuR)kB3oJm9B?O8GuH;C207?&y#b2E=Tq8w$7@d=B_)z~_RW1wK#vGW!><`2QZNY(5llEfg>t z{B__9z+Vr(ki4A!W(pKQpa=@M0sJEHbHEpap9{W3c{Tq{0}3Hf3I)sqzZCp@@a5na zfL|^>Eq|r~MG&Zf0v3W_1^!0xmEdmzzq-c0c-S;x5d^BBfSbXu1%C_pYVeD}*VNdT zl9&b*L%pME+1U{-W=*c42e736!>?GNr^X~5N zTXeumDiryUG`l6HjkCcA_Mq%Q0rYagBkb2C@N#b^^id`}5a=DM<@;pt=6f6yO93Big#DTdz7O+qMM{G} zUlx$R>IFUyd^&iw-y8av06rm9suei{0*Rpk`92f88+jWB;dT9@7y^R=0sW=~ zd@}e_@PolG1wX`iyM9p)fuSZ~zb*%#0=@$LFz~Cu54XLkUsOV11Qf6ud@A@V@FT&m z1wTsrK>eZ`0;8dT8t`f0{ou!d-wJ*#`C$E`76Rj-fI9Hw!PkR72mCJZ=PExtB-8+b z^Fo|@+Xy}#d=q$n>jtwm!*4e5n~_@BHmc+H&_^Bkl29q%*K^eNHoWV;VaEJxE-nyX~{5>%IVem`Aw}g7NfZYp$BcTE5?gQTnK4c-KWE=Q$ z7`{Ek)3D1Ra6BZS)cxQ)z&`-~B>3gvJA+=g{Xqzv3I@om053M^ZIIc)TM0fKd<@NwW*fp;5k=G~(ZNHPJNKL*|d{&Dch;48tW*gkNi_yhz} zp@1jBr-Ao^9}j*t_;l@qM~Y8DU;-4d27Cs1ANWl0Rp7J9%js_xqNgE{0|h(-J{SD6 z;Pb$*1)s0Hn*XK&&p}`|6tE6_0r=;^7lN+_UnD&(f2IL1KwuFR@FMtP@aw^sfUf~x zx&e+9H$Y%16tEF|Irx{rF9+`jU$Fs>6gNR&6%?=;d?k2(>jv|!t2r59rt=mkWi^zt z6?_%=m%*>)zU!?Us(GNxts82<*Fpt;@UMd3%Dh~rwn3nl1-jn4p$>c<3||laHSoJa zy;`PsK%gNsAl**zjo|CSH-Ud0{Gkv}>(m<%XbuS|^(Oek;CF#<0sj{GBSBB~ce!;# zYcN2A?tuc@z&C(z2fr8m@qpL$ce!;#M!6#`S+@^j20S^@LA^2qQ&EQkO ze*`|2yqx}~r~4QJX;8o?;KzeM3_czFr{E_jujapLz-JK1fC4@Tp9#JNd^Y$mz~@L$ z%b#h$mk`K>0=@#D2mT27eDGg`pS=;bsYfAD00n#lz7TvX_#*J%f?u=|wyEDipco4H z3-}W7{MHTT_|#mtWu|i*lu`<%d=Gvp_#eQRbKmvW4a<3;%dHzK!2bjltO9=wd?oX8 znfe(5t68AytsAPq{{q9W1^+Af>QJwi=@SsBu>%GhZ)VAvZR(kAYIx{EBo~x3+tf4L m)HB=ER?}tLum3aK)c@uxv`JyL8~gi~)LCzZf4$|Z3;qjPIC_2n delta 3284 zcmZ`*e^AuN6@TBm``$0UaEAiN-JVGBI1u3gJ&d7~A}1gnP!YA!qM-3kVQK`FB$`+u zlT1irD(GUqDw@Q^s$i!In%;DpAVgwoog`z6wKT1bDU+F|Gl`9@#rBW%?f1J2q$4}? zZojwhvu{89-tK#Qthe#~;U9K-*M>gx0Bd>LFTjR5>ntnEOR5do`|taz_MwRWq9IMg z7-Ko~{GdMd@&&#B8}ID~kbu@KNKbM)8MN>nUEjp#9XJwL>DIr^*cnikl;z|H3QQaT z-KRc^*GS~B#it5LQUXI6BuTo%;aX2``{ub1_UbDgUzC<~CB z5oEA3)$-zHr7#mAj3^_rrEy1O^QN6Gbxn~S4UJ*YcTk&zDr|M2Hg6=v)ALGu!h`wI zd_$$-y(D_l_&B|7sPX*s`52f6`F&%hVTB8k1Ky#s%OwX!ropwEe9|3S>~_%|#;gMC z4CD_+Y0RZl8r-jMS!~q?p{x2xfRi5Y_o&;qS%|^uRIjfm!8Q7xyr*?oIhtebt*wS# zWvf<|LTE@=NCN;NrhIF|%^$2-QeGjJ?uzPwHj6e9lnek9U1{;Edpe0CE9n3-J>c_u zf1jy5!=lBE1X!rCC8FF}C=bv)Yn7BnpR?xTZ=coWE;VX;;6^b%dP}3XtuCpQ-nZtu z(^!bG#pvAxrj!=ja&5vw94*Ef4*HPY70=?zny^+%e_>lntE^Tjf$80h4r@Zcn6pS? zHfWXIsbQ1ZznT2ZPGGw(c`i`<@z`K!b$JPnBDr`p70w*t5wL+xGY)Veyzr|2**xR( zKQm|2TvLvwRZ`J=>J0~Z$W>kP)H$Zp1LQlaySZ2`vXO?l_5615aG|xU-|Nksgh##w zcMCZ87{VKNr8M@r`k!wKwl@VO6a$J8#e^cGn4?WW%j*30lD!rbD~b)pjuL~SL`(L@ zUR>V3Xfzfj4#k1uM2Sb45nVJo6SdjF1e8RSBor4)GKxEz9rRpW9!>4=R5jBl=GVJh zo(xBXuT%QdO>7^^e7u9~9MOF_E?M#7F;_cAh}Tk~_|Zo|2Xg#or2tPC(6O8Xvr+_L z2D&~_WLAQJ2ZweCf}Yds;_{}sNq_JKjUzoC>PWw)6wX<9xNH^14em;6p_BONqZ#&; zXL$-EdaywBUp@NV=~=&_o9DehOWcAh+~0%z#4z@44db?j%cNinKSW($?*rXC$$^;g!bR&OmxDfWrGhhYAvpuA)VqlIw4g{ zT`bbBXRE$e%by{cQRmy_s8RiT7dhx<<^<~(&*=uc^w~gL&2^^3PuEfPh23QMThJq% z*|31lw$yn#l+Jq(_O9yhCW~!(+<_J9rV)~*9_c0>6|a5rVyAs4|2mPDJ9N}<; z@I$ZPqyH~I7#6*q{4m#|@sb%YmQ%u+T|D{O!Ol@v9{Evao3OUnaeh*EMa+)x&uq9K}$3ZTN)za)Mbs~ z$hOplfxLo=kH2Q@R3B_QZ@S{B=5~{Z-3ga;(i{2CM9DB~)lqGAkW#xy(@CJdYb6WJ zESVcXeX)b2yk$HRlIKWAM!HwekXA*61j%kg0q8_T(WgW7MBb7py1CDo-?_l-V8r^h>u9Q<%XE!lBr#=eUH%-qkF`mET zq1*De7*C$^&>Q*Hrk|ECeOb*uM*K!BYU-L}q$V-+E5T(N zGY}xrA3gW>aAR4aYqQo@zVuI5>J*uzBbyv`=fGLg`es0lTSZc zt5rNN9ILIW1-35&V^4m%8623$T@YVwJPhyPXcqi;_|+lsC*UwXT5ewZE&K(C-LSM~ zEFA5l;uG!%5oxRA-NnN1UOLxr_05u|~hJq|vqoW!d5+ zCgUKkF!6xUF2}^a0Df6p@v)j^BYTU*F=J@f_CFQ)S6O*U2*g8`4cr%3gxe1iC5!e?||lF9sona59`T@1Lcq? zl?`A8w9+p!=9H)tnKp^bBwr?a_iwW(7lr=g&0V~8G*u~wIMp1y?CC8Eqi%04&wg<9*b&G2MN&C7B@Z$I%#d-6={z8 SpKg+(r0QNyWnR)4NdE)R46Y^s diff --git a/sample_data/hybrid_height.nc b/sample_data/hybrid_height.nc index a35b5b13435cee31f7b39d32a22e0b2a08a90956..740980bb128abb9f6aa21511ea0e975ecb4c6b1b 100644 GIT binary patch delta 8060 zcmbsuX?RrC^}d;x$=qQ`X0lB-$UqCkjQ>u!##+R1hjupvB?>f(kAzme9`@p%i+z^T<3C{n7hN^PV`nzw;|Rt`7c6~do?vp8k;x8l=kYBgOfWxdrZAKjX6@1E^os8 z>+k>pXsr@`n+Yil5}k-t`Z6ozWB|vnS(g@rG6bcCzD)EPQ;?oDXxN&zm>fpPTilg% z!}$E^9p`4>NKv8~iAhnEa6yc2$wuC&#&M|w85cVBGKYiV-{>!>H1X2qv2%H}?dGi5 zbVhGRA;`)^oj7_zK4PnjBd7J}}q5n5bHejdPC+VPmF+9&2gsX~;cp-cvZ zo)<4UBNyHkiawfZR`0;WLGT4VpVxZ#%UsHaD*u% zR!Hh=a5_1ehm-LAWw+eFX9#j44?i&Ej9-jiDTl!(1xf{cwNxk<5uxMVwp=9psCj6^ z)xO$Ou3=ibIiAY1_~_d|lN_{}mSQDYd(o0Xe?PUj;tI?HBh7)DuV~3ySg1=W&{9G( z;*zDJx=9?+6lk@f8KwfWwkb45`MdXWi&PA6pi+Xa86hgNHG zX??MS%@&!|31u;!X$p9?mvBrx2;A2|MrS^U`OLfji)Y*8r%uVwS5UYfTfj>^>oun> z#aTa8jQ6Y;Z-qjWoNXt9p(YM(4~VS>S~G`cVzDy-RXMc&*KcceUTDNs7ALX}bfH6X7;BL3sSMM}#9QzXuPE95Mu6m3*Fy2ZkKT_1H^@nuq)IzB$wHA=C5lec4yHMJkMySZi;(Al9~s7qUi* z=xO?=%%03r*#6lQ=_<(P8r(C)jZi2e*GjOv1%*-7$f58m!V?M##*WKHkcb(=zgXxr zdAE3j^}_r46D@N@^7Z=WH28dt2<@>|E##5V$|50icNz&bFAh3sTD1JWA$ zjO6qojGl!TZ>qmp&^r)My%cuPvuhiw`PI>+uKY=*8|rGwltEFLt_a&~5@K-?I&Pme@KruK`c(szal}wn`EJ6tV+z-+mVEz9Gy} zEyD56)!_I^*UP zBtwy?08Z;aHzj;yK}ktIijvl)P~~4|bjlG0jbNs*s@tMegdUEJT@e|p(8oj{GksM0 zgwV%ApHTW(_eI9y1O3R1Vw}A%tq8AH_H~Q>5S=DcW5L$$XoP**#rx~bXg`=<2k16<)Q>HbQZ37ulXpm$DZ zA19tV>~E=-E(D%SuF+oP<7s$W1Q?couHfJ{lbm~Tt*+sC%Ali$6Uc#LoM~s0 zIKDb+;|4gl24`6>Ca-A)os%2Ic$9mNLDkgTb!~5WYyW9-qX74GQ?2m^YUF{PL8R)m zGFJ@E#POu25Xbb)cv>)}^?KWc<%|qe_kO^=f8Xn~u-iQ7#D9d~-Rsbxs@$nX1t{cCBeo)RnGWcMQm_6e61_)1`=ZU4D*K{Z z{#o$xeS^?&y0G5pgK3SKD4MD$8I9jmKO8wIGzI0vK0FGIrnKp3z-u0#)EVvTyD5gKOEgf|v0 zU8wNiXH)31ElPC3*u3D%A?q*Y3&>ZcxQ|jz>@#s%3zS}p3W|uBSLgKw$^yOx^tZNo zUWHFT8!0NNr@M$L>RUtDaNUwD}=G!~Qq}EnG7KPnY0rT=+{Mf?W$~te!KkyqW+Uw$H}l zd{Ei50x34G^H7(hDBN1zY%}~KSl@TB}S`^E@jd<^s2XTsa!vmi?8Df{Xj;j z!$E7rDD5ALMf*T;XOH_yKa*n`ByuQdyQEOb5!)fj_YMVUa0P` z5E4U-pc)Tq2Mlg-SdXIE9?3#Q7Oewe_TMhKD}C zoapG*2o1U%v7J=R!%Ka7oTCF+)1#Q^B0Y-8J_HD`!TE1}K91EXQ1~RjE@pz`xP1&2 zlAVl49sj06C=!MmS!4bU(K#Z?lD$ov3y@q3P!LS3Sif=Uq7t&l5p z(v~-69SN6P@b0#;wTKr~R`eJfGLW5`=w0VFw*&EG6(wYWW5pY}cc3whqJW;n=Ww=QyE{ci z><6JrZOu1m52J94I(K<5=he?0J(klFVx~gU?!opJNeZkv=b(v$1P_uic=C?R30i%S z1Y7HlonfLu+ug&~tPJivbB+6n%$1~PEg8j%&t3I>1nmQe)a2qE+WYp;!!|M4;rVww&1*5pnVhiCBhBTc zogRY8_!=D54P>Aj#qOQc$h}i=hL)?wZjCF!3f+w~Q(RYGEACn9rIA8fVGHconhbad zjX<-=fH$VEhYWa_M~?TRWd;Ic`g+8GH>R&g4R~YvdYp!Vj`u99fnKJsC-pv0U#-xU z0vglTdIR2=zJ6!G8`IYo1O5l;Ynvnh1^?ga>lad^KbgLUZM|*dhV(T846Qd6iOWuoVsC!xIY9ag{q>q=~gFKV^2l;m*y{Z)LcEFhw z*n5gw?tGeF$(3+hd}%SIVZNZBk;Gxqt@neEhVG6BK{r3Sxf`1s{N?oMSRN23UpEsu XWZ8hrWSYY^j0&Eq?rW;Yp7s9&HYRS* delta 5160 zcmahN3v^V)^}g&&vcqo3CSmj2kZi(k*gzl~2sYLTA)if9G=YF%3tf|JK8qoV`Je%6 z_*7G&^ngI2_=uuOu;`iVY8?_3(^XARE8yCIj%+Bx5 zoqO-x-`n#~lVu!maBIH+!=r20R8^KOjL7Z$%2C~cP4lq7yY;hP=Jl>C)wL}=K3raj`XWYacF_)08+#t3cCnwrJiZ|(`Q|gu<)kA zG>waC{Py4miZX#8jZhRNmJbkO*_CFJPaifI=UU(!kKp3B6{0_+_r`eS3Jl5xoy#I| z7GN#lMqrgK|0S6i$(c6L5tA!Hu)jP}U~l0a>+dCr6W#CDyi!$Kd84FjS}N|d8xV!| zN4hidI%r&R>Z@R^;a`lwn!u&MExtO6V~NWam9UN`GNUMd^&zvSoTd)H&)ZVlv{9D#}sP(>T)c*8XMbhLBQXDh0%(bu6MJVn~8I;h>2 z%$^<~ZgM(Sn94oupFWZkEtow|r$%}lRmDQA4J(UFwynnl4vDy(Q2bHWxQAn9J!!dbu680Oj0Zy zSBc}}6lwpGxHOIgjx$fQPNX$<*E)v2Jf&M!0AJplq}`Q#l~XBvx0maL!`QpHZaBzw zE4N96G|_!;ed1Nt)5ok!-ZI2GtuV$ZcFqYcHWkm4yN)rZq@>QPtY?o|w+{Zz)MRGf z{Ujm^p_4*JsJU@x5db8%*?FjFw+ z0Yp-Jft7w5ZA;`d4De;Oz}cj8$KGiqF4LTCiBoXFvg+*5Hntb^wHGMZBCr{-sn{a1 zMPZA^W*llSAjc=OC+o>w#)2ak@8CR8le1sN;1fx|bY7oofAZSA_-jwTs>2U4lpROd8|NU}OhLHhx-0-rt$vJx!mTd|A1f4!B(GQAXk^ zbtS5aCFUwha&qK~)vzwVl-`zeYijB6Q!&9exVsU`Y4%iFhv(Mx`hfAnO@ zqRKLe;t`7nW}oY9qri;S=A1?4z5jmgsX;6lVYY6P=RC7SpiB}1fb7p??s*q}|HKnz zMb^0)>}VZ1Wb@S2l*1N1cq=ggpYz9ir0}!LpU}Ib^NaXwQ(cbC2^vAq!pdsR*dsv8 zLgE&gWLE31f3AXkyn>{&p-m*2z0*Y!*sR-0LhbVl1aAP>JveW!CY+HA7l%}QuG~to zYXXy?3igaQZ_=;`$3g}<)(vB^?i)X#3Kgv8Dsc9?5K%BPS|^aDpUB=`^Zpkli{^V! z5)%(yQ!#X5H~*mB0-D(9?mXc29V__>XyF7gIK1%wgGKOHPTT{)WA}@9CId|{yo z4j0Q?UZ@f2VH%nwZp?V^(H@wK1DGuM@@`(g37Yu=Qsvh-m}gUXoD;!tqH*lsJpi}x zttZ9PCtmtw7_Q+9b^tSW8#3cYN&mZd4vv_vQ<(n&V(!&h47n=gnh>r+K5z0PfK!dvagSd%i;vW#q1DzrC@xP=hjbqjeS zX88w565H(~u86vo>_s2BlzJ}e++V_7@6V&NyES&W}d(o2uW^f;tdPQ0qR zJw>KOZe#o%o$&I&PV&87H1mkTR|qqwF$%A$?-VkVQDl|eLdsDBdooiT%Y>Pv2}x<%CeF`@ugwy z|I=HK({0LBdQyp3V%W<&(6e6lW4b*Zsj{OxNFE#NBR8-I8;K=DRs`BX{3)suSL$>6 zkEk@J?j%>(JqO6L(jUQw&Hj+n@a>Awxu=ETi<=ZJ$HaXv$FMh=$zomoHg3iUVbsqv zda(KvHTB2XMs1~<4&73RIZebCUp!WcZd77HO5YnfT8Sjf$PIebDPnhSCpAC)=#`FB z;Uk~&F9mWra}(p!A6+cu6VD9rz(X5ThQ3)&o>6nPXOWdU7i3@m_OBIHFS0rpm(t2Xjlqs)iPCbA!LBx$BBp14Wg!R$TD5m(!JxsHn0K zi3urXWnrMG@_FI|`?Gkht5JDg3a9cA#$mj`YHMxMDq>mPD5{C|tR^? zg0XaVcQn@4H+XAXM8WQ8O_1VgsN!d+@pfW2Bwi&PtnDNv2CoY2L?=0xQ#hD?L^FbO z`nQa-7T~zJ_c8`$t-5npKUB`D_Y><6+jGm+?U`4B!Smowx4ThGr2L`7Cek6@ zLbjB%C%Q;FJ>!@_FPbK>i(MqKcZ%+d3i#_fp4akNg?`AW)^Dg0sZqIHu0rF!1PQd`n{%5|>XQ z{-lT%A1Hpk?rmWg2FR39X6&39YAeV^cH%bg|4sNSexPDF+bo_WOv=jArICV95*5R$ zSLc7kXdk(|DG;)S(%{FQED-;6{PXFp2f{OqUv3VL6oF6ZqnWlO>b4oh9}CCym{h_} ztwC*M^pWyFjpM{4+WK)W?ywXNB1eOp_H&)-_x u)gUZ*HeSxfY;5ajsITqp^0qg0_}hH;?`J`_`!GpU4TiRO_UFT7*8c&@S-?5~ diff --git a/sample_data/ostia_monthly.nc b/sample_data/ostia_monthly.nc index f405edaaf04fee0547e03f5a356f7b71e242ed3a..f217936937280476c72326e8f958a886021c3850 100644 GIT binary patch delta 10128 zcmai44^UNA8b9~F3%-Y7Tm&u$P(TntLH^2`Yb`1wmZPRr?5^v-o+YKNJ7JsUrju4S zS)*Os>ulE2u1(gI=}19p=FXs^yKUJ}p`?<_D3T==+t_Nk=eyrI_i;aPx--Xj@9+1! z=bZb!`_6ZMAO7R)J>5r_Zc3fmot4bsD0Kd1?Bi3V-_5G3nB&WPdrtFheOrTXgro7t zNNLR=JZZ|s&k~(&zISt3eFt0T^pM~8Za8O}77tRNrfFq(54$E0t4PaOGbSKW#$KA1 zkjiv%>5pdFZQQ|k<{1;)*$4{E-1ej@VcLy5+2Ilq{frJhwj5 z2B%}El5{`9H(3Q{O3L1qfhPxk#;ev3zNvA6IA>;B;eT;r2P3VS{1`w%lBgkq} zaZdjmjmqpiYluHRk@S^}0cNfd*VX|&(Q4I_tX*2d;N^3^|FClI+^MXo&3ZZ_ z*}P+Z{7ti_RNabSO#LYAZ>F_BVtNHT^8I%Oi2{L%6w3F)B-b$4yMJ?H+Q0z_TGJD6 z4Mv^9ui!cB#l&^iaesakhTx~&E!p;>lRVpBp3Cg}n)xvM)sy_@jGjUGc#OsTsWUb_ zy|&M=pPBHrg)rrnN4zV0xzb=HQ6L4}2NLd8K^AQs%DxxY>I60KcjO^V2lT9JoitX*Z#YUhbohZbw}yrMAt z^jXKXG-E{?#(M24MEwrhv(=00sExzeDQyI5pZ21WsExpQD(xWDKI%mSsJ-k(gMdc$ z`~lxo@eD7EcfeTjrMkwM4#pls5zeI=CmMS^?L}j-$2;*G5J4(3xd=*(vEsObu+tl` z(_8Q@P@TqNr_El}K~rO<){rtd*+ag+I3zu1+`bbXew$S#Pp&@8;06TUUWEWYlnD7OGE8hFPe_p zd9b1~ruxl5?I&Jz1ZrnFnZng(qPE40W}$Yn;pulIYWI23Y}Ag!WmwH&6l!;R(Hzv~ z;4fKeb5Z-U7afh-6eq2^^BsfQEnYMawL=m-^UX)?lU{TzYBgAZ!!avB?ZaNQ5VgHt zvkC1UbF|sv&6bZ{YwC#+|gYGiA4;_8X8BFdwx_a^E zDikv7O<$4z`i9@GwSU#l%ht5!Ag(f<`g}s)IKi1m5at0tIWdoRI*+gGkNw3_MSu#H z0kKi>V;|v-%HOwcx-f9Kn0@DAKE+r0nZ4~Wzb80x3{Ixojy}R?R6bBMq?X2_TW&eR zvAFl0rgA41=RU>hX}dKxOH7Bc_J@y+D7Rla!bA4%vwZO_^`A%u4#Ps^NXg779m=sXoExxh!zwrO0>bEMT-_A+7Quh z5G__TL$o;2hKd$1T7qbaqCwwDuXdUIWq!RR3|Y+#*2`T@!2V(%H-kIe%grZw;u^L~ zUCU>d+YfogSQFtf`MESZNO(+mVZMd%YQKYu`xyI(2-L<$buq8#R?)w3E%6SJw*6a!mlO#Fv4#j{O~?JjL$X^L0TUH zoNgt2I^pXHpF#K?gdY*+o$=W&BFGF2oU_w@ z%beq?ItVY5n#*?)UfNB@gfup=g zpOW(D`|Y;(c$q(a46}^`JfWmMC9_YiN?P1?N^Xn9O@2*HNiFkinLw>?0?W-uhPJN) zC2ry_$#?ZO>mlDbk=s(QO<8Q)=lJj_!`Whp1lnf5e~!-^(sfZST`xZU<8b?~W?q*1 zxSWCdMXAq{FL#`M_3>j|9~vG>Rih?Zi_F?&RT_fF#iPZz^dp=#mLK4|wCipP91j`X zZ~-S13!4J4ppn6?E(99`>x>pOb|^>=0jAtEng|M}+%#GU3Z~pN+6W4!+%&of3Z~pN zE)&$*z8HaSIa4_YFy&@2NuiL`(~KN-DMjGPO=CPk!IPUt6+yw1o531|inN_<9n>O> zr3#UKf+aVNhX@McG>r`e1xs!kTL}u5+%$F&6fC(2DnY@Li=Yw|EV&3OLBW!Xpb``; zxdn!IF#E5)>@Ch%G_ElAA^iLBW!n#!`ZUB^U7} zC|Gh4UxI=q7x5)1SaK0xf`TO%@g*o&auHvGf+ZL6B`8>O5nqCWB^U7}C|Gh4U-?^C z+c~i0BEAF#OD^I|P_X18z61qJF5*j2u;e1X1O-bj;!9Al;h}Ge~f+j-~3#hS_oBV^`F=WKIGq)^vBimi%R`l@3ruHB2)cMOB^^;fpIVM_gvxj z9f$bQvc?}`V+Iyh{LEi2Qa$MnlgL`>!57&LC;0E9Rr-PnI_b;KKFV)~edn8x?MWYa zSUa_ar+O3VXwOJ)B2ACOEwYix>lrJllqh`{Oz!@JH<^VbDbr2=caZX&-IPd>?@~uu^#H&duL*W?%)0G|87vCRxK$O*GM{24;=zjrlqBdbK@8D^^t{(f@4xZ!FzqOk>c)nk}xME-O z+I15?Rsnji;Oj1{09V~htODp)C5;nMf98@`j71aP4X{+i6JEtfb&yJU_1MUZDnJ$y zr~s%A@(8aYt@u*HyWgf1OjJBp0V>p`q-hmk2H~UCIP$Za@U#jrPw`F_Af&G3-Lwj@ zjPNRLs6p!pPpbf{3E!^@u$~BL6<{ObX%%1#;b|3M8{zv^0k#tXtpdD3cv=N$AiQ_& zW^6y<`&9v&iGWrC{zZ6N1!yHatpa>b_3Ujk`AH?Ppbg&gda)nPbIuu;)F`N&@3VtMLNhMymVHss+90@v2^)~gzr}cs2~Db z1(-p2$P0asYQobhz&yhDs{(|GfK~yP5uR27>IhG(0ILb#uL`i92xt{xBjITkU<=`C z6=0j>Rn$lo;JP^qP3d-dcOO&Q;6kn`ZB~#R1WjoxLD7_U5)@5oFG10i>ObM{LoL2& zO5+KNrZkJ7Xi7^Jil($eUCN1~DXk_bn$i$K(UjIH)HS8+6(ajYQ@VwqXiB#e6isOZ zLD7^p6BJEpD?!neb`lg#X)i(1lKK}AyT1rqfr4G^HVe zqA9H-D4Nps1VvN2g`j9kw-XdiX#+vglr|F-O=&AZ(Uf))6isO_LD7`z3^Sd2Y0;F% z6BJEp7D3UJmJ$?AX$3*ilvWcIO=*arXiDn{il%fuLD7_MAt;*C?F2Giu&}VKv>ECK8Hf)Au^KZ{1T5$!rB1Iu%93=n zskYN}EE}_jnw14nvt?pdYC@J-lh0vPEG>(Iu%;K;cka2@aHW~f%y-Z0JKy7XzVAEt zbVtJn-8Vh%nNi$_5o>OOLn5jnab$*2k@n0Og1?yC6k;iSK{s)-&r=ahu4yTSVA6C6lsQr$ZBOi zIU!z8Y-%#uE9K@cLWPHKU67l?Tfx`>DEM^O(TPLLBKOPbLXg~%rjc?nQ1GqB%|Xv| zmo1PSh@Rvtl_pr#AM-h?Hr-apFHS$WuOwE#%eNjj@j7wE&67%J%_>24JIEt8U6a24 z>q&-z%f*!a{T}BoC+)T=zNu>t6e0#3I~!0K147KG!ShitZwquyC4RM1NFz(tJovV% z8ZMC-J3@BeA?v==$Z=H@O33eO9!@h69~Q+^LW&)erJ`M!taA8(_gzOG4wzSX(FMGM zGtXh|CXixPBb83qIeum!YnJ|Qa*uN!liha2c}*-`_Z3#$7sM@@s!Vg#nbzG;)+!Ih zafD*pIYN8LUqYgm_r~}`|DyM{8LL$^sS7t=+`vK_nkfZF4U_hsVJnfB+nf(I2}lOA0VzOsAQi{~P(08`padWdXw=509Q-_K zO201M3Xm(AcGh?_EcY(TW#-2~Zi+PAhdpF;#u|6NKhHldFF%lyKeo^xEG)<;dox<~ z(3H-mq2reR-HUS_9G@9mG1)iflccVSX)g}Wt?((1n~}F<%O!sM3vVvl_Ia>~^j@+} zW}GC|o=?-eaIM1vP9OvJO!?G(x6h;}_v4`S;K!IIAHXB^Z^V<2A){g-lw#n^Fi!GK zmLu!BaElOG@g1&_#$S1L*;R1uu;+&D`pn|*xD5{EU>HD)9Zr`37+)gfV~A9I*jQFt zRaS{qGS@!J={9e!ZU=q9qgU59f5y)vw6ncqPPFvq&S>e&{S#*xTn=UfE28z8&3mG` zrj&=4md*l8ky|87?MbxjAa)DO=uZc+mNL#MY#Iz>^B4U1b~78vfE3WD%T<^E_ss;Q zV;Z8D_2G^Q`si4%k^EO2;bMfB5UOUL|FG(YFav&>h3^-_}keJ;F))i(P2eq5j)uWManoM_#T^W|?whgU;H$ER4;h8ZQH z&E5RmN<#^Ju7yH9M~*-G#f=-ChpOi?viM34oj(sB$8CMi4ML>g7ra%;9vbHVezk`JJj=}N8oKCbyf9>8%k=Ua z*<#+v<)%@5HO*WtwM<37X%)|q+4eCKL?e~UE#vpOmuDE!Wn-DdUW@w7uo=Tn za+qo1N27>+>` zI${?B>?;J?Yt!Rge9?Rz5UQeaik_6qWm&JKDZ~mRQ>j(A&`vr>(F5|v&_~rExCR76 z0Hs?0XP^|VKeg2_$Tl{FCNFNhZNXVY%T+yHXrOmOr+07R{tL2MOpsiAd;KDK8EK2E z<0NIbwQpIW-a3DA(Z`4f+l^}PE$-ob2u*9h)5pZ|z};kGT4S6>D5j7U zWoWeR49-$BvbGz=EYkZYeoKg)_UJa*ZT+no{B=LPG3lvPy(ruA<$wHN5x0~{7xrL< TKAoz+Fy!+bimXi2Ys7y5#a)R% diff --git a/sample_data/rotated_pole.nc b/sample_data/rotated_pole.nc index 67b549d8dc3ef9f4fd59900e396e25f93229494f..8aed4667af6297090a297bd220be6ecd6a40b6dd 100644 GIT binary patch delta 2741 zcmZ`*YfPI}7(Um zUS{LKRU4R}Z5-+~P{3}-_O%k)OH+d)@4TNS-lK({aT`gJ6N_X?l3YxEhdt0D)p)}S zLDyfRo*$=O1PHDeLa-#az*kb6`8PJvnW~kuT!;5ku9ZJaRvCL1R{owEK!gty+q-mSJ;cCO7HJ zhM0#6unxPxkyTa^V=zHF?1MiIJ&70#5F4!IeY4g^nWN5 zf5W(&zcD7oreyfV$Ka-OG{rd^Q;s%0NSMo87gwK*Pg{cs*G`kPzF z=t)I%p6A$d2e^AO&xy<()(@2J8ijMtJ<5@fe_nyK>V&kC^w!CgW~XBFG}OU&O0!hA zd|Nqb^A4qYdL-mR{C#sfr@7f(KGsNG%(00i4Rx321ym;JOjxfq5Y!XVqY6QD0;&?U zGpg4b3A!qvB?KLA(rdLtrA4+Tpk{)eY1V5?32I_*jjwG7LC*`Qg`mL}y$35nKN8R~ zg7&xSwKjrY6Hq%r$J+GTa)Ras)IrdhcD=TOpq@Q?RD1Q>KI8TdliRNaCPNANmpEJweXgVzZoCVp8~wF^ zuv$8^c;r^mVS8;fBXg`k9Du{&C#1b#ZQR}E)&jsOH##wuIi7xYdSogynwiQxeLORQ z7bq47r}s(&UH~tPT5g5CNh8TENny7Dau@N|%1X~{$F(zs+pxyAf$uFm8Hp;n*^9T} zWMq2zY~_fhXv&k+%;xW<`^c*XFPs!b=s$sHMTXM|u>kmZcAIYUiokOn6J(J|5U^`{ z(}6|c*_E`(+XbHQBo{aZp6?Xpvj}+vf&L;ws7l~TG3MSiTP^TA3OL;K3%sR}&uX2Z zAg~q+w01<`%LG0q@HT;O7kGOCpT+G?K~P>Gz|DOE?-2Mdfv*twZhH(tr4f)(aN6{6@xLRn`@y7DF%0upiuvq z%M3^ROFTRyg~6=C5*aT#NFL^)tYTZD3}s@2kZI+!Q$mnPzLa>}pi?d~H3%u7R#l(( z$^Pz0YGDwDkYP55pwmqRO;a*|HCVaB%-yC#Hn^%PtGH80&Pnv867n@H2cx8_a4h3d z{2rFV9o!CrQPAu~MCE*>W07tYfybD`Xq$K^gpB4?UL75>+Tn|_Q)hoQr~Z4#SQ$xp z1uigtZ9>LG5eTXFz0tr$Q{^6&e8vp<&-pgxY6pien~=40$_HC@s|6A(n<+gyjoai# zn~-6=*eDFFnT#+~bU=Ud`s#dY3E}tACaRi-hT;Xr=FmHU5{zoMnXvG7v85oUY}PBE z&aV82z?&@+XyJk2hNz=rv`>D>s>Hn%i4B^5zP_GZ_l@1e<>BR&B1Va5v2j^R@}58- z*y|tg4*3RAmKY)ZK+^l32?TvT-oc?$=Y2tcpf_I