From d8a0cdc1c67a45345c6c6751e0349189f0fdcee9 Mon Sep 17 00:00:00 2001
From: Sander Roosendaal <sanderroosendaal@proton.me>
Date: Wed, 6 Mar 2024 13:10:59 +0100
Subject: [PATCH] starting map compare

---
 rowers/interactiveplots.py            | 330 +++-----------------------
 rowers/tests/testdata/testdata.tcx.gz | Bin 4000 -> 4001 bytes
 2 files changed, 28 insertions(+), 302 deletions(-)

diff --git a/rowers/interactiveplots.py b/rowers/interactiveplots.py
index 67c1996b..73b041c5 100644
--- a/rowers/interactiveplots.py
+++ b/rowers/interactiveplots.py
@@ -92,7 +92,13 @@ def get_chart(end_point, chart_data):
 
     url = settings.ROWSANDALL_CHARTS_URL+end_point
     headers = {'authorization':"Bearer {token}".format(token=settings.ROWSANDALL_CHARTS_TOKEN)}
-    response = requests.post(url, json=chart_data, headers=headers)
+    try:
+        response = requests.post(url, json=chart_data, headers=headers)
+    except:
+        script = ''
+        div = 'Chart Server Error'
+        return script, div
+        
     if response.status_code == 200:
         script = response.json()['script']
         div = response.json()['div']
@@ -2202,273 +2208,6 @@ def course_map(course):
 
 
 
-def get_map_script_course(
-    latmean,
-    lonmean,
-    latbegin,
-    latend,
-    longbegin,
-    longend,
-    scoordinates,
-    course,
-):  # pragma: no cover
-    latmean, lonmean, coordinates = course_coord_center(course)
-    lat_min, lat_max, long_min, long_max = course_coord_maxmin(course)
-
-    coordinates = course_spline(coordinates)
-
-    scoordinates = "["
-
-    for index, row in coordinates.iterrows():
-        scoordinates += """[{x},{y}],
-        """.format(
-            x=row['latitude'],
-            y=row['longitude']
-        )
-
-    scoordinates += "]"
-
-    polygons = GeoPolygon.objects.filter(
-        course=course).order_by("order_in_course")
-
-    plabels = ''
-
-    for p in polygons:
-        coords = polygon_coord_center(p)
-
-        plabels += """
-        var marker = L.marker([{latbegin}, {longbegin}]).addTo(mymap);
-        marker.bindPopup("<b>{name}</b>");
-
-        """.format(
-            latbegin=coords[0],
-            longbegin=coords[1],
-            name=p.name
-        )
-
-    pcoordinates = """[
-    """
-
-    for p in polygons:
-        pcoordinates += """[
-        ["""
-
-        points = GeoPoint.objects.filter(polygon=p).order_by("order_in_poly")
-
-        for pt in points:
-            pcoordinates += "[{x},{y}],".format(
-                x=pt.latitude,
-                y=pt.longitude
-            )
-
-        # remove last comma
-        pcoordinates = pcoordinates[:-1]
-        pcoordinates += """]
-        ],
-        """
-
-    pcoordinates += """
-    ]"""
-
-    script = """
-    <script>
-
-
-    var streets = L.tileLayer(
-    'https://api.mapbox.com/styles/v1/{{id}}/tiles/{{z}}/{{x}}/{{y}}?access_token={{accessToken}}', {{
-    attribution: '© <a href="https://www.mapbox.com/about/maps/">Mapbox</a> © <a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a> <strong><a href="https://www.mapbox.com/map-feedback/" target="_blank">Improve this map</a></strong>',
-    tileSize: 512,
-    maxZoom: 18,
-    zoomOffset: -1,
-    id: 'mapbox/streets-v11',
-    accessToken: 'pk.eyJ1Ijoic2FuZGVycm9vc2VuZGFhbCIsImEiOiJjajY3aTRkeWQwNmx6MzJvMTN3andlcnBlIn0.MFG8Xt0kDeSA9j7puZQ9hA'
-}}
-    ),
-
-    satellite = L.tileLayer(
-    'https://api.mapbox.com/styles/v1/{{id}}/tiles/{{z}}/{{x}}/{{y}}?access_token={{accessToken}}', {{
-    attribution: '© <a href="https://www.mapbox.com/about/maps/">Mapbox</a> © <a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a> <strong><a href="https://www.mapbox.com/map-feedback/" target="_blank">Improve this map</a></strong>',
-    tileSize: 512,
-    maxZoom: 18,
-    zoomOffset: -1,
-    id: 'mapbox/satellite-v9',
-    accessToken: 'pk.eyJ1Ijoic2FuZGVycm9vc2VuZGFhbCIsImEiOiJjajY3aTRkeWQwNmx6MzJvMTN3andlcnBlIn0.MFG8Xt0kDeSA9j7puZQ9hA'
-}}
-    ),
-
-    outdoors = L.tileLayer(
-    'https://api.mapbox.com/styles/v1/{{id}}/tiles/{{z}}/{{x}}/{{y}}?access_token={{accessToken}}', {{
-    attribution: '© <a href="https://www.mapbox.com/about/maps/">Mapbox</a> © <a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a> <strong><a href="https://www.mapbox.com/map-feedback/" target="_blank">Improve this map</a></strong>',
-    tileSize: 512,
-    maxZoom: 18,
-    zoomOffset: -1,
-    id: 'mapbox/outdoors-v11',
-    accessToken: 'pk.eyJ1Ijoic2FuZGVycm9vc2VuZGFhbCIsImEiOiJjajY3aTRkeWQwNmx6MzJvMTN3andlcnBlIn0.MFG8Xt0kDeSA9j7puZQ9hA'
-}}
-    );
-
-
-
-    var mymap = L.map('map_canvas', {{
-              center: [{latmean}, {lonmean}],
-              zoom: 13,
-              layers: [streets, satellite]
-              }}).setView([{latmean},{lonmean}], 13);
-
-    var navionics = new JNC.Leaflet.NavionicsOverlay({{
-         navKey: 'Navionics_webapi_03205',
-         chartType: JNC.NAVIONICS_CHARTS.NAUTICAL,
-         isTransparent: true,
-         zIndex: 1
-    }});
-
-
-   var osmUrl2='http://tiles.openseamap.org/seamark/{{z}}/{{x}}/{{y}}.png';
-   var osmUrl='http://{{s}}.tile.openstreetmap.org/{{z}}/{{x}}/{{y}}.png';
-
-
-    //create two TileLayer
-    var nautical=new L.TileLayer(osmUrl,{{
-             maxZoom:18}});
-
-
-    L.control.layers({{
-        "Streets": streets,
-        "Satellite": satellite,
-        "Outdoors": outdoors,
-        "Nautical": nautical,
-    }},{{
-        "Navionics":navionics,
-    }}).addTo(mymap);
-
-    var latlongs = {scoordinates}
-    var polyline = L.polyline(latlongs, {{color:'red'}}).addTo(mymap)
-    mymap.fitBounds(polyline.getBounds())
-
-    var platlongs = {pcoordinates}
-    var polygons = L.polygon(platlongs, {{color:'blue'}}).addTo(mymap)
-
-    {plabels}
-
-
-
-    var latlongs = {scoordinates}
-    var polyline = L.polyline(latlongs, {{color:'red'}}).addTo(mymap)
-    mymap.fitBounds(polyline.getBounds())
-
-    </script>
-      """.format(
-        latmean=latmean,
-        lonmean=lonmean,
-        scoordinates=scoordinates,
-        pcoordinates=pcoordinates,
-        plabels=plabels
-    )
-
-    return script
-
-
-def get_map_script(
-    latmean,
-    lonmean,
-    latbegin,
-    latend,
-    longbegin,
-    longend,
-    scoordinates,
-):
-    script = """
-    <script>
-
-
-    var streets = L.tileLayer(
-    'https://api.mapbox.com/styles/v1/{{id}}/tiles/{{z}}/{{x}}/{{y}}?access_token={{accessToken}}', {{
-    attribution: '© <a href="https://www.mapbox.com/about/maps/">Mapbox</a> © <a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a> <strong><a href="https://www.mapbox.com/map-feedback/" target="_blank">Improve this map</a></strong>',
-    tileSize: 512,
-    maxZoom: 18,
-    zoomOffset: -1,
-    id: 'mapbox/streets-v11',
-    accessToken: 'pk.eyJ1Ijoic2FuZGVycm9vc2VuZGFhbCIsImEiOiJjajY3aTRkeWQwNmx6MzJvMTN3andlcnBlIn0.MFG8Xt0kDeSA9j7puZQ9hA'
-}}
-    ),
-
-    satellite = L.tileLayer(
-    'https://api.mapbox.com/styles/v1/{{id}}/tiles/{{z}}/{{x}}/{{y}}?access_token={{accessToken}}', {{
-    attribution: '© <a href="https://www.mapbox.com/about/maps/">Mapbox</a> © <a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a> <strong><a href="https://www.mapbox.com/map-feedback/" target="_blank">Improve this map</a></strong>',
-    tileSize: 512,
-    maxZoom: 18,
-    zoomOffset: -1,
-    id: 'mapbox/satellite-v9',
-    accessToken: 'pk.eyJ1Ijoic2FuZGVycm9vc2VuZGFhbCIsImEiOiJjajY3aTRkeWQwNmx6MzJvMTN3andlcnBlIn0.MFG8Xt0kDeSA9j7puZQ9hA'
-}}
-    ),
-
-    outdoors = L.tileLayer(
-    'https://api.mapbox.com/styles/v1/{{id}}/tiles/{{z}}/{{x}}/{{y}}?access_token={{accessToken}}', {{
-    attribution: '© <a href="https://www.mapbox.com/about/maps/">Mapbox</a> © <a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a> <strong><a href="https://www.mapbox.com/map-feedback/" target="_blank">Improve this map</a></strong>',
-    tileSize: 512,
-    maxZoom: 18,
-    zoomOffset: -1,
-    id: 'mapbox/outdoors-v11',
-    accessToken: 'pk.eyJ1Ijoic2FuZGVycm9vc2VuZGFhbCIsImEiOiJjajY3aTRkeWQwNmx6MzJvMTN3andlcnBlIn0.MFG8Xt0kDeSA9j7puZQ9hA'
-}}
-    );
-
-
-
-    var mymap = L.map('map_canvas', {{
-              center: [{latmean}, {lonmean}],
-              zoom: 13,
-              layers: [streets, satellite]
-              }}).setView([{latmean},{lonmean}], 13);
-
-    var navionics = new JNC.Leaflet.NavionicsOverlay({{
-         navKey: 'Navionics_webapi_03205',
-         chartType: JNC.NAVIONICS_CHARTS.NAUTICAL,
-         isTransparent: true,
-         zIndex: 1
-    }});
-
-
-   var osmUrl2='http://tiles.openseamap.org/seamark/{{z}}/{{x}}/{{y}}.png';
-   var osmUrl='http://{{s}}.tile.openstreetmap.org/{{z}}/{{x}}/{{y}}.png';
-
-
-    //create two TileLayer
-    var nautical=new L.TileLayer(osmUrl,{{
-             maxZoom:18}});
-
-
-    L.control.layers({{
-        "Streets": streets,
-        "Satellite": satellite,
-        "Outdoors": outdoors,
-        "Nautical": nautical,
-    }},{{
-        "Navionics":navionics,
-    }}).addTo(mymap);
-
-    var marker = L.marker([{latbegin}, {longbegin}]).addTo(mymap);
-    marker.bindPopup("<b>Start</b>");
-    var emarker = new L.marker([{latend}, {longend}]).addTo(mymap);
-    emarker.bindPopup("<b>End</b>");
-
-    var latlongs = {scoordinates}
-    var polyline = L.polyline(latlongs, {{color:'red'}}).addTo(mymap)
-    mymap.fitBounds(polyline.getBounds())
-
-    </script>
-      """.format(
-        latmean=latmean,
-        lonmean=lonmean,
-        latbegin=latbegin,
-        latend=latend,
-        longbegin=longbegin,
-        longend=longend,
-        scoordinates=scoordinates,
-    )
-
-    return script
 
 
 def leaflet_chart(lat, lon, name="", raceresult=0):
@@ -2500,48 +2239,31 @@ def leaflet_chart(lat, lon, name="", raceresult=0):
 
     coordinates = zip(lat, lon)
 
-    scoordinates = "["
+    data = {
+        'coordinates': [{'latitude': c[0], 'longitude': c[1]} for c in list(coordinates)],
+        'latmean': latmean,
+        'lonmean': lonmean,
+        'latbegin': latbegin,
+        'latend': latend,
+        'longbegin': longbegin,
+        'longend': longend,
+    }
 
-    for x, y in coordinates:
-        scoordinates += """[{x},{y}],
-        """.format(
-            x=x,
-            y=y
-        )
-
-    scoordinates += "]"
-
-    if raceresult == 0:
-        script = get_map_script(
-            latmean,
-            lonmean,
-            latbegin,
-            latend,
-            longbegin,
-            longend,
-            scoordinates,
-        )
-    else:  # pragma: no cover
+    if raceresult != 0:
         record = VirtualRaceResult.objects.get(id=raceresult)
         course = record.course
-        script = get_map_script_course(
-            latmean,
-            lonmean,
-            latbegin,
-            latend,
-            longbegin,
-            longend,
-            scoordinates,
-            course,
-        )
+        course_dict = GeoCourseSerializer(course).data
+        data['course'] = course_dict
 
-    div = """
-    <div id="map_canvas" style="width: 100%; height: 400px;"><p>&nbsp;</p></div>
-    """
+
+    coordinates = zip(lat, lon)
+
+    script, div = get_chart("/workoutmap", data)
 
     return script, div
 
 
+
 def leaflet_chart_compare(course, workoutids, labeldict={}, startenddict={}):
     data = []
     for id in workoutids:
@@ -2564,6 +2286,10 @@ def leaflet_chart_compare(course, workoutids, labeldict={}, startenddict={}):
     except ValueError: # pragma: no cover
         df = pd.DataFrame()
 
+    mapdata = {
+        'trajectories': df.reset_index().to_json()
+    }
+    
     latmean, lonmean, coordinates = course_coord_center(course)
     lat_min, lat_max, long_min, long_max = course_coord_maxmin(course)
 
diff --git a/rowers/tests/testdata/testdata.tcx.gz b/rowers/tests/testdata/testdata.tcx.gz
index e380513e867d867d4c7b9150dbb1e0a971ee706e..fad5c58be2eec72a8a29b0f1eab8f4138b3dab35 100644
GIT binary patch
literal 4001
zcmV;S4_@#eiwFquQs`v@|8!+@bYx+4VJ>uIcmVC4TW=dT7J%RLD-1uh4~tUg!o!Q}
zx@eOwHbBxX+MwIF#ne_E?b@;-O)mZSOG>HZBp$TK-XS;x<{_{~bJKit$aClM*KaQ`
z4qh)e+tt<jt3xz!c<}YfqvIEwZna*m&mS+>{c`i9>$_jO?eZY)yxx9w__FV>j~0tJ
zZ{C>mZgaU>o71byMf!4k`f_>MZ5MYRzIaWC`%@fsZ~wwszB%nzuUGx=-@WaZ>wJT^
zPVjcS`rr(2*j#PS7X*NdU!FhvX?Kb*R_krwtxuQfsO{=#f8J+Tr(M7MwfDM#1N~92
z$-Ld3{rQ1^@A}8x#nooD+%8_@pYQ)$*Prh9U*22_|JtPn0p0h|{ck)xd311a{O$py
zpOHVl-|+jvPuEwQ{;R`(_Ai&4!~K`(>!)WY1n9m1`GVjD@{z?Oawa<B2!8`O0ytiz
zo&MPWS$BQ#Q{Qd+7pu$VSBHOg?9GMJp)aob?jj%d)AICceYQOj_Q#7)w!S~$$!h-w
zJzu6DyG{3*@8zS7@Av)#-w9$kUfg_ozr#04)VlNKKbGn4f9(3@x7U}uSF`{9I{lKH
z%}>|gF1vpF!*cU{wSLtv|J~pE5$QBPcNedgCuq|}_J6*=!s3og-1_=?_jYyp>hk||
z1LE$ufjcg7>+7fM-K*OFezJQjmWk}mmfH{VY?&VL31NPW`@e5d+^o*e)2%K0mse*e
z&%5=j?jlkAVC(I@yVL)EaQC|{*H^1`e|vNOBE18D;a%MA_TiBqy%J*V@N2R=;DH2g
zQy~8dg*}P85q$7~;g09K+u>&P+J17`7eD;+muG@I$MBEOn9ljHuJ4oB#N);OuaB>}
zNEiM1=iPma`9-)3ua-K4NzyJ>K^@uPG_dCfT@iP36>(1x5nNDkr=&L(<R0$No$8Wz
zD;05fh8R_&d>SISAon;M?#@=kJtld#;x*i3NqJ|O4R>D=cT2l*-YK}-V5lH>d;oW>
zihF>HxJ%-mgi52$d*dd;y)JnVToHF~f(NJ3=bbBQ#xc!@yRV14fbjvXM&7*`nF#mF
z#`%iuIAIgPMh*8MCGGb)$$O}XyE6_5G<sc>qWb)N<Q<_R#ZIDO*12gB%97^?kf*xH
zQ#DgO9YHeXL}$t&N=e@h+xf^mOGV_ZiAkwf$vlY(1<A{d$j6GvCv#4agDO1ZirVpR
zKIE;dhdhFD7D*-Z$@nknx%2ax=LISv&&DzeYseS&-NhM^r;5lsgHGqfLFY@(^A_4@
zKMMI^B8k<6=g^0e{<{Z}Z=?Mnc|zj|ylS)$Wy#McLWfe3RY%qcAgVOH7)#pob0Tj;
zMdXv7&V;JuD~Khdbe3i`+G9oJqX`J!s?zQn`5BOBGA?3Nk<ZDw1<B9IZ(q}B&tlkt
z*Ua-|%aXS<B41I&hJyhBFB<X!C8P89A%1&qR^1T!B)`_7Dh-Pymh|DxihM;rJj;d=
zttLuKzG#MhBR>*(XTU1S2bhe^*EY{r6s5Jn`1C5O?D+`BHw|y<^Sq!*BCoRN8K7vy
zW=5lZO)@Xmq-P7HA&)p2-M&R>hsitR1%oC^%Y;Sk`8koV$ezy@8e7#o&x|Hcggmx+
z{%-P|b8i`tG@hLFqLj?Kwyyc@M<Fi;7=lXXg^EUO8u@X^dlQnxqlUcWCiC|s&-IW;
zFcKiBnq+t`D)YYS_M?zz!xlv~&j&85;rS->qmXw7BZ^Aqqb=&aZ8AR&`Caa<ss)-Y
zTT%HS^4NTML*#Srtyq<AZ_8Fx+R$Wv6!O-jp+0JKdrPjU%!g*xjY8fVis&@nTg$F!
zzI|w;{TSpKb}g|qy1liosBRxxe0~)2WE?|OHNlgz<mXgT=_(>G#tO4)o_9rSX=g=V
zDk2|D;;zcTv}G$Q1)9u{LEfTC8mA*w%hRI0Ze>PQ)M6_l&&E4ll-9;rGACzF<m;)i
z;hcNxJfdc_7lg^ke6#0=$wxy0Vo;HfO}C$Z#6~E0OHoB>g`w&8)0cU%CM4ZnCG)my
z)^%*Pw4;#sMpzfI5h+W)`S3;{Zxd{T-m0>$#WnKxC2uPtPbPR>b(MIWOnkn{{4jYj
zmQYn)C82CxUF(544*6gRQPnv?V%ZuS-*o#i$Rn7D#2W9d18ljs)0cT<BY{X0rFA$R
z^4vT*gXCS#x>j{Ubc989`#F)XX|#_95Q8Q@@0f~KQTrC3AA@|(y$u*NeQ%wmqPe$o
zB43d`Pu7q_R?YK5lgYj1RzE*T-W#+?s?pw2QG0$y<Q-K+K6g$a*7Ub`o{P@&b0S}n
zZlAO@j;LC-%mr>y+8K}+10bp{=Pq)K(oTQ0XKzGwoe(`V^7kcQk#29nuosml$D=K&
z;XNo1uW5M0<T-agRe5iHvIJwxz??jJG2S}V#Oi#~4GPBU1ZQN=R}`PO!B}L}<mCAD
zw%Cb~ug!-S>LD+fYiU($T6ig1KR+j#uSvrbHeLudK0NPSNx%K9=J|^3`6Ry_NmXOx
zy=xVf_at9YlvbhvCQyy`!A(ZDZ{FJh@^S=dywZo~W6@Q}XC?De5qWDYh)U)I6kUsY
zRx<A@BJT|XfySx}7$zh0jr=hAXj~fXHKTo?P}H8E(`X;+A&+FN#i+@;4jh}z-;;c-
zh`bmG%F*5yU8!(Zo}7wCd-TSOsJf*DS9EpKIeBkA)kB^j_r_w4CnuC$qf)4mAB8*{
zBwY{85Q?r^HYf5m&GRHC2u4i?W{6lcGlOO{&-;qV2O~-5HI)_-sLA})V|8-GXe^T|
zGb185@>7r4B=Xi!`lljKMc1*(^`|W}V+8U(=i#aP;YBXGj*ZwEbw90%d^EvfRLOkd
zY%JUQ-i|@ul5wJ|u8M`TvF5ZcZADSqB)`H+-QKyPR^6P)SLDgD!B8a7<l9H@i;niL
ztwlWs`F-xKZk`XJs8u&7@)db<Bw1)opdueb(VQID=J`>`I}<FiCQ6F{Wy7>{LSK<9
zM`BL6hoC}F(9qwr|7kVjy~~|%UDS>D2#iHj?jM34Dnl<7p_kk+1#|<xrOBYj*6(@{
zeK3NAq(YxmL!O9H-=uyJeQsw;$U5}S7Pal9&Gn<8x5f&pe7Fc+n$+JDJynF>=N4HW
zRp^78jMP_#o?SiYlbvP*fmNZ8zNnfvrvV=;LeFWgCs57x2w^hlYa8%ul6q$tkW_O$
zV%c(=Iicr@(8sLdN1fD{UA3^G9|wKzgz8w8nT2V_Y8w6w&`TchRUunMp`y|HIjMOS
z8Fw5Ec?>EIA6@h2-PgDa^`K8<H@Ms%PlMhU?SuQEj7(_8-5`2v0El%`AM9k%V>9js
z(ffp6*F+;{a}`C;4gDzSqls4O(jjm~jXP@IyfM%xWj&&<>;@@2{+rZK9=#anwdkpA
ztF$>uJye7~7!gu6z(z_Ip>OD??b7kNGp2PSD&vllf*=!VpMz~#`a$&8u-CQE$qQpl
zrgKuW>xa>MLqURW!1qwJupXP5HwJo40<0{!Lkc*Vj=0#;G6&JKL8t4Qg%pb}Tso@(
zpDIG{3|P`-<TCJN*0^bN{V@7qqA06x7;RCPPMhmTK~HER5vUThn5MI)S@Yo!qqioo
zvigSEO{QsD%hDM}?~S*@Dh;1)*>!HH)!dDNK4(}8>gwxR%C=LXrsj=;o{a!Toz%N9
znE@YLg4Q5<XDsU0GDP;VXyY_$cKs;mqoIgar{RaP;Tvw=yfM&Q%Kfu+-aL-6sd>{M
z@U0jSW#N@Ym`uFBwq)#zg0?nkY)CPx&=a+}e){MIOjIshjFzcr<gWSf$3f2~IwT!>
z%ad8&tT}Rr(eGL&dX<`I#Wt7D^c}gomWk+fbG?(vBxtp?%mMUa*D?`wOWvUM(u}*Q
z+w~#0Otd1Xa^%|3s&DR#K2(H0<PKW8t{OI$4c~AZ@W()pXqd&Rq+X!xa`ikTH4iF6
zZ+EQ~RBE0SZlg4(`g&UtdT+=ogSNs@)SGA9fIkL$0)xcJD)bg6v$kRL<_)9Y6<L!i
zUY{(2R%AT`sdu@nM%1BqRMfa@T{Xr*AF|YY-GJ|TGE(0L{6X}LCJ5`|^%87R4Zj6#
z$3bt6Lqt_q4T-jS^QP~}Wp6x_RRwJwxXGkyG&OG!eX`T6LtU1RL)QlUspsY;^lXH(
zlSYG{rc;8^+9(a9cg8|awA7&I&<6bJJ92GI=!sN?SI)-CRNA&Iok8?MhP6I?=b+8?
zGe9qfNUP>~Xc@UvFT4u5Wg@0+I`jc2(@>?Uc?0NUZkY&(QD@gjYH67>K+h(C(wpa@
z>~hA?dg+XU-evSrW!HP;qBVEWjJr|LN8?eMhwrIu2?n&F?HK5hazhoR)N|Pq3}`{y
zQP4|ns6slaw^B562Tjcz1${6q7<5uEu4t-8Q}aeaPsz9=WorfRTyx}3-K7(A%S6zw
zY3se8%v?_`Xgi2L7g;l?!Z$wnX4lUEeKeA+JDt==KN-88nwmF+p2>){y6{ba(Cqpd
zpqE^D1v>Ot6n$PQs15idpbrKZb^Y)|PTrpgdTQRhA@q4AYrvpG&v7#7+xqpx=-C7z
zUF0skC?z%g+=-&09|65Ho^4dQbb`Z@E}dDSmx|EGyrwPaq~2rE0pGV0j4{v)8D}x5
zvUEZyT50<r^sVM@7`^1`o1n_b4Y6pY?SlsV)=Os?eaLIt>LPa$iY|6GD-FMvEFFQ}
zO4fj0<;aaFlX2-l>zg--ez%gfZk@YGWrMcR=K4|4`&@maM9+RQt6R3Yeh__L$(oT=
z-n?ji(drv&ee*^^&$~s;RmNSE*jBfkxHm7Oce(INiQdJcQr{Y-j)Okt!YiF!9|N@J
zX%qM6CG@+Mtd)_w7_oK6o&kDm(oC;f0EGZrkKFsB_hX@dxAl(=$BSp(_05OxulnWo
z<g<(Z%kM6hm&<j3@cC}j=LcPX0B`+ozXn{MSrQAMo%dfJFZKsK+Mnr<-<Ut&d|K_l
zKYqMP9IrO-_Wiy~Z_noH&9+;gbr%=Ar<Fdg*5_wk-=*_*`?o7md4JsTH?R7aSDRa}
z$v4;67pv2*Pe1+U#ke{5Cnrnaefw&4adz`_e@^dXI`r1VdAH?x_n&mh>@2wb`P1E2
z>$KU+`~7wsZXWGN=RbatZu;cM<;AkwE)$3R>D{BdbzjHNy7l?fCnw*1pV{pG+;9IO
zy~W?Zy8LyS?hN1yfc^i+i(6mppZxAIojls1Uaiw}c$`o!H&4>j_%&Vn<k9~CEl<2T
H#DD<+Yt1?T

literal 4000
zcmV;R4`1*fiwFpRxZ!02|8!+@bYx+4VJ>uIcmVC4NpBoC7J%>m6@m}RVG#CP__!#F
zf^35^U}u8F$mBMtNENwfq?VvAFZuUFc000ctB^ctieMG!gQ2djt?R2tzMV&3zPUI*
zc-3ET*O!|w9;1QB2VWjPI(mNHtvBn<**E=W=&!%+hVIvH+aIKzH`^~BzZizA!`15b
z>(}P2yS`X&%*o}&Dt);<dC^~V+tuBNuU^sP{V5K+H-F(QU!4r=SL@;T@7@gkCg0$#
z6TI23-#f!=Hka3DD+0jPFVCL-v^&LT>&<rPHYa^LYP&w%pZDqINjL0%?Y(Z`K!4OL
zGH<r0e}3TKyWyMe{PKF;Z&$DI&-eeW8&3B7_cxcqzjmonK=(a#{~I44KRP%#diwy<
z&&VI&ZTS7*r>o2B;fu%r9A5O-kN01uub-SA6QKJH<THZj$cGjW$(iVgBm52E5a4K)
zcKT!gr`^@TPeXS-JYQe*Up)S^V{a~$4t;((bm#f7pZb%_&FS`7*dMJv+WPK*Z`b=b
z=vkkB>^9wJzLyU+zT5i`e8-64Xm#`H-40(RQR~k7fAs0@f9!_->#K|1t2z9Bm43<1
z<|mu4`)=6&&|g1WZ(a`lzlU2tBAw>v?)+tcj5b|l|L40atnRqPt*@VTZ`KztFaA$A
zAnuMExZ@JHzJ9XVy{g0S$Gf+pPh@Yl+<uU!eR{yhg!wV<|Gq`>dVO}5Zml0)T%I03
z>ozaD^F;Bzt+)5?PX7De-S4(sU9LC7?ald%^bY)mcX7Ad`$vB8N{F%JugUI!2NJkV
zf&51l_9X5`@ZJN4JD%@uhnvxB`_W~e|M1ISo(b+8!#_M@I_JN-VMt;Vk5>D?KD_2C
zUG&4BclRyk7vV0vTIviYNxN7Db!3C{z@8s;Mcm0%#63Yoa6!SHlHOF1d$>P$s!QIj
zRK(pGVpNUtX^7y0+~Z=nJ6jR=nB?7x*Km&|<(*+M+<isdE$zm6r{Hdbp@Q7;0o<`F
z?g1*|E{S^*DvdVpjhhMgy5v1@Mclm!9-Ky>cdn!v$Fv;oz8>xZ#s{<-dG}&uCfq9<
z=PR<~giQn+HQa-gwBMH`@1Y{@&Nv{@=yg$w>hsHycZ7-*JBfx_=cYv{OP(J<p6Vh`
z)lBhp2+5QaohgSXC4Dz+mm~8m6_K|lCZ%2_^CTt|BrgjhA1flC%sD|0s_={}YR9|f
zkhiWL@(9LRB$do3<G-Zm&M#-47pRCl8_OuHAz#>c7Z*gHDkAR;I-L^-oi91hTWF*G
zB;<pMBvuoiLmx`|?;b?HjrODD35_H0s?k1_CBK{q9ZE%39a$rQsM7FaENRa#iM$OJ
zkxzO$6RMK0AeM~MSz6F&j}?)RCLnmLO1o?17eJoLxQJ0jJ}2iEB)=fPeNCf1i(v;|
zGtZMPOWrPsd_@r(4h8_cXvhncjLzGK`0cq_bz|g{{91>qG%S)>(ucPw@)h~;EE`6&
znkX&#q8av${6yrP0jne*U^X&e+dN-Ul-360)2pbm=OY;3G`zXb^MWRcyvm+ufT9tb
z1&#JK$-G#To-L4uJmPG0`xd1gC-00G44Nn{6Bf1SmqflIdp=ueY*q6-GnzOP^4RA2
zyUBCTy=6escyiK<QZnn>y5_f^guEDF2r8KuDjKnA<fkF;O-K@t8uE^t%-@qd*Fzq`
zNPwVflHs|i%=@O>PePszTNKqiAGoN7=bOw=Lf#pSC@Ptcwy5{E$^10rce%H!7HGC?
zMdgFYWAovSk<YoeVpY1mEn87(LzDSQ$Xk<!`l!+EExDpHADUG+33+cQqSJV9ExV%m
z_MwgTQ;=uawZzis_SU+hx_xNz`ANu=aSTz_1W(G6Us6S-tBAZ9E6l2S-W9E-T@-n#
zh<q@KyD9_ImaV81Xfi(qd5b1#oQ_m2PmA`tl?7E%i>-(}8}D>cS{q}@oSY?*ucyX_
zbMCG4h?>z}5N0Fu&7L189}NYFK}9|`-G2TN8=>4SMHQtLhNj!kU*^S{kaT;M%-gbA
z*Rj>oPD0)rVO_*Vq%8U7!<&G-O|T7mtIE0-*T~<Oysd~lnc#KRRpN0r@%bk6<K)Fy
zLREE@gtB#Ytq0~b<bxqZRp$hWWovAF)9t4qk6<DaYrMA(u;t#)U*?gG1R_n8*5Q1}
zbMxekl6N`lTGa{B5f;_$mqfm%(LNeL44U}7V=7uj?OS|)3i3JkHek^7y>*s~=H4!e
zd`0#=Swjw4HO~vpCij+G{ro6-Z_pyCMtes^?fC_fcT^Gi+&O_*)8F2CE;`RIiF`%6
zebUxAqH5JL7q~@f7eHPNfT+5hyT~m{JO9z1y%EuMLiEtc-<Nzvy1fO%UR0hOkG7<S
z_n<tyrs0i~=iL2N<-PUE5{xYabN1xLc<WFTtMf@WC>W~~T#!9qQGDJ8W06slljGCd
zVrN3WHXmN7hrD2}rB$tI;iYK({E}q8CJj&6cp=pI@Vs*+{q~ER=PR=3ll*cdRgI1J
zu2odtlYB){T8RdjKsDM2Hyho)d2dI^%ORlgN*|t&MOPtTl*~&-<gKwFDwz*ZbS>&d
z$-JwGyf+908mlf~n2pRg^5f*AacQvEjP`*-QG0$#qkX7{Jd&{%qbBP*aBMPvPx7%M
z@?sz;M|)dzrNTvdaw;0_(Hk$K>Xs5*(bY+p<h}J&4|#&z8;doboKSX+N})!667p=2
zbUiRbD7tFdlE~LI&y$!S7&RG~A!5<Y3|i1U?<*o7j3k-YR9Zx!Ci8QT)yW~Fu}rGW
zjELOG&pl$3$Xi3{pNc#cUB@QZpSH}53CR1Lho|a?7rE#<Hewgl{j?_X(FBK4CG&-|
zv25#mI|X@5#)+=FDi+ShTGG0-6-8;2{0b{|d*_N;bxR^&ktfFnLy<s}Zy&ucI@-Io
z7WEY5_qn&ac|L@qR^5`wSLDf&WT7#EihK-3b8=jp=O-cWOt8e7C@lh%4bv_OeMPPu
zi8<jOf(kuBLx0cyr`3%2E_c3lQ8(TrFcwX@e+YW0482r@UUI_}&<*&OW`iDEzw1%-
z!3YwP3Vl)yc_v1ElloEgxt%E?>(Dz})V7m0*H41p8Y`&s;UaixQh!hMR1tcgTV#1u
zp$~30QePQ*cJ-i7cA5<YR)s$LqH5lf27IgtJ*T;zKsDDRgxR34ZNRTd>YZUgQqA><
zWy@`rgq|xxAG3xZby8n;)xw5;8uYmns$*4V7N!}iY4{63FL}UMg=`Upibm^~q~=v*
z+;K4EF{m_rbj_Q0U*j&+gFcPj;BtRF4SHX+5AK69GNBoFqv)*xAl6BJu(Lsr&A1yy
z?-P1m6OEkBRTMop^pl{ECR(LShrksz?x=b5ra+&R^@zH%8>H;`Z&E*d^kSUXqNlR0
z(v~FkP!alIL`c;D8!1_YzM-GDOULKVnAVA?j5|sSg3P3S4z^|KN6}luUe`V+FN`sp
z&PmO#A4l&E1qr$V-$T*DdTeUm6zDMtu(IF|Dd226;$ln797WFtovv#ZQY^Y~>7oXF
zstCO^U`dye%fPc)<EG8^<LHBlqO87Qv_)MyZLXgLJ)wz2pi0nUn$MbM&4)jZ-kQM5
z>KkS^o2F?kOJ^LtH{J@XG<>#Y*SVoqb2kP0oM9=btFLD%+fId=nl}l0HUb!RQt!fS
z27GJ@TBGQlv8Y?i5ZT9~jnk;v^^>5Fh9X*>h9Am?Z@78$ra*5g_s`OK^Ek$)=FNY=
zw_-q)g;y3~Hu3t}lCdia+S;VCA;qXdPt@l6`J)#wQMqt2TBf3ryXM2620feTkaXxR
z&t`eE=ExmKziXN3Rcf9U+gv*HcjWF`CZgBP^-gA!pw-ecN6?2|%S6yEd4twVGw$YY
z*N5CP(TbqTk!wS%zPT^@P!al&J80>;YS>sde8X+Pp8`FiVHT&7dV#Xb)$@YXJg5l0
z-L+Owsd-Ykjnb0p>up8oy&<a%+6qHaZ=P)f{uJm53=$)&&|8?z+J?=WH;#T+WKF7g
zeX<Bzk@W(k-sP?uQHS19QRA+4)tCl-$Wrfh1HR|kNPQddN6|ByAgqhmORz;X{1&vG
z2E8>75mj9^B--Z9o4+HMz41&|6|{BWW|OMX)Vxvj$xgElby+$NT^sP{o|~7@vk}Tp
z8V!1yPYFhAqco1*84EelQiGmD8}R4v$h9${CsGw&IU8qFY1^`NM$roy*81?BgErSM
z0KFI@t(xbdW#rDi@G9h%iI}$O&<C7NLzSlHji8UYWg;L(on0TPrDZMvJ(~bZZ=Q#;
z%NaxKr85b7m(fR+UGI^L*4#le?j}JWjYnl3zNfM!7|?>YQ=mu64ONs<&t*$6papFw
zK`*(X3hAWYO3}z2G&OG$^ue%T&`G_xqNy59&6@;0CF72itrfg;&5=8Imrl$r6G6MC
zt@nO5b3L`7?I`+OWX+%o-}vC0UB3YI(MYoHbW$JvZ0ve!YTg)nCL`AB!Z!gzv+Eat
zUUK0T=+I+P^m(bEHsDWyJ{Vxs^}`Q2d4DG8sd@9p(C3w`0fP=b$JwB7>(`H?XA^{U
zk-PMwl+^HZCyIuC0`$&!wo&EM2@XrTbQXnPDncLgnzo>mdXGg1eBVkira&)boW-EZ
z(g~qxrR{^zx0<_g^pdM@f+{07#G;k94;t`WFP(ApA+Kqxi`+#hy4cyGH2hk!bOd%Q
zSp#~NBR8VV#-#(TZ{8^S-AdNFb?zdS4cbDR>nB0)bM=i9J^R_LZrSGgQS^BwYerIe
z^P=@dt8b|F&6@;0?-nsv8Fx`)TitTz-n@+7<-#i^dKZgIeQT6D4f>c1uXJ{O4A7dV
z&D@)p(C=2VRz~h(#MT*m0qCtsGrej76ar{Ha_@`YPlf*N);~5Jt)6yQHy^&g9Qy6?
zC+EZG-<|gt{bo4$bhqi#gKjv0H@^R^KNWCY_!T}m8$Lf;?GJdgKhqz-F@L`KwBCPz
z^v!kRczOMH-|xHh_FP}S-gcYQ?)-fBw9?1*=Ipc^x^&)d|8^xR?~XhA>gDj_^7__m
z^3~PV`TC?A(oes6F>cQN(aF+xU%y<RpWgi3pVRx84!!ko-fnr;{U=>AI}2`q{&=_b
zCT;feZol1zn@9V>`H!Bbn?C-rKkvJ3pE%@CZy(*Q`#O5sZO)#2d;Hz^na%Fc{q`Ty
zTm1dYi(mV6X8@l8?EgPn-TGqx<adwh_|XpadXt{RHwk5b{cU<0zottcKl&dk%KkUR
GfB^tO#W3Fh