From d2d6ec6afb630636db4b2fed4f0918e17320c490 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Thu, 12 Aug 2021 20:58:36 -0400
Subject: [PATCH 01/17] Updates for MOM vertical and horizontal grids

 - Added MOM6 35-level z-coordinate
 - Renamed to mom6_z_coord
 - MOM6 rho2 coordinate
 - Moved construct_rect_grid() to horizontal submodule
 - Added horizontal submodule
 - Added tripolar grid generator
 - Added MOM6 35-level z-coordinate
 - Renamed to mom6_z_coord
---
 mdtf_test_data/resources/ocean_static_5deg.nc | Bin 0 -> 233984 bytes
 mdtf_test_data/synthetic/__init__.py          |   1 +
 .../synthetic/horizontal/__init__.py          |   2 +
 .../horizontal/construct_rect_grid.py         | 145 ++++++++++++++++++
 .../horizontal/construct_tripolar_grid.py     |  87 +++++++++++
 mdtf_test_data/synthetic/synthetic_data.py    |   2 +-
 mdtf_test_data/synthetic/vertical/__init__.py |   2 +
 .../vertical/gfdl_plev19_vertical_coord.py    |   1 -
 .../synthetic/vertical/gfdl_vertical_coord.py |   1 -
 .../synthetic/vertical/mom6_rho2_coord.py     | 126 +++++++++++++++
 .../synthetic/vertical/mom6_z_coord.py        | 126 +++++++++++++++
 .../synthetic/vertical/ncar_hybrid_coord.py   |   1 -
 mdtf_test_data/tests/test_synthetic_data.py   |  17 ++
 mdtf_test_data/util/__init__.py               |   4 +-
 mdtf_test_data/util/rectilinear.py            | 136 +---------------
 setup.cfg                                     |   1 +
 16 files changed, 512 insertions(+), 140 deletions(-)
 create mode 100644 mdtf_test_data/resources/ocean_static_5deg.nc
 create mode 100644 mdtf_test_data/synthetic/horizontal/__init__.py
 create mode 100755 mdtf_test_data/synthetic/horizontal/construct_rect_grid.py
 create mode 100755 mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py
 create mode 100755 mdtf_test_data/synthetic/vertical/mom6_rho2_coord.py
 create mode 100755 mdtf_test_data/synthetic/vertical/mom6_z_coord.py

diff --git a/mdtf_test_data/resources/ocean_static_5deg.nc b/mdtf_test_data/resources/ocean_static_5deg.nc
new file mode 100644
index 0000000000000000000000000000000000000000..3d05dc8eecd2615a23001bac2f4b0a4224ebe500
GIT binary patch
literal 233984
zcmeFa2VB%j*Z2>L4IB2Z*ijJ+ii%gw4to?4rHZH^AVqoyv4Zp}h=3FoRP0@`w~$z{
z_uhN&<=XFmW|R2byU$&?!t=h5@BcUZ@pyJ7GiT2EW;U5jCc#KPw33>dnmQIuELE^n
z3UkN*tEm~2?^A^D9q@OJU=JTxvH^c8g#{?|aPobM@Vz7cUJHu`mg;UE-mbn*it@yM
zT|Z}6Ctru4V5eXYXFXqMoTXMJ*dy3mG)T%c8s_ac(aBp|W(wL6=I<JaoAV7OCDqhw
zObYaHaR{E`|EKcR0$nGCcsm9DX*+IIQCyAKt30ZQjHGH}X@o_Ji^NJk#MjDLYIytk
zPI3^P#c5KvnY*^XpNDU-wy&R$hc9-Dl<^1-aS`oT3Gwv^4#EZThm^17;yNkN)iub$
z)hUR$T*uidFxWN71IO$X<{=a$<<xEQfA~`sS7WfopV5EjOnf6!nr}`2%{L5HbWPc>
z+LB+seu2Tlu2k*YSmKVtH+60DMNMr?JpNDG@5C=ftXNmU@^{>%vdVt6c5#{#q#fkp
z>+GuCvuDp<-4*}vS3|`=Y^@Y#){4z)wVj=CpK@^uRO$y+O-7pzmyJ>ssxFnmDdOv&
zw_)aL)fJ^G&-qBUi+z~HMfyQ$%$o)L&-Vw?5Aer!1^!R^O%dP7@2NE3n*N(_|M~u)
zgMGlCaj)=yRVTUn;h%#dL(!)d-vl33^^4Y~((gyx&rLfx(8<{~P+P~ki$GQQeJJ-C
z<=<3oACDm1gC;qII(dhPzeBQ5)#Z)Ww+h>m%;4{u4nsY>y={IJP~Kk4+11<I!N)b&
z-Oq)zO8l=$j%$X3Cpbxdn-)alk6lyzpCVo-;Xw3H<85Us?gQ1~swv0&cfZyDLd<JP
z%+-Z`>_86GBxceFrQZl)`ggzmHh<a2&VGTuWFMoDKksDanE%N>CVfiXac$T$td;lg
zpSS&!eOntJ(_I7o#U6kS;*TOeUq4^crLeuOeH6KB-+yp$3jQe$t^d65{~klR_x+#a
z{oiAuu49UsNPQwiWcKIdrm8;C`7eE9svkL(hPrzDIeP?8(bh5gzt%UpXes)Jx@CmK
z1NV)e;nC^OeWP?7|5y5m#AI5U#Kh0nArya?jvM%^;{2tmkLdoTkNDtmpNofbM+yDY
zLFE6yF{DcY30h6PkZ8s42KSSnVbbl-{iJLp|7ZFN;eolFVDSG{odR8*9DeS+s%;e#
z=!ORm+N7@1xs{yS$vIuQe^m16jdRu1{-@5R?T9aiK_NccQ`{9QOi)V0ia(7j)PKK)
z{i*5V>L2Xxptv$3{wa?ws^Pr9JI}V4lX{Ba-Cea6)wP|1r9qYAvPChVQs$uULt6QJ
z4i!1CbX3IegRk!h0ENZ`6j~Vp9>WXcFQ=f8Kv!Ydr--|j!iHk6Dtbh1e5FW+`YuY>
zWi5hmmvQs<^9ytk=uN~YlgSQvwBq4HV#9IhM366Pm5guU|0_>+Rdg_AwyKUroqB3}
zxM_Riwsp08k*})p{#!rO5XTVzQ~gYNocd2$;#kG;Y^4B4%%lotq8QNN@6s{4tg&4;
zm0bRQfTP^^N@G)zzV~lT{_*&(68D1zIf3Jk)JMu<p|6rf=_mAGn2C4}uk<Vfb5QgT
za<al7seh2~wXmrA-Sk(<r0kQ|ugs-@1g=)OQ6c`nx{+U?ho84%Hb>%8R$tMP>S+5r
z1>)&2eDOh!)qne=QRYxNsCzHc`rjQ-{yV>!%0uMXpBm&?O^)M!#Jhj}Zpio2`ba}F
zG6$$V43F=+8{z8OGHeIeQI5N`^L0hy{-@_@Rp*`>igOnUo<#cRb4_KBFd>WnPYF;0
zlmI0_2~Yx*03|>PPy&<yB|r&K0+awHKnYL+lmI0_2~Yx*03|>PPy&<yB|r&K0+awH
zKnYL+lmI0_2~Yx*03|>PPy&<yB|r&K0+awHKnYL+lmI0_2~Yx*03|>PPy&<yB|r&K
z0{=|{Y=Etzc{V_J#U$&$T}b=8i<-PnPEC#VBd^KA+V_>H%?0hVLDX7;_VyRF*C|o=
z6triWpgmrR+Fa1?8$_)o=t=&9b~`2No`QBw6ST`KQJV|ed4s661U=DT&`zgB-BZwx
zX@YjRCFltSg0^=Qw4Jt~$K4Y2*aAV@ItqG>wxCDf60~)JpsgGQZK*A2^IL*8V+3u=
z3fe?W4-?ZMrjeLt#Pq-jHS%E%MBPi&T}Up@Z!KyKK{JvT^BJ)mBi3U?JB(<T5!+$J
zb{WwRM)Zpj{bWSH88Hq<jEfQDWTdz$53(cFMIP#m&@a>(iAz9)Vn?w7Vn6i1C8Yhu
z{^-vLX+JUFPfYt32>HVMO4V2&M<MOaidtLHUJ-)!yd`Qz&>jV%Rwy;r-BCzSVnwYj
zXtxMKyWSEtBWRZbQ7e=h>+C3`C$gf}7PM1@pdD|Cnh~^vs3$0rYOEbA<d2UK^f*S)
z;>&y4F|3fbi4gQ?M$py@rN&yaLfRri(B_PwM=6vVJ0e0zlQ%o#4?8qM&~k-RV`XA`
zu$Ueo>OP{@6LkksHx#t=KO?3YaXTZHXT<u9Xpa%?Gh%y;*ghlr!-)PdqQ8vjKO@G&
zi19I&tfFbgTJV6STm&8h7yLU@B=PTz5%MKI^n9k^!P4`Y&@X=Wd`9~Z?LU}6xyM7a
z|Iq$J`wt#3Y6<;-UXLj_{I6Y)(c=$2{?OwOJ^uW+=7Z?*haP|E@rNFN=<!E54$%7<
z0tbT#L5ue#*kAR220ecX`-#pUDmeV~{V>{pX#c?k==>p_{~<qw|MmOv^!P)MKlJ!R
zufOQ^mv9`Q*Ixn$di^DE@DkTy(e;P|kN>}|N2KE)VZYMrcLfJ}{q8UH2f7|w=nvv~
z@PBMQG(G;%;}0f4=YQz@5BVX``5*E8L!ajnj*ImEm%xEue+eAu^_RfG*Idx_c|w81
z|Etdv(($vf|LOI+f&;yNr`PYKBhdBpr-b8zTTemL=ivno;(abwLy=VbuRjk@&tJlR
zq316J2YUXZ=P%L`==qDDzeqn|>HH7DpvKbqAAtj%{}DLQ`5%FUm7}2l-<|&%cz*5k
z68t{0jhmOudd>RpXlL+hU?A)7k?=WVgB$DL?Ecl6@oQPX3^(n2E!MEU?{0)^J$udi
zj_%dgsD3o-vt@jk$M!(hr;(j*rFVs_cZhx`SEnJY*SXvC>>l1^z4~}OJJQ^i^_+XM
z*l0y#*5hrHVyz3|tcT^WSiKG2toue+Z?6&eS$D0l3GZ*3u#^1a7HMj4VckyT&dn%l
z!Mb(dKIODo2J4#qXz<G}5v<Fzo`ZLfdBVDk$ZF8{{&LoNb!+z}O|@BP&108y7qaX`
z&n0KS)f~?{9Vv{<8kWL3bv<!o_w#A2V^Z6QmksK$4j1+J_A6P&PRMpouX5dlwI6$_
z(b~aotX+%A9fs8J$&R}|pjW_0Lw0Pw)=FdJd#tVFOPg)mwAnH34xX4@Z9hBuUi?+|
z`6<?VS?4E~#Tl&Cq<KD{+ZC~vI>%%s{cKtD$Gg9`y?l!`v)&(5%S4McJ@0dS+3i!T
z$!Px#HP586!_FPc$&68B!7A21C_jQlr<ZTm@vvnXvmsUL6xgx@v-ehy-TaC*NdK)@
zNR&C-E2*3FI0tREOSPDn_U~A>)vI@H<&9ak^0LDF-K$108gp(Zt!ciDX+7~{yXGff
zGhKfB+R$*8E7NO<nVb2Ut&G9)T00*OQD+9O%+c>w;}3>W_>s&Yg<t0wQ25z^8FqHc
zSEJ4Q7!!-^s5Ia8jOj_etM9;zF*7wgw=_zRF;~QWjImU3IEdfBMYyzO1IER)37E&T
z0d;4jqOI1f|9<Al44Db*AGYTA8F!Yl{<@1~!Fw*Ue*6B|>^#?=_3LOS51E?I`X+Q^
zENAa!eVz_qJjr1`>oX$%<B-XItoN$+w`%1VvEEfL!;FrtSTCou{oP#iS<hX4`Wb4t
zu$~RNAGwh1%6fSGF}`m6FxGwFJAGzOAJ)B*Bll`;Av?)4e^6aZJJyY@d2FBWQP!=o
zTH&k7+gVrdGkv=R<+3jOkCq(xlE=EV);*My>&QA!>9jNQ${lv%`D<<7I*wo`_HnGn
z#;37P8C&hM3!GWU59ju+4Un;p<Cg47jPYe1BGhvFUzW2IMo#_Qz4~I-zLjaa>LH%&
z_!nPoG(VqV$8j~d_nlU;V>6#d4k&oX+IlveanPX_J7&;_L7JDXS(|#5SH7>J%Ua*r
zeEs)C+<!L3R5|mz9cvl>hufh}yIBj9?eS0EcVNxy+e9pFHJu%~wo}4_nf~kuz3^M#
zu5V|JSENSto_UHL+QPNbyNoc_D9-eDPF5JpJUeymQN1+ww;W%aBIllLKeNx5rlX#0
zuRj*$Y9th}U2JOg(7eL1t*jm{dNLw{Rjc&2$JI>{jAr-+{!YDErfu(QYuyh{VY+QL
z3Mv?C%IMe0;q`UDGyM&YbbG9Gml<TyM$P5hQbwlmV+UiT@GFEF+F`A_iQyN<I5VK$
z`6tGV$sga3&mK9M8R@E7bLoq%jM;%nyIwZZWh@kNcV#RU94r|t1(!N(fa8EUJz7?0
z1GHzpK8*@l|67Jnx~{y=`sdf1dPgse^|xDQ{Cj#9>)&+#yM=RavwlnUU!+W_&-yit
zpVmUXHS0Sw)hGW#ZPw?`ruHM86ImZ7BuoEPL)Lp<T+^iIGgz;OM{5~|v}C<{c{jQI
zu`TNv{K)dU*%#Jhb2E+n-kPlY^IneK7wWO@9X1vYSlf}EG=9|hgO*oVxAYVFXPi#6
zuE!pHaXkKkb*<TCZMzLsSeH?AA9bpBgLTgC^Y%dEU3TK5u!q6*qu7Z<5<QOfc*;60
z7<sK|ekkjt?w{`1<qhlTH97j;_uH()mh5)-YoB8$EJ*GivtcT0pW=3{?UrS%UDVEo
zlZR@v<HH+{8Eu=$j+?2`f9HE|c5HaFT}G+5SlfsP^4<RU_%rLvA<)dg?<IS>V%73s
zcJy?38o%QfYaN)={)UMyYc(<G%v7^z*3x9s_2zGMSc@)N9}Ft}#+rQ|v^CteKRYt9
zB&GXo9d?BJ1>2L1Icq$j<-{IaUa^p`+w<|cFcuxy;4tY{1UqCo%6bxhiXD($(&1AY
z!|Hn^se8YUV7rcJZj^h*k!|5|H0;QDX5Srx=ad@L*jan-16dQM^Qf5a=l<Bj^h$Ev
zd$PC=({He?)7cdtnBR8y-|yh{l#wagHD!zhKNwK>)r&E<Z2QbxV+LcAadKzhd*+Pk
z{)dI{@@p}p?)M(HJEMp(7vjY4X(MrmFqQ%b#!A5jAAdXtjUDNEoDCQ}y>6}t?oahK
z4%}XHjP<`ahPmjU$odyODrmpVo%N4uoph_u8`gi+@W<bV8?gS37U^y7ca8Pi5Mpy2
z@vL9hy1PQ(Eo6N&mmk@aSBLd^mG0jw^e*euzh;A0V;-^I9!Xzpi{z|Vly^EiES>d?
zIXwO9WqckBitP}#z7Fed`5=1nu6%Y<_h!@V<&LadttnUI%wMxEpADiqbc$e|zs~g;
zvfYPuZXB$!e8etx;;?tm?yP9SI;Hxgo=9lMIzF-Oi(Fo_j-y}66KXoK4*88NU8jbz
z6IM<GeeE>X{@^8B&8wYQySFi?!gUQ;J6-$m0JjP3c;Cfe7Z19?j@$O^(G>qv?6|6R
zi?haa?AS4yxks#zv9>GTyY^k(haL0j)nJdp-`Fw64VuKJ_F--2*jKjablA~n7e~y_
zDPTv}n%X1cYzb>Up!4$g9WJp}Hqkqe-SlQH9iObtjfrM0tlNY>fBYM3-pwp?nn6o;
zBsZKra5arJxut6w;*iY_Th`<G#M4@=y!NIM>aC4grb)-Qey$(ce)D^^KM>~6cGq(^
z>4y7jqa{7ORwS!2wQeW0?#7>FbnNCG-!tzf)2D??R;Ax^n1QAaR}L(0&M*qQd5n>w
zoh)XkO@7x8-Uf{EiowsXo|(>=6rUg8w`N_&bZw20oV_=gQGsS=PloSg%oTAQX3Q0F
zEoLkfadu%Wg}50@1qVaMO2MTs8xYpp)w3js4RD{b<Id9;`2DPIZ#rF6X9I?9ubQgY
zhYc`j+xPdo5o|#FVUAug?rcDv?VS(2na2A6e#7=fJvG+<`m6e{hpuD&_nm21uZ11!
zpZBPRjo~5I-;cc<at5*f{bnY=epAT$UG|$-#bhAsSAD6T-(+*v=lt}*-y^28-fEvL
zMnub5535&I^fK{zOSkUl)1FIM=kTiQil*&iomNGqR1c409Vd@Foff~Ibr`9A`n^{?
zJE2{bWBvR3v-Zr_T%E1&Si3Zb^yLk2v36?PChtm5X2-`5DKfk-XUA*D&L6hxGdu1;
z8^=cm_1STOcV=2Vv17+|o~EN6sK$=H-R3?&tS38mp-VrPC12UG6Ph&*HkY$w+ZcFS
z-{{HO-a4GQx-ITUiyk~oTECmMb;>Krv(sR0Ti<H_+xAoJm`l+YZ_ZX@$E38HwXRVi
zJI3%<ev`w`Seu(uoVUw0Sewc0hy%@cv!nkoj+guRv!neAd~Bb#WUcpAcD*>f7Hf6=
zv;LtYX{_a@n|hiV!K_6+6T2HP3fPgOUmv}p{)!!1`9A2?(_#(sBf12ydc|mbRUff*
zj2Y9v+L1)NcAXd_MY*?(al4Eid#(>=M!4>^Yozmn8M*7-q|{7p#%yNQqwCV^Gv<nR
zTQTMeKY|zwg<rcE3x%JR7)yoUFBnT94#rXu7k<yE5GP}$h`S17rQooWu~Kk3!H&Kq
zTt84vIbZ@m-5;d)8dELwAIjn9I7sssR14fFho9ph&3jQTaHAZ4j)VA~dX{Q|8|CnG
z9K^qOz)~%6qa1#YgY@@4s1~?U4nN01`a3aH3*0D&pW`5{C!kv3MmhW(2kGxpQ7v$z
z9Da@iVukWlQx2E_`?LLkVFYeeQx2s6{A@pv{(d6WLjR#0evX6mo;0cjZj{5%agf&c
zQ!Q|#9Da_2^j>_b1#XnX&vB5RhoxHJMmhW(2k|*->3M3Z75LOBhkwR_5&I9-0x!zp
zCpn16YldnCILblnKmW8ph~vc{?_H#6p&wBW|BQopz53&Q4Kyuqqa6Mj2XVfPp<3WY
zIs6<4alZb?`^jlq;72+9GY;bYGKOk_7v=DC9K`$e4AlZR%HiiYi0h>oss(P8!_RRL
z*K7Xx`x!JX@S`0583%E__>cE)(zL*ha`<N)#OGBQss&z@!_RRLpO<5(7PwIkKgYp*
z*zn$`2DmWhqqo-B=CPMCcfNA4u<j(rJYd212i==6<})A0ecdvMF;6N;*>Ua?V?J+l
zq5rx=jQNsU7oUFBVa(U&Y`e5Olri5ne4@doQTWYTwUXuC-ZJKg&L<99ZN`|N$hOfl
zS;Cl~Gs%2+Je4uObhjXRQWwVjnq6depCHEk#)e8A*PAfrH?J(pcy9YIE#*t)q<WGa
z$*$Cn)UM=*<d@{9<hK-u6qgjI6t~1d;v#7&Un(cnlk7-#rFNutB|jv;BtIp;r8uOx
zq&TIxB@VwTpD~k1pMExjF^?JPk^OoPW4<x`Ug)(OjQQ2;V{g~cVJtKTjBIjnFJsZw
z3+>d=WGqIF{8)Q#GGpPlcYm$ws~L;g3p(U%z~xpQ3EMvQDPys}%Y%;vhZ&121wIRI
zM=%yI8a26=i|c(4uo$~(6=PX*(WsNw#f)X6jEG}qVT@(V!5XGER*YqvOpW7ioBpMx
ze5ssRPqOpFuGEgyuH=X0m*l79w-kpMmlWrZxFrq}7fDO`l#kF)BrfjR@mCwRWz4sp
z*<)7SiLp>`p{7^Qld%{UqSt54ea0f;-T-d&cgCV*(yDFI7a7Y+p9*^3`Oa7lNI}WB
z>=;Y$+KbkWkTaIK_BB>JSTL6Rt?z%>;>uV)K4=qS9Kcv<PFwh?_9MoslVSSeyB8R%
z-%2VU&(3D7Mm}p`IyscF8XNkh5E?O74zM+|W9@%wDPJll)syUqcBOWHv@7```6c-&
z`YpxrBQ7aUF>Z;2xx_`%l#ilc{2QOiXWWO3t-+Yz={8}}0S(3izQrK3Z;Zu~%k@vM
zU%^;>?KgGW)KJFKY+9o0g4>Mcg4cCtwRy@|K0M0Tdg{nnb?-RhlFdrS%I9;(%r2>n
z)vAtT`&Eu$tZt9jpOyQWv2M_Q_rREN#(MCDlLoro7;EPRQ%vnAG1d_?TX?K*&RA!J
zJvun~0%N^!i|NwVn*Y*LzEn=CC)pA0O6~k;SMo#jOY-xF-(nn6TtDIz<CZv>OI##P
z`6&9uzwwz>ul?7>c>MCUuczt0_4vHhWbFfuOZfbA+7X??#}&rB$DT$EV_9Gyb&?sz
zSbiz&G{1@kV`W)n-Ll_E#%ftrgD+oCGFI<DXO%oy%2=bzRT_M3!B{6Beijk9iLpLg
za%cQteP(o%*gLEC*JehK(aHF*e+n}?ZR45n30}<TeTRMht|l;}ueEwPrn||%w3IKE
zlj@0fB)dP_5!;pg_~Do6r{wpKIK;T5IHkBH4(1XUNmD+GexbyNPgYao6C%|3IIN?w
zj>I|~%Tz2OSp2bgVR6IagvAca7%Y}pMq)9>VuWQdmH}A$VCjLSE0&H}w6V0n(ilrU
zEVZyy!%|62jr)w{6_y8Bu3<TYr3A}fEJavWVJX0pg$4V`Vc$6H4~Ko=aJw9Ci^KLg
zY_k~KDaLh+art6grWoh5_&fP91!`(CZ7eunM!w^6GF%?vdI;wuYzJX`2-`)t48rXo
z+#bU1B3zy_!W<|g+K=#Wh5pts1A087$87vtpyxHrfS&v5_y~`c=ok`@k?0tbj<e~t
z2pu2cYfE}fOUFlaZUpxMJfA}6kLdV_&L7cvNIF-G#|3ykna&^4@exnwkLbN4dVh`H
zU!(UR>3v!}Z%FTDV+QnoKV1_+*Nyx?xV9$EhL6+k%}1|q$w&TH$cLxx=cg|Dz=v!C
z-hbC1-fQn{-fedy-f8o7-fme7e$3pyyydhV{76T8-dOg9H)<QBSbx<2#x!1kZxvoI
ze+l0udKllq`zYUL?0vquX?wnr(NkVa7Rc9>^L*8jyLfdwP4453$=u7@e%yVdRPJhc
zf9~|&<-+<S*;uY9qJ&$uJAf;AGl!ekHi1hS_Kb^iN#X*Zz2zo+?#G$Dvf&I)2664O
zy}3H3XNuomKTv$M(zD{lb*~l&Hnl9q<9jymq}|}C52u+mH8;pQ?@EwO(KnGT?Ri;t
zeBlGx`&+9~ozMMHyJv5a!QsZpB=s_yU{oIk+;4{>cL$*4h!T_~8~d;IN5^;YaXB;j
z7^ew*RQnJ<;`&Q|M(g%`m}e3nw9t(A-PVrx*zCl+r0enyR$2VmkGpy6sGIz#kDmPS
zaW{FR4I2bTX6AhFrEU4{)w1}`)7$dxs>%75YdrZTKIi!Qrulqr;~RW+r;U8&%n<Gy
zzl(bpFXx^HcH{0$jpHuoR_9I~ZlPeb&)`0{)%7E{Hf9P}m~x)Wo4lJ#{|&hK^G&#!
zQwDOL*Ah9K<53(k-pF+yWy>|W^OCDncY5(f-mG}lg?Gi%G(RaAB|ND+WUZeObF=m!
znem|4veW}NWIN}tm))D=jH(`wLd~xwBE2aCk!+P7vJOi`ZZD^xDJ@o`n3ySOw(CbU
z&*1(KMu{e4`1m=k_}EkH`IsuEd{j4ee%7F+{LFw1erk_2K4f(|?|=0^@3qaIcQY8k
zJ6Z+t<LZp%t-Q^6)2aP=Fny$8WU#p@-y^1k*O|MKZ~J;V-(pP)-{|lmUdwn1U$a9C
zzN)n^uYSs(`_l3f_qK*F_tfeOclXI@?%EzB?)<J~?!<*93P#0U?75xtU~V(o!>wu5
zgDYH9g`58&oy)kZ!zDYk<)X%P<bwRNxCz&TIO9tXxLzYqa7_oe6n~i*T+E*?E>4_e
zq+q0b;4zartdh)N|07viV2O;&%ags{cnsCK9)#NaWg&xg!;$e`gvK5;M;^~RqA9kT
zD7sM(l!W%6ba^n!>0<wbQBtjAd}3o)KB1WpA78&8ANw_lk2&kjN2e9@QMz~dSqm%i
z;p?{Z)6(R8Xi*sNe@Di9oPESQ4YA;DUM2G;-`gk{^-{mhcdQ%CH)q%L^>v%_HCyfD
ztBi=})%<pI9~aNzUdMgpo|;bK?(hE2-P&Z%U8{YLySS|*cQ$tfcXF8zcl6jv1*84?
zz3?A9pUdrb-NJ1*zRGQWqsy%yxq(|{REaCx9LFsTtIo}{_v2Ev6Sx@nS6om+6>f}K
zKG)|Pa5cGJ#oJ!JCybaOaVuqh)ep)Je@;Nv%<@tDgeu7J^i4EMRu4_sc>{UB(n3=v
z)j^R)l~CNs4Jhf2E=n`LiZWccpzN`8(R}@0KNykw*%cH;Y<QjFucS$Q{|Qm7%`C#P
zQBlo+jAf~m&H(pEj!Nka$XY`xr8B_CYQ_3+s#SBqYuXj-=c!iB0pF`uJYPVyY7XSt
z3o4~EAoF=tN@swdCsE9KQmvW;zOGT+GpAZLhgz5sl~Nffo{ht65~(hQ0lrpOJU>Xa
zY7Y2XP4Rps)v7t*v7_SoQ>s;Sz~{4KDy1?Y&n{Cboq-ztlu3X*`%a~F1~P5@J^(7E
zGr;FF#rqSeR?VRRGon&D13bP~ynlmg)f`@7MpQ~=pm=u+zV4^G6bAU5rg*;%)v7t*
z`A5aGNmQ%mV2l}2DV+g%cNCS<8Q`&_;{94wtLA|3&nn(OMs*n+Ue`oX*ZQMqwlRvy
zvqG`nTTz^TYZUi962<4fL<xqmDB<W#lsHlwCGI_plG;2&N&a<E(z;_P>FNcP^x>`I
z9dAS{(j=d3$K^<QQjgRpc8EP<m)IxmkoHKsr2R<jhhks-YX0DU6!m2%in%ix#ckF^
z3EsXasrELMoK+dk)(AqWwyjZGniiU~+8Cv;*n#H8nV}4LhcfOKpiBoxl)2p%WxlAV
zcn2QQiZsqA+etZ6p421ti5+5(*d_K!JET3*E@_|mK^XmN{t%y{|9FVv&LyD4h-N57
zV-`wvHbZk3tVDB*2B6I4ozT2c4V2U7Fv?wAjPkzDMGLwPM+<cV(89OlQGQ$y%D=r%
zv8Ik_k|z0NJ1IxXlX|2+u|w<;yTm?ehqOo9CG8VGeno$XPca9ppo9_leZ|?GQR=3C
zC_TRg%5-K>cBQ^3cjhuQ|B4mLZ}b=~>RF7I==h<+kJ@P2ytQcg=Pb0cZX8<4RYiC$
zoI;Z{$tT-MIZ~d~BlU?LVvpD*_DMUWJ<={|pZM`B`a^t*tv3`UT-HG;6V*^!Q3RTM
zNegA2RYy6ij-mN3y-|MUiD*%BJyfWkhL*YAL@W00LaS@8Lu<_&qV==8qm3o4Q4zL7
zG)a?uvYnJ8<w-qKpV%Sxh+Sfzv_sk>?UMG1AHSeK^c|Xb-;!$#GL9zWL^1{<`w6dE
zQygRQyb)!DSx`pwJ)?NOg1%=IGobGorSH+gzXkf9QOtn8XOzB23y+oPdq(jXiH;%F
zWcav9-{&u`Su7COoM9$p%@FyJH6Jo<%1E*1wFL9QJQVBlFcZoMk1g?a7JVNZW`O$z
z9cN<(^lz}xaW)=T&~Y|qK>ywm9cR;dNIHhZV?{a-NykU@Z`9#)BK|GVc}O}wqVMy^
zeSpqG()+Y(bdHVAAK|glFEf7>wUdvlHinN*+{H)MKF)`0Kjf!&?8k@bxbgm7?Rc+l
z)p@tBlXxecX}n!WD}GGdYP@9&{0_}V>AZ0rjyI|nF1#P|<1AkP^<ZA_i7(&f9)5@B
zjZb`=D{J}Y7i0KF=W}?ib07Jd=L`9&7jt>_D~GvHdwX*)OOA8*&o|(%-dW0>e%(y*
z&bz&J9&$yk#&fH>KIaPh>2q1K6I`;14;N{El=ENd%h|1I#~H88<N7Rm!f9ukaJ8aa
z7QecFq_|{FdT~K${bK*gkBjk`gU#3KG03|^4zqaEa9P`38L~jF0NJ9hnX<!1=VY&~
z`l4EE+9B<gzNpWldB`~93bKn@fc$;kP-NHVD7n>mlvU>;S)V5yvsdWxab^r3{jLok
znaT1I`a1lytjB!lu6=yK*>v9fdMH2X(p`Syf!(})p#yIlvWT}bJj0KCzk@ek5y=~k
z7T=+{%AD69Gn3cT9LDP$(Bs=DI`FN>FXEdS%K3(!Tk&<<CGwiOzI+wAj92q$!u_uA
z$UR^22Y0VWBzJ}TgFEdoOkh;7gWEDlom+isCbxKuI+uN{fJ^Ounv0n{lbbwhIOiO4
znKN@v;s!Q;!*!Z7ovU{ts`%4|^x~6iONy7u>lKHL`>J4+X;N*7PGm=B@3uBF{k837
zQERKn*0%d7yU_if?E9fcs6ohO)Opl!WccD68Z~|ia>|jR;C1RK>eNh>GDsa|73}!I
zD6#TjKE853AM+jfs8?qEtP7X<@CjA;sZU?<p%&lxfFwQMdu~sDlG|Q>qDEcb&I`X6
zGQS0HUU-N%iJr~Nhg&H6QNOM0`Ciik`EKI}@EslR^KCNL^Uc1z<{Rc9zOMHvUNiU%
zUu9ioUcF&`?vv36?p0wP_h|5b?pD=j+{M>1+{w3^3Px<L58U=T)427Y8gk2vf8+AU
zyK-5(HMx|FSGnjVpSV!n9o$4yjvJZ#iR-U5h-<r6&S_5fDZbz7x8kA>LB)|Ku?j{$
ztB*5n+h@r%O()AdUKh!hj7yN6(w-*!()SIjZ_1%gx;xN-VQ<lh_1lnrSVQEy^)L#L
z3`23B8ltp1(@=J;5C3AcWCEXnKJsx_9`mu2H2LT|NBOAXwfI?cU-IEo#_`ia&+(J%
z=<<PGKl9$(5bw6Ym3KILhqoOV!<+r?#1CtKU%{x)wHCZ?^?bhl{%L&6KCSsi6Z81G
z#(Vi1=d1D>ug~#n3A)^;+9SBv8-ROq^)z=c_$+tB`2u$-V?1}}oi=woFhRkHOPS2=
zsd<ar=CzdDIKL~mI;lIiv_lVWq1i4j>&+H!_VeXjY)e}%Z1qEK;&n~VBqW3D(IbJ=
zdeNx(=B}E>Ic>TSM*HlVpJ%>5>n59Zt($DGLp|BMr=3y#>zYVs^+WVq{|#v5RWoFN
z-xGP~hoY&*_9!Z;J4%?}6{UJDMVU2kq1=?oKNuxNKI4-rf8`Tr{l+JJoXjUoYR1PO
z{=mn#{LaUXJIu#Ud&I|t?Bb*Qhw+iC*YXj|IeyyNI(*Rk8@#)AB5%{JiGoqbbC3B3
zH$(aAn(cTs-8bBO$4}gIixJ$z)@`{vT^e#Xv?g#@GOlqKlQp?>mpXH2ob$L-tvhhX
zztrcBzN*6=erKUzbg<@nZhzf39H;Su+jsjKw`aa9x2yAIZhO!NZi`n3ZsW_v+}aDz
zxD~f=a*Nk(;d1k<ajDB%aN!rHa^{g;x$0}i6GpO$8$FR)WLMPV;#6d2YJgk^*G3`s
zo+#q>O%!wC8A^D$7^Qf1KxsiA(A>_KQRaMCH1GB`l&kRq&9D0gEv&hoFyhuzDU-oc
z%!o?q4Dv7|Dy1`+gBejNok1LCM5S~FQ!yhdr898DjHr~(z!Ec}QaXddm=TrI8Q{N_
z!ci%mL1WB_O6d$LVMbI+Wl($#uOFmRI)ft2h)U@U@V%pADy1{P>(Gj+l*)j`?}wsN
z1_L}_$5JVk0fYHaDT4w2tua(eWgx@HYbs?hz~_A#mC_mD@uQ4N=?w67h>S|93=qDa
zrBVh1JTHQ%l+FOpuOTX>Gr;q-h)U@U@O(d_QaS^C-w07DodLf8h^Um#AQCf*494$f
ze2|TzY%EX|mxrRe?ncqG)lu}*11QGG55;5zqL>>aP^@+WiXF$J*yKehcFP46d!Yl0
zeLfDwei@DM{cwd=q)9&6PRfz;xE`smu!HRpy9)cb9nv0Ym$a|&<5%<t`xF&CA4Ok#
zi(>mZp}5EZ6u<i<O1M81C4Sz5lD^hQ$&V(Wl)Y=w?3ry)YR7vhb;}@>)-3|1#kW9d
zM{-cwo5cvPt59e~n&ji{q#P+v>XG`y4zWk<68oeb(jIA-v`_pX`|MZqhximd3-9l{
z%TYq7?kH)vIZ7GX3Z?d%j?zA~N9noAD5KSLlo{R^W$mer=AFNdvJX8+IdjIM+@97b
zcZng&eY6|l^*IWyNaK96os=WxNj*}Z*dg|aU1FcKL)s(llJ<!ozoI|Hrx@SkD1O-u
zl(be0&5lV#bM$wkxjQGKtQONzwyg~12Gl|GU5BED`rT0e#Uiw5lrvg9#{ezKzK9C#
z4N>8RV+gPFQfQJU`D8mON6M3Wq&~4j>=C=fK52)vN7^Or6F+`Me-u7_`-BpTZle@E
zGn5u+gWsQe0%ZlyLOFwMQQm`Kv@obYDtIs)Es+JFr7>1$`GzBC<-QcOdU*r1)@&nM
z7k2>RHG>LG(j=d3C*??aQjgRpc8EP<m)IxmkoHKsr2QD|$1ms)pUA55@r)WDgLM?v
z5m--;P~#_K3B=-q#T|<?mI+wKVzI_zhQ$O6U}3Nf#A1M@7ZzPCow2mT(h5scEDf;K
z!BPW@MuZyo2fp|D9?MfKx3OHravTefWd|0rPLr%bBx~%*8Zxp*imU;`;}pdj4zfl8
zUmqy$F;^1r5t4ggc%Fyc!%)o8OLMGv4pA{jM&^LX90$IJS6ri#YcO(+LB`Nzj7Y{n
z<Qz}VKjip|+opUl2g--`AN-q9++(J{FM9l;$7g!}qUU!y{-NVTIzFW1XT|+_dVN91
zhxGasUysx4cRIg8=YQz@7@a?)^M`bPRaoOo=bwe=SLpmay?;T+hje^M$A@%$NXLhC
zd`QQKbbN^W0C@&p@r*YeAO07{hb8)ayzO{C=HYQZDzpJVs}jf0X#ScH>&Edx{Q`O4
zK?`_~!8Li8{_A;%?ko6lt>paZYJc!%FQ)P)r^9*J_?BnV&hZ2Nba{hOwfSCpvv}Ql
zv3#eG>U`TfyZM%9Z22aKEcp7x@A=xCGhe-AHedPNVD8(I*4*27AGjx-H*>cfPjD9&
zW^l)D_T_l(qudTBZ*F~lZ*J+Ok=*>6np}F{J6!yj7;gIY63+cY2xtAVJICDF#p!O%
z=bDDp=QP@UDZUwWx_E2ShvJmftm5Id4cX`SQ`iPkBbh!9YnX2*<udc|XR@>>9b`Lv
zhRALo2t^tnyQ8KbLXhtD62z>pfvmmuBlpG)(X<V@C@#M@n&ae+^0bfgi4VT=@tZ&K
zu`_G%(MB>pQk~<&{gU{phY%lHZ51D&*O&M1|B#>5X)5pZA(FRSnZb{dP2erqt-NWs
zOnzwWH=emv%n$55gzr10GT-Bz7vFV=3EyGwV!n-yH{ZO*OupgeSibI|X1wO<qkNT*
zHoSWOz1%15C)~@AzjF`XeB*9ZyULyK*@8RfqQ&ufuehBj?{gbJcH~ypuFMrwtHEX6
zY|N$1j^d*0T64kM&vA}-7jULI54nC<s&nm@S#Y(hgcZN2vbMO`Zh3J=zk9_alk2kg
z9^Tovd)96yCNNjlvgSBhaHE^DrKXEyCyVOIKIxx8T6Y(q4%^S6{;ijz5v{5r`@Ooz
zzfOHLt9E6S@UbIGKY1VJ<-OvQI=b=+zHWTnrs{mmvs^y9Iq;DNpZSQv{`_>?zWn45
zK77!~BfPKme%_<fe%{%lFK<6cgSRQM;?2Hx;*B3W^0F0EdBZWWynf>XUbpfUzGJ_&
zd~5D2-*nwkzW%r0_}Z&~<Ew8n;5BN@<^C|s<bJO@gnQYcG50X*1$S#?Hg~z-E$+-v
zckZaCFSkE4o!he`o7-BV&8^?$&aLpM&lOyK!exId=F&EJadFSCa5MVpa2{6mIqRBr
zxk2*tT!(jBT+Ozlif`WTS-jljTk*KK<?NAHn}&EeEMoqM=49S6{bVa2Y0Iu8tU?+!
z>!M~>^-=eJI%vq{B4pus0Xe;mL;>5*q6m-rC~lKGnq8ueGIwO7`I+f_Qfv;Nn6A$!
zWWVL(=Z)jz5~lL8p8NTjzT5byyD@yES0W$5<?u6#TJuu|8S_CejCh|o3*L3_Qr^zo
zjJN2xm6y*+;QRSb=DS`o;oGd=$Tyaa=j(*6=d0J-!&l<Ab6@B2+`IP+xEB_exkmv{
zxVv52aW|}AbC<6)<<2k8<xbD<%^lyS&mF!$pF7yxh2xBmb9<dna69eYxUC%@aT{0H
z<<?xd!7ZQNm|L{2Gnc37%*}Nl#3fFh!-emt#ksW_$_@RT$+eXk7C-L3xi~sm&W>Mj
zT~>d*t88&(i0t|MuBdUVp{U1>T1fuBEwVW^6}dG$fP(BsqVSg8Q1pcxC}DM7G`r&?
zly2{avYbz#T;t<tVRIM7`(g2#qKaz8yCzHUx#SlSRDjpb{Rf1@=R_){Fj$4RQBlnR
z_d$+|Y6e-D4;9r6@UfYrqM89-)6Y>+&A=b?p`x0B9p*ztH3MVJhl*+jeJ~#?su^fw
zK2%gQsD=4ZQOQ8@cdYUKd#Y71z}MTwR8%v-*ItUhfll?`bHHOm#ovdg`tLd5b5}7H
zl?+(?yP~3+fg1gkN`S%rhl*+jGHv|4JQdXp@Oeu{MKyy0%!i6<26*f&qoSI@E6j(A
zN(KmDyHhEJ0X}yjDykXaxkE%nH3MVJhl*+j_}&GgqM8978zL&I8Q^=Tiua;Xt%Adp
zE+}eJ2#PlAgktorqu45zDE9O?6gTG<iiZa%{#hhSh`xgo>&Q`J+&7f?ItV41bU;aS
z@=?;EWR&!DDN6dbTM+V{<dJQ-3@NLqgX<C-#1^qhY?C%fTcl0WHt``6`$FVb@+V5o
z9Yw!BhGLJ;L-FYwQ6ht5KVK82n8%>m8#kcT-xs1eZLLvy_oHZTi`QuGlMN^%`xDCS
zcob!(0m{6%4rP9c7KD5!d3YNsL&}mmq%N^RY!REpHfe*jMcO276Ca2#znniYkto)0
zI!fqah4+0EG`oK+N{jx6($~76j7^<U*4z##+iVfadBC7N>n~`2PA9ZreF9pTk%02$
zp(uaXJCy&?LlE*E=aFrs3@J<Mkh;VMu|;eW+oTQB7HN~TO?-&Lz7YAP{CS#<;?itU
z;>RR3yITh|N3Rmfs9FtWZR?40<nvM9hA6b)YkO4CISwuEH2@Vh+KiT-$U@64)}s|^
zwb9B@d$jVElOXsz$s^lH8B&(iA$5rjVvE=$wn-bLEz%}woA?lgeIfEI`4e3~9>smz
zjFJj(qttpSD1F2(lsOLXZ}bJ_)vSsZY(0w#Osvr2!&T5y11+>X!xF80;Dy%coj~jS
zO3;SY324*3o(dt~Ngmlo%8;_84yj9Q5L?71u}#__ZIL!f+r$Ur%dh57jTZW?;*2>|
zE<IeceD#bUlu@su+ogCaYUiAipA}6X6YZ6y4~pfED(S;w`bZgCDyLjuvRhXB(sqfX
z>bTpLbSrt>^&Tph?$*3~^}2p9qh3X~OYv0H&Ko5^E1G^P+AB*x6w5tP(vOtV%G;&1
zs&bOuvf7unOB_|l-J`>s57kXqp>pYw8s)2Jda;ap72PhyQ&BsdlAje#Zx!v8r8lT3
zS0v{DQ?1;ND!a<%q_o6Qb=-Xy$99SRdIOb9e`{I3dSAbnQLm!grFbf8=Y^7=6-_@7
z?UkkPs3><&MZU5<RpnIKRjx0kC621&9vJUkC;E?WsC;^AIqO-LQLm!grFbf8XNZ!Y
z6;1aQ?Uki_iscNHw7!_`Q-+qxDc6_mmeszrUE-)ZZf3qsQI8CBR6ZS1&U)@;)T`)r
zDV~bju~YK1qG=n^URl~&MLFgC(PBACOLmpZm6b1TSJo0o)o~kHH=phtISG|dFDz%h
zs50tRbh{K!MeT$r`B~AluV}9<?Wv-ix0vr!hL*}Hx1-9gaycn2aa0{QJm`C3L?1m=
zKJ8o1dZWsySJCZKJQcMwSjo?drVT`UW$E53%Joso?<=MyE!kBrC*`ZMD{WWS5=Yf>
z8;@KT`ngDk%BN$?S<j`6dKKL+#Zyr`mP&qBG;Jc<D@zYmQO;P*m-KKkEomuV(#q|q
zva4K9N=qD7$88d`$7SR)HB>%rTF!c%%cxh;?NU4ywbM|^&x)osMSErGDk{oV7xQa~
zT1rb=IbXRQRd$ujNok3r>bOl4J0Gau*cz2jFD+-i>1EWb=yoZdirR5k^0T7p38KBS
z^f(pe?8JO~72BnJshn~<s_ZJ4lhP7L)p3uq`8_1G&U{op{kEL-O3J8L(d|+^6}7WT
z$<K<Wmx=bu(u-A;D^!uMY)@4=Rd$u@OKFLt>bT7wty(@Rsy!;7o>|U%=4I5Y=yoZd
zirO($^0T68J<(oST1Q2>u42A!8Coi*+>R={%H^cA#8GwJ<_nvRU)D1gl}~RiXT61G
z)T`)rDV~bjNmuf-qUmJOURgR(MLAXZ%Jx*1Q<bk=UrI|HRmW}7d4|@*YuTuL`batJ
zZ7QQ)MYl`wRMbwPlAje#FBI*SrRR&~@|AR1^^~_O*O%-nms6FNIF=H3YVd)Z&CHQS
zxpZVX>v@(@ucF(fcq(egPRY-Trfo!fWoc^_<y7S>+f!9eRlah4DJ^j<CGPWmpB%h=
z2U(OyUuj&vdY8YJQLm!grFbf8=cSUL6-_@9?Ukh;s3`YPMZU5<RpnIKRjx0kC61-U
z{jq=c@FaU=Q7)ZU&U&F`)T`)rDV~bj!N;$^aH(k8MYLCzc2rSrqL}YohL*}Hx1-9g
zaycn2aV#rtsoc+M+JAnoe@XkLEUjukOa0;(>ObW<Z%OAndj9&!^P6;DC`+q4&!{>t
zDce(3PL*Bd`cnEIoWII7ZkEPXkxJuBDZeaT(Q$2A^`-X9Dp%I_QpSx6Za=%-M#Z>T
zIbYQ{S5^Lhy1ryj;#f+Y^!mM&`HMn&{jTUAsyUQz-cXuXqVtD;<)?a{O(-YL%ax^7
z&GV_sSGK3BoT_}~`chij&#L1p-+c>dUYp*(_$xov^ZY_N>ApZ&TGf39>Au9TbN`}z
z_f4hyYV`i;U-_xN&nT3W?n{=XRo&-QbzfB3o~m-H>?+rn($aob9as6*ZAka^>H3Yo
z@>9LeK`1A!izrL0T4y1x%aGP-DBF?prE<#KRoPW8C#9wRtU9jpts9fp718x$f90oo
zot02dT9;LpR<+Jc)jBa%`O596va4K9N=y4$bzJ3Jw<@h`r0ZAz%1`w=PobQ&F0?GI
zYMrUHE>+Ufx>ZR_`I1&{N0nXWa#C8_&#L1p-@18eUA3>cep^caQ$4DTdKKL+`B70j
zs@IJR+<tbu^7?dDX;tgomDj~f>*ggb*;6hj<*TwQZCBP3N7Zqb@3|f6xeogL&R_Yd
ze$GWGCp{NamR9wgjr3g3uk-xQKYMO)lK7mOvi4AsuUt-g{!G&LVm;~kG)V`F`9oFM
zmFg=$7dTVN-p|@mt}mq-CA(7llD|?~ig$`w-cZ!icB#Db^LhSaIZ6A7>3{nCpY+_I
zrC4v2s3(Y8YTsK-yNTLWX}jc?)b0pzJBT_+)RKP^e>>TSFZ=hO2yjGNuZ)I#otKVU
zESmgy$K(E}-Q8%rNl8(t^Xh|x2i{+ebUSXnJQ?g!&(2*c<-|Cn-u32h9rvUu(my>S
zsJUw$)Yo-Zo8t{Uk-_mTR#o3#M*VBdI+eS#6B^JocemAJga&*rak$rE9Wq>h$M*1+
zm1tnUqu*vW3PS_a0}h>6tB(fmJCtd$6F(2F(2K<U0b;obu^unlnJU_K5Zf6owmV<+
zBU<pwWvl4taM5o!F^)hnt}HRmcrk8Wk%O<mC9sbQJ_TaE0MX7V(XO8opRS@GpGCir
z=;u4p?*U>Qwqjg)Vw~Ax+;JiY{2a1^%R*7l7V|a5a{gky0@2P^(JmvlGf`~!rRYab
z(XSICpW33|j$$0^#JIx5IM<1MOhi7`B9}L^&nKTA|Dx6%HRzcAV%xbu)Vi1T#l(?K
zQ2WnAGP}!PAf2JDjSD+tyP8%@8#Ek%dRVuZz*}BGy_#I?zNfzy>TO-WfO$Fr={pw{
zhYblueMU45_@mM))VGeBVW%@qQQ!5s^S<6nL<YU%LvsVL&b~it^VzS+;GVgI+0_-O
zU!xED>ECTpKfN4RgNWIJo>fWA?<tl$DAYTL(+WF5qTSU(JEz}@?G}lCj2HcyBKX<w
zo8Y&Vwh%{yPhwo|Vw@kuxW@<_&Mp?YX#JJX+4Ew#Q)0a{qMeZ<pN1lzuS$HTihf0i
zein#)@b6B+XPFq+bdk>`G46W;hi|V%E@~pDXJY<+vD_<>kEUp+NVL1)$9|CZi?E-3
zL_Qg!pVdUazleMuh;bQ;ab6YUz9jY&7m<sL?EB_vaXoDdP$Pp_%XiQBBW>OK_T#@G
z)M3w7<Bz8*A)V<n?Fu{3M%_{z+ea1XBE7ox9@e~b5B01a)o*o{G3t4@-PaNIeNeAC
z_10>wzJhwUb~~y*U_a`ec5}kh_j6J2=T*ig_~GM}fl1FsS5_c>zZOL<b0Uy__V>YR
z8X8EyXxm!0+jpeT&GzjXiqi^>e;<ncyJ9)HST9kulPTEseJ8fFU2Hd5@MBE}(XVvT
z&kjHQ7UHNX#<f|DQ(cU^zsTW)$YsA0pV?x$9zW`p&S$*nhp`Gi%S1j^#kl;%IOmFS
z*HPkAUEmY)M9A0eE|jbMS*+*sS3Z8tl=$G^y<$H-7yPRCUF2gR^4TlyXLA*N{6#)*
zB|b>K+NV0kHH}f@D?@V)T6aU*uj-iBJ~{?<>^``DyEZOJ$8Xi?-xqB|U0)TfwVmD>
z>0WDT-r+Qd^sMWS?$Ko&>JfNZ{bTAk)U#982aC5wsOOa9Zx`N(M?H7kzq8}US=95{
zp{z~J3*66-9DT1@9rYSAN3C9~IMmB+nTutxHtLl)y_RPE9jMo$%&J-MT|{ju=HC#?
z1>v|A_0|Y>YT*7W+7;SK5ZfK{!w-dDX=$RLT7utaY{WQrigEb}akhT{BW{61oXBOI
z5}y#U9Ef~|h<4HiyVGun?bK4i=d$Rhi^yl47)Oc7NAF*85<Wp+RPgcsSH8qi;;P`>
zU9@}nU;9C^UnD-M!hTB|D*E~2hu@0*jN=gYd%75Bb0t25MK0c`lBVzQ_qF<?CI{ZP
z8kAIm+Kvt}+ate$I@+bI=rKDR>Da7^Q5(?{b&cP+{>Y(Re104-&F$?&)ZJob{Q4a>
zNUzJb)QN|-Bfa9A8~2@zKs~xNncnv$e!r2&YmY{oW}_ZSc1<tOOG7>84D^*XYlV8u
z&<?1%*BtdQ9cZ#<)eWR~wo3b@u`iJRpzfy(4^9;I;2-&la%**kdV9M3%MP)d<sr0F
zG+b;~`-dNrUxJ^te)uiLA;l%cxmSoADvKOuid_7a_{<Ucyi&pEp4iR;kq_=if6GT(
z@OyK+^6@$FuAF>|1U@I8{qXCrd=`lNIq%1Q7x)DK=odmi87FY5`NDB+*fT1W!GF6C
z7{$DPWE2*?&}icyOXcdlwB)@Gn8;o3&XT*gERi=HZ1A7-58dcm-Kb90I3wMUnnnjU
zME@sk|4-SHPk%o|zSX#){LrQ$U>ougg3jo{)CI{f%hw!c&2I(~(+7fURw9fT&==Ya
zY%YJ*+)}<Lrl~w``Db}l<zo5l$JOQiJ8v*rVLRPuDECmd^>uTj_TdGnm+t95`)u@T
z!!UX3b6xrVGgm-!pbzvN_7>#aWiWfV3A*`a$dA71_UH2djqgS|gYU|FH?Ak&(kd6K
ztL*?S?*x!lo(p5%&4!VB_rRgaYH;+r4<4KsOda$cqEpX6)SC?uwWvKLSS3MX<7N<f
z>j6aNMnUYbMG*h+E+pF8K<t{ikg9PWqEFO@(8Pf-aj7Qs^&BBz-!estQ{M3W3F!93
z6MTIiK+LQQklfiH(uQq@jKt0`uk$L%9(n`jnh%2LF9tCE@m8>{y$&2k41jhUkAicu
z2~2w#1CHNo$=}&71054azM<h8qdG(HNcN4Ab~(#+#vGF;8dWkXto07Mwweo`TeQKg
z?O~8<MS;cW>)_E?A6!DufxDR{gzV`6Q@*SL|3zUCc<=>;PIv^~H<!Ths5{W@y)Ec}
zv;)IoBVl5|JMgWX4U;BKgwaD>0c9~@l2ionqzH`r4FIRNCNMeoHN?JTAuTx&GAB7h
z-s9>p|HW!3n56{^*OtKiYb{{D_6o>a*&H%cPC@$DevsVY2+XMA2=JblzYBaX&nniH
z-`iyaz1r;srj-$RwOj*Jy}coJu_4T984l?OqhZ0Zo3JSMEiAs*2j-sK4|xY$K+dWn
z$jcZH1^rsWB6B-f5aJ07KNmu7&;^*=;4@@&_Jp)gP7t@C0A|Jog0FLZaEhG|Ciw>;
zWyJ!BS=j)h6R$&*OLd4odJPf>ErXPm^C0zoWk|a>7m`+{LgeAf;PUzi3~%fN7SkJm
z%L7A*Z+aWjmKj0eU_DsCJcsNgOUSyo0n%1?gm|Zh;5%jvjPSnzgC=-`Wqc}Fj=u&j
zv(z9g#s$LlJRoXY93<7*0x3N*A$`tWn0u-c<gD8X3!9CBoa?J0?dDMM?{Wh=Tp9yD
z0dv8-(M1SctOFAt+QF#dTF^_QDH#1$6-@0NVUAZmL`Thl)M4LXa;qSivbhGtuGWF{
z@6TXa!(_<&b`7%3qF`Bq7p&;g81hfQfV?|>V9C}4uw;8>n5TaW(%!a%NlhET&{a{O
zuDcC_61G5`og2(BPXn*kEI8d^!O@^CIG=EZ5P2<#tl<GOhNghyVJjGwIt504j0AtQ
z3?epYK%~n|@UdzqPv7Y&Kh<0pQVlF%{+Gp2Sa%96d>#TB=Sm>CvN^=wp9Hg9W#FxO
z9^9XNhgpVeAgFjbSnf3g!!My=U$hoNG_FHTkqIPkiHDrZ_aRH~0VKCy2sZUnm}cv=
zN&h#>2-lJuMc#rYJ4OR@sv5Ly*#ZpXPlMNG14xW?hByshh~D}ILj8|J&j(wfW!v@8
zJ;?^fp9lxX)jPm>U32j0jDO2MwV}ZJ4~TxS78u=h2u3oPX&eq|Gd4n^*)kZ*8_3sd
zHI^T!mkMrH(cpPyJlMJ(1o<TwFx6fMwhbGD*YJlBa`p*?H~b2-u9m>G?q1-tDimUx
z9)=-<UPGKk4rDH=1n~_W<yN}m<<nLK$q!GBf`D}cpg<=a3ja6``Dbk)JFYX#jgNrL
zaCOKD=><9N*C5@g3dA<M3I5GT!Kh}QVDM!L*p~!CM4u%PdwVa$Uk`?;@H7ZK&=d4*
zpZqvBG^v#it$H+s!eRAc_2TZZ_)!~JVYmkhDm{SN2lqmBUU%@!GKNvvmtaH<b?|do
z0uk?Z!KGgqB&M%`6b8XOjiFFb%^s%rsRZ_qDud~IFBqPC9!6W7hIkte+--NuH(lNW
z)t4^<Gj|;rni&hsZ%4q$JRck{SA$UY4$Qhy4DnUl!R&k^nCm_ea>q@9yoY~4mhn*t
zit?8)iroo`qnE>6ukSGD^(WAn6o8t943MvvpOe2HdLD+>dJc0niXpGLHRMd%0t>rE
z!g6*SEDs+Ci?5D?yw43G-8>6swXp(^d4@2ssw-G;O#{C}!yv-H2L$h^0j?`Tz&5}U
zMqdMo5Au&)0#!S0g4Mq6P`G^>6vUUn;wi5nH?IdI>&=1b+m3<pv)y2J>v!->O@a^)
ze~8{c5Yjrng5>AgkUgM3%x$p+l6reX=HxFhw#i&*xAZIYu?z**o_N3cT!Ora?GV$b
zANYN21ZFW!!LFwXOnQF|{5EMp%>8{Z#}DuGI^i%ct0|=5-*mpgPFU{t7IHp+hJk7A
zjJD()BK|^5ojgdN9}DJtljWP=?~>b(dnKRK(+yI)FM|BL>M$qI2o^7B00kqBV7Z0`
zEMEw)I2s@q|8C~IPKJcZZy~bDKJabl1pU9n!nj^);9K1Vob}s-^?L{K3o(VsA)R4n
zQZ1v%s-1r@>SBRl(90#Tpx0xVfAKZs)(eDW_9#S#AAw-C(|~$xhdx*5!{`HpAxZx`
zEO@;T^6eYJybD<n)p{)0E$sw;_!t^l+Y9myTY~4pi!k)%NEqMA1Tv?jLe7~k5N|RZ
zB9Hci0SQ&Wc~U)qVI9EJa2JgEeJrFrehCHZ3t-%c#l$C+H{4KeQP0z;|H^rik3`EK
zd)0=32A}1}{HMsDe%A)AmvaCd*Fz43!|Jp+SkZ76EaO|l;^2m`h`R~devcvBL<f=v
z%!lx4iQxHq3^*@|2Aln6U~qOR7`MuSaZ3k-^P>WA*l-!#=Q)GlH~hQaz{sC(GyQYh
zM%^6ug3;a&Fzbj1csD!-qhF@L%nuuI%$LFGR0x>l%!V1a3c)Kf6rv*!!NS_7U{<H5
z^1jt86FWw$Hyww5^EDvaa4p!rx(Fi-Pe5GQ4468qI?Szm7?OW4f|)nB!*3o&5WS)h
z97i33S*NGK^t_djo;w|%8_z)1iR1q4lk#_Y<A{gg`1FAM;jRx*{r3}KwR-~u>N)|S
zP*~ENht-=rpwRgQ<ex8rdDE{!bhq~~*+LG{0~f%I?snkU#TM+<kAR8CMd(-UJjj-}
z0n0Z}!MY=Y(VnMZe04YQdcGS%HtE6u`+|RKzeTU>Fl5g~h-jG!kxw&V-dP0czi)<V
zOO}I-!RN)2{V?a92SnpBLW=8EFzCMS&o+$GZl8j&c~ikBXCL&DWkde1eK7w+4kVnN
z4pYoFgUz1v;8dd@q<Y|f+4>mFe`*11_G`mJdrL?x`7HlxbN0{nl)uY+zkUMM#s$m&
zxEe2C+;ISmoBsja-rWTCFBY)=Q~<28c@DV_jUYXF38b_-3W+EO5;dEH_sd$~Qn^34
zXP$*|v)+Pf>?N@H<_z8~bimJ34nePOL%@4u2<!6|oQ8jg*&jl|clQqE_DR~vaC3dI
z8mtegw%cIY-lkBn%oY;t<H6y;U6_4#CFBQ&!U8T4#=rHIZ_~T0Yy&j6HIx7TC{gac
zMF#0L^dS0d4VYQ<1k9`Tg}?=UAtBHW@~6**$(uc3(ZNozZbATTy>}MYF`K~oVh5;d
zo`h=M&r!B(^r2oo`8&4==(8*b?5FF)O!Lmr|63qbdS)hnH*qvf>AoKp-6(?j>FSVf
z$ij^DjWD|8JPhfi4|WHBhw1$_AmoNIxagb(8-o-WF{(cJn5RP6)Gjc!=2(dQwjL5Y
z89}nvT?l5<<nbGlm3@%XP&2p)QrZI)zP%3%2Hb(kF_tjo)M|*hkp=Vnr@_k1t&ly-
z)M(Pd_EH^58=Z*$C_gz-F?NPH&xT-oHWJX&7#Q}|3z9dwK<3OzkW<YHqP0)M(lZ5+
zeP}u?=(G_Q4w(RHzG;w?Q&T>>&0Tcz>Q-VCX*s+!%4=!}?GFH$kGBE+2I=6~@eQ=U
z8we)&81tLEHyE^X!DGI|kaxT*gl)}$MJww;wofyNIXo0%CO?KLyW2v<XeXH5C=xvS
zjR8mF%`nVx2{=4>22o8Sz<*B)1omnVF<bFIUF{(EF)fnf`8zGYS8Eo`U0Mh$>xRN2
zwHGimJQz^IG?=z*1LP#vg0*I=Am_Y~e37f>-_?gk*Ec}z!Fc?0eJh0R*bA*Txx?r!
z2O)_mhNXStVL|<-c-)!@`F+2_q8T+HzJI8^#J7dKgY7lBb6Ao5c>;oVkEY7gch{5M
zjMpM<$m=teU`Rw~sN8lPbTX<A6IOMF`9{Mb=gk*b=4=Vs3tb@vbDlL}DJ=d{A9AAB
zL7wgbi2mpSvp(Z%EX!_?)UFT2X*)r*a~e##{QtFcCGa&|?Hi>^LP&^MqpGdg8e3)T
zxie2|$wwykv|_2+m)JESLXd=P1QldM7GhsRQTjzoYpGThwb%NgwN$@S%YRP%&iBoK
zzL`67@12|6v@^fo@t*UZ_kEt{JehlQm$W6(ZA1{lQ(wn~{tFQD^af@Q^1|Tu6Xnh%
z&0$WEFV>a6ft3#fF~hGh-W6Uq&3qAwxLVtg+AavORa<d=W`;}q$Qs@1=3!P_Z;bsm
z4PLdK(fhC;`oI1R?>);v`mxur@J=;E4>Vyx;B-vgUXg9qy0&lnV0Y5XemU|Gw|r9*
zl>F@%?yVcs*vC%oxjlzJ;r`fq0$#x@xChanqU({x_=sPPIX*E+Y5N^=@`40UW(v>E
zjgc}o8A%^RAo0!<#Mkq}ytvDVKKBeU(+?re&lAyCHzFdZ7G`uFfg#O3(Zj`rzGteU
zQ-TkgzE_od5nfvEOf<*!7?X$ET}NTgiz8S(trkA5uo$za%?5k^cg&l95{sS+eDjjf
zyhBymBaLwVaVi#-T!MKI2Vq=DG0a|n1HB?%#n8z1NV>ilQLP3eRCxVA^n4iC+?^-=
zrcQnHAFEu6a!1hn;qP!R^_lQG!3&XDQ_=5YEzAr40^vV}!h4|5D?D#R^^U^q0sW9U
z{zI(JoP_LCjgXZVfhAjh#Pr#dv7pvgEbMp?^A1)=R3|4yeI0@MCpKe#y+H_%i^GKO
z6)^17Jp>5v(fkIlK)rj<QN`_5G!OKao`<fvt=XxVKKv<WPx=D=nohy+pQ~VIbT=$+
z@gCyquSRIgMwsJ#5Yyt*xGksV(LT9|+7XMfJarG^%eKInfN#)u+&r`x`6<HEDxtyh
zQy6^uXViG~D))+?FX?9+pGxKSZ#cjW%Nd3ed8r5+d;<$7xg%t5TVQ(!u#2KG>w^?T
z*Ls44djqie&KHQSlY@18o+A7GJJ|U9QY2q{hOka|kT6wv4Bz?<$rFbodDmksPV0b#
zTM>xqzaNoYHo_5ti96%)aTyNp?MuO+KaOBvR1R9iOp!Ya_u;c@a1V;X?0@$|R7M{}
zc$7uN>%oY>{ZHZdL<pug?S#?|$IxDOfM;boHy4H937=t7jXyD^L_Fr*aDsQ&LENUZ
zzk%`U3OuaI-9O!g{pI#mih=f_tMKOdC^m2XRHS7r#G1aTNL&$u$e;Eh^VE0PdGa-^
zomd|6Er%kuTL>}|-$d4z8<F{MCuDp$2TQ+Of{X^9$h`9y*^OTyW7BaYuSmqwkSj>K
z7=^?aKA3m<FlKd}fN>9^5Y*>5xDB;2df8~Ifkfa?CU@p+AZkyX#O>|mkAV^KSl*!l
zqTU{j@G-9-edz#f@-2yso5E|i$6nmk?j6PbTITGX(^t6rhpu2mxjf<h?fZy$w+nYJ
zr9J88USE~~=Q0<$-Ou}qvACRe*U%^I8+gnbfV30!urZv)(jz&THsRlhJn#{c<KDyi
zl~=GLAq%U93_$eJ`NC%;;}H9}6k==i#?o&mBfHnPSg>gtR!z@E);|YgjrTLGoU<3n
z*F&)2+!w<8i!m5?=P7!A9SQ#dTi|gg|9LN8*V6Uh!72R(DkNs1;#eL-f4+)tWkQiv
zY%r31f5q^pH4q#2Ph>Q^i`3Ca5tI5V_~_EyiS5h8=g>Mg1>b6gax-gT*vPu@k6_s=
zqZg|6a-{>0pxJ>gNcsL#q&CmP;#a=M$G4{l&%I$-{mn+CmJ}X4&5mHBe-*6q{S$GM
zS0Xd8xbQyaA!eT)gs5_7F|ocs(x)aM+wBoDx?VzZIpMv~#8Q~QV+>~A>Wzu-)I`(9
zCAf^zm-2map7u3r9CqRsUAWGjKfQ;mzIiMLb$f#BZ;NBgx<9btsw>ie8i^?-lMuN-
z5^IvOkm>#^GK$wky!&`=_l47H=Th2dled2Z@5yJlIp@v3rV#hDl_#&@)|RP*&@RcC
z>(g0yt#=0Nw|HV($-CG!W;}NESdCq;CS!ZEomkU77)$OqL+tW$h#PSk$(@HG>#eVl
z)3+(o+X?Rhp7g_lAIc(TeGulaZGl+>^Uz~z9rnVcl9bc*!$FQq9rYvk@MtH~#;Fxh
zqVG<mU@S5!r(^xX>ew+j6qz@_#B4qU12dZ-v8V7Dd%g!Lo#!H{ZUb~mAH%L+ahT%S
zJNo{1)RnE<qye`eb_uGt`~j)e+F@7w6l|Z{A8`Xdz?O=y;Ir^zSXzA$VwoOTayt_1
zh1cX;8qC1DPnu%eo-xQ7ITbmf4X|=SDJ)qKhDDFJVt!}~3?CZ7l|KKU=RCQyaPc26
z$FfSTk=p%BEIjr*GFlEsR-<BAxNQW2z79mEFE3)mxJd|WwifA^h4*7?mV>$3Miy)F
zNfS4=ZyI;H{dz7W;(uOG%k`+`eD_^=%d~NRF>dETC$6lg8{XLwhfNiGV4MF5tm2)K
zRDM6Qeq4zamj)uS^?ao6J%F6jQ?Th|3vBfK0qfgz!n%4FvAW(PWG0nI+8*I~tm|&X
ze8+MzyY{Q`lYLY#7K!L*Um_`XG}3#;Vf7cn{nAnR+&AhWMt<fE*KJ!-<8~Y(ei(|x
zQ!@~r%y7jw4EJq(`5rf+$qBAk`E+hfmAQDW-W)D)cqH2|?EifpphbgP4Y&HzaMb$r
zIQItI6Vs*(XO20Ioa8CUS?-R^>d%q7V*)m(jKa26ld#Rj8Jjv+Lr$;L$a?K_qzC+o
zqy*tPwEqCiyg40%zX(NZ|B7my;(hkw+MX!md=!KId%^dQpV0r6%~(8onebgl3v`;`
zhPv?;@X?WGNG$D&WiRG)O>^p!ezwc%YFx;CH#V(a7m8<!JiCN@yTP|yr;Pmf12nHk
zWcPhN8^ym~jpk|ISo%vIHhI>@j?ex?PG)OlhbLlf@;>C0c@0?)Mj^H3bfo@k7E%%d
zu{c+FEjlt0VILht=RRY(7r!^MRF5g<LV~HE@9U=9W%qKcb6>-V9<hl0{QyRvyN{NS
zdSc3qI>PU<hL|;E8!E(YVP|?*vy{Wudegddy|_gI2RLR&6RuR7?-BKRI%eMk#&?MZ
z$_n2N=K2YrmHvs&%Raye*JX(I8;sa5BQfh>8I1bz7<#t4$DIt&^u9_-D>rIn7<#qd
zjj_j8qT96&@ZFz@?oEKehF<VLvy$7EdX)Bd82*Q*$1tt8HJU=-o69!5*^eEzA&X#f
zwWH3U!GILBDZLj%_P>EfXQy-V<({$~1Bd#yKm1#V3SW40^EQ;iyK}!st2yc1(#;mX
zeiLiv9wiB{bIAM0&3)j#Fc#iV_Mz0=d#F8UGHTxF%_Vo9%Uw+$$nN&czjsXEUc4{v
zFU&btiR;Fe7+Kp?efVSFDXg>JTG$`oisEh$o5>v>^AC#6{d%$^JPwrMGHN{I=C$+U
zmb@Rxg~X<DAwkjH{(INCLBF=8wT={hombShhaS!~K2Zsm-5)lUE9rNg8{McDN}U`B
zeorm99vH%X)@Uw!#I2Q5pS|+F@XYhggOBgCTOKB_JDjPjzo#iPIn$KQj4~yMC$PJ_
z4&j1xo^xqF9^5Q;G*_qCG1K)1`M+1J(ZsF(@Ddj|rJcE6&a2KfE+?cZs+KK{CQgxP
zkr;*7n^ombq>Ny<ciC;)vWGKmX;YW2dh!L^_;fYX>WAx0)em$qhpD#bTeIeA)0?k$
zG+l`7XzEgBrm1S(U{l@Dou+BO_GSn6Okh9n^EMZMZZwzV*OxoHuQE3$vnscCVhl=j
z>44zZFJk2W90Wb+2cKpmxwUN<nSNzjo35-1Fs(eM`FWc4=(m=;a{Mc9Rmf3x?IVA4
zn%Ts?tGSi^HleLw7Ah<~&PA2H!S*VDz;t>1zu5Shq3kD5g4o~%|6~{cGKFn5s3*ZP
zee;p<H|ib#G>5CSXByIJ1h;O+L9|T!2fE@9;l1P-RQ>%ulyXWzwL4*G((x_4x%Un#
zmh*s#uMhua)zI^jH->wz!HBI1=r6o~ZJ$z|i|g5*9Z)sf)Y<i(soWdamH3qM?2^N6
z*j~5(^L?sZU5Dtw!CcJ=PH0r7Eb0e)pxE7;+{4LE+@xb!?5vC-Y|y3}tneo{?9}XQ
z?4X<d35IFe<<h22eW#hqX3UqvRLYq;KDuPuKI%i?oUd1~8F|a${%tmT*Yw5kBqxNQ
z=!=MTXVK)L4{&WPx&<~tn+^R?DQ_QlZ$eqPT}uH|rY#uHOzwGyQ{2dr*H}M~cy_=C
zr3fGE&t#h#HQs3|{mWP-el2<S$zh(Gb9)ZA^?VNZxLP18zasq2#r!Ua*w-IZm(9S(
zhZ-Z;WiE#HTLG6Py|{Phb!W#e%VXm{@60Aw_<(i)<gn>Nn78kVva>pTmG?i7vtq_&
zIAAznIAA!C??C%|FO9G^I;xFo=W$fsMNOBQ&yKpU<g;FSmGJ*AdA)pebR6l=)6w-x
z^%Zsf((@SWmp@U%2g3*AfN>uf_kp?l;H5jpePTRLEPI^L`-VotKg<04yZc4&8yf%a
ze&~Jvzl*Q*oHXj~<6i|&%^$iijenK%E1g4nUb<GQw|Bm<`IXK`by*sf@ct@!H9rdK
zoJu+hieHUKNuOF?4Y#Q8tDT?j+xwi9M=9Ey-%@-foiuMLUQw;5=hYhRbxt*YDKD*;
z<_mi5C^`!3-1PifqrG$%HosoyqUX>H&(Uz``D~4PpG&JAy?NE}s4lg-zso-L!_sIk
ze;h>zJs;Jhw5IvLYoz{YjrR6Qs}5<My>w9?d*K&VT<VKj)cbr=K8i2RTf3%sdPS)Y
zy>N>fj?@RbE~s-*K6^z=&p~x4_0fH4RI2~0<{jnyQXTeMYyCOe%CD8DsMb+`97R9X
zQ&3Uq+<NmX@f6j()E`GXhg4rtUptz=*6OkL=U^?LwS8K3s^M$Jv3H%CPDkI@%O6KS
zC)KCddF_Qmbvk-1>gT8Wm3*fAG-|D*^^OvibQgBs-u@KS`Rsj;!p1KsAGPunmLFR6
zsntpK(|T$Cud=4+lSX<zY2IFI%2&|Q-shqC1${0@=cB$T#iFWTs$1!Nmgc3n##-J1
z>Z|;&<7KVYOZ%*ij^~8Fimdg;+CF{NtD{MEs&~j&NOA1FmRDi_ti8`=FMPe7QPXKH
zpBk>M`&5?{hvw}arM!;5E@(fj<+HbsRF@JD-Pb!R@i}Uqo?l5H-7n~<q{mV61?`8y
zXE<OuV1EbNhYSCY4uzpNyO>9BY<uE~b?Z5oYF5LR))&@VO{d=b;&bYaP4L8*7@y?D
zYhzxt{h&W*<L4;<bEN(l_r1J7<h~JKMdd!CnhQ!V<)g8nc(w3IHPJPVS}^pgp`0{Y
zi)(AYl$YWuMQOb~*K!&x<&$Dq+Amix#gwk)dhC^xa?`aXe6&X#E$Nr56=TWG<$COu
zljEm3t(>%897!*|wBpFsQ4DdEtFu>5%qut7$}86^&S}L_suN?8xl*0pc`2{bS_`i&
zy>j>zPr9d+m)7aErhIBqjn7=aTCa3piY47I=(SXjbgfm7v`@O$ibLw8mnFS=*GjcW
z*XH=7wOZCnc%)anC&m)z3vw>jAU9XjBG;oOC&d-_&^>XE%xS(Lqf~=6I^;04<fP}+
z(qoN2DUT&<sU}N$#kEo#OV$*}k{(h^`|KITnr!6}V-?igR?TYtYUfbvm)lpv)O%k}
zgI+nUS}o;M!lgJ$`?lsOzFyH%J(luPeA-9X1raHyt<f5-mT+2wuUC&HI+bc|)uMz;
z^S0vK+Aq~X*IJ{k{G?xvOAiivVTiTclS7Ox<`mDxHOBnE4v0OX^Z)J|GwJ~KMjGjU
zK}V?$Eo)me5_~O}79F%k3!au5il=AP!Yi#I{8}!pI>lI8vBf>&xmHdsbz)q5%(ZBs
zHI!evrg<%q_F5aYaM)U7tv1?cE5FfiI8e|I(0f-d>a9i0p*OZY@x;3IoJ%#UVN2@^
zYpteJ?|tz(^~NT6;!BKA^5V5IH~jdk-vjiq@K-;lwdbS$7|#vL&jm%%BdXcj`@)(&
z%IRqL64gb|MAtOZd_hFYNu#x#w)RVTDV|c4*4uL}r@>M_DTbx}a`jS7=~}MGUO6c@
zU0cFOd&JR_esi^AOgcALZ;w67D;Kr$DfQC4Rve`|y>OK9&E+Y#QZ&~q-m|q|jAd?4
zdCb-5xkvfbq83gi3@yEKH55;7uOM?&hf>t4MX8_UwPMgZDHhEaR+Q=yuPxCe?vtBa
zgJnsNy=vvoLULO8lzPqOm2jn;Vhp7od*{U(r24dQOZ&B~<#5G)bYGk!bDA&6DAr&t
zuNcFUxs*?9ue8seYpTgs4vJ$<q}-PF(ms1eu_iru#CQcU*F&>847D?u>s8uQ!W7pl
z)mW1k<LW)vs;97cwd$jFg{{+4eyMJ{CXr%^b9&AxzOB(3O_p+6!=Xoyr8?#6ZPg)%
zDbCr7Z)-o*LUHI?YovYF#)9IpR-am*Sif4WHTx7_+)HyxQCw@x|NC4l=T~9*Bd0|z
zXNy*Xt#+yL%k8USYT1{=)gz~co7RZ=Xih1LYt`lygVfn_sqt8|Z;c);=deTz)j)gb
z+LB0n^^BCCV$!v>k>cnPt<_-k84ef@6s`kmkAFSYNxAgIwilLEvz}|QW-Ykl`ofrN
z(WrL~J*VEd6wluHl*fqX-}l(dAM*3H7j{v{HST@$_eN3kr=ah3>Q6!Qsy*(cdJ0pg
zUfj|%sd->-U#v+wH`iy6Jw5ov7<y<F<4Naw;9G(r<tJ-P`bn+yvb0}XFU1qDrM326
zi@C*fOL)b7WKMg<bDFaxD&bm!VN0!)U%a*iSKLSEq@P@t;K<cd3>xKX?2#jUYL`}C
zIUKFM(mH}A&5<>^<ns2)y<{~1uNH141TDSd8iGUC<Pzr#YyMJ+{J)lHl1F%%wG@Zc
zlgpAmd(^&EA^)!$AMK)Np=+h6hDY~lKV26_q#ER+7Cmx(YB@Pv(j(4`F~s?ToQpL`
z^=Z{1#n7^r!&d8)>#=7}tj87}F^(m3F{jo(agROcRF57U6ssVj9(vU<%$-B6S83lI
zrk;CB{CefJXtjn*4p)pLw{J^MjIGz)5<S*%5`1wV<s&(JU1C04=ay=;l+RLpJ?bsh
zs8nx@79~t^-WGgYdZik~Ypruz_(`u89?C&#9p#ehwudgtL9yvtE|MCOcf`wHy3FCq
zd1S7>FnzO^E;U>wkJS3?;gvZp=Ju@BLNV3)X`R-nhNrf#mB-RLHD2jH#h2!#Yg(rz
z(q39i*Oo-ut5+o41j`<mr5cF_agQaM#C_JzE#b4K)>>WSK6AXgJnq&vhMEU`ymxsR
z!>kAH%ZBVShItR%c6#qJhItQMwl&^m4D%k8O=_~s80I}F8QgG}G0b|9`($1lW2kuW
zFTY~B#xU=}X}q6n4D%l3B{j`8hItP*j&aL1hItQ`-mRZ&4D%k$G<D52hItS8mha{o
z!@LKz0-EL;!>k88PmY^l3>6Q0`Zn2V46`1X`;5}qd+9(+;T|%EiU-2;QA=Z(^}t8?
zon{Ob4}{l7KE^Qbf$$pJ#~9{45Z<f#7{j~=!h3!nW0?0q`0U8X80I|?K7;fzhItQ!
z&$4}tVcrAbdl4UFnDu}WzFRYfiU-1Xw2U#-J#g~8&p4%5VVrA}Wt<Z_G0xA%GUYz<
zV#*!xVait-#gzXzl_{UIiz$CAi*fPn$GG^NWL(CVVO-*#FfOYy7?(}XjPN%DFC$r#
zJlQ98LOtor?<Y6{9>FE}geRX*;3fP-2l;t9nx1@3TNuW9`y!^?s!L4y=|3?pTyMtZ
zkA_TzWlfoi4H%~4y3&km%>c%A#^;Rd_m>#AYAqSJ?%s^sjJk|lN*~5;(=$f++mn}(
ztVy2i3w5MEzfb5VI0TR25`4l#_y{lIcM^0sq8_5jxm{(Z+_r6uOY!4O1+SG%#dmHq
zu2mi~u2(;0+?IA{Dz)%qD(&sdRBr0bR9>`$seFAL<6iGg#yv2SaSwZgaZe0nguh{W
z8OfUDg?&<&UoZ5Leu6{r2rj`VJcN(%5`Ln?5%uJ2%G}J9AAN&yd0dC7=zorJo#w*0
zh0kIt4SK><emj_{e4-KKKIksve&Ie-<^9S`l@;fiDwo3<kID}ik5;!B5C3b7@V9U;
zBUuZ1vY%fk)RR8aPjCnx!6o>Fhwu?z!tX5Ta6~=%nvT9;%5Uk!ROs%;RNQ@(alL<@
zaeGmPsdTdpQ+ab+#yw~Z<NhFrsS-SisdDBt<KaDx@tCxg@mM>D@i?B#cw9Tlc-*c3
zGLkjPlYOC%)RR8aPjCnx!6o>Fhwu?z!tX5Ta6mo&9nT~gL)in~-8aY>W<2=h?WMdi
zls#BqsR?fkWe?7=19)R7dJs_Z$_`_g@t|7T$pB+0d+<(^4FSeb_P{i{d4MsLJ?PQs
zc7QRIJs8Zl3owSV2jh3{2{4AT2Q$z!z!=IN%)a+cfH9Okh;Dr=z!=IN#Dr!B7(>~E
zxHTICjG^p7{DFf3#xOjv_8u@ETgs0u;*Ig#to+<eyfJ>aEB|gM-Wab9m0ueYZ;baA
z%I__RH^zHY<@ct<8{@ML<<B;VH^ygU%Abu9Z;a1Yl|Nf0-WZ?FD}OdmyfMDpQT}d+
zcw>A&X!iTjZaXi$<sUh#n*VvH&-p4f-r(B?HRT63+R4vkI`RoI*ZGurAM!a3%kVo!
z=kWXeT=+wIr}*QiSpJ6@Zv2^j!uPU?as2t`ZSx12lN{M2HKdmGkY0j8um~o>CLDxI
z;3V8c1JUAedQQy<@=t2`mH+*DS-j_oXS~V#0Y4xhnGZ?(iI3X1i(i<wm|tEqlV8~<
zkIxP+!LOh81HX9@$8Y<pKEJb01fRQc3%_gPH~E9iNsjE18d6JoNH4)4SOk+`6Ar>9
za1w5!p_`zE42RP5X!$w+HQSo;#RnwvZ{FR)vp+xNKO9+{51ty&hmWepM^^RVW1?dD
zgop3)i{}2xe>$uIpR^@_Uq1B$zv8P0d~(KKKBdxx{6Xd<NA^e!sU<z6mtYVqf=RFm
z2jL=|0yojnP0&JyBkA#ft^FbY9pAO!ONUS3y@ucC+3CIc-g8RvLoWTyk53EbLzf)l
zXEx8}BaQ|0bB@2mM-FK#+ylM%=%-Qq{1(mlnC_SOSm*5gLFOb!_DBt>B|W5<U=S>V
zNw5hA;Ub&@H_?!<#nJRs=sw<m-#csgvN>z{`d3HurtO3I9*E@!dL7|M#QO2${;`Lj
zT>d%#Nyj~WXt(+Nw7|>!wCJDs&|4+>kb!6U$xFTXq+4V22bq%`*&{Wimh_Nbf<dqd
YCc!2ggo|(r+(d&v(c)lw=zkCXf5`HMW&i*H

literal 0
HcmV?d00001

diff --git a/mdtf_test_data/synthetic/__init__.py b/mdtf_test_data/synthetic/__init__.py
index 37f6993..2e732ed 100755
--- a/mdtf_test_data/synthetic/__init__.py
+++ b/mdtf_test_data/synthetic/__init__.py
@@ -8,3 +8,4 @@
 
 from . import time
 from . import vertical
+from . import horizontal
diff --git a/mdtf_test_data/synthetic/horizontal/__init__.py b/mdtf_test_data/synthetic/horizontal/__init__.py
new file mode 100644
index 0000000..04662f9
--- /dev/null
+++ b/mdtf_test_data/synthetic/horizontal/__init__.py
@@ -0,0 +1,2 @@
+from .construct_rect_grid import construct_rect_grid
+from .construct_tripolar_grid import construct_tripolar_grid
diff --git a/mdtf_test_data/synthetic/horizontal/construct_rect_grid.py b/mdtf_test_data/synthetic/horizontal/construct_rect_grid.py
new file mode 100755
index 0000000..29649b8
--- /dev/null
+++ b/mdtf_test_data/synthetic/horizontal/construct_rect_grid.py
@@ -0,0 +1,145 @@
+""" Collection of tools to coarsen model data """
+
+__all__ = ["construct_rect_grid"]
+
+import warnings
+
+import numpy as np
+import xarray as xr
+
+try:
+    import xesmf as xe
+except:
+    warnings.warn("Unable to load `xesmf`. Regridding functionality will not work.")
+
+
+def construct_rect_grid(dlon, dlat, add_attrs=False, attr_fmt="ncar", bounds=False):
+    """Generate a rectilinear grid based on values of dx and dy
+
+    Parameters
+    ----------
+    dlon : float
+        Grid spacing in the x-dimension (longitude)
+    dlat : float
+        Grid spacing in the y-dimension (latitude)
+    add_attrs : bool, optional
+        Include lat and lon variable attributes, by default False
+
+    Returns
+    -------
+    xarray.Dataset
+        Empty shell dataset with lat and lon dimensions
+    """
+
+    if 180.0 % dlat != 0:
+        dlat = 180.0 / np.floor(180.0 / dlat)
+        warnings.warn(
+            f"180 degrees does not divide evenly by dlat. Adjusting dlat to {dlat}"
+        )
+
+    if 360.0 % dlon != 0:
+        dlon = 360.0 / np.floor(360.0 / dlon)
+        warnings.warn(
+            f"360 degrees does not divide evenly by dlon. Adjusting dlon to {dlon}"
+        )
+
+    lat = np.arange(-90.0 + (dlat / 2.0), 90.0, dlat)
+    lon = np.arange(0.0 + (dlon / 2.0), 360.0, dlon)
+
+    dset = xr.Dataset({"lat": (["lat"], lat), "lon": (["lon"], lon)})
+
+    if bounds:
+        lat_bnds = np.arange(-90.0, 90.0 + (dlat / 2.0), dlat)
+        lon_bnds = np.arange(0.0, 360.0 + (dlon / 2.0), dlon)
+
+        lat_bnds = np.array(list(zip(lat_bnds[0:-1], lat_bnds[1::])))
+        lon_bnds = np.array(list(zip(lon_bnds[0:-1], lon_bnds[1::])))
+
+        bnds = np.array([0.0, 1.0])
+
+        dset["lat_bnds"] = xr.DataArray(lat_bnds, coords=(dset.lat, ("bnds", bnds)))
+        dset["lon_bnds"] = xr.DataArray(lon_bnds, coords=(dset.lon, ("bnds", bnds)))
+
+        dset = dset.drop("bnds")
+
+    if attr_fmt == "ncar":
+        dset["lat"].attrs = (
+            {"long_name": "latitude", "units": "degrees_north"} if add_attrs else {}
+        )
+        dset["lon"].attrs = (
+            {"long_name": "longitude", "units": "degrees_east"} if add_attrs else {}
+        )
+
+    elif attr_fmt == "gfdl":
+        dset["lat"].attrs = (
+            {"long_name": "latitude", "units": "degrees_N", "cartesian_axis": "Y"}
+            if add_attrs
+            else {}
+        )
+        dset["lon"].attrs = (
+            {"long_name": "longitude", "units": "degrees_E", "cartesian_axis": "X"}
+            if add_attrs
+            else {}
+        )
+
+        if bounds:
+            dset["lat"].attrs["bounds"] = "lat_bnds"
+            dset["lon"].attrs["bounds"] = "lon_bnds"
+
+            dset["lat_bnds"].attrs = (
+                {"long_name": "latitude bounds", "cartesian_axis": "Y"}
+                if add_attrs
+                else {}
+            )
+            dset["lon_bnds"].attrs = (
+                {"long_name": "longitude bounds", "cartesian_axis": "X"}
+                if add_attrs
+                else {}
+            )
+
+    elif attr_fmt == "cmip":
+        dset["lat"].attrs = (
+            {
+                "long_name": "latitude",
+                "units": "degrees_north",
+                "axis": "Y",
+                "standard_name": "latitude",
+                "cell_methods": "time: point",
+            }
+            if add_attrs
+            else {}
+        )
+        dset["lon"].attrs = (
+            {
+                "long_name": "longitude",
+                "units": "degrees_east",
+                "axis": "X",
+                "standard_name": "longitude",
+                "cell_methods": "time: point",
+            }
+            if add_attrs
+            else {}
+        )
+
+        if bounds:
+            dset["bnds"] = xr.DataArray(bnds, dims={"bnds": bnds})
+            dset["bnds"].attrs["long_name"] = "vertex number"
+
+            dset["lat"].attrs["bounds"] = "lat_bnds"
+            dset["lon"].attrs["bounds"] = "lon_bnds"
+
+            dset["lat_bnds"].attrs = (
+                {"long_name": "latitude bounds", "axis": "Y", "units": "degrees_north"}
+                if add_attrs
+                else {}
+            )
+            dset["lon_bnds"].attrs = (
+                {"long_name": "longitude bounds", "axis": "X", "units": "degrees_east"}
+                if add_attrs
+                else {}
+            )
+
+    else:
+        raise ValueError("Unknown model attribute format")
+
+    return dset
diff --git a/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py b/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py
new file mode 100755
index 0000000..ab6c757
--- /dev/null
+++ b/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py
@@ -0,0 +1,87 @@
+""" Collection of tools to coarsen model data """
+
+__all__ = ["construct_tripolar_grid"]
+
+import warnings
+
+import numpy as np
+import xarray as xr
+import pkg_resources as pkgr
+
+try:
+    import xesmf as xe
+except:
+    warnings.warn("Unable to load `xesmf`. Regridding functionality will not work.")
+
+
+def construct_tripolar_grid(
+    point_type="t", add_attrs=False, attr_fmt="gfdl", retain_coords=False
+):
+    """Generate a tripolar grid based on a real 5-degree MOM6 configuration
+
+    Parameters
+    ----------
+    point_type : str, optional
+        Grid type. Options are t, u, v, c. By default "t" (tracer)
+    add_attrs : bool, optional
+        Include lat and lon variable attributes, by default False
+    attr_fmt : str, optional
+        Modeling center attribute format, by default "gfdl"
+    retain_coords : bool, optional
+        Keep geolon, geolat, and wet in the dataset, by default False
+
+    Returns
+    -------
+    xarray.Dataset
+        Shell dataset with masked variable and ocean depth field
+    """
+
+    ds_in = pkgr.resource_filename("mdtf_test_data", "resources/ocean_static_5deg.nc")
+    ds_in = xr.open_dataset(ds_in)
+
+    if point_type == "t":
+        lat = ds_in["geolat"]
+        lon = ds_in["geolon"]
+        wet = ds_in["wet"]
+    elif point_type == "u":
+        lat = ds_in["geolat_u"]
+        lon = ds_in["geolon_u"]
+        wet = ds_in["wet_u"]
+    elif point_type == "v":
+        lat = ds_in["geolat_v"]
+        lon = ds_in["geolon_v"]
+        wet = ds_in["wet_v"]
+    elif point_type == "c":
+        lat = ds_in["geolat_c"]
+        lon = ds_in["geolon_c"]
+        wet = ds_in["wet_c"]
+
+    dset = xr.Dataset()
+    dset["mask"] = xr.where(wet == 0.0, np.nan, 1.0)
+
+    if point_type == "t":
+        dset["depth"] = xr.DataArray(ds_in["depth_ocean"])
+
+    if retain_coords is True:
+        dset[lat.name] = lat
+        dset[lon.name] = lat
+        dset[wet.name] = wet
+
+        if add_attrs is True:
+            if attr_fmt == "gfdl":
+                dset[lat.name].attrs = {}
+                dset[lon.name].attrs = {}
+                dset[wet.name].attrs = {}
+            elif attr_fmt == "ncar":
+                dset[lat.name].attrs = {}
+                dset[lon.name].attrs = {}
+                dset[wet.name].attrs = {}
+            else:
+                raise ValueError("Unknown attribute format")
+
+        else:
+            dset[lat.name].attrs = {}
+            dset[lon.name].attrs = {}
+            dset[wet.name].attrs = {}
+
+    return dset
diff --git a/mdtf_test_data/synthetic/synthetic_data.py b/mdtf_test_data/synthetic/synthetic_data.py
index 93c7117..3f853c0 100755
--- a/mdtf_test_data/synthetic/synthetic_data.py
+++ b/mdtf_test_data/synthetic/synthetic_data.py
@@ -11,7 +11,7 @@
 import cftime
 import xarray as xr
 import numpy as np
-from mdtf_test_data.util.rectilinear import construct_rect_grid
+from mdtf_test_data.synthetic.horizontal import construct_rect_grid
 import mdtf_test_data.generators as generators
 
 from mdtf_test_data.synthetic.time import generate_monthly_time_axis
diff --git a/mdtf_test_data/synthetic/vertical/__init__.py b/mdtf_test_data/synthetic/vertical/__init__.py
index fca56a5..47a9b2a 100644
--- a/mdtf_test_data/synthetic/vertical/__init__.py
+++ b/mdtf_test_data/synthetic/vertical/__init__.py
@@ -1,3 +1,5 @@
 from .gfdl_plev19_vertical_coord import gfdl_plev19_vertical_coord
 from .gfdl_vertical_coord import gfdl_vertical_coord
+from .mom6_z_coord import mom6_z_coord
+from .mom6_rho2_coord import mom6_rho2_coord
 from .ncar_hybrid_coord import ncar_hybrid_coord
diff --git a/mdtf_test_data/synthetic/vertical/gfdl_plev19_vertical_coord.py b/mdtf_test_data/synthetic/vertical/gfdl_plev19_vertical_coord.py
index 4d2e605..50fc800 100755
--- a/mdtf_test_data/synthetic/vertical/gfdl_plev19_vertical_coord.py
+++ b/mdtf_test_data/synthetic/vertical/gfdl_plev19_vertical_coord.py
@@ -7,7 +7,6 @@
 import cftime
 import xarray as xr
 import numpy as np
-from mdtf_test_data.util.rectilinear import construct_rect_grid
 import mdtf_test_data.generators as generators
 
 
diff --git a/mdtf_test_data/synthetic/vertical/gfdl_vertical_coord.py b/mdtf_test_data/synthetic/vertical/gfdl_vertical_coord.py
index 747bc6f..ebcd5d1 100755
--- a/mdtf_test_data/synthetic/vertical/gfdl_vertical_coord.py
+++ b/mdtf_test_data/synthetic/vertical/gfdl_vertical_coord.py
@@ -7,7 +7,6 @@
 import cftime
 import xarray as xr
 import numpy as np
-from mdtf_test_data.util.rectilinear import construct_rect_grid
 import mdtf_test_data.generators as generators
 
 
diff --git a/mdtf_test_data/synthetic/vertical/mom6_rho2_coord.py b/mdtf_test_data/synthetic/vertical/mom6_rho2_coord.py
new file mode 100755
index 0000000..f5379e7
--- /dev/null
+++ b/mdtf_test_data/synthetic/vertical/mom6_rho2_coord.py
@@ -0,0 +1,126 @@
+""" Module for generating synthetic datasets """
+
+___all__ = [
+    "mom6_rho2_coord",
+]
+
+import cftime
+import xarray as xr
+import numpy as np
+import mdtf_test_data.generators as generators
+
+from mdtf_test_data.synthetic.time import generate_monthly_time_axis
+from mdtf_test_data.synthetic.time import generate_daily_time_axis
+from mdtf_test_data.synthetic.time import generate_hourly_time_axis
+
+
+def mom6_rho2_coord():
+    """Generates GFDL AM4 pressure coordinate
+
+    Returns
+    -------
+    xarray.DataArray
+        GFDL AM4 pressure levels and half levels
+    """
+
+    rho2_l = [
+        1013.75,
+        1028.5,
+        1029.2421875,
+        1029.71875,
+        1030.1796875,
+        1030.625,
+        1031.0546875,
+        1031.46875,
+        1031.8671875,
+        1032.25,
+        1032.6171875,
+        1032.96875,
+        1033.3046875,
+        1033.625,
+        1033.9296875,
+        1034.21875,
+        1034.4921875,
+        1034.75,
+        1034.9921875,
+        1035.21875,
+        1035.4296875,
+        1035.625,
+        1035.8046875,
+        1035.96875,
+        1036.1171875,
+        1036.25,
+        1036.375,
+        1036.5,
+        1036.625,
+        1036.75,
+        1036.875,
+        1037,
+        1037.125,
+        1037.25,
+        1037.65625,
+    ]
+
+    rho2_i = [
+        999.5,
+        1028,
+        1029,
+        1029.484375,
+        1029.953125,
+        1030.40625,
+        1030.84375,
+        1031.265625,
+        1031.671875,
+        1032.0625,
+        1032.4375,
+        1032.796875,
+        1033.140625,
+        1033.46875,
+        1033.78125,
+        1034.078125,
+        1034.359375,
+        1034.625,
+        1034.875,
+        1035.109375,
+        1035.328125,
+        1035.53125,
+        1035.71875,
+        1035.890625,
+        1036.046875,
+        1036.1875,
+        1036.3125,
+        1036.4375,
+        1036.5625,
+        1036.6875,
+        1036.8125,
+        1036.9375,
+        1037.0625,
+        1037.1875,
+        1037.3125,
+        1038,
+    ]
+
+    rho2_i_attrs = {
+        "long_name": "Target Potential Density at interface",
+        "units": "kg m-3",
+        "axis": "Z",
+        "positive": "down",
+    }
+
+    rho2_l_attrs = {
+        "long_name": "Target Potential Density at cell center",
+        "units": "kg m-3",
+        "axis": "Z",
+        "positive": "down",
+        "edges": "rho2_i",
+    }
+
+    dset_out = xr.Dataset()
+    dset_out["rho2_l"] = xr.DataArray(
+        rho2_l, dims={"rho2_l": rho2_l}, coords={"rho2_l": rho2_l}, attrs=rho2_l_attrs
+    )
+    dset_out["rho2_i"] = xr.DataArray(
+        rho2_i, dims={"rho2_i": rho2_i}, coords={"rho2_i": rho2_i}, attrs=rho2_i_attrs
+    )
+
+    return dset_out
diff --git a/mdtf_test_data/synthetic/vertical/mom6_z_coord.py b/mdtf_test_data/synthetic/vertical/mom6_z_coord.py
new file mode 100755
index 0000000..34d7cd4
--- /dev/null
+++ b/mdtf_test_data/synthetic/vertical/mom6_z_coord.py
@@ -0,0 +1,126 @@
+""" Module for generating synthetic datasets """
+
+___all__ = [
+    "mom6_z_coord",
+]
+
+import cftime
+import xarray as xr
+import numpy as np
+import mdtf_test_data.generators as generators
+
+from mdtf_test_data.synthetic.time import generate_monthly_time_axis
+from mdtf_test_data.synthetic.time import generate_daily_time_axis
+from mdtf_test_data.synthetic.time import generate_hourly_time_axis
+
+
+def mom6_z_coord():
+    """Generates GFDL AM4 pressure coordinate
+
+    Returns
+    -------
+    xarray.DataArray
+        GFDL AM4 pressure levels and half levels
+    """
+
+    z_l = [
+        2.5,
+        10,
+        20,
+        32.5,
+        51.25,
+        75,
+        100,
+        125,
+        156.25,
+        200,
+        250,
+        312.5,
+        400,
+        500,
+        600,
+        700,
+        800,
+        900,
+        1000,
+        1100,
+        1200,
+        1300,
+        1400,
+        1537.5,
+        1750,
+        2062.5,
+        2500,
+        3000,
+        3500,
+        4000,
+        4500,
+        5000,
+        5500,
+        6000,
+        6500,
+    ]
+
+    z_i = [
+        0,
+        5,
+        15,
+        25,
+        40,
+        62.5,
+        87.5,
+        112.5,
+        137.5,
+        175,
+        225,
+        275,
+        350,
+        450,
+        550,
+        650,
+        750,
+        850,
+        950,
+        1050,
+        1150,
+        1250,
+        1350,
+        1450,
+        1625,
+        1875,
+        2250,
+        2750,
+        3250,
+        3750,
+        4250,
+        4750,
+        5250,
+        5750,
+        6250,
+        6750,
+    ]
+
+    z_i_attrs = {
+        "long_name": "Depth at interface",
+        "units": "meters",
+        "axis": "Z",
+        "positive": "down",
+    }
+
+    z_l_attrs = {
+        "long_name": "Depth at cell center",
+        "units": "meters",
+        "axis": "Z",
+        "positive": "down",
+        "edges": "z_i",
+    }
+
+    dset_out = xr.Dataset()
+    dset_out["z_l"] = xr.DataArray(
+        z_l, dims={"z_l": z_l}, coords={"z_l": z_l}, attrs=z_l_attrs
+    )
+    dset_out["z_i"] = xr.DataArray(
+        z_i, dims={"z_i": z_i}, coords={"z_i": z_i}, attrs=z_i_attrs
+    )
+
+    return dset_out
diff --git a/mdtf_test_data/synthetic/vertical/ncar_hybrid_coord.py b/mdtf_test_data/synthetic/vertical/ncar_hybrid_coord.py
index 35d866f..00b9d08 100755
--- a/mdtf_test_data/synthetic/vertical/ncar_hybrid_coord.py
+++ b/mdtf_test_data/synthetic/vertical/ncar_hybrid_coord.py
@@ -7,7 +7,6 @@
 import cftime
 import xarray as xr
 import numpy as np
-from mdtf_test_data.util.rectilinear import construct_rect_grid
 import mdtf_test_data.generators as generators
 
 
diff --git a/mdtf_test_data/tests/test_synthetic_data.py b/mdtf_test_data/tests/test_synthetic_data.py
index 0ed9a63..74b24ae 100644
--- a/mdtf_test_data/tests/test_synthetic_data.py
+++ b/mdtf_test_data/tests/test_synthetic_data.py
@@ -18,6 +18,8 @@
 from mdtf_test_data.synthetic.vertical import gfdl_plev19_vertical_coord
 from mdtf_test_data.synthetic.vertical import gfdl_vertical_coord
 from mdtf_test_data.synthetic.vertical import ncar_hybrid_coord
+from mdtf_test_data.synthetic.vertical import mom6_z_coord
+from mdtf_test_data.synthetic.vertical import mom6_rho2_coord
 
 __all__ = [
     "test_xr_times_from_tuples_ncar",
@@ -25,6 +27,7 @@
     "test_ncar_hybrid_coord",
     "test_gfdl_plev19_vertical_coord",
     "test_gfdl_vertical_coord",
+    "test_mom6_z_coord",
     "test_generate_random_array",
     "test_generate_daily_time_axis",
     "test_generate_hourly_time_axis",
@@ -104,6 +107,20 @@ def test_gfdl_vertical_coord():
     pytest.gfdl_vert = result
 
 
+def test_mom6_z_coord():
+    result = mom6_z_coord()
+    assert isinstance(result, xr.Dataset)
+    assert np.allclose(float(result.z_l.sum()), 57085.0)
+    assert np.allclose(float(result.z_i.sum()), 60460.0)
+
+
+def test_mom6_rho2_coord():
+    result = mom6_rho2_coord()
+    assert isinstance(result, xr.Dataset)
+    assert np.allclose(float(result.rho2_l.sum()), 36176.125)
+    assert np.allclose(float(result.rho2_i.sum()), 37194.875)
+
+
 def test_generate_daily_time_axis():
     result = generate_daily_time_axis(1850, 2)
     assert isinstance(result, xr.Dataset)
diff --git a/mdtf_test_data/util/__init__.py b/mdtf_test_data/util/__init__.py
index 048bdc9..2c904c9 100755
--- a/mdtf_test_data/util/__init__.py
+++ b/mdtf_test_data/util/__init__.py
@@ -1,4 +1,2 @@
 """ init file for util directory """
-from .rectilinear import (regrid_lat_lon_dataset,
-    construct_rect_grid
-)
\ No newline at end of file
+from .rectilinear import regrid_lat_lon_dataset
diff --git a/mdtf_test_data/util/rectilinear.py b/mdtf_test_data/util/rectilinear.py
index 4634d6c..be1d028 100755
--- a/mdtf_test_data/util/rectilinear.py
+++ b/mdtf_test_data/util/rectilinear.py
@@ -1,151 +1,21 @@
 #!/usr/bin/env python
 """ Collection of tools to coarsen model data """
 
-__all__ = ["construct_rect_grid", "regrid_lat_lon_dataset"]
+__all__ = ["regrid_lat_lon_dataset"]
 
 import warnings
 
 import numpy as np
 import xarray as xr
 
+from mdtf_test_data.synthetic.horizontal import construct_rect_grid
+
 try:
     import xesmf as xe
 except:
     warnings.warn("Unable to load `xesmf`. Regridding functionality will not work.")
 
 
-def construct_rect_grid(dlon, dlat, add_attrs=False, attr_fmt="ncar", bounds=False):
-    """Generate a rectilinear grid based on values of dx and dy
-
-    Parameters
-    ----------
-    dlon : float
-        Grid spacing in the x-dimension (longitude)
-    dlat : float
-        Grid spacing in the y-dimension (latitude)
-    add_attrs : bool, optional
-        Include lat and lon variable attributes, by default False
-
-    Returns
-    -------
-    xarray.Dataset
-        Empty shell dataset with lat and lon dimensions
-    """
-
-    if 180.0 % dlat != 0:
-        dlat = 180.0 / np.floor(180.0 / dlat)
-        warnings.warn(
-            f"180 degrees does not divide evenly by dlat. Adjusting dlat to {dlat}"
-        )
-
-    if 360.0 % dlon != 0:
-        dlon = 360.0 / np.floor(360.0 / dlon)
-        warnings.warn(
-            f"360 degrees does not divide evenly by dlon. Adjusting dlon to {dlon}"
-        )
-
-    lat = np.arange(-90.0 + (dlat / 2.0), 90.0, dlat)
-    lon = np.arange(0.0 + (dlon / 2.0), 360.0, dlon)
-
-    dset = xr.Dataset({"lat": (["lat"], lat), "lon": (["lon"], lon)})
-
-    if bounds:
-        lat_bnds = np.arange(-90.0, 90.0 + (dlat / 2.0), dlat)
-        lon_bnds = np.arange(0.0, 360.0 + (dlon / 2.0), dlon)
-
-        lat_bnds = np.array(list(zip(lat_bnds[0:-1], lat_bnds[1::])))
-        lon_bnds = np.array(list(zip(lon_bnds[0:-1], lon_bnds[1::])))
-
-        bnds = np.array([0.0, 1.0])
-
-        dset["lat_bnds"] = xr.DataArray(lat_bnds, coords=(dset.lat, ("bnds", bnds)))
-        dset["lon_bnds"] = xr.DataArray(lon_bnds, coords=(dset.lon, ("bnds", bnds)))
-
-        dset = dset.drop("bnds")
-
-    if attr_fmt == "ncar":
-        dset["lat"].attrs = (
-            {"long_name": "latitude", "units": "degrees_north"} if add_attrs else {}
-        )
-        dset["lon"].attrs = (
-            {"long_name": "longitude", "units": "degrees_east"} if add_attrs else {}
-        )
-
-    elif attr_fmt == "gfdl":
-        dset["lat"].attrs = (
-            {"long_name": "latitude", "units": "degrees_N", "cartesian_axis": "Y"}
-            if add_attrs
-            else {}
-        )
-        dset["lon"].attrs = (
-            {"long_name": "longitude", "units": "degrees_E", "cartesian_axis": "X"}
-            if add_attrs
-            else {}
-        )
-
-        if bounds:
-            dset["lat"].attrs["bounds"] = "lat_bnds"
-            dset["lon"].attrs["bounds"] = "lon_bnds"
-
-            dset["lat_bnds"].attrs = (
-                {"long_name": "latitude bounds", "cartesian_axis": "Y"}
-                if add_attrs
-                else {}
-            )
-            dset["lon_bnds"].attrs = (
-                {"long_name": "longitude bounds", "cartesian_axis": "X"}
-                if add_attrs
-                else {}
-            )
-
-    elif attr_fmt == "cmip":
-        dset["lat"].attrs = (
-            {
-                "long_name": "latitude",
-                "units": "degrees_north",
-                "axis": "Y",
-                "standard_name": "latitude",
-                "cell_methods": "time: point",
-            }
-            if add_attrs
-            else {}
-        )
-        dset["lon"].attrs = (
-            {
-                "long_name": "longitude",
-                "units": "degrees_east",
-                "axis": "X",
-                "standard_name": "longitude",
-                "cell_methods": "time: point",
-            }
-            if add_attrs
-            else {}
-        )
-
-        if bounds:
-            dset["bnds"] = xr.DataArray(bnds, dims={"bnds": bnds})
-            dset["bnds"].attrs["long_name"] = "vertex number"
-
-            dset["lat"].attrs["bounds"] = "lat_bnds"
-            dset["lon"].attrs["bounds"] = "lon_bnds"
-
-            dset["lat_bnds"].attrs = (
-                {"long_name": "latitude bounds", "axis": "Y", "units": "degrees_north"}
-                if add_attrs
-                else {}
-            )
-            dset["lon_bnds"].attrs = (
-                {"long_name": "longitude bounds", "axis": "X", "units": "degrees_east"}
-                if add_attrs
-                else {}
-            )
-
-    else:
-        raise ValueError("Unknown model attribute format")
-
-    return dset
-
-
 def regrid_lat_lon_dataset(dset, dlon=10.0, dlat=10.0, method="bilinear"):
     """Regrids xarray dataset to a standard lat-lon grid
 
diff --git a/setup.cfg b/setup.cfg
index f0bc74d..195c291 100644
--- a/setup.cfg
+++ b/setup.cfg
@@ -49,3 +49,4 @@ install_requires =
 mdtf_test_data =
     config/*
     tests/*
+    resources/*

From 0420fc343552a665a3b251a97bdc02773a5819e3 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Sun, 10 Oct 2021 16:09:48 -0400
Subject: [PATCH 02/17] NCAR tripolar ocean lat/lon bounds

---
 .../horizontal/construct_tripolar_grid.py     | 41 +++++++++++++++++--
 1 file changed, 38 insertions(+), 3 deletions(-)

diff --git a/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py b/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py
index ab6c757..ea5a956 100755
--- a/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py
+++ b/mdtf_test_data/synthetic/horizontal/construct_tripolar_grid.py
@@ -64,7 +64,7 @@ def construct_tripolar_grid(
 
     if retain_coords is True:
         dset[lat.name] = lat
-        dset[lon.name] = lat
+        dset[lon.name] = lon
         dset[wet.name] = wet
 
         if add_attrs is True:
@@ -73,8 +73,24 @@ def construct_tripolar_grid(
                 dset[lon.name].attrs = {}
                 dset[wet.name].attrs = {}
             elif attr_fmt == "ncar":
-                dset[lat.name].attrs = {}
-                dset[lon.name].attrs = {}
+                dset[lat.name].attrs = {
+                    "axis": "Y",
+                    "standard_name": "latitude",
+                    "title": "Latitude",
+                    "type": "double",
+                    "units": "degrees_north",
+                    "valid_max": 90.0,
+                    "valid_min": -90.0,
+                }
+                dset[lon.name].attrs = {
+                    "axis": "X",
+                    "standard_name": "longitude",
+                    "title": "Longitude",
+                    "type": "double",
+                    "units": "degrees_east",
+                    "valid_max": 360.0,
+                    "valid_min": 0.0,
+                }
                 dset[wet.name].attrs = {}
             else:
                 raise ValueError("Unknown attribute format")
@@ -84,4 +100,23 @@ def construct_tripolar_grid(
             dset[lon.name].attrs = {}
             dset[wet.name].attrs = {}
 
+    if attr_fmt == "ncar":
+        dset = dset.rename({"xh": "nlon", "yh": "nlat"})
+
+        lat_range = np.array(np.arange(1, len(dset["nlat"]) + 1, dtype=np.intc))
+        dset["nlat"] = xr.DataArray(lat_range, dims=("nlat"))
+        dset["nlat"].attrs = {
+            "long_name": "cell index along second dimension",
+            "units": "1",
+        }
+
+        lon_range = np.array(np.arange(1, len(dset["nlon"]) + 1, dtype=np.intc))
+        dset["nlon"] = xr.DataArray(lon_range, dims=("nlon"))
+        dset["nlon"].attrs = {
+            "long_name": "cell index along first dimension",
+            "units": "1",
+        }
+
+        dset = dset.rename({lat.name: "lat", lon.name: "lon"})
+
     return dset

From f0b9de7273df04ca4dc5f913eb77227dec523809 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Sun, 10 Oct 2021 16:12:07 -0400
Subject: [PATCH 03/17] Specify accompanying static file for synthetic data

---
 mdtf_test_data/synthetic/synthetic_setup.py | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/mdtf_test_data/synthetic/synthetic_setup.py b/mdtf_test_data/synthetic/synthetic_setup.py
index 1693c0f..fda6aac 100755
--- a/mdtf_test_data/synthetic/synthetic_setup.py
+++ b/mdtf_test_data/synthetic/synthetic_setup.py
@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 
+import xarray as xr
+
 __all__ = ["create_output_dirs", "synthetic_main"]
 """ Script to generate synthetic GFDL CM4 output """
 import os
@@ -67,6 +69,11 @@ def synthetic_main(
     # -- Create Data
     print("Generating data with time resolution of ", TIME_RES)
     for v in var_names:
+        static = (
+            yaml_dict[v + ".static"]
+            if str(v + ".static") in list(yaml_dict.keys())
+            else False
+        )
         stats = (
             yaml_dict[v + ".stats"]
             if str(v + ".stats") in list(yaml_dict.keys())
@@ -91,6 +98,12 @@ def synthetic_main(
             else None
         )
 
+        if str(v + ".source") in list(yaml_dict.keys()):
+            _ds = xr.open_dataset(yaml_dict[v + ".source.filename"])
+            data = _ds[yaml_dict[v + ".source.variable"]].values
+        else:
+            data = None
+
         dset_out = generate_synthetic_dataset(
             DLON,
             DLAT,
@@ -102,7 +115,9 @@ def synthetic_main(
             fmt=DATA_FORMAT,
             generator=generator,
             stats=stats,
+            static=static,
             coords=coords,
+            data=data,
             generator_kwargs=generator_kwargs,
         )
 

From 88e362883be41ac9c66a29abec277cdf0d1c02f1 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Mon, 11 Oct 2021 11:35:19 -0400
Subject: [PATCH 04/17] Use default static file source

- If not specified in the yaml file, use the default
  ocean static file and the areacello variable
---
 mdtf_test_data/synthetic/synthetic_setup.py | 26 ++++++++++++++++++---
 1 file changed, 23 insertions(+), 3 deletions(-)

diff --git a/mdtf_test_data/synthetic/synthetic_setup.py b/mdtf_test_data/synthetic/synthetic_setup.py
index fda6aac..5d27618 100755
--- a/mdtf_test_data/synthetic/synthetic_setup.py
+++ b/mdtf_test_data/synthetic/synthetic_setup.py
@@ -1,6 +1,8 @@
 #!/usr/bin/env python
 
+import warnings
 import xarray as xr
+import pkg_resources as pkgr
 
 __all__ = ["create_output_dirs", "synthetic_main"]
 """ Script to generate synthetic GFDL CM4 output """
@@ -98,9 +100,27 @@ def synthetic_main(
             else None
         )
 
-        if str(v + ".source") in list(yaml_dict.keys()):
-            _ds = xr.open_dataset(yaml_dict[v + ".source.filename"])
-            data = _ds[yaml_dict[v + ".source.variable"]].values
+        def _load_default_static():
+            """Function to read packaged static file"""
+            _ds = pkgr.resource_filename(
+                "mdtf_test_data", f"resources/ocean_static_5deg.nc"
+            )
+            return xr.open_dataset(_ds)["areacello"].values
+
+        # Load the ocean static file
+        if static:
+            if str(v + ".source") in list(yaml_dict.keys()):
+                staticfilepath = yaml_dict[v + ".source.filename"]
+                if os.path.exists(staticfilepath):
+                    _ds = xr.open_dataset(staticfilepath)
+                    data = _ds[yaml_dict[v + ".source.variable"]].values
+                else:
+                    raise ValueError(
+                        f"Specified ocean static file does not exist: {staticfilepath}"
+                    )
+            else:
+                warnings.warn("Using default 5-degree ocean static file for grid")
+                data = _load_default_static()
         else:
             data = None
 

From 6a18bb308e6839afb4f2e087fefac1ec07252be1 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Mon, 11 Oct 2021 12:32:28 -0400
Subject: [PATCH 05/17] Monthly ocean data in CMIP format

- used by tropical sea level POD
---
 mdtf_test_data/config/cmip_mon.yml | 110 +++++++++++++++++++++++++++++
 1 file changed, 110 insertions(+)
 create mode 100644 mdtf_test_data/config/cmip_mon.yml

diff --git a/mdtf_test_data/config/cmip_mon.yml b/mdtf_test_data/config/cmip_mon.yml
new file mode 100644
index 0000000..a7dca8f
--- /dev/null
+++ b/mdtf_test_data/config/cmip_mon.yml
@@ -0,0 +1,110 @@
+variables :
+  name :
+    - "areacello"
+    - "zos"
+    - "tauuo"
+    - "tauvo"
+
+areacello :
+  atts :
+    cell_methods : "area: sum"
+    comment : "TAREA"
+    coordinates : "lat lon"
+    description : "Cell areas for any grid used to report ocean variables and variables which are requested as used on the model ocean grid (e.g. hfsso, which is a downward heat flux from the atmosphere interpolated onto the ocean grid). These cell areas should be defined to enable exact calculation of global integrals (e.g., of vertical fluxes of energy at the surface and top of the atmosphere)."
+    frequency : "fx"
+    id : "areacello"
+    long_name : "Grid-Cell Area for Ocean Variables"
+    mipTable : "Ofx"
+    missing_value : 1.e+20
+    out_name : "areacello"
+    prov : "Ofx ((isd.003))"
+    realm : "ocean"
+    standard_name : "cell_area"
+    time_label : "None"
+    time_title : "No temporal dimensions ... fixed field"
+    title : "Grid-Cell Area for Ocean Variables"
+    type : "real"
+    units : "m2"
+    variable_id : "areacello"
+  grid : "tripolar"
+  static : true
+
+zos :
+  atts :
+    cell_measures : 'area: areacello'
+    cell_methods : 'area: mean where sea time: mean'
+    comment : 'SSH'
+    coordinates : "lat lon"
+    description : 'This is the dynamic sea level, so should have zero global area mean. It should not include inverse barometer depressions from sea ice.'
+    frequency : 'mon'
+    id : 'zos'
+    long_name : 'Sea Surface Height Above Geoid'
+    mipTable : 'Omon'
+    out_name : 'zos'
+    prov : 'Omon ((isd.003))'
+    realm : 'ocean'
+    standard_name : 'sea_surface_height_above_geoid'
+    time : 'time'
+    time_label : 'time-mean'
+    time_title : 'Temporal mean'
+    title : 'Sea Surface Height Above Geoid'
+    type : 'real'
+    units : 'm'
+    variable_id : 'zos'
+  grid : "tripolar"
+  stats :
+    - [-0.09352626651525497, 0.7232627272605896]
+
+tauuo :
+  atts :
+    cell_measures : '--OPT'
+    cell_methods : 'time: mean'
+    comment : 'TAUX'
+    coordinates : "lat lon"
+    description : 'This is the stress on the liquid ocean from overlying atmosphere, sea ice, ice shelf, etc.'
+    frequency : 'mon'
+    id : 'tauuo'
+    long_name : 'Surface Downward X Stress'
+    mipTable : 'Omon'
+    out_name : 'tauuo'
+    positive : 'down'
+    prov : 'Omon ((isd.003))'
+    realm : 'ocean'
+    standard_name : 'surface_downward_x_stress'
+    time : 'time'
+    time_label : 'time-mean'
+    time_title : 'Temporal mean'
+    title : 'Surface Downward X Stress'
+    type : 'real'
+    units : 'N m-2'
+    variable_id : 'tauuo'
+  grid : "tripolar"
+  stats :
+    - [0.001569337327964604, 0.07642165571451187]
+
+tauvo :
+  atts :
+    cell_measures : '--OPT'
+    cell_methods : 'time: mean'
+    comment : 'TAUY'
+    coordinates : "lat lon"
+    description : 'This is the stress on the liquid ocean from overlying atmosphere, sea ice, ice shelf, etc.'
+    frequency : 'mon'
+    id : 'tauuo'
+    long_name : 'Surface Downward Y Stress'
+    mipTable : 'Omon'
+    out_name : 'tauvo'
+    positive : 'down'
+    prov : 'Omon ((isd.003))'
+    realm : 'ocean'
+    standard_name : 'surface_downward_y_stress'
+    time : 'time'
+    time_label : 'time-mean'
+    time_title : 'Temporal mean'
+    title : 'Surface Downward Y Stress'
+    type : 'real'
+    units : 'N m-2'
+    variable_id : 'tauuo'
+  grid : "tripolar"
+  stats :
+    - [0.0020488635636866093, 0.04154461994767189]

From bec696f52b81c201a37e0cd7d2dbb035d71a83bd Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Mon, 11 Oct 2021 12:37:06 -0400
Subject: [PATCH 06/17] Generate tripolar synthetic data in CMIP format

- Generates data on a 5-degree tripolar grid
- Sample entries in the `cmip_mon.yml` file
- Use the `CMIP` convention to generate sample data
  for the tropical sea level MDTF POD
---
 mdtf_test_data/synthetic/synthetic_data.py  | 113 +++++++++++++++-----
 mdtf_test_data/synthetic/synthetic_setup.py |  23 +++-
 scripts/mdtf_synthetic.py                   |  30 +++---
 3 files changed, 123 insertions(+), 43 deletions(-)

diff --git a/mdtf_test_data/synthetic/synthetic_data.py b/mdtf_test_data/synthetic/synthetic_data.py
index 3f853c0..3823c24 100755
--- a/mdtf_test_data/synthetic/synthetic_data.py
+++ b/mdtf_test_data/synthetic/synthetic_data.py
@@ -12,6 +12,7 @@
 import xarray as xr
 import numpy as np
 from mdtf_test_data.synthetic.horizontal import construct_rect_grid
+from mdtf_test_data.synthetic.horizontal import construct_tripolar_grid
 import mdtf_test_data.generators as generators
 
 from mdtf_test_data.synthetic.time import generate_monthly_time_axis
@@ -68,6 +69,9 @@ def generate_synthetic_dataset(
     generator="normal",
     generator_kwargs=None,
     stats=None,
+    static=False,
+    data=None,
+    grid="standard",
 ):
     """Generates xarray dataset of syntheic data in NCAR format
 
@@ -91,6 +95,11 @@ def generate_synthetic_dataset(
         Variable attributes, by default None
     stats : tuple or list of tuples
         Array statistics in the format of [(mean,stddev)]
+    static : bool
+        Flag denoting if variable is static
+    grid : str
+        Type of output grid, either "standard" or "tripolar",
+        by default "standard"
 
     Returns
     -------
@@ -100,32 +109,46 @@ def generate_synthetic_dataset(
 
     attrs = {} if attrs is None else attrs
 
+    # some logical control flags
     do_bounds = True if fmt == "cmip" else False
 
-    dset = construct_rect_grid(
-        dlon, dlat, add_attrs=True, attr_fmt=fmt, bounds=do_bounds
-    )
-    lat = dset.lat
-    lon = dset.lon
-    xyshape = (len(dset["lat"]), len(dset["lon"]))
-
-    if timeres == "mon":
-        ds_time = generate_monthly_time_axis(startyear, nyears, timefmt=fmt)
-    elif timeres == "day":
-        ds_time = generate_daily_time_axis(startyear, nyears, timefmt=fmt)
-    elif timeres == "3hr":
-        ds_time = generate_hourly_time_axis(startyear, nyears, 3, timefmt=fmt)
-    elif timeres == "1hr":
-        ds_time = generate_hourly_time_axis(startyear, nyears, 1, timefmt=fmt)
+    # Step 1: set up the horizontal grid
+    if grid == "tripolar":
+        dset = construct_tripolar_grid(
+            attr_fmt="ncar", retain_coords=True, add_attrs=True
+        )
+        xyshape = dset["mask"].shape
+        lat = dset.nlat
+        lon = dset.nlon
     else:
-        print(timeres)
-        raise ValueError("Unknown time resolution requested")
-
-    dset = ds_time.merge(dset)
-    time = dset["time"]
+        dset = construct_rect_grid(
+            dlon, dlat, add_attrs=True, attr_fmt=fmt, bounds=do_bounds
+        )
+        lat = dset.lat
+        lon = dset.lon
+        xyshape = (len(dset["lat"]), len(dset["lon"]))
+
+    # Step 2: set up the time axis
+    if static is False:
+        if timeres == "mon":
+            ds_time = generate_monthly_time_axis(startyear, nyears, timefmt=fmt)
+        elif timeres == "day":
+            ds_time = generate_daily_time_axis(startyear, nyears, timefmt=fmt)
+        elif timeres == "3hr":
+            ds_time = generate_hourly_time_axis(startyear, nyears, 3, timefmt=fmt)
+        elif timeres == "1hr":
+            ds_time = generate_hourly_time_axis(startyear, nyears, 1, timefmt=fmt)
+        else:
+            print(timeres)
+            raise ValueError("Unknown time resolution requested")
 
-    generator_kwargs = {} if generator_kwargs is None else generator_kwargs
+        dset = ds_time.merge(dset)
+        time = dset["time"]
+        ntimes = len(time)
+    else:
+        ntimes = 1
 
+    # Step 3: generate the vertical coordinate
     if stats is not None:
         stats = [stats] if not isinstance(stats, list) else stats
         if len(stats) > 1:
@@ -139,6 +162,10 @@ def generate_synthetic_dataset(
                 else:
                     dset = dset.merge(gfdl_vertical_coord())
                     lev = dset.pfull
+
+    # Step 4: define the synthetic data generator kernel
+    generator_kwargs = {} if generator_kwargs is None else generator_kwargs
+    if stats is not None:
         generator_kwargs["stats"] = stats
 
     assert generator in list(
@@ -146,25 +173,51 @@ def generate_synthetic_dataset(
     ), f"Unknown generator method: {generator}"
     generator = generators.__dict__[generator]
 
-    data = generators.generate_random_array(
-        xyshape, len(time), generator=generator, generator_kwargs=generator_kwargs
+    # Step 5: generate the synthetic data array
+    data = (
+        generators.generate_random_array(
+            xyshape, ntimes, generator=generator, generator_kwargs=generator_kwargs
+        )
+        if data is None
+        else data
     )
     data = data.squeeze()
 
-    if len(data.shape) == 4:
-        assert data.shape[1] == len(lev), "Length of stats must match number of levels"
-        dset[varname] = xr.DataArray(data, coords=(time, lev, lat, lon), attrs=attrs)
+    # Step 6: convert to Xarray DataArray by assigning coords
+    mask = dset["mask"].values if "mask" in dset.variables else 1.0
+    data = np.array(data * mask, dtype=np.float32)
+
+    if static is True:
+        if len(data.shape) == 4:
+            assert data.shape[1] == len(
+                lev
+            ), "Length of stats must match number of levels"
+            dset[varname] = xr.DataArray(data, coords=(lev, lat, lon), attrs=attrs)
+        else:
+            dset[varname] = xr.DataArray(data, coords=(lat, lon), attrs=attrs)
     else:
-        dset[varname] = xr.DataArray(data, coords=(time, lat, lon), attrs=attrs)
+        if len(data.shape) == 4:
+            print(varname)
+            assert data.shape[1] == len(
+                lev
+            ), "Length of stats must match number of levels"
+            dset[varname] = xr.DataArray(
+                data, coords=(time, lev, lat, lon), attrs=attrs
+            )
+        else:
+            dset[varname] = xr.DataArray(data, coords=(time, lat, lon), attrs=attrs)
+        dset.set_coords(("lat", "lon"))
 
     if coords is not None:
+        print(coords)
         dset[coords["name"]] = xr.DataArray(coords["value"], attrs=coords["atts"])
         dset[varname].attrs = {**dset[varname].attrs, "coordinates": coords["name"]}
 
     dset.attrs["convention"] = fmt
 
     if fmt == "cmip":
-        dset["bnds"].attrs = {"long_name": "vertex number"}
+        if "bnds" in dset.variables:
+            dset["bnds"].attrs = {"long_name": "vertex number"}
         cmip_global_atts = [
             "external_variables",
             "history",
@@ -215,6 +268,10 @@ def generate_synthetic_dataset(
         cmip_global_atts = {x: "" for x in cmip_global_atts}
         dset.attrs = {**dset.attrs, **cmip_global_atts}
 
+    # remove unused fields
+    if grid == "tripolar":
+        dset = dset.drop(["mask", "wet", "depth"])
+
     return dset
 
 
diff --git a/mdtf_test_data/synthetic/synthetic_setup.py b/mdtf_test_data/synthetic/synthetic_setup.py
index 5d27618..190a6e1 100755
--- a/mdtf_test_data/synthetic/synthetic_setup.py
+++ b/mdtf_test_data/synthetic/synthetic_setup.py
@@ -44,6 +44,16 @@ def create_output_dirs(CASENAME="",
                                            )
         out_dir_root= f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}"
 
+    if not os.path.exists(f"{CASENAME}/day"):
+        os.makedirs(f"{CASENAME}/day")
+    if not os.path.exists(f"{CASENAME}/mon"):
+        os.makedirs(f"{CASENAME}/mon")
+    if "NCAR" in CASENAME:
+        if not os.path.exists(f"{CASENAME}/3hr"):
+            os.makedirs(f"{CASENAME}/3hr")
+        if not os.path.exists(f"{CASENAME}/1hr"):
+            os.makedirs(f"{CASENAME}/1hr")
+
     if not os.path.exists(f"{out_dir_root}/day"):
         os.makedirs(f"{out_dir_root}/day")
     if "NCAR" in CASENAME:
@@ -91,8 +101,16 @@ def synthetic_main(
             if str(v + ".generator.args") in list(yaml_dict.keys())
             else {}
         )
-        # vinfo = yaml_dict[v]
-        # print(vinfo)
+        grid = (
+            yaml_dict[v + ".grid"]
+            if str(v + ".grid") in list(yaml_dict.keys())
+            else "standard"
+        )
+
+        assert grid in [
+            "tripolar",
+            "standard",
+        ], f"Unknown grid `{grid}` specified for variable `{v}`"
 
         coords = (
             yaml_dict[v]["coordinates"]
@@ -139,6 +157,7 @@ def _load_default_static():
             coords=coords,
             data=data,
             generator_kwargs=generator_kwargs,
+            grid=grid,
         )
 
         if DATA_FORMAT == "cmip":
diff --git a/scripts/mdtf_synthetic.py b/scripts/mdtf_synthetic.py
index 3a3cb55..ee1fa70 100755
--- a/scripts/mdtf_synthetic.py
+++ b/scripts/mdtf_synthetic.py
@@ -143,20 +143,24 @@ def main():
             )
     if cli_info.convention == "CMIP":
         print("Importing CMIP variable information")
-        input_data = pkgr.resource_filename("mdtf_test_data", "config/cmip_day.yml")
-        input_data = read_yaml(input_data)
+        time_res = ["mon", "day"]
+        for t in time_res:
+            input_data = pkgr.resource_filename(
+                "mdtf_test_data", f"config/cmip_{t}.yml"
+            )
+            input_data = read_yaml(input_data)
 
-        print("Calling Synthetic Data Generator for CMIP data")
-        synthetic_main(
-            input_data,
-            DLAT=cli_info.dlat,
-            DLON=cli_info.dlon,
-            STARTYEAR=cli_info.startyear,
-            NYEARS=cli_info.nyears,
-            CASENAME="CMIP.Synthetic",
-            TIME_RES="day",
-            DATA_FORMAT="cmip",
-        )
+            print("Calling Synthetic Data Generator for CMIP data")
+            synthetic_main(
+                input_data,
+                DLAT=cli_info.dlat,
+                DLON=cli_info.dlon,
+                STARTYEAR=cli_info.startyear,
+                NYEARS=cli_info.nyears,
+                CASENAME="CMIP.Synthetic",
+                TIME_RES=t,
+                DATA_FORMAT="cmip",
+            )
 
 
 if __name__ == "__main__":

From 6274d5b9e0c515cee1c5cb7203a9aa35a273fe23 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Mon, 11 Oct 2021 13:04:28 -0400
Subject: [PATCH 07/17] Issue with CMIP date string bug

- Current implementation only worked with daily data
- More flexible solution works with monthly and daily
---
 mdtf_test_data/synthetic/synthetic_data.py  |  1 -
 mdtf_test_data/synthetic/synthetic_setup.py | 33 +++++++--------------
 2 files changed, 11 insertions(+), 23 deletions(-)

diff --git a/mdtf_test_data/synthetic/synthetic_data.py b/mdtf_test_data/synthetic/synthetic_data.py
index 3823c24..ecc5ef3 100755
--- a/mdtf_test_data/synthetic/synthetic_data.py
+++ b/mdtf_test_data/synthetic/synthetic_data.py
@@ -209,7 +209,6 @@ def generate_synthetic_dataset(
         dset.set_coords(("lat", "lon"))
 
     if coords is not None:
-        print(coords)
         dset[coords["name"]] = xr.DataArray(coords["value"], attrs=coords["atts"])
         dset[varname].attrs = {**dset[varname].attrs, "coordinates": coords["name"]}
 
diff --git a/mdtf_test_data/synthetic/synthetic_setup.py b/mdtf_test_data/synthetic/synthetic_setup.py
index 190a6e1..5a7f8df 100755
--- a/mdtf_test_data/synthetic/synthetic_setup.py
+++ b/mdtf_test_data/synthetic/synthetic_setup.py
@@ -11,11 +11,8 @@
 from .synthetic_data import write_to_netcdf
 
 
-def generate_date_string(STARTYEAR=1,
-                         NYEARS=1,
-                         TIME_RES=""
-                        ):
-    """ formulate the date string in the file name"""
+def generate_date_string(STARTYEAR=1, NYEARS=1, TIME_RES=""):
+    """formulate the date string in the file name"""
     date_string = (
         str(STARTYEAR).zfill(4),
         str(STARTYEAR + NYEARS - 1).zfill(4),
@@ -24,25 +21,16 @@ def generate_date_string(STARTYEAR=1,
         date_string = (date_string[0] + "01", date_string[1] + "12")
     elif TIME_RES == "day":
         date_string = (date_string[0] + "0101", date_string[1] + "1231")
-        date_string = ("-").join(list(date_string))
+    date_string = ("-").join(list(date_string))
 
     return date_string
 
-def create_output_dirs(CASENAME="",
-                       STARTYEAR=1,
-                       NYEARS=10,
-                       TIME_RES="day"
-                      ):
+
+def create_output_dirs(CASENAME="", STARTYEAR=1, NYEARS=10, TIME_RES="day"):
     """Create output data directories"""
 
     out_dir_root = CASENAME
     print("Creating output data directories")
-    if "CMIP" in CASENAME:
-        date_string = generate_date_string(STARTYEAR=STARTYEAR,
-                                           NYEARS=NYEARS,
-                                           TIME_RES=TIME_RES
-                                           )
-        out_dir_root= f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}"
 
     if not os.path.exists(f"{CASENAME}/day"):
         os.makedirs(f"{CASENAME}/day")
@@ -64,6 +52,7 @@ def create_output_dirs(CASENAME="",
         if not os.path.exists(f"{out_dir_root}/1hr"):
             os.makedirs(f"{out_dir_root}/1hr")
 
+
 def synthetic_main(
     yaml_dict={},
     DLAT=20.0,
@@ -162,14 +151,14 @@ def _load_default_static():
 
         if DATA_FORMAT == "cmip":
             # formulate the date string in the file name
-            date_string = generate_date_string(STARTYEAR=STARTYEAR,
-                                               NYEARS=NYEARS,
-                                               TIME_RES=TIME_RES
-                                               )
+            date_string = generate_date_string(
+                STARTYEAR=STARTYEAR, NYEARS=NYEARS, TIME_RES=TIME_RES
+            )
 
             outname = f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}.{v}.{TIME_RES}.nc"
             # output root directory and file name base must match
-            out_dir_root = f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}"
+            # out_dir_root = f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}"
+            out_dir_root = CASENAME
         else:
             outname = f"{CASENAME}.{v}.{TIME_RES}.nc"
             out_dir_root = CASENAME

From ada08cbe11e32764d8a17a7fd72c084fa88cf37e Mon Sep 17 00:00:00 2001
From: Jess <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Tue, 12 Oct 2021 13:19:27 -0400
Subject: [PATCH 08/17] Update README.md to reflect PyPI installation

---
 README.md | 20 ++++++++++----------
 1 file changed, 10 insertions(+), 10 deletions(-)

diff --git a/README.md b/README.md
index 399b890..b5aae73 100644
--- a/README.md
+++ b/README.md
@@ -16,14 +16,13 @@ cftime
 * xESMF
 
 ## Getting the code
-```
-git clone https://github.com/jkrasting/mdtf_test_data.git
-```
+
+The mdtf_test_data package is available on [PyPI](https://pypi.org/project/mdtf-test-data/). 
 
 ## Installation
 ```
-cd mdtf_test_data
-pip install .
+mkdir mdtf_test_data && cd mdtf_test_data
+pip install mdtf-test-data
 ```
 
 ## Usage
@@ -31,7 +30,7 @@ This package can be used to generate fully-synthetic datasets based on NCAR's
 CESM2 and GFDL's CM4 model output that can be used to test the MDTF-Diagnostics package.
 
 ```
-usage: ./scripts/mdtf_synthetic.py [-h] [-c CONVENTION] [--startyear year] [--nyears years]
+usage: mdtf_synthetic.py [-h] [-c CONVENTION] [--startyear year] [--nyears years]
 [--dlat latitude resolution in degrees] [--dlon longitude resolution in degrees] [--unittest]
 
 Required arguments:
@@ -48,23 +47,24 @@ Optional arguments:
 To generate NCAR CESM output in a directory called `NCAR.Synthetic`:
 
 ```
-./scripts/mdtf_synthetic.py -c [ NCAR | CESM ] --nyears 7 
+mdtf_synthetic.py -c [ NCAR | CESM ] --nyears 7 
 ```
 Note that "NCAR" and "CESM" are both valid input values that generate output files appended with "NCAR_Synthetic".
 
 To generate GFDL CM4 output in a directory called `GFDL.Synthetic`:
 ```
-./scripts/mdtf_synthetic.py -c GFDL --nyears 10
+mdtf_synthetic.py -c GFDL --nyears 10
 ```
 
 To generate CMIP output in a directory called `CMIP.Synthetic`:
 ```
-./scripts/mdtf_synthetic.py -c CMIP --nyears 10
+mdtf_synthetic.py -c CMIP --nyears 10
 ```
 
 To coarsen an existing NetCDF file:
 ```
-mdtf_synthetic/util/mdtf-coarsen.py
+git clone https://github.com/jkrasting/mdtf_test_data.git
+cd mdtf_test_data
 usage: mdtf_synthetic/util/mdtf-coarsen.py [-h] [-r REGRID_METHOD] [-o OUTFILE] [-O] infile
 
 Coarsen a NetCDF file.

From 526183e02207f19f5c21fa1efb94f53acbe2b1c0 Mon Sep 17 00:00:00 2001
From: wrongkindofdoctor <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Fri, 15 Oct 2021 17:56:08 -0400
Subject: [PATCH 09/17] added thetao and so fields to cmip_mon.yml

---
 mdtf_test_data/config/cmip_mon.yml | 168 +++++++++++++++++++++++++++++
 1 file changed, 168 insertions(+)

diff --git a/mdtf_test_data/config/cmip_mon.yml b/mdtf_test_data/config/cmip_mon.yml
index a7dca8f..12cdfd7 100644
--- a/mdtf_test_data/config/cmip_mon.yml
+++ b/mdtf_test_data/config/cmip_mon.yml
@@ -4,6 +4,8 @@ variables :
     - "zos"
     - "tauuo"
     - "tauvo"
+    - "so"
+    - "thetao"
 
 areacello :
   atts :
@@ -108,3 +110,169 @@ tauvo :
   grid : "tripolar"
   stats :
     - [0.0020488635636866093, 0.04154461994767189]
+  so:
+    atts:
+      cell_measures: '--OPT'
+      cell_methods: 'time: mean'
+      coordinates: "lat lon"
+      description: 'sea water salinity'
+      frequency: 'mon'
+      id: 'so'
+      long_name: 'sea water salinity'
+      out_name: 'so'
+      positive: 'down'
+      realm: 'ocean'
+      standard_name: 'sea_water_salinity'
+      time: 'time'
+      time_label: 'time-mean'
+      time_title: 'Temporal mean'
+      title: 'Sea Water Salinity'
+      type: 'real'
+      units: 'psu'
+      variable_id: 'so'
+    grid: "tripolar"
+    stats:
+      - [ 37.091705, 37.169098 ]
+      - [ 37.090427, 37.16815 ]
+      - [ 37.086548, 37.163586 ]
+      - [ 37.054672, 37.13838 ]
+      - [ 36.948807, 37.03889 ]
+      - [ 36.872665, 36.96365 ]
+      - [ 36.825863, 36.92407 ]
+      - [ 36.793346, 36.899162 ]
+      - [ 36.767387, 36.878147 ]
+      - [ 36.744244, 36.85687 ]
+      - [ 36.721256, 36.8339 ]
+      - [ 36.698128, 36.80898 ]
+      - [ 36.673008, 36.782074 ]
+      - [ 36.644700, 36.753128 ]
+      - [ 36.61258,  36.72212 ]
+      - [ 36.576668, 36.689003 ]
+      - [ 36.537296, 36.65366 ]
+      - [ 36.495052, 36.615482 ]
+      - [ 36.45111,  36.574173 ]
+      - [ 36.406807, 36.529804 ]
+      - [ 36.36356, 36.483067 ]
+      - [ 36.322468, 36.43506 ]
+      - [ 36.28378, 36.387257 ]
+      - [ 36.247128, 36.340984 ]
+      - [ 36.21161, 36.297478 ]
+      - [ 36.176785, 36.257397 ]
+      - [ 36.14187, 36.22075 ]
+      - [ 36.106518, 36.186737 ]
+      - [ 36.07005, 36.154045 ]
+      - [ 36.03362, 36.121113 ]
+      - [ 36.000076, 36.087433 ]
+      - [ 35.974377, 36.05332 ]
+      - [ 35.96569, 36.020916 ]
+      - [ 35.982994, 35.994526 ]
+      - [ 36.03397, 35.97875 ]
+      - [ 36.11469, 35.44826 ]
+      - [ 36.207542,  35.27616 ]
+      - [ 36.277607, 35.158367 ]
+      - [ 36.288494, 35.102154 ]
+      - [ 36.21689, 35.088215 ]
+      - [ 36.06323, 35.089135 ]
+      - [ 35.852272,  35.091137 ]
+      - [ 35.620174, 35.089703 ]
+      - [ 35.406887, 35.08471 ]
+      - [ 35.24506,  35.078117 ]
+      - [ 35.145275,  35.071728 ]
+      - [ 35.097454,  35.06888 ]
+      - [ 35.080467, 35.06888 ]
+      - [ 35.077976, 35.072856 ]
+      - [ 35.079536, 35.078241 ]
+      - [ 35.081894, 35.0827532 ]
+      - [ 35.084568, 35.083969 ]
+      - [ 35.087505, 35.086932 ]
+      - [ 35.09258, 35.092854 ]
+      - [ 35.09881, 35.099214 ]
+      - [ 35.10581, 35.11024 ]
+      - [ 35.11267, 35.11221 ]
+      - [ 35.11715, 35.116893 ]
+      - [ 35.118637, 35.117593]
+      - [ 35.1188980, 45.118823 ]
+
+    thetao:
+      atts:
+        cell_measures: '--OPT'
+        cell_methods: 'time: mean'
+        coordinates: "lat lon"
+        description: 'sea water potential temperature'
+        frequency: 'mon'
+        id: 'thetao'
+        long_name: 'sea water potential temperature'
+        out_name: 'thetao'
+        positive: 'down'
+        realm: 'ocean'
+        standard_name: 'sea_water_potential_temperature'
+        time: 'time'
+        time_label: 'time-mean'
+        time_title: 'Temporal mean'
+        title: 'Sea Water Potential Temperature'
+        type: 'real'
+        units: 'degC'
+        variable_id: 'thetao'
+      grid: "tripolar"
+      stats:
+        - [ 23.438206, 23.501827 ]
+        - [ 23.393858, 23.457973 ]
+        - [ 23.332743, 23.399124 ]
+        - [ 23.145956, 23.248714 ]
+        - [ 22.192282, 22.395288 ]
+        - [ 21.162468, 21.3618 ]
+        - [ 20.480297, 20.662258 ]
+        - [ 20.052425, 20.224604 ]
+        - [ 19.769545, 19.938387 ]
+        - [ 19.559475, 19.728683 ]
+        - [ 19.38311, 19.554249 ]
+        - [ 19.219982, 19.393152 ]
+        - [ 19.060669, 19.235872 ]
+        - [ 18.89744, 19.073635 ]
+        - [ 18.726213, 18.900887 ]
+        - [ 18.545187, 18.714998 ]
+        - [ 18.352053, 18.513506 ]
+        - [ 18.143354, 18.293566 ]
+        - [ 17.917408, 18.054794 ]
+        - [ 17.674921, 17.799715 ]
+        - [ 17.41897, 17.532667 ]
+        - [ 17.154121, 17.25874 ]
+        - [ 16.885664, 16.982744 ]
+        - [ 16.61834, 16.708748 ]
+        - [ 16.356245, 16.439966 ]
+        - [ 16.10178, 16.178211 ]
+        - [ 15.854486, 15.922733 ]
+        - [ 15.610342, 15.66957 ]
+        - [ 15.362126, 15.412205 ]
+        - [ 15.101099, 15.14296 ]
+        - [ 14.81898, 14.854497 ]
+        - [ 14.509788, 14.541605 ]
+        - [ 14.171214, 14.202383 ]
+        - [ 13.806399, 13.840008 ]
+        - [ 13.422476, 13.461085 ]
+        - [ 13.022508, 13.067213 ]
+        - [ 12.592179, 12.6415415 ]
+        - [ 12.087998, 12.138475 ]
+        - [ 11.438775, 11.487448 ]
+        - [ 10.568429, 10.617012 ]
+        - [ 9.438543, 9.494733 ]
+        - [ 8.101263, 8.172442 ]
+        - [ 6.705785, 6.790224 ]
+        - [ 5.4549055, 5.536766 ]
+        - [ 4.5219636, 4.5789433 ]
+        - [ 3.9483228, 3.969062 ]
+        - [ 3.6311398, 3.623945 ]
+        - [ 3.431602, 3.4139373 ]
+        - [ 3.279515, 3.261523 ]
+        - [ 3.1576657, 3.1415262 ]
+        - [ 3.0560117, 3.044555 ]
+        - [ 2.9734735, 2.9688902 ]
+        - [ 2.9163926, 2.9166791 ]
+        - [ 2.8895073, 2.890756 ]
+        - [ 2.885688, 2.8863788 ]
+        - [ 2.886533, 2.8871552 ]
+        - [ 2.88546, 2.8871353 ]
+        - [ 2.8842352, 2.886324]
+        - [ 2.883122, 2.885201 ]
+        - [ 2.8819532, 2.886345]
+

From 088246b5e859b44482f514694a1ee9e66b9e88a0 Mon Sep 17 00:00:00 2001
From: wrongkindofdoctor <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Mon, 18 Oct 2021 12:00:43 -0400
Subject: [PATCH 10/17] changed lengths of so and thetao stats to match mom6
 z_l size

---
 mdtf_test_data/config/cmip_mon.yml | 329 +++++++++++++++--------------
 1 file changed, 165 insertions(+), 164 deletions(-)

diff --git a/mdtf_test_data/config/cmip_mon.yml b/mdtf_test_data/config/cmip_mon.yml
index 12cdfd7..51ec3b8 100644
--- a/mdtf_test_data/config/cmip_mon.yml
+++ b/mdtf_test_data/config/cmip_mon.yml
@@ -110,169 +110,170 @@ tauvo :
   grid : "tripolar"
   stats :
     - [0.0020488635636866093, 0.04154461994767189]
-  so:
-    atts:
-      cell_measures: '--OPT'
-      cell_methods: 'time: mean'
-      coordinates: "lat lon"
-      description: 'sea water salinity'
-      frequency: 'mon'
-      id: 'so'
-      long_name: 'sea water salinity'
-      out_name: 'so'
-      positive: 'down'
-      realm: 'ocean'
-      standard_name: 'sea_water_salinity'
-      time: 'time'
-      time_label: 'time-mean'
-      time_title: 'Temporal mean'
-      title: 'Sea Water Salinity'
-      type: 'real'
-      units: 'psu'
-      variable_id: 'so'
-    grid: "tripolar"
-    stats:
-      - [ 37.091705, 37.169098 ]
-      - [ 37.090427, 37.16815 ]
-      - [ 37.086548, 37.163586 ]
-      - [ 37.054672, 37.13838 ]
-      - [ 36.948807, 37.03889 ]
-      - [ 36.872665, 36.96365 ]
-      - [ 36.825863, 36.92407 ]
-      - [ 36.793346, 36.899162 ]
-      - [ 36.767387, 36.878147 ]
-      - [ 36.744244, 36.85687 ]
-      - [ 36.721256, 36.8339 ]
-      - [ 36.698128, 36.80898 ]
-      - [ 36.673008, 36.782074 ]
-      - [ 36.644700, 36.753128 ]
-      - [ 36.61258,  36.72212 ]
-      - [ 36.576668, 36.689003 ]
-      - [ 36.537296, 36.65366 ]
-      - [ 36.495052, 36.615482 ]
-      - [ 36.45111,  36.574173 ]
-      - [ 36.406807, 36.529804 ]
-      - [ 36.36356, 36.483067 ]
-      - [ 36.322468, 36.43506 ]
-      - [ 36.28378, 36.387257 ]
-      - [ 36.247128, 36.340984 ]
-      - [ 36.21161, 36.297478 ]
-      - [ 36.176785, 36.257397 ]
-      - [ 36.14187, 36.22075 ]
-      - [ 36.106518, 36.186737 ]
-      - [ 36.07005, 36.154045 ]
-      - [ 36.03362, 36.121113 ]
-      - [ 36.000076, 36.087433 ]
-      - [ 35.974377, 36.05332 ]
-      - [ 35.96569, 36.020916 ]
-      - [ 35.982994, 35.994526 ]
-      - [ 36.03397, 35.97875 ]
-      - [ 36.11469, 35.44826 ]
-      - [ 36.207542,  35.27616 ]
-      - [ 36.277607, 35.158367 ]
-      - [ 36.288494, 35.102154 ]
-      - [ 36.21689, 35.088215 ]
-      - [ 36.06323, 35.089135 ]
-      - [ 35.852272,  35.091137 ]
-      - [ 35.620174, 35.089703 ]
-      - [ 35.406887, 35.08471 ]
-      - [ 35.24506,  35.078117 ]
-      - [ 35.145275,  35.071728 ]
-      - [ 35.097454,  35.06888 ]
-      - [ 35.080467, 35.06888 ]
-      - [ 35.077976, 35.072856 ]
-      - [ 35.079536, 35.078241 ]
-      - [ 35.081894, 35.0827532 ]
-      - [ 35.084568, 35.083969 ]
-      - [ 35.087505, 35.086932 ]
-      - [ 35.09258, 35.092854 ]
-      - [ 35.09881, 35.099214 ]
-      - [ 35.10581, 35.11024 ]
-      - [ 35.11267, 35.11221 ]
-      - [ 35.11715, 35.116893 ]
-      - [ 35.118637, 35.117593]
-      - [ 35.1188980, 45.118823 ]
 
-    thetao:
-      atts:
-        cell_measures: '--OPT'
-        cell_methods: 'time: mean'
-        coordinates: "lat lon"
-        description: 'sea water potential temperature'
-        frequency: 'mon'
-        id: 'thetao'
-        long_name: 'sea water potential temperature'
-        out_name: 'thetao'
-        positive: 'down'
-        realm: 'ocean'
-        standard_name: 'sea_water_potential_temperature'
-        time: 'time'
-        time_label: 'time-mean'
-        time_title: 'Temporal mean'
-        title: 'Sea Water Potential Temperature'
-        type: 'real'
-        units: 'degC'
-        variable_id: 'thetao'
-      grid: "tripolar"
-      stats:
-        - [ 23.438206, 23.501827 ]
-        - [ 23.393858, 23.457973 ]
-        - [ 23.332743, 23.399124 ]
-        - [ 23.145956, 23.248714 ]
-        - [ 22.192282, 22.395288 ]
-        - [ 21.162468, 21.3618 ]
-        - [ 20.480297, 20.662258 ]
-        - [ 20.052425, 20.224604 ]
-        - [ 19.769545, 19.938387 ]
-        - [ 19.559475, 19.728683 ]
-        - [ 19.38311, 19.554249 ]
-        - [ 19.219982, 19.393152 ]
-        - [ 19.060669, 19.235872 ]
-        - [ 18.89744, 19.073635 ]
-        - [ 18.726213, 18.900887 ]
-        - [ 18.545187, 18.714998 ]
-        - [ 18.352053, 18.513506 ]
-        - [ 18.143354, 18.293566 ]
-        - [ 17.917408, 18.054794 ]
-        - [ 17.674921, 17.799715 ]
-        - [ 17.41897, 17.532667 ]
-        - [ 17.154121, 17.25874 ]
-        - [ 16.885664, 16.982744 ]
-        - [ 16.61834, 16.708748 ]
-        - [ 16.356245, 16.439966 ]
-        - [ 16.10178, 16.178211 ]
-        - [ 15.854486, 15.922733 ]
-        - [ 15.610342, 15.66957 ]
-        - [ 15.362126, 15.412205 ]
-        - [ 15.101099, 15.14296 ]
-        - [ 14.81898, 14.854497 ]
-        - [ 14.509788, 14.541605 ]
-        - [ 14.171214, 14.202383 ]
-        - [ 13.806399, 13.840008 ]
-        - [ 13.422476, 13.461085 ]
-        - [ 13.022508, 13.067213 ]
-        - [ 12.592179, 12.6415415 ]
-        - [ 12.087998, 12.138475 ]
-        - [ 11.438775, 11.487448 ]
-        - [ 10.568429, 10.617012 ]
-        - [ 9.438543, 9.494733 ]
-        - [ 8.101263, 8.172442 ]
-        - [ 6.705785, 6.790224 ]
-        - [ 5.4549055, 5.536766 ]
-        - [ 4.5219636, 4.5789433 ]
-        - [ 3.9483228, 3.969062 ]
-        - [ 3.6311398, 3.623945 ]
-        - [ 3.431602, 3.4139373 ]
-        - [ 3.279515, 3.261523 ]
-        - [ 3.1576657, 3.1415262 ]
-        - [ 3.0560117, 3.044555 ]
-        - [ 2.9734735, 2.9688902 ]
-        - [ 2.9163926, 2.9166791 ]
-        - [ 2.8895073, 2.890756 ]
-        - [ 2.885688, 2.8863788 ]
-        - [ 2.886533, 2.8871552 ]
-        - [ 2.88546, 2.8871353 ]
-        - [ 2.8842352, 2.886324]
-        - [ 2.883122, 2.885201 ]
-        - [ 2.8819532, 2.886345]
+so :
+  atts :
+    cell_measures : '--OPT'
+    cell_methods : 'time: mean'
+    coordinates : "lat lon"
+    description : 'sea water salinity'
+    frequency : 'mon'
+    id : 'so'
+    long_name : 'sea water salinity'
+    out_name : 'so'
+    positive : 'down'
+    realm : 'ocean'
+    standard_name : 'sea_water_salinity'
+    time : 'time'
+    time_label : 'time-mean'
+    time_title : 'Temporal mean'
+    title : 'Sea Water Salinity'
+    type : 'real'
+    units : 'psu'
+    variable_id : 'so'
+  grid : "tripolar"
+  stats :
+    - [ 37.091705, 37.169098 ]
+    - [ 37.090427, 37.16815 ]
+    - [ 37.086548, 37.163586 ]
+    - [ 37.054672, 37.13838 ]
+    - [ 36.948807, 37.03889 ]
+    - [ 36.872665, 36.96365 ]
+    - [ 36.825863, 36.92407 ]
+    - [ 36.793346, 36.899162 ]
+    - [ 36.767387, 36.878147 ]
+    - [ 36.744244, 36.85687 ]
+    - [ 36.721256, 36.8339 ]
+    - [ 36.698128, 36.80898 ]
+    - [ 36.673008, 36.782074 ]
+    - [ 36.644700, 36.753128 ]
+    - [ 36.61258,  36.72212 ]
+    - [ 36.576668, 36.689003 ]
+    - [ 36.537296, 36.65366 ]
+    - [ 36.495052, 36.615482 ]
+    - [ 36.45111, 36.574173 ]
+    - [ 36.406807, 36.529804 ]
+    - [ 36.36356, 36.483067 ]
+    - [ 36.322468, 36.43506 ]
+    - [ 36.28378, 36.387257 ]
+    - [ 36.247128, 36.340984 ]
+    - [ 36.21161, 36.297478 ]
+    - [ 36.176785, 36.257397 ]
+    - [ 36.14187, 36.22075 ]
+    - [ 36.106518, 36.186737 ]
+    - [ 36.07005, 36.154045 ]
+    - [ 36.03362, 36.121113 ]
+    - [ 36.000076, 36.087433 ]
+    - [ 35.974377, 36.05332 ]
+    - [ 35.96569, 36.020916 ]
+    - [ 35.982994, 35.994526 ]
+    - [ 36.03397, 35.97875 ]
+#    - [ 36.11469, 35.44826 ]
+#    - [ 36.207542, 35.27616 ]
+#    - [ 36.277607, 35.158367 ]
+#    - [ 36.288494, 35.102154 ]
+#    - [ 36.21689, 35.088215 ]
+#    - [ 36.06323, 35.089135 ]
+#    - [ 35.852272, 35.091137 ]
+#    - [ 35.620174, 35.089703 ]
+#    - [ 35.406887, 35.08471 ]
+#    - [ 35.24506, 35.078117 ]
+#    - [ 35.145275, 35.071728 ]
+#    - [ 35.097454, 35.06888 ]
+#    - [ 35.080467, 35.06888 ]
+#    - [ 35.077976, 35.072856 ]
+#    - [ 35.079536, 35.078241 ]
+#    - [ 35.081894, 35.0827532 ]
+#    - [ 35.084568, 35.083969 ]
+#    - [ 35.087505, 35.086932 ]
+#    - [ 35.09258, 35.092854 ]
+#    - [ 35.09881, 35.099214 ]
+#    - [ 35.10581, 35.11024 ]
+#    - [ 35.11267, 35.11221 ]
+#    - [ 35.11715, 35.116893 ]
+#    - [ 35.118637, 35.117593]
+#    - [ 35.1188980, 45.118823 ]
+
+thetao:
+  atts:
+    cell_measures: '--OPT'
+    cell_methods: 'time: mean'
+    coordinates: "lat lon"
+    description: 'sea water potential temperature'
+    frequency: 'mon'
+    id: 'thetao'
+    long_name: 'sea water potential temperature'
+    out_name: 'thetao'
+    positive: 'down'
+    realm: 'ocean'
+    standard_name: 'sea_water_potential_temperature'
+    time: 'time'
+    time_label: 'time-mean'
+    time_title: 'Temporal mean'
+    title: 'Sea Water Potential Temperature'
+    type: 'real'
+    units: 'degC'
+    variable_id: 'thetao'
+  grid: "tripolar"
+  stats:
+    - [ 23.438206, 23.501827 ]
+    - [ 23.393858, 23.457973 ]
+    - [ 23.332743, 23.399124 ]
+    - [ 23.145956, 23.248714 ]
+    - [ 22.192282, 22.395288 ]
+    - [ 21.162468, 21.3618 ]
+    - [ 20.480297, 20.662258 ]
+    - [ 20.052425, 20.224604 ]
+    - [ 19.769545, 19.938387 ]
+    - [ 19.559475, 19.728683 ]
+    - [ 19.38311, 19.554249 ]
+    - [ 19.219982, 19.393152 ]
+    - [ 19.060669, 19.235872 ]
+    - [ 18.89744, 19.073635 ]
+    - [ 18.726213, 18.900887 ]
+    - [ 18.545187, 18.714998 ]
+    - [ 18.352053, 18.513506 ]
+    - [ 18.143354, 18.293566 ]
+    - [ 17.917408, 18.054794 ]
+    - [ 17.674921, 17.799715 ]
+    - [ 17.41897, 17.532667 ]
+    - [ 17.154121, 17.25874 ]
+    - [ 16.885664, 16.982744 ]
+    - [ 16.61834, 16.708748 ]
+    - [ 16.356245, 16.439966 ]
+    - [ 16.10178, 16.178211 ]
+    - [ 15.854486, 15.922733 ]
+    - [ 15.610342, 15.66957 ]
+    - [ 15.362126, 15.412205 ]
+    - [ 15.101099, 15.14296 ]
+    - [ 14.81898, 14.854497 ]
+    - [ 14.509788, 14.541605 ]
+    - [ 14.171214, 14.202383 ]
+    - [ 13.806399, 13.840008 ]
+    - [ 13.422476, 13.461085 ]
+#    - [ 13.022508, 13.067213 ]
+#    - [ 12.592179, 12.6415415 ]
+#    - [ 12.087998, 12.138475 ]
+#    - [ 11.438775, 11.487448 ]
+#    - [ 10.568429, 10.617012 ]
+#    - [ 9.438543, 9.494733 ]
+#    - [ 8.101263, 8.172442 ]
+#    - [ 6.705785, 6.790224 ]
+#    - [ 5.4549055, 5.536766 ]
+#    - [ 4.5219636, 4.5789433 ]
+#    - [ 3.9483228, 3.969062 ]
+#    - [ 3.6311398, 3.623945 ]
+#    - [ 3.431602, 3.4139373 ]
+#    - [ 3.279515, 3.261523 ]
+#    - [ 3.1576657, 3.1415262 ]
+#    - [ 3.0560117, 3.044555 ]
+#    - [ 2.9734735, 2.9688902 ]
+#    - [ 2.9163926, 2.9166791 ]
+#    - [ 2.8895073, 2.890756 ]
+#    - [ 2.885688, 2.8863788 ]
+#    - [ 2.886533, 2.8871552 ]
+#    - [ 2.88546, 2.8871353 ]
+#    - [ 2.8842352, 2.886324]
+#    - [ 2.883122, 2.885201 ]
+#    - [ 2.8819532, 2.886345 ]
 

From 533c969b321d295495d0c4963d174ee52fa2c284 Mon Sep 17 00:00:00 2001
From: wrongkindofdoctor <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Mon, 18 Oct 2021 12:02:28 -0400
Subject: [PATCH 11/17] added vertical level generation for cmip tripolar
 variables that uses mom6 vertical coordinates

---
 mdtf_test_data/synthetic/synthetic_data.py | 15 +++++++++++----
 scripts/mdtf_synthetic.py                  |  3 ++-
 2 files changed, 13 insertions(+), 5 deletions(-)

diff --git a/mdtf_test_data/synthetic/synthetic_data.py b/mdtf_test_data/synthetic/synthetic_data.py
index ecc5ef3..a26e2a2 100755
--- a/mdtf_test_data/synthetic/synthetic_data.py
+++ b/mdtf_test_data/synthetic/synthetic_data.py
@@ -22,6 +22,7 @@
 from mdtf_test_data.synthetic.vertical import gfdl_plev19_vertical_coord
 from mdtf_test_data.synthetic.vertical import gfdl_vertical_coord
 from mdtf_test_data.synthetic.vertical import ncar_hybrid_coord
+from mdtf_test_data.synthetic.vertical import mom6_z_coord
 
 
 def dataset_stats(filename, var=None, limit=None):
@@ -155,14 +156,20 @@ def generate_synthetic_dataset(
             if fmt == "ncar":
                 dset = dset.merge(ncar_hybrid_coord())
                 lev = dset.lev
-            elif fmt == "gfdl":
+            elif fmt == "gfdl" :
                 if len(stats) == 19:
                     dset = dset.merge(gfdl_plev19_vertical_coord())
                     lev = dset.plev19
                 else:
                     dset = dset.merge(gfdl_vertical_coord())
                     lev = dset.pfull
-
+            elif fmt == "cmip" and grid == "tripolar":
+                if len(stats) == 35:
+                    dset = dset.merge(mom6_z_coord())
+                    lev = dset.z_l
+                elif len(stats) == 36:
+                    dset = dset.merge(mom6_z_coord())
+                    lev = dset.z_i
     # Step 4: define the synthetic data generator kernel
     generator_kwargs = {} if generator_kwargs is None else generator_kwargs
     if stats is not None:
@@ -191,7 +198,7 @@ def generate_synthetic_dataset(
         if len(data.shape) == 4:
             assert data.shape[1] == len(
                 lev
-            ), "Length of stats must match number of levels"
+            ),  f' Length of stats {data.shape[1]} must match number of levels {len(lev)}.'
             dset[varname] = xr.DataArray(data, coords=(lev, lat, lon), attrs=attrs)
         else:
             dset[varname] = xr.DataArray(data, coords=(lat, lon), attrs=attrs)
@@ -200,7 +207,7 @@ def generate_synthetic_dataset(
             print(varname)
             assert data.shape[1] == len(
                 lev
-            ), "Length of stats must match number of levels"
+            ), f' Length of stats {data.shape[1]} must match number of levels {len(lev)}.'
             dset[varname] = xr.DataArray(
                 data, coords=(time, lev, lat, lon), attrs=attrs
             )
diff --git a/scripts/mdtf_synthetic.py b/scripts/mdtf_synthetic.py
index ee1fa70..a317fd2 100755
--- a/scripts/mdtf_synthetic.py
+++ b/scripts/mdtf_synthetic.py
@@ -1,6 +1,8 @@
 #!/usr/bin/env python
 """ mdtf_test_data driver program """
 import sys
+sys.path.insert(0, '../mdtf_test_data')
+import sys
 import mdtf_test_data
 from mdtf_test_data.synthetic.synthetic_setup import synthetic_main
 from mdtf_test_data.util.cli import cli_holder
@@ -10,7 +12,6 @@
 
 MDTF_PACKAGE_PATH = mdtf_test_data.__path__[0]
 
-
 def read_yaml(file_name):
     """A function to read YAML files"""
     config = EnvYAML(file_name)

From 0f1bedf644100b462bef74bb62be71c0fbb95236 Mon Sep 17 00:00:00 2001
From: wrongkindofdoctor <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Mon, 18 Oct 2021 12:56:01 -0400
Subject: [PATCH 12/17] changed the root directory name for the cmip output to
 match the output data file base name and the file file base names to have the
 format yyyymmdd for framework functionality

---
 mdtf_test_data/synthetic/synthetic_setup.py | 31 ++++++++++-----------
 1 file changed, 14 insertions(+), 17 deletions(-)

diff --git a/mdtf_test_data/synthetic/synthetic_setup.py b/mdtf_test_data/synthetic/synthetic_setup.py
index 5a7f8df..ac4c220 100755
--- a/mdtf_test_data/synthetic/synthetic_setup.py
+++ b/mdtf_test_data/synthetic/synthetic_setup.py
@@ -28,25 +28,23 @@ def generate_date_string(STARTYEAR=1, NYEARS=1, TIME_RES=""):
 
 def create_output_dirs(CASENAME="", STARTYEAR=1, NYEARS=10, TIME_RES="day"):
     """Create output data directories"""
+    if "cmip" in str.lower(CASENAME):
+        # formulate the date string in the file name
+        date_string = generate_date_string(
+            STARTYEAR=STARTYEAR, NYEARS=NYEARS, TIME_RES="day"
+        )
+        # output root directory and file name base must match
+        out_dir_root = f"{CASENAME.replace('.', '_')}_r1i1p1f1_gr1_{date_string}"
+    else:
+        out_dir_root = CASENAME
 
-    out_dir_root = CASENAME
     print("Creating output data directories")
 
-    if not os.path.exists(f"{CASENAME}/day"):
-        os.makedirs(f"{CASENAME}/day")
-    if not os.path.exists(f"{CASENAME}/mon"):
-        os.makedirs(f"{CASENAME}/mon")
-    if "NCAR" in CASENAME:
-        if not os.path.exists(f"{CASENAME}/3hr"):
-            os.makedirs(f"{CASENAME}/3hr")
-        if not os.path.exists(f"{CASENAME}/1hr"):
-            os.makedirs(f"{CASENAME}/1hr")
-
     if not os.path.exists(f"{out_dir_root}/day"):
         os.makedirs(f"{out_dir_root}/day")
-    if "NCAR" in CASENAME:
-        if not os.path.exists(f"{out_dir_root}/mon"):
-            os.makedirs(f"{out_dir_root}/mon")
+    if not os.path.exists(f"{out_dir_root}/mon"):
+        os.makedirs(f"{out_dir_root}/mon")
+    if "ncar" in str.lower(out_dir_root):
         if not os.path.exists(f"{out_dir_root}/3hr"):
             os.makedirs(f"{out_dir_root}/3hr")
         if not os.path.exists(f"{out_dir_root}/1hr"):
@@ -152,13 +150,12 @@ def _load_default_static():
         if DATA_FORMAT == "cmip":
             # formulate the date string in the file name
             date_string = generate_date_string(
-                STARTYEAR=STARTYEAR, NYEARS=NYEARS, TIME_RES=TIME_RES
+                STARTYEAR=STARTYEAR, NYEARS=NYEARS, TIME_RES="day"
             )
 
             outname = f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}.{v}.{TIME_RES}.nc"
             # output root directory and file name base must match
-            # out_dir_root = f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}"
-            out_dir_root = CASENAME
+            out_dir_root = f"{CASENAME.replace('.','_')}_r1i1p1f1_gr1_{date_string}"
         else:
             outname = f"{CASENAME}.{v}.{TIME_RES}.nc"
             out_dir_root = CASENAME

From af54d026c42b27be4dc62e4ed41c404436cb18c2 Mon Sep 17 00:00:00 2001
From: wrongkindofdoctor <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Mon, 18 Oct 2021 15:48:49 -0400
Subject: [PATCH 13/17] added STARTYEAR and NYEARS arguments to
 create_output_dirs call in synthetic_setup

---
 mdtf_test_data/synthetic/synthetic_setup.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/mdtf_test_data/synthetic/synthetic_setup.py b/mdtf_test_data/synthetic/synthetic_setup.py
index ac4c220..8512b50 100755
--- a/mdtf_test_data/synthetic/synthetic_setup.py
+++ b/mdtf_test_data/synthetic/synthetic_setup.py
@@ -62,7 +62,7 @@ def synthetic_main(
     DATA_FORMAT="",
 ):
     """Main script to generate synthetic data using GFDL naming conventions"""
-    create_output_dirs(CASENAME)
+    create_output_dirs(CASENAME, STARTYEAR=STARTYEAR, NYEARS=NYEARS)
     # parse the yaml dictionary
     var_names = yaml_dict["variables.name"]
     # -- Create Data

From 23cca37cf5124a02e16b984505eba2b1861eb966 Mon Sep 17 00:00:00 2001
From: wrongkindofdoctor <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Mon, 18 Oct 2021 16:34:32 -0400
Subject: [PATCH 14/17] cleaned up cmip_mon.yml added cmip standard lev
 variable to mom6_z_coord.py changed data check procedure and lev definition
 for cmip tripolar level generation

---
 mdtf_test_data/config/cmip_mon.yml            | 51 +------------------
 mdtf_test_data/synthetic/synthetic_data.py    | 11 ++--
 .../synthetic/vertical/mom6_z_coord.py        | 14 +++++
 3 files changed, 20 insertions(+), 56 deletions(-)
 mode change 100644 => 100755 mdtf_test_data/config/cmip_mon.yml

diff --git a/mdtf_test_data/config/cmip_mon.yml b/mdtf_test_data/config/cmip_mon.yml
old mode 100644
new mode 100755
index 51ec3b8..aa3f278
--- a/mdtf_test_data/config/cmip_mon.yml
+++ b/mdtf_test_data/config/cmip_mon.yml
@@ -168,31 +168,6 @@ so :
     - [ 35.96569, 36.020916 ]
     - [ 35.982994, 35.994526 ]
     - [ 36.03397, 35.97875 ]
-#    - [ 36.11469, 35.44826 ]
-#    - [ 36.207542, 35.27616 ]
-#    - [ 36.277607, 35.158367 ]
-#    - [ 36.288494, 35.102154 ]
-#    - [ 36.21689, 35.088215 ]
-#    - [ 36.06323, 35.089135 ]
-#    - [ 35.852272, 35.091137 ]
-#    - [ 35.620174, 35.089703 ]
-#    - [ 35.406887, 35.08471 ]
-#    - [ 35.24506, 35.078117 ]
-#    - [ 35.145275, 35.071728 ]
-#    - [ 35.097454, 35.06888 ]
-#    - [ 35.080467, 35.06888 ]
-#    - [ 35.077976, 35.072856 ]
-#    - [ 35.079536, 35.078241 ]
-#    - [ 35.081894, 35.0827532 ]
-#    - [ 35.084568, 35.083969 ]
-#    - [ 35.087505, 35.086932 ]
-#    - [ 35.09258, 35.092854 ]
-#    - [ 35.09881, 35.099214 ]
-#    - [ 35.10581, 35.11024 ]
-#    - [ 35.11267, 35.11221 ]
-#    - [ 35.11715, 35.116893 ]
-#    - [ 35.118637, 35.117593]
-#    - [ 35.1188980, 45.118823 ]
 
 thetao:
   atts:
@@ -251,29 +226,5 @@ thetao:
     - [ 14.171214, 14.202383 ]
     - [ 13.806399, 13.840008 ]
     - [ 13.422476, 13.461085 ]
-#    - [ 13.022508, 13.067213 ]
-#    - [ 12.592179, 12.6415415 ]
-#    - [ 12.087998, 12.138475 ]
-#    - [ 11.438775, 11.487448 ]
-#    - [ 10.568429, 10.617012 ]
-#    - [ 9.438543, 9.494733 ]
-#    - [ 8.101263, 8.172442 ]
-#    - [ 6.705785, 6.790224 ]
-#    - [ 5.4549055, 5.536766 ]
-#    - [ 4.5219636, 4.5789433 ]
-#    - [ 3.9483228, 3.969062 ]
-#    - [ 3.6311398, 3.623945 ]
-#    - [ 3.431602, 3.4139373 ]
-#    - [ 3.279515, 3.261523 ]
-#    - [ 3.1576657, 3.1415262 ]
-#    - [ 3.0560117, 3.044555 ]
-#    - [ 2.9734735, 2.9688902 ]
-#    - [ 2.9163926, 2.9166791 ]
-#    - [ 2.8895073, 2.890756 ]
-#    - [ 2.885688, 2.8863788 ]
-#    - [ 2.886533, 2.8871552 ]
-#    - [ 2.88546, 2.8871353 ]
-#    - [ 2.8842352, 2.886324]
-#    - [ 2.883122, 2.885201 ]
-#    - [ 2.8819532, 2.886345 ]
+
 
diff --git a/mdtf_test_data/synthetic/synthetic_data.py b/mdtf_test_data/synthetic/synthetic_data.py
index a26e2a2..6f99bad 100755
--- a/mdtf_test_data/synthetic/synthetic_data.py
+++ b/mdtf_test_data/synthetic/synthetic_data.py
@@ -164,12 +164,11 @@ def generate_synthetic_dataset(
                     dset = dset.merge(gfdl_vertical_coord())
                     lev = dset.pfull
             elif fmt == "cmip" and grid == "tripolar":
-                if len(stats) == 35:
-                    dset = dset.merge(mom6_z_coord())
-                    lev = dset.z_l
-                elif len(stats) == 36:
-                    dset = dset.merge(mom6_z_coord())
-                    lev = dset.z_i
+                dset = dset.merge(mom6_z_coord())
+                lev = dset.lev
+                assert len(stats) == len(lev),\
+                    f' Length of stats {data.shape[1]} must match number of levels {len(lev)}.'
+
     # Step 4: define the synthetic data generator kernel
     generator_kwargs = {} if generator_kwargs is None else generator_kwargs
     if stats is not None:
diff --git a/mdtf_test_data/synthetic/vertical/mom6_z_coord.py b/mdtf_test_data/synthetic/vertical/mom6_z_coord.py
index 34d7cd4..846b762 100755
--- a/mdtf_test_data/synthetic/vertical/mom6_z_coord.py
+++ b/mdtf_test_data/synthetic/vertical/mom6_z_coord.py
@@ -101,6 +101,7 @@ def mom6_z_coord():
     ]
 
     z_i_attrs = {
+        "standard_name": "depth_at_cell_interface",
         "long_name": "Depth at interface",
         "units": "meters",
         "axis": "Z",
@@ -108,6 +109,7 @@ def mom6_z_coord():
     }
 
     z_l_attrs = {
+        "standard_name": "depth_at_cell_center",
         "long_name": "Depth at cell center",
         "units": "meters",
         "axis": "Z",
@@ -115,6 +117,15 @@ def mom6_z_coord():
         "edges": "z_i",
     }
 
+    # duplicate of z_l with CMIP standard attributes
+    lev_attrs = {
+        "standard_name": "depth",
+        "long_name": "depth",
+        "units": "meters",
+        "axis": "Z",
+        "positive": "down"
+    }
+
     dset_out = xr.Dataset()
     dset_out["z_l"] = xr.DataArray(
         z_l, dims={"z_l": z_l}, coords={"z_l": z_l}, attrs=z_l_attrs
@@ -122,5 +133,8 @@ def mom6_z_coord():
     dset_out["z_i"] = xr.DataArray(
         z_i, dims={"z_i": z_i}, coords={"z_i": z_i}, attrs=z_i_attrs
     )
+    dset_out["lev"] = xr.DataArray(
+        z_l, dims={"lev": z_l}, coords={"lev": z_l}, attrs=lev_attrs
+    )
 
     return dset_out

From d6a7ff86e14381bfabdb9851874d06a7b5a7ebad Mon Sep 17 00:00:00 2001
From: Jess <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Wed, 20 Oct 2021 09:36:32 -0400
Subject: [PATCH 15/17] revert mdtf_synthetic.py mods

---
 scripts/mdtf_synthetic.py | 2 --
 1 file changed, 2 deletions(-)

diff --git a/scripts/mdtf_synthetic.py b/scripts/mdtf_synthetic.py
index a317fd2..be2ceb6 100755
--- a/scripts/mdtf_synthetic.py
+++ b/scripts/mdtf_synthetic.py
@@ -1,8 +1,6 @@
 #!/usr/bin/env python
 """ mdtf_test_data driver program """
 import sys
-sys.path.insert(0, '../mdtf_test_data')
-import sys
 import mdtf_test_data
 from mdtf_test_data.synthetic.synthetic_setup import synthetic_main
 from mdtf_test_data.util.cli import cli_holder

From 5654440ba2b277c94e8ad0eabcc84797fdfea454 Mon Sep 17 00:00:00 2001
From: Jess <20195932+wrongkindofdoctor@users.noreply.github.com>
Date: Wed, 20 Oct 2021 09:37:48 -0400
Subject: [PATCH 16/17] return space to mdtf_synthetic.py

---
 scripts/mdtf_synthetic.py | 1 +
 1 file changed, 1 insertion(+)

diff --git a/scripts/mdtf_synthetic.py b/scripts/mdtf_synthetic.py
index be2ceb6..ee1fa70 100755
--- a/scripts/mdtf_synthetic.py
+++ b/scripts/mdtf_synthetic.py
@@ -10,6 +10,7 @@
 
 MDTF_PACKAGE_PATH = mdtf_test_data.__path__[0]
 
+
 def read_yaml(file_name):
     """A function to read YAML files"""
     config = EnvYAML(file_name)

From 63676f367f1e8e51657f7d8809041f5d35b085f4 Mon Sep 17 00:00:00 2001
From: John Krasting <John.Krasting@noaa.gov>
Date: Thu, 21 Oct 2021 10:26:15 -0400
Subject: [PATCH 17/17] Bumped version number to v1.0.4rc1

---
 setup.cfg | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/setup.cfg b/setup.cfg
index 195c291..d8e5a51 100644
--- a/setup.cfg
+++ b/setup.cfg
@@ -1,6 +1,6 @@
 [metadata]
 name = mdtf_test_data
-version = 1.0.3.post2
+version = 1.0.4.rc1
 description = Tools for working with MDTF Diagnostics test data sets
 url = https://github.com/jkrasting/mdtf_test_data
 author = MDTF Development Team