From 38925348778a9194684c6f0316ec5811a9025933 Mon Sep 17 00:00:00 2001
From: Anna Petrasova <kratochanna@gmail.com>
Date: Fri, 27 Sep 2024 16:13:42 -0400
Subject: [PATCH 1/2] Add tutorial comparing R and Python APIs

---
 .../execute-results/html.json                 |  14 +
 .../images/R_Python_compare.png               | Bin 0 -> 30715 bytes
 ...omparison_r_vs_python_grass_interfaces.qmd | 289 ++++++++++++++++++
 3 files changed, 303 insertions(+)
 create mode 100644 _freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json
 create mode 100644 content/tutorials/r_python_interfaces_comparison/images/R_Python_compare.png
 create mode 100644 content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd

diff --git a/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json b/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json
new file mode 100644
index 0000000..085866a
--- /dev/null
+++ b/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json
@@ -0,0 +1,14 @@
+{
+  "hash": "02b5f6931312f20d479fb4fbc02fb26f",
+  "result": {
+    "markdown": "---\ntitle: \"Quick comparison: R and Python GRASS interfaces\"\nauthor: \"Veronica Andreo\"\ndate: 2024-04-01\ndate-modified: today\nformat:\n  html:\n    toc: true\n    code-tools: true\n    code-copy: true\n    code-fold: false\ncategories: [GRASS GIS, Python, R, Intermediate]\nengine: knitr\nexecute:\n  eval: false\n---\n\n\n![](images/R_Python_compare.png){.preview-image width=50%}\n\nIn this short tutorial we will highlight the similarities of R and Python GRASS interfaces\nin order to streamline the use of GRASS GIS within R and Python communities.\nAs you may know, there's\nan R package called [rgrass](https://github.com/rsbivand/rgrass/) that provides\nbasic functionality to read and write data from and into GRASS database as well\nas to execute GRASS tools in either existing or temporary GRASS projects.\nThe [GRASS Python API](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html),\non the other hand, is composed of various packages that provide classes and\nfunctions for low and high level tasks, including those that can be executed\nwith rgrass.\n\n\n\nThere are some parallelisms between the\n**rgrass** and **grass.script**/**grass.jupyter** packages, i.e.,\nR and Python interfaces to GRASS GIS.\nLet's review them and go through some examples.\n\n\n| Task                                                       | rgrass function                | GRASS Python API function                               |\n|------------------------------------------------------------|--------------------------------|---------------------------------------------------------|\n| Load library                                               | library(rgrass)                | import grass.script as gs<br>import grass.jupyter as gj |\n| Start GRASS and set all needed <br>environmental variables | initGRASS()                    | gs.setup.init() for scripts,<br>gj.init() for notebooks |\n| Execute GRASS commands                                     | execGRASS()                    | gs.run_command(),<br>gs.read_command(),<br>gs.parse_command()  |\n| Read raster and vector data <br>from GRASS                 | read_RAST(),<br>read_VECT()    | gs.array.array(),<br>n/a                                |\n| Write raster and vector data<br>into GRASS                 | write_RAST(),<br>write_VECT()  | gs.array.write(),<br>n/a                                |\n| Get raster and vector info                                 | n/a,<br>vInfo()                | gs.raster_info(),<br>gs.vector_info()                   |\n| Close GRASS session                                        | unlink_.gislock()              | gs.setup.finish(),<br>gj.finish()                       |\n\n: R and Python GRASS interfaces compared {.striped .hover}\n\n## Comparison examples\n\nLet's see how usage examples would look like.\n\n1. **Load the library**: We need to\nload the libraries that allow us to interface with GRASS GIS\nfunctionality and (optionally) data. For the Python case, we first need to add\nthe GRASS python package path to our system's path.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nlibrary(rgrass)\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\nimport sys\nimport subprocess\n\nsys.path.append(\n    subprocess.check_output([\"grass\", \"--config\", \"python_path\"], text=True).strip()\n)\n\nimport grass.script as gs\nimport grass.jupyter as gj\n```\n:::\n\n\n:::\n\n2. **Start a GRASS session**: Once we loaded or imported the packages, we\nstart a GRASS session. We need to pass the path to a\ntemporary or existing GRASS project.\nIn the case of R, `initGRASS` will automatically look for GRASS binaries, alternatively we can\nspecify the path to the binaries ourselves.\nIn the case of Python, it is worth noting that while grass.script and grass.jupyter init functions\ntake the same arguments, `gj.init` also sets other environmental variables to\nstreamline work within Jupyter Notebooks, e.g., overwrite is set to true so cells\ncan be executed multiple times.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nsession <- initGRASS(gisBase = \"/usr/lib/grass84\", # where grass binaries live, `grass --config path`\n                     gisDbase = \"/home/user/grassdata\", # path to grass database or folder where your project lives\n                     location = \"nc_basic_spm_grass7\", # existing project name\n                     mapset = \"PERMANENT\" # mapset name\n                     )\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# With grass.script for scripts\nsession = gs.setup.init(path=\"/home/user/grassdata\",\n              location=\"nc_basic_spm_grass7\",\n              mapset=\"PERMANENT\")\n# Optionally, the path to a mapset\nsession = gs.setup.init(\"/home/user/grassdata/nc_basic_spm_grass7/PERMANENT\")\n\n# With grass.jupyter for notebooks\nsession = gj.init(path=\"/home/user/grassdata\",\n                  location=\"nc_basic_spm_grass7\",\n                  mapset=\"PERMANENT\")\n# Optionally, the path to a mapset\nsession = gj.init(\"~/grassdata/nc_basic_spm_grass7/PERMANENT\")\n```\n:::\n\n\n:::\n\n3. **Execute GRASS commands**: Both interfaces work pretty similarly, the\nfirst argument is always the GRASS tool name and then we pass the parameters\nand flags. While in R we basically use `execGRASS()` for all GRASS commands, in\nthe Python API, we have different wrappers to execute GRASS commands depending\non the nature of their output.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Map output\nexecGRASS(\"r.slope.aspect\",\n          elevation = \"elevation\",\n          slope = \"slope\",\n          aspect = \"aspect\")\n\n# Text output\nexecGRASS(\"g.region\",\n          raster = \"elevation\",\n          flags = \"p\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Map output\ngs.run_command(\"r.slope.aspect\",\n               elevation=\"elevation\",\n               slope=\"slope\",\n               aspect=\"aspect\")\n# Text output\nprint(gs.read_command(\"g.region\",\n                      raster=\"elevation\",\n                      flags=\"p\"))\n# Text output - dictionary\nregion = gs.parse_command(\"g.region\",\n                          raster=\"elevation\",\n                          flags=\"g\")\nregion\n```\n:::\n\n\n:::\n\n4. **Read raster and vector data into other R or Python formats**:\n*rgrass* functions `read_RAST()` and `read_VECT()` convert GRASS raster and\nvector maps into terra's SpatRaster and SpatVector objects within R.\nIn the case of Python, GRASS\nraster maps that can be converted into numpy arrays through\n`gs.array.array()`. Vector attribute data can be converted into\nPandas data frames in various ways.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Raster\nelevr <- read_RAST(\"elevation\")\n\n# Vector\nschoolsr <- read_VECT(\"schools\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Raster into numpy array\nelev = gs.array.array(\"elevation\")\n\n# Vector attributes\nimport pandas as pd\nschools = gs.parse_command(\"v.db.select\", map=\"schools\", format=\"json\")\npd.DataFrame(schools[\"records\"])\n\n# Vector geometry and attributes to GeoJSON\ngs.run_command(\"v.out.ogr\", input=\"schools\", output=\"schools.geojson\", format=\"GeoJSON\")\n```\n:::\n\n\n:::\n\n\n5. **Write R or Python objects into GRASS raster and vector maps**: R terra's\nSpatRaster and SpatVector objects can be written (back) into GRASS format with\n`write_RAST()` and `write_VECT()` functions. Within the Python environment,\nnumpy arrays can also be written (back) into GRASS raster maps with the\n`write()` method.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Raster\nwrite_RAST(elevr, \"elevation_r\")\n\n# Vector\nwrite_VECT(schoolsr, \"schools_r\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Raster\nelev.write(mapname=\"elev_np\", overwrite=True)\n\n# GeoJSON into GRASS vector\ngs.run_command(\"v.in.ogr\", input=\"schools.geojson\", output=\"schools2\")\n```\n:::\n\n\n:::\n\n\n6. **Close GRASS GIS session**: In general, just closing R or Rstudio, as well\nas shutting down Jupyter notebook, will clean up and close the GRASS session\nproperly. Sometimes, however, especially if the user changed mapset within the\nworkflow, it is better to clean up explicitly before closing.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nunlink_.gislock()\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\nsession.finish()\n```\n:::\n\n\n:::\n\n## Final remarks\n\nThe examples and comparisons presented here are intended to facilitate the\ncombination of tools and languages as well as the exchange of data and format\nconversions. We hope that's useful as a starting point for the implementation\nof different use cases and workflows that suit the needs of users.\nSee R and Python tutorials for more examples:\n\n* [GRASS and Python tutorial for beginners](../get_started/fast_track_grass_and_python.qmd)\n* [GRASS and R tutorial for beginners](../get_started/fast_track_grass_and_R.qmd)\n\n## References\n\n* [GRASS Python API docs](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html)\n* [rgrass docs](https://rsbivand.github.io/rgrass/)\n\n\n***\n\n:::{.smaller}\nThe development of this tutorial was funded by the US\n[National Science Foundation (NSF)](https://www.nsf.gov/),\naward [2303651](https://www.nsf.gov/awardsearch/showAward?AWD_ID=2303651).\n:::\n",
+    "supporting": [],
+    "filters": [
+      "rmarkdown/pagebreak.lua"
+    ],
+    "includes": {},
+    "engineDependencies": {},
+    "preserve": {},
+    "postProcess": true
+  }
+}
\ No newline at end of file
diff --git a/content/tutorials/r_python_interfaces_comparison/images/R_Python_compare.png b/content/tutorials/r_python_interfaces_comparison/images/R_Python_compare.png
new file mode 100644
index 0000000000000000000000000000000000000000..438104c32fa44b1d8fe6e396766492a16110c2b7
GIT binary patch
literal 30715
zcmZ6ybx_;E_CB110KpSzaZ4yx+)Hs99E!BX-Cc?ncPUQM;ssh-ptxI*0>x=@x8e@P
z<(J-jzwaOK&Sa9wOg=gL>^aXq=j_=iwU=@PxRkg60DwS29;pEUfG_|6AP|Cu{^rlk
z=MMB6j<dXuD*(XB@b3lmDw6O(e@WpcqwV&}$<oc!<gEq3)6<jF#?j8z%*5G()9I~M
z=7A_B06-5=KuW##%G&dLmqq=i_Tl37h_YsC1q1>Hv4{jOhhO6pR}q@HN<)CjjG;6F
zkpqVmn_^PQey}H1A))5lATTj>Fd+;;+|B$6qNt4lhQJdnEU$Bx=vo(Qo4z%BXErXW
z@*21wv~FhU8MwDEwbjD-y8b%~4iUxfAlm<q+AuPV|9<kLHpKV(&p31O;#pq#H^zEb
zegsDR8*+4VugXY{{*7WI>hZ;ycB9j9_;#KZfKQ;v+lv9P{5enJ%a1^PF)8i2Cl&c3
z%{hqD$K7;J)}(Y?_hXj=0t$-^ynRc$TmK$g?~IDE7}E2-yeE`;I-QDHXGJ<k{JW0Y
z*oSyJsUwq@8?5{y5PyXa(C6C`(@D8`|7QS~G3zcQ9$-kYY-Z_?jS|6cQ30`mqaHsp
z#zYN}DNCFT;d(B0fw-_l(ZidS(Dra|$tiG}nZA+?)y!+H1=+(odI|W5Cz>r2e(Xxu
zC?OLnna}=^=6dibjQ-7#RszF1sLZTq0AL6|8lZ1MgnAwS-nbu;5NUOH_(|2p@aVtj
zoDs*guaBA7jPX+g+i;G^bc|I2Dm;<XP*8>kuPXQW1Trkp&jVgP{{pp?amf>8eE;q9
zStr@&diTc*EVPNs|AA93#8lcxZeX@UkGeaC0YFR8LnuiOBb!Xo30g(k8^v>giDMtc
z(+?4Zj^3+Ua(9@t79G(*A|j%AtgeOu5u!YPRUmgzkV*~7`38T$Pb7P-6hAz#XXHEc
zqb2ujP{_k0$gSmy3<+(JoWhYQ3cc(FKuG|Nx6P(W@G&2hE>F2yD>$;&{Q@(j^rSXN
zHS}`lQwuGhxEU~VSnsT)&f+K%^S?muP)8WFEbk)nN+_neco10QMi+mGT@;Q_%M@eE
zCkJ{|F&ez@0;MEQ$*5%`{_9EZ??YxqZipzA(?IFC&d?ypsa@DAE^%dzM)qYGuR)CV
zYI~=24ae0`;pU3^|4<!=6Smv0oJ&XgV6Ftl$S{ffiQ&RZb`vADuiWFGA@#XDIZH9N
z#u?jlyhKE$SLo#S&|ai&_3XAPZ0-v@hXn$BEeq*8m+oq56-7*{i?+Snl$<%n2Keqy
zjkKSV)y0V7{tt(^QSY%;xXn>+IuE;3=_6RgVyD&jt>edVlkWN~T1afU1@-YTCQa&P
zcg(ht$l>s|tVO(N3;BQio&+j6p41kW;rWwlfT8mwiu~Ym7C+#fCuiuO+VtA46?rF@
z_LF5Y>>SVDjH8N9Bfft@q|R*JTI$Nn<EUkwQ~-6E*Ms2w<&Ay{Ba3$$R^;!?#*>Ej
zi2ltqqYI1ppY!ejuxmZf=~`f|-GeXqQoRwTwT3DO`tqFjK3`wkw*P3)xfB#VIV$Hq
z{}+rH_>?gP`_h1XX-Q7t5s-X<{FX4#A5e$ArUdw?1UDV3=WDwQmT6S%79rb%e~5d?
z$Z;ES0#XCxMQ;cH#kL&H@~DzX8t>5?fU^Pf@)PME1T6Al!F+m0V#v`>X~18P%iCIB
zNJtfY8i3KF(q`e;;>(hMnd&XG8sakjONm(mKm_=2aEwzIlxx?hyn22|3V5#&y3oBk
z7!e@)R<@EHJtIO#8^Zv;z{Bz5fr!pECI5kFBEap%7oqzRaZ$gfz=-zD|2XQTeAF3B
zHyOwVD=rcN1{imot>lVc{~jYj-TKVrEOYpxgAGO}vzqxel{^<{JW}p2S~|3)P4_<@
ztHoMlbBRh%4-kB~g%SP0!?-sd7|D`~DS9XO_A5i1FOy{UuXyS1bI9bSeRRCn@xD~q
zF~u_)TEPF-SQv6EL;oB`e00)8FqD=`ygS2s<H>hCV|*p_jxkF6({i6DVIQY{XCG|o
zStgV?_jm{+z*%l3gV2xse;82480;YjH4&CaMR=cy4p{}}3nZv+*dg(MS|&mx|H$^Q
zN9ddOsIHUtNV5@ViBr@5?{~yz<vH~{7@6DWPy}ffsf7_H7kGq6zn96h)n>aMx?@^g
zdg<LOeV}i=4E6diL8~kAJCN7Bvt;FkeNTV*t<KgTn+f75ndS(6iI{vR?C47qcX_<n
zcqD*G81gX+{%?E=Ahn*0aTB7a_>$Nq49Zprd4n*0^<?X~cIoD$(!h4m=w%a~dpP}#
zX93sU^2UBlo*(}IOoWa<!Fv7@+x-4x_k0p^{t;41AgUqfP@g_%BtO3EzYIkgjO6C7
zess(cI5@`y6UxS|`>ooNCfz;y?Gfoj9LpEz`fvQAo_eelL~nzulruN@|Ms16VjvON
z6bd(~<{=wi`Z8ZJ=$6o5MExPM=ofF~^kbuF;~pV?g|o{h64T1cu}X2Q|1y+qB)!Uz
zOH@Ks94E7-kB<+)`F`7{2<Bqr1n2h^S=IVGjH$RWuEg+v<{=1@B}1S2bce+|y|Eqq
z(Q*HwV_IyDlLC=SxxprjSnWHW_zI9Q*IwJ*?lY{pyx#}^Z}Qltu$5(|okn+QIk~5f
z$uDL%2`zelyi?UgXe!rTdWW&pAN_!Wo>L#QId4iE{m&f*<a`n48|o!Si(B8G1sj{+
zVB}hKh5*FPnxP{NFD2v16}xqN1lqzXYDCDc1^?%R7|HUXU&`;NWL<mp`2%`Q@8xb!
zR>p#ZB`~JtGOM|m+n#^OPic(Omgq{@;|T2^#FUFV93kPKcn*wTADx*I)`nC*cHQAa
zPV*Ds6Gny?c7(_qbE85~9Zbk9vLNZB9!SEKMMn!E%A4NR8XWT%-w+#I5_rapf6Z@v
zA4)+7?C7W8(zzOwqibE7oSb~~fzZe64EHQh_yW+96zS~tJT@S%$#~IWu%Wp~;r87@
zd~cbplM_*n{jBF)j+-qQ7KgD$Uu5taQHPzfrW9U`9JJ_3Loj=9N-v>dAqeiJ3-CRq
z->NsGG>r(fVF~dBdtpQUQiB{&(;&3C0;^eGFae}1uCr8#jpchuq=Xv)>x=K1bii+4
z%1lb2<Xr<}1!wS2>E8K1#yI3TG6Wb7qjr&J&%y1fFW3IkZ{3fi^Oy;~etYBP#V_8E
z==2}-#l`qyzXhj08HkvKKg-L&{1}<+4M%Q8mfI3*-GU>ZU0sHM7`^nPGQ=Mz+fohi
ze#`vn%n&sF7<j0dw~$`RrooYkHw+bd<%rGnVuelVdk8$q;OXypwy*4SmXJU=c&{HN
z1PqkfKovehS(6ft1ndb&wjg5yK4%yrNV$MwBk@#1j&ClRj*qtR<g*~y<S&hxJ42U1
z#S!fAw#AwgPtLO|G7?mD(8v6O0*2ZNO}%}1!u2O;3ms?=PwEB%&0P1d2w5;ELS^al
z?Trm&IpD+9ko4O75R@o>o(YIt1sn(haCliS77Rok`+mIQdNgYE5yRhOy!gbZb9H6R
zIZ{4UJ?hz0m@L9{rjnzr>ggW}*0uw~H<ERiDXp&gPMy&~PxGs*UAaeFxIp(!#=d6J
zBlfPwA1ihIwaBg{-f2%%n0(`%N&6*BkriX^`3Mz<bbCg=^^#;m{HzZ*kTm2S4r*Ys
zOz$+NowK7riy<KW5pyS&7haxN8)?D8f&SLok!hav*501|>gUf~2`C5?CD8Vo6g9Nm
znf&^lIh_hPG(|>&LyZ*u7nW35^+BUp&}*Xk%ioQzbUOV{Xz>!f?dKa3M|2jvz4-YR
z;j5vM2_*EUkD{C=_=>tcaTv6GU6ubzqe96ra8{3hL`&on^V5c^Lq4_>1_1}TQ(CH)
zi9tmWaBXRv6-Jgv0sHS`Df}4Y&oCQ^D7WW5uBE7KOhlNr#KTb+(L<knM^}M*g6zuk
zJaVM|Z0p&?J3M1{geBuAj$*zxUfXADu49*-M*Pb^B!#!Ow&Wcu|9ULZla_&FRD;Bc
z<a;A6XupeA#WR6A`**wqQKnqlwFmYume(Yo#eKYXUwT+>3Bw<J9OSj!=dyNh7J7bU
z4c=OlMj%;)O1|h_QXjuj(fDy`NZ`c-rIXrBV#L<?I9oTrQ%!<0wacr>3`_FU*7*sN
z;@Fj`2F3~MQWQ|2I=Y`dOCjWJe+Y^&GS6d>+=)A0lWR3C++I=NEEb1RZcHB&FvnfG
z`)x|op@kZ((t7Nb(WB*`gQ7s>+)nuk`<^SIEl0*)K~!Kxl_4(`)Kshg9Kdbhe7)&<
z5a7j%QSIGiZK*exo+9f_OZzvt7@5ujOM@;+bzjZ0PdpW9n@syQ-<Tkqi+(s9y~^|a
z);mH$8LAe`(G(>=lk?5#G&QRv)dsHT?E_a(&n00fY+`5QS9nIgOBk+nG9MkOc+E%s
zx*~`Ud74&N)L)%d)dRv|m9-ePZLS?9KrSqH?F#GddeG5p#aim&!8->LR>1CXNpgzv
zxdI;eYw0g@qa;ZdZ3Y~$Umnk1#^axhpwF66T`j+}pP8;3gRPCtI8*htb!`@fhb!vt
zg&m``axuJ-6{bl0_h7L$cjV}RHuLK2QseT<4@uh7jiFs3>jIT{%MCg(iY$<~PrS6$
zJSZ<!VA5W)gEDF_y6z`?G7GfMY`Oi_!sU$E#gJcm`ghTjK9~6I6*19oVf=52D5B*<
z1?X4BUC@Zj<&t5<OWd7ls)%h@V-qf#<MKHrRyT)A{sh4}(GK4;T%5vuV+#r0GN-Fm
zk(6YUuUYh23K>Yhd_N?#aLE~f+fguUO}bhx_9*M=Tx8(w<Ad#*lq!}BAuWPZv3ho`
z)ZsS?vyXSy{a(QpEqS!~Vr%B*Gt@^$to8d*%Sq_1;()<``zszWjr(nPS`T{w7f|y0
zwVF{!ZjSHUh-8RAfRq`7JiR(i`P{@p0S6qhpl~irg+}$Cd1MYBbXb(=wL!K<LvOi_
z;N9UA4h{~xXS3+M(ZhF*xmq(g@s1}&jP_5+fMRzM`hUtd#kSX^VlKfsO?&mPzQ<2&
z8-|Mu@l|*o^W`3{oT<JFRWKf7Piy@D=KopK+=w6unzpB8Q^3^70MSzWwshsWbb~6o
z(+l3!)zh2WSwM^3QMKgAzfCyv(&S2zXRDcUH|*Cj<Xd~b-XgZd4+P6dkPro0yQ=<o
z4pU>x2%m!3*P1r}9{g{`p2(kv@urDYHxL9WO*%n0m;5e5(g-QX-NFJet0UxiLSdDn
z>EoN(cz!gkhWMiQgM67^n!4YOcUxhIN-_hdUaD6~Mk<{b)qH(CqxFKmoAIWistOZS
zXHHD!cVaWOO+(|#321V%x3fFFpnvB@&JrJKhmJR0KRNweq}nd{X_&=8ssNXM%URav
zgR^!VUvcfnpqTl#4i7l0Vm2BNsM-HoNH?SE+iaOb%xgRg{a~nCwK|6O&&5gEoxFwN
zb3vKVK)>IQZWT&2sG`Olq*1m+nXG+IE&23KwBsk6&4+CdncgbsH1StdgZP2|k!_yV
zFh09(&6%-O12HOxG_~d+;hTZZZ)MKrqQPOeqT0%YHwO&gA>0L|{ie0s+0{Vfn{WW?
z2B)VpE>76>w@sFS%j(L2L#?Pcoi3NIXlSvI_46MuqAG<$^ysb4pc<01{fwP5+|u&y
z^A^JEUp|Fas##lVCXu*UF(gAiV#7KUZG(*Q`}7lyaA5=Cl<lrkz0DG$<kZy5-=VM|
zd5kljbaCL=4Rqn5b277yP`RAy@vo9@$qhna6IQ4u0B-J7Bq}r$k;2Fb)mPs9&vKGV
z?+CzeJ~{B>S>S~BnH#l98s`bfuguZVTo^RiPV_G+XYdxy*{_WCJv*zyG%}Vdz&s5?
zP7|>%<^Z8eHxM3PWO8$uDl_6OG#{DSlXCj2E265IT83V_s%He@GeB<>D{LgI;maz*
z#Qldx@tDs-@iw^kSJbn;VIV_W5$drK?b0d)2msLFuDZT|{gg1ZT?>l^v!!{1WVA*G
zAQAp(q;n%};tX!x`e#yFx4~<h6<;4@%s{n^cK3<Cwq(J?pKT5Qf!r1WUHwPpGaxBC
z;~L^#c#^{@D0+E4iiF5#)^MOeo%{aR@B;ZWsla%~p7Sd)zR<|`_lT7Ew4^9CPM)ZW
z)%v3emNpHMc3Q@^CR#>uK2^W_rmq`iVTUi9L6TRdc`VNejA`ZdBS-YxYEF)o`ELhN
zNsLP8%C0+DSZ%l{W3@R#gHq*%DRzr|52cKdG4CJtlB^Ev9pfhTr5_y-AN*u6I}IO2
zGFuHbU4#QMK7Nj+lkkr%3zNGuy)DXKGG6w@>H#X7pQw0@sY4Y*H_~A#kiX9!Gnj6=
zXe;`*ex)$7C8v~7hv#{ge)k9S0SZy`xSus!PVcbkgV8ya(o(S{X{;|0&halKXbDq!
zE=w}R;2)vJ@;8=SKSFjKyA@x2jn$a=B^zt2pirylW8%iO8&Bg<h<9(CyggG4RBjs_
z7o3%+sVTCoQ5S0YabS&6JLt%T01pyc&1&xftiD^;#gJVIJ_SKlbeg|!J{&adItB!p
zgyj8A8g3<QJ&WqatcVQK3{+Su$5r|XHF=HT%!Y6r*S3)#r@XoHiX%4iG|GzH-Ot9W
zec^~d65D#!VJ#Ybwp6it)2qps8J<O<U1k1Lt)rt1A0OLz?{OVXo_0o4u1>%7ek$)&
zo+&(MK7xTkf@=8@P$Of2wj@@3%i1DN3R!gs!6ls+xLH)vNWOokr>B>;@Ij7)=X~eg
zUGRZadm&(x*VP1{Ju@<zB%kMCJRQQ&J7r4#MQ+_Y|F|}0W?s<Sch~xJnv2b`U4;Jk
z((hz@nE+pG&8#~M`|7KU5|@0Ft|5HUO0R7mG3whlfUeVTiqhmWoAFCeEHTE7EvW>)
zF&0T@<Puni&H!#$Oe?>jp0zV)>+;)8W@VSfirR)xvXLI-?9g8F2YumDk_fkz?yhC*
zRBXc>rCHixdD134*%R#)uO{-S$GHT`y>GR)pCSa>f4sR94#UXK&Q`%AlW>OxiU94u
zRhtC4;$uyCDA(1f+XBN4?7HH3u?sR-6+|kt)w?i3Zq6nWynFUeU2%fHhOHe6NIct2
z^OaR+tP9vfalm9>BWAk5!R~7vJb!apU2(u35b!WT)vZOx_NL`<%O*>+RlKd)?H@cj
ztl%`zrF2=;GfLRb&D&sPaj%u0ub*{{b4MY?IE4|#o9yGp!Ap(xkE(d<Z-)u=X9L0{
z*Hy1M7j@s+6jj?nct&Tg&ODNL>JvIP^~lG)kN)go{;H_?Ij7@Lz+Y3NUSNA55P&7}
zcF}P^Y^-%YZq@p?d^J#wjlv->b!OjPm@1*(u62PtG*0BOWwW1x&yrt{>fLM1HLiWp
zj0`$U;}44Wh>3Yh7;m(C6t}Feof<YQRlZx6TCA;-b+u=zSuOukK#EgXbLupt#(pIB
zs<HATu(MIyrH3_Tu;SncibN@u8mY1>3~k%6!l>b$tWE??&yITwq9()t3dV|38QWoi
zj%Px!odnAZ3e5O?knxJEqXW5@v?K=cgOBGJVJ0q!RK>2y;z^`+l+%lv-6v<4`l4;F
zEoxeyjK08Vjcm0U%(2-srcWCyQ|;2Tfg2^qlN2|aPk%R&VyLqFpZF{ri=*_eX8>Q>
z&+R4Qtj!A}SdV6>vE#_Afmh-cS%8)0gPi4z48Fsep`E@0gmvp)lpwbw0Qx#Lk^ySM
z0%KMUmG$j)$tWMvDAQQ8^W#}%D>=C&W<+Np2Z3V#4u{lze|mSk{e3q(tMq%k7g}2B
z=pBDVljuH8bZA{vfwH4n0u<;^{RW*We7%V)8hpFly!YA1W4fy!QB!&HIAV5484tK;
zqHhuIhY^@!>=txohET*?eg*XeP6MRufQ2}h+5ilk@pu44Rsuvc=xu9i%mBh3<Oidm
zfJ_$b%Gic7kFRg$1>-8==uHtKWwWLj=t;W#!*fyircdDm6MK!wxjgRL_pNy$zydzK
zGc3PNqtU2PfyBiH!Qa#IByxHdIOR^P&^y)k<<qglJ1q=b$M+mrdoT9h&52K&_wfVB
zn6lj$qJqebd@C*fM3fu?8dT_19gG@mwI=yL4A8c+*TmX6S5{YqjT?AOY(1@}EAOvZ
zscI)&W&}9Fe&YiD)5s5jw;s=-hkyA6TIU39$0GzU{(Nd}6`Quo^$Rb2^%ee=mJ$J?
zB!=LxNZJ$hs9}W@W6xob;~ojltSiyqB)kGmKp{JTwnl)MckkfLQl{@<XoiOz?Vr_1
zw&y|w-~}S%gyXxZWip%nsOd*IXr9?J-F9U~Qtu9jb+1lI6Rf|jPG~Zn>Nmv!--4a8
z-0<q!kN3NO$!13dq2&$1t((DAt1&I(oB@Kql>H-G@)I-OXv~H}g0w*$kAS<t7Njy4
zPq&itEfMjVMmwpu&3a#AF6H-<(Cyz}@Bnn?4-uw|E${BcuL^y9_{Y=~qq8Hexh(ex
zIu`*>6h`!D5eK30Wp{vVVM}?=yIA$geqgR0s=z*F!M2goLQk}tepQ;Wcz>PelQN<d
z2%CfSks8rGtp2E%tLQ;fg6nN8u@RnaypaIfBzbeAZn+L|a``)I)OaB;_RiJ<SQV|)
zdv=?DyLPDMENjxv$dLd|ETVA1A!+ZU0erUkgV#eA!%A)7WgCR(4;%lt3r+Dx)^>JE
zbKfqo2D7eEnqU1}gvfoc2Yq<9lmLf6?c13?V^ss+w;pxicg8cResqpnXq02MZna<W
z0xc13EoDYW_A^S~O>0G|QK)1g`z*S-#?3f#-*=~(Zkah~9OF6a9~U`5NhJ{dK0-m!
zsM9DV!8jYw>qA*ZU&6WhBV8}YF?BM&f4rNiF?u2gJt==*Lu37;-yS<+p<8DEv(H^%
zyas80Y)YEG`N-R9YC-M_@SrO3=_>;R6;&FS?hcs%n+qmt*TG9qU<=-IW!C<%2pwgT
z?e_I^oNFe?64)6KxS>0Q^=MC#3beK=ahq?pT-jT97cT~kf0jOaZ;IzIm=>5%<7R>c
z_nAqmBlsx?qTkS)QZT|QbLtD|lOqRUAa62THgJE4666UBbgwOQETkHEThVfhSKNqJ
z@O#wCQUfh4()5)7g*)G&pXmmHcA!$@0GgP0kRxYTRQ%N2Z|#OWHoqF)05&Wi%ns%o
zF{ZxnU*o2Shj7zM-jilO9l;{E0L}Z~>{n0GdCrF}g<kX*kWUgAzT+g%EwMwei}<fC
zHs*<x&1d7KgjX^Y=Q7T4Sgat{z%6gjRVxSHQ1x%g(+WJsP)g&ZkyPOz3gmO5e%<cI
z78X8DKr{hIyYVAtpM@N(m3JMYY2-w4R6-}h3@-6%HFb60F)QeA0<N<!7Movu%OoPp
z)9i0qlYS@bhiQdDU%oK_X>Z$)(-G@c6~C?`UpPtN_V1|*qR&cWhowq(|A^$)M9Zai
z*Nz^{39$u(m|~8)en{xh*vZd)zyW6^1FkI-Rw{?y<oraekK{85OWJ?)^ijbSCLhmC
z3BKR^S&koPW}n~s6-*bHD?dN!n*!+LEceejwIA}SaEOWn+8D?q0QQD<YMe82)-Pg|
z9Ge}^v3uzxRwk!uP6D-0e<~nbWU7!XS~f90+jc750`10CpCB+$VqaTiyu8OxSvO)z
z1AO16y4|PHcy*#@X;2m!`eKWyBtKS%Pk5pAfHrK9P&|!K->ED~e7DQIi(Cu#G8RC{
z&CShzF8lLO@N*2*d7vyYZxS&CQzWSS&<SmT-ai9)ypUKxH<lG@byq|TxlyzUG+&<!
z3pSY*$q{N1T^`@d*z0-5{H`JLc{8;F<Ta8$raRyz`TJ}~KQT>D_Trx;5He2x6JIRV
zH*Na+c*9O?nb5wuT0z?v)da>)7)n$jdwAGoZy@C7Sv+$+TbKjAa%7U8;AL2VzxU74
z$^s(cUc8{$yt1<WG6j<M(_I0&fC@gm98@q}{eCDrk!F&ZGQ|A-#N8kQ`cZ1$1?^iN
zs^3)ofDwjm%k5y-6^=v!`cCl&QdA$Aa!;fa(bFrQ_K|g%Y+qxL(S0t5S?l+K94a+8
ziAJfe6$Ue+(U`+0V*}qt@Y%J`*42T|`3{>{>weu9T9fQdoVG#>%|y<gI-!E4avbD*
z9=e!Ef5lkayl>$Qi&%kU1w=f$Ua#q1eOBJa8U^N=O_u89s=EFr>bN!F)+Np%sF;5D
z3X8VsY9p@P_t3|IZQwweF3|Y-_lgasRYl@_4P62RdsO6qRM0mAZq^8o?s9T4>a~2~
zld2I8-JMkPskJe3fQ6CW>xqmYnP}2|DHKPR-%Bq<HP99ljb@$to<;_7otkFQOrRF<
zxs~z5I2u5?Tk~6Z#FOhMtnub9XJc{?cE7p{URuER22pSNM7vo-I;U{&s77WO2C>LJ
zu);CF9%Rta(0qMz^NY$w&Buy_T+VXxsY18AVyHX$lJjb$Q9tGf+1n?Q;}*4i^^Y1P
zugmw8e|`Dwc$^Q7Jg2Y~^7SQQW5$o0G9h_Gf)Ou%y+t4NnGg`tWrdNv&_4}J7PL>T
zT$MfjQu^)Oy8RlD5TSfmVv(`B3HIvzVyEfnxVIHjPATqt$t?ROzUfD0FYVPv<%c#`
zLU=7z`T^jDd7GxY#H${>2m@9#_e$fukf|~y)pG3!ZyYTA^>P6DuS>MFyzt~KW@zT{
z=Y>G*N--bVmqM$gxpCok;19dlFBbd-n}^~jndqiWu`jy5U=VT^E%AP1!qUaUzo)65
zE;z3S>Sma5?G^tV`mVm)^?$tpSfg*5s6{N#A~AYm#S06BBnT?FjDHm}7>S%uA5!<S
zw;`8{6ip5qCeBqvW=iOKjQaq54|}P4mT}*a65e>h^yJ&d)q5JeNb;0OC1z<Q-!d*)
zUyC=3A)mUZ7RZ=ay_;R(tdP;(Nf0U3=7*}prU_11GcV!#wbil7w1Gdh3kDC9saAIJ
zd-5&)hIKHvYO2EySs6?y&%yXBL8e7%!`Llr$xeda)eJzF5j>3Hxi9E39w`z__>P~S
z2~BH_X7Sla?%cVc9_+9?N!pggC5dESO)`$~EJnbwgU~#;Yl-4{QDw&xXYRkO<~XFN
z2T3m6M%wjE3^9d=K6oA6s2f|;lap<YQHsueqOo-q--@(<oRJWX171T&eRV9I0Vs*b
zz>a;&KCPNA$1A1!lA``Jx~EjD*iNJ`*cW^dTYpm)yJ|T1HCZq$<VpTapd+_sV6z8}
zkYzWcSx=0@I*vN6J5G-1)fQEH=7sSQb&BXZ(ImkbCdh2yHST?wEGad&nGRq21@<?U
zf9Xz7gKK^;nLShV_gIPQw~e)-+oy>*6bTS1JF{$pimN`l)6eJ>=rHCFNnPBc1xX4i
z#yGn(1MlJYDG%P*24Go17;l+$J^|HZ_Jav{{c1HZFAVt0Cx$|7b)3hn=aQ}C1`;Q7
z5P1=o>3aqPXF+;C&fn+bo0f9w;-3owBE4>hghpOlPKgFOi0ih}1p+6YtbQcPr_FUZ
z{tK@r_2f*ZGJ3E7(*}5CM5dg(XqLY8Y_}TV^lgJwv5#9`%&ER6$}=%f?7OMcy<*Hu
zs^VgD?F#4P3lw?T+d73Z%IeMiQcdnTfF;$BRu4%whi6Wz)`H8gAX#^VucC=41#HLJ
zyp`hHf_{F8ZxhZYH6let_xrQ?RV?J;g6m$xT_$2EJndg?i}OFy=$Ptf`jn2;e}#yl
z&NbOjU~6k)Rg5h6%$Gi91MSp&%nb2Ukxq^3mhI#G7a6i}?xF8Dfmupz2SF}AoCT&M
z6kg`rzlN@Z{0*k*l&{&@wfosck_H@`JZs1L5i>=+lLZ^G*$5TpSSrt>fSwrEn{gHO
zJ@3gR9sr$VYGjpPbDaXO^ColhC$Bqb;y0^n=}=l|c{r*W|6Wy*|H0~R`TpVET2*;j
z^kaVT(%k{i!`{}`Wq*_VHkDMmaZ+niQd|A(H~QvW0vks{PLQ*-Mx)^KU+is)k}R;t
z%%(Vhl|xOAQ&=p0PoLw6`G?{$>CD*<^9Z3!95%@!WKY3C61w|?Gi+xx=1$W^f08UQ
zj8kRaXFlA1tfhmg^YD!R;<%3p;xPTts*~@PsraF_d5F2&?S5+v+sNvPz7+V$x!LF=
zKDzhLdkt*KOkn-FFXIgW$wV_gjMp_yq=Bp?!8qIWu^Y9yreD6sb76l;;~C4P=;$<=
zE4*(pzIJxq;YQ5P&WZ!MkL%_Ryi<ZsQ$+D;ooA|*b}rt|s|XF<E<`cY7A*c4m6OWD
zNMXH8ihRZi)14WlXn%y@UoC*Zs$Z9kUpB3cUY_JoF}CG=YZ;MicFmXea{BGCsAJ3l
zk8!@=@q|@Lox|PD<{a)*8nBmt;BZWv)!k*_$VjG4FUND~)~KR7T#Yqp(>26r>>RSN
z3k}~ac6^#H;=NmDB{x#WUM-$XWO%T7q($DG+q(4yze(~ZwgU0naX-X%!?CH{ocn$X
zqT7JLnVDp$<9b85d{OkaATo-g;qRn~D#n-q0QBL4AmoGFlbdR%mixnU<y`UB+pl?8
zJzH#~4j*>5{YAuKx{>l!{AYM);}Z=?a(LVNR7`nU8nXK>VJy_h_LTy+!F}59=UY@o
zyU!$d;jxoBtwiyL#!0#Vf@5V08<gCr=cVV}%$6$n_4>ykNidmS7~sWjf}?4{AS3jt
zzk}rWSQD$j<A~SvZTfg{V`+)4&RrcP>ba^Wm$e_?E1-Zt6~OZta-mo42c8}|*VBGo
zpT*yj*VgxnTfTf|^2>b&4^IZJ>u;pdgH<rtmjLVmzKA&BWp$)&C6H?JS+Yig3kK;C
zU!bFR$<m?kg?w&#qo(VK(&W11E{ip9i_N5Twy?V2OeXQuVZVT^_jR7MKM68E?ucu8
z*?fR)^nTKvBq}-x6OYu7W;wAcHbugImcvNAx*)d7mi1Mu$ll4Y4XatJN*gcV8nc;f
ze7fk{^*+_VloN8(XdzLZ`m9dJeXS?j027{9K=*>0@UyT;+g)v1go@w#t;y2vJilTW
zb)3X~x*qVJ){54{F;p`S*PSc{%P`6}{h_&B;@dOmG>IAu=0+Q<(+qD^1s2l&6@T63
z4wUjDSiwO#e8%c$6l0rss9bb8W^HddyXF0#rmq7kd@xAw0EG;L?RNXg?yGC_hW#d@
zXTGQv3_!wSvvoO*yuT%Vp3yq#>D^?z?~iW9jH)tw5YRH?y2OY^q0173onrgOtAoyr
zCY>pUMr;)^j(H}<)sj=&T`5*VC1kiq2?2^7yfkita2S;1XP)%B5<YX+aYZseTM1RD
z5BforZSF*6?5gxJJ7}a_r+Q_&xYCCKJZ9J*K8UC^Dol)Q;78ODxM0EF)*k35+$9tL
zV>iCGoci(89erj1U>7C2pnGyovDVR2t=Pw_`WSBhgY(WRI{vl}ue$lx+C2j}nghFW
zqE|<kmm1fpBsUb~(VOa?hGP&aPO$O5%aG}v+8j#5vewm~iVCe#CW5WU(|NF$9JW1%
zmct~^rXNQ(3rCh;0_bl2dRYTR5uf;cFPeJ|rIO!4(!@vmebvd`I4~p<D2i_`_W(SH
zDJ-y5WVgPhFE;G?+Ab$G=J(5A;+k)L8Q)VF!QM<>s&l57az}DijNk%{?{z$+)>xxr
zUf@v3>~TsoI#@>UR#A!Xm|F&3VqbB7m18xYPBAF!tOnXsX-x=*iv<9D`=nzig?5%@
zEF<i&xxk-#zdO&5z>ZtHUiz))?5p~lSF;$VaRt1i<lFdplaY}jd`;?car=wYZnBh$
zj?vQT1-d~4wY*UBGW5NOoD}9N)l2GQ6~ttwU?(~w!+H>&Ai3?)J)6uZqxWS^PRO4l
zkWSqrrOWroq|?$b@q5+scHG0H)f3Lcb2OdA_yW_}I%sC1yai~50a(T$?5^*Lqq(YS
zovM{OR_a^8P9*~=^bXyKR<VX{MyO<5l$q9kQO)ISYTowZHWX^?oM*+0fKbs8=kXO-
z&SW0#p)fqL%#1CEWRTZADuIwcx;W6vac+Kplch*1BM)e7<b<0!m-_0ibGO#Z5v2HN
z66z^8lc>*KvSYz|O7%@pH>0CzymDXC!fFVOd$va!osvG5WQ*2%860h}q?qYbauv<z
z!w>47Kc43h05D2q=w)N6gikPi;&&e(U<}{|6|PZ=*JESG6N&sCJ0!Ccu;}MG<ijQb
z7UKC-Rld^29?8lEN*+lGE!{?leX^X>qlfNg>AY9y)}Z=ovX}#b?X2m_`vz<Va~!_U
zj|SBUs}kYw-Sdz)hm*J?r9EayNGRfbdEIT9|C358$9zc1rhe;~N@F;~jp<Wng;8}*
z05!dBNMuM8ZqonaOsA9b7zC${sT(sBN{9CSP~kbX==v<GaVGtnL{-$1m4vf;KRN70
zXzZ|c$kU?V+h5iB4#e{IlbEGlf1tZs@aN(w12&~YFXQfCLl(M4w;SzV`pxTJqic%6
z$7^65ojddB_)nmqp{B9!-r?agXfU!HwWB&252&zk%ubQmNi*8422wB_vzqnB>02K<
zw08Qby)lGLv;CEV)2HRMud&AK>+LmysR(Q{2z!(?N50PhUeF(fwt2wec#E|1agdu^
zJwRiRC+BNYJDd)a5d)#Ww07J5d{`K!Oo8`5VO8s0=^GJ9jgMR<3o95e>1vxT1GTBM
zl^D_G&F}cByx381tOe<RlXZrf$<e%RZvFBX_d@CqIw*n6tndWvV;hl{zkk@eL~Cnn
z(T$^x!sS;|yuN~xm3-DE6D2~vBDx2ym@NiL$l_$p5Qv?0w}$z@<AfC!=O0?_xSqHD
zs}q3L_MGPP<ag($2)8II;2vpif(2*{B0TPRieZy!ovZ01FD~iOEoo(C6@W#k)$VLD
zg_}!Ikvw56NUwF?%}iJR@nFy101eDLIMfdn<q!Ly!{zpo(ru&S)kZee#&O{XrV%US
zwjmvVI(MAVJzx(zldOE(RNyIwVNqm)0NG<ir!Gk)W7|qM%qU}_La=O~BNRNS_+qbp
zZ}#{67oVk%-aIIh*TI;tS7l($e<Vb9p6?UVjlIOPKBpYQhuDfoQe+lq4}ShqQC2lD
zF<~z$i1J}2?yae<EgCYMWShrMAp?zJXZ7O*?NN^X5XpM}KxLBi3V0{4TId|9ec|k9
zOtECxS&&3Uer1RCRJxn^sRx%>$w|>KB)dSe7!~|<rs_4DhOc{92}2u)|Bs{cK<<r3
zJJsyVfT^6iex2R@ufGBJCgSD3x#rTA$?|=l=~_9d#g~r~FV-UUw3n!zUXm~G)noi2
zqLrpb5U4hecm6b-t9*N|T@7p$933zgE*5rrUCjzM{`6Zmma-O6Q|^xWzcv2yOPW5K
z$#ii3YP?s_<=w-9XZ#%%ifBW0J&M6xP*P;0ZrU<4|0et*Unqz#J|s>)^lG2=jT;?P
z5i}oN(tgI@ru+p>7GJAg7x5hE7X$8J4)z8d+zp*;pU=J8jF|L2Jv;4ac$TR>NiI8H
zi8EvBc{R2pOH7wU#78GR)ihm=fz#5NPNY6bEtg)E?e^TG)`j5Y8n{(cBKzCGf2!SX
zehYR%sj^E7`aX5*{cridCPDj5nrRQ(ASfHSJEyU63rC>$mUcfmAQ%S7!BhHF9~cz$
zC<GF_uVaQ_muOq1m-gi_6ReM=mGEs(>z%(mCh2ZZ_Yw!+AKkjY%4IcdZN+i=-T(^a
zQwEK2k}>{e@t9T+OKpwu^|odMa4}%@6aT9@>hT4OoqkQ0IWJ{@9^R=437r^X74I$n
z6?i0_=NJY61WLZDlekLETxu{$b{~DdUAUIpxDq~>13G%s<ys8sV~3i60XH8eyeS|%
z7!j*57l)_{MJJ#>fiYv?1bYP~kSym7llC=hqD>;#qefzrJqLsn)bhP_w9{F|UcG~7
zh&2J9IJkm(@IeuKb*BC9{7eYp<F&cRk$BKX5(+zxnl&vE3ZkDP@Pe9tV&@s{>Ied>
zg$!evOSMxApCnx1lM#gf3JokPtJdtUGlYb86M$BufxZI<5srw}@~4I`!UBEwMyF=y
z=0eaF>NK&FNj`LQx&6)(5yxL2u2!*N0r9!jOf##*OFN(;{*f7|rBvAZesk4G)t8DE
zvqjhfKX{o5pqHlF2M6UdAzmzaR?|&9#sFV~^iq^wS7}}Nv&~RB6jnjuPz8E!LxSQP
z+XCB=d^fn-;l=62Q+R^jZ8?pq_g3B?p*TI`k~Bs4B}U|L3w(S+GGHAqZ3|i;@5QA5
z4v7s1Y742+-)M)WBM0Bw5F_mgp2}=!tm<KT^7#`5(OPy6Wqx0BTkO(`odFP$42(=u
z4vDn`F?%0>7Su<UrKBkJG2(@SPx9vnb=;CYhXAMkwx%*zCGsGtr>;1m#pdNH6z73l
zlzdH<L<lZ0r6R<9lJyx`;2Zobuhl3!|9~!s(~iqHJxBv4sCt^HRUf|)yXg9rs$-cx
zV$uMKFwJZ{yHtF{Cth83%c`wE_ktr$>@>5^lbur(usrlZ2t(ijFofn$s)^790-wYN
zg-e#{5)nSsG8fBkVxOT{LaDh^O=7=vs3dz3ss$l=ty!Re(SaAVRZIp)*HCH>y!s2e
zQXl-+OX>u{^CQcLEQzIHizFy)AJ5pMaGmlIg2rB6Dp{sGVd_+A1!?c3TPs554bYz$
zA9ow;#GBLX?27aJd=SK+4YB<QI+A_K%Mm?3;WtoBHIRc8Cw-U@RpEA8Z1Equ!*a}3
z?Gn(IA%<vVDE(^#F$D!3&+(X1d1|82_Eheop@BjQnSl)`&FO`OGqqQ*Jc$CCvaOTx
z_j>0!6O-(W`HA)nxx^JZWV-<fm`hGw0j%=%n>Y0D`zm^-ekZmfqtLz4)bMBS$K}l8
zctiw9uhW`V*r<HPGu*uL$?ya_cBt~w9h8JII>`U-+JU)~lh!RUC3-JbP(uguwU-zY
zhM``V>5@5Cqn@WvyAdW?+<Q#|LSJBs3x9#1<uP0MhR;7>*QjH9iL^RWa#~y?m>?!%
zgtFS?teh5>74Y)CX!sMA?)O651y4OWxdWH%Z){ITEJrw3Oso%cgdksYVmtUDVnL@A
zWVK;eOeR9VLg*`aBI`jFR1lv))f+@)H5n$8Lf?*Oc0+v954rV*UD!Wc`NupkHE@Li
zN{E|P1S&?Cg>5oCNREg~nD?54)g6xuEG;(q720QdnEle<E8=0E`&>XyQ3<fP6nQ`o
zWi>A+!vLT%Ivn%R@*?{rl3cAPpG6wAEgxs!J=H#CUwr{R;!*N<AIyybJ9t&LVrG`z
zQ?j^2P)5KCkwAyUVMA*jvXRd<iCW7QRuX>?gA1D&VYNd6KOl`BSCCQ}<>a<#$gH6X
z;^uBTG;m@#OwhT8Cz!Nv(YlKej!G4bY3^#8BtMd-A2;*E0|F#f14HyJklLomdMNn|
zX0F91^eJ3WMAU@XWG2!+DZoSKZa_Cmp&Qez7+N9M_d~0VG464|yM#|pvc#anoO-gw
zyOoq_-yvMJQJkX@Kn$S5VXW@W3y-rm5NU>*QM{d+2ljF$2L!BYD4xv`^MvR6n-f(u
zk)HToAM+;jvU_Ay5;g(CzQ3m<-N+T4Xa7xyRF=qU>BI;Q+PrP)_$Jhe7xb*Mgb@tD
z;vv9xj_#Uz_`@(ajXi9}S^JYo1Rv{2J;t1pC}sj6IKWPNZ%L}u9LDt;Le5juM~&XV
zI(3gcL9(1pM_Wgm;L_W%HmK?Iz%D?QL|hFgG@`gjPzFbMFefiwyHP7&RekO4;~#*)
zd&OJP_X&Qbf|zFPL3zqJwkGMh=#CL{x_9()+l3G>xLQU7X0Gy{ZyM_T<MP(3zk(ue
zL1Ujb)s^Yb4(afvxyUfU{?S1uI3zA6I%J>|ZIV|Ia&(G-xi<0nbIuD=L++Woxa7zC
z7>U=Ho5mlp)@ke1%30C?6Z^K~u@;AzSKW<F@a5?*IEthggQ5vKq=iKZelD;)Z_NQk
z7=&Kbf`2!qk5E%a_?u{laY?UoJjrF(s{0Tff2ZvvZE^3sG(^zyXp_FRxj8Wu0DziZ
zV7xS#yxseIL#3szf-X30Uhg`+qgWPOcsCm_3~w@PuPhxyA)P8+46E{&K$o<^If?nw
z%h@(Y(t(8B`^%-T$XJ3}Y^D3)SNuIP!6sLwvdYPvSm0|itbK%e{4-F#F?-nFzDXhP
zDuk=r+AWz2tY1jN%6arVz{u12P2Fh)-}Cy<ZmOBaE}}r}M^5T=()3jIe=tg?DNqtL
ze~mDpHw$z~>YHt|SS(O`k;IyjK}v7PCnS5bI+-~<Pm{athIjP~Mf%x#Ia9mowk_6S
zjwdms9TR$!TDkKHe4-cwa4<6Gz-bJ>g1~M-SHDzAphZdy=-8C}W|2J8>yU<S?Eqp2
zR5xyk6<0zdyJB&cCzc2Uuu&6$?~5)ww+vAsc0hZJasw@CGz0%qnd@Z5vJR7a`<z1u
zEDbI69vy&(dMV{dKQ)5%%}v0;SyS>hH^B^Tb`K0aBWbsmt820*gojz4;*kNEQzeSK
zpcHz^=yF;5_1O#7*xC}aTwJazM^F-{8aH-HXAr3t<dN=xM+CI*B?+ZNLqhW6=dV#l
z?ib(2#_B@@<NmZ$`8lbM)<^^fq!%4uYt(3hr0Fq^b*$J{JzrP->g$s=$3HuZR5g1t
zl|VnL*ZAjYE2f&ARZ>9ZrEPUq-&&XwR;yd$Oc(p2*da508-#!ntbJF$S6T$EI*=>b
zMVSAKya<{r7LT;t(3|Awk{j+lLWyoCvj?j+Q4~;;!c-DkF7R^ZCw9uu;V=hC*`Un%
z%b-ACK=T4CbTWog;Av;$o-ESCC!K0ccSP6phEpk0NowD`gKEm}rVP+wjn!^cG5gzP
z*pd;~!D&jGd(0|J1yg*~kHsuYg&rVZ)MRF)Fz*a)p$=5&zH&j*sO=%yDTIY=$x}F@
ztMW_rAk7*!E|krSyjR0DB!D{c^l$!ZoN%QZ=O@zPJ8~E5f^fU5E%B`L(348X%x81t
zfYV(TV9^LZB<Z^*yIUKr6~<dCEr21UK`|ya!+d^cb!jTzD_in6uig22KMeFozu4M^
z)T46Adq!^BJ}@f@ia1Z{jt#%R^EL#t2?<3&dGzWtp9|{9F<YR9+56)Lnq7m^x%9k#
zR&wczTrR%52|N-R-;a?bPwat9hF(g6LaCEfi!WI|B)o_cd-A=yoTsib@J`B7My-l`
z8F2Vi!5Ske%Q@Bx<0Kln754kf%zXycuYbtd*F0I9#2S!jKJy5Ca9}G*fmR3>ar4j%
zu8Tg-BuMBN7rs0DEfJM&-nDv97i-idd5&P6(gv}?*nZz;Q3<a<g&B%>pNn9hr6gRI
zu>NST^l<z`5z^r+R8{XZMW;w#uHAm^ieBYsXanQ<E+YMK8P-3Y<F0?0O9i%v602R4
zKz`qbx44#IfAY=SW}$A7&BPtVVsRiRxG_bDtalW7iP=T@AeHDH_Y*R_|A^o>Hx?i-
zGHRfJ5Y&*PL=Ak$KG1Ejm4fGII5y8cA&UCZNoJ=)xA&W`ipgn2+tQM>=rbE^h}*;`
zBNbhXeQj?P7r=K|SMyb}R9_S<??>9(gx66Cr+&4tD#l!S)4q+~$3lgw(zg_hd1rQ5
z$WyBR^8x1Gur`#sr5D@uP^GT;sGARdP{P+~jT$oo40mT%O*D#ep+CD-qKk*0;$n&A
zmu8BnyBR%(f&i%CEUrf;Of~^2Ed3sX7W<#NPs3s%Of71-FgwG-_+l<PUai$1n|~L+
zerGGo_TI$y33jflE&W`~1h3W+BIeP2pTKd2+_96-jd9qnYQ7Eac7QRqE>xNGa9a)v
zEm=u6seN@;(+ThDtP?=DeSH|Wmuy2|EKeVLeOERAF7!<aE#|)FaSpqF{(WiZ*hq`c
zAm|~PWvl)tL3F(6CM#^<-(WIl84Mmn@QI$>Y^ljoi0I2o^QSAMozm+K$A>=Goof;m
zZa9+PXe9gW)j9r!0K?nT-<Uj*Q5(;zs26=Wn`@$&tC71kdZ#*Z2BhDSI4qk6n)pXc
ze?oWN1SK7-Z3=S)0N>E%t)4QvmUnB#Z}88?%b)<j=ZfgNlg$`(*>N3V*`oGj%n$?g
zz!5O6WL^@}9(i+;-#N)I`^W%N@j4>E4lkB&-gJNTy-jn_CB#lg1u9$U0)fK)69Djt
zb&{vzPi6Y3X{TnUkaaE!)gZ}6t~ztka@t3>`LF1^vkCfAg!l8SfmMFLg<UuJTj(p6
zBA#^J3Z5`-$i9y}Srk&k;|}Q<<@_9-V<bT9G}vblX%`(o^l1>?qk&I<r+@70gDpO7
zyZeCQWq`gADGM&B+>;YbK0jgT5?<?aNW(tI#>ad{<9<@LlGMjGS$>}-h9gb*vp2)1
zTRa8X=X|L^pKChTplKQyxLJa}46x}Y)%MBd*w~}UqSy4!);?WN%*fQfcuU+2C+K;j
zWgrA>?A$>tR`#uRGZ5MLq0`H*R&4VwvzHox<>`5GBT?VJ7-sclZ;7C%YTjGGNZkBJ
zkO6%!?Q(C%!KV7}*CBvPjYrM51O;X$xYm1g;_*$n@2YJUGLxkGgvb1J?LE~O)tLv@
zV{%2Lj%gspY&n$kY^*|r@KgT?Rx!3U$|vmNTcQcorz*yC0^*-m1w)W^?z5T*v_y%Y
z|GtA4?|qco4Ao=XpE(yJm5&Gt$n`n!^!Yo7>kj?>3V<N41P|A1U^#R04{Ev5J!fx>
zC(r#U&UUoOgbi#dD#&~j981n7(`P>!i;PJ@m@A1VUQ_;JHq&cD7fmZq-(fd%aK5vw
zVCSj$lbR`xXRMOXFu78Bluaq~3s^k&A)NqF7@NY%@QF)i*Px(d{^W}I!di3OssI-V
z$UE4$!?hcGTY8d7B-GNlwyb3khBPf!X8*<<Wqz`T=N^u^n?is5+6|7<NMWgm+LrWB
zcO{8v`O@Hju34+&a;aD_<VaTP`_)Ia>n1e=CSPSHgfnNUntqWqby$1%U5pDCoul6z
zwg(8IRM%ZK74c*ln$gEoya6$SR|KZ&N2!ABM=92=c264GAKIRR0fCY2DLU^{&Hfq%
zDRUj4txTrjy1Tg*rrT&`?BR*yg))%Aij;B?_OGF*zlxFimyzCh@By~B;4PjB56i)M
zYUG3J&E3E(tS3Nc|FPWvkLZh^!^yVoFzb1|X!|8+$V#{`c?MS0rt4%AebsEvM`r@9
z&5fMV9lX@so(J<Bojw!t6>+{9#7dh5*Y;lU(M4KrXZvyL5-zlv535hF(D!saO=SZn
zqLG}b7?F*CRAIsAPeq^m-OgnN5W_nowJEr2`ZvzUUp{6DC?fCoDx(&fye#LsbjbX(
zoq+6X`_|MXpV0F3WnY}HX;e(Q-+PrloIfOXmW)=Lv~B`-FpE5U9;oy0+fLnlhSW;j
z+K&XIni^ta{7u%ql%O6`w&Zjw_NtM$p#WD0qu0q@aZUmI&jLN;H0lMV0mh`!`WD-6
zwoVgtn`>VePeW#^s)w%)6)wnxo<nMo3r=N7e2bO@!B0Hb&N^IlRlOsZwU0oq693l=
zuq1raaUp?c%wUC05~2M3FY-U$H*T>`GRdqbo)|g%kMdvUQeT;qlA0`GS{QJ>t|KQG
zyUt6K+~ElCQPaMe2`m+R#S|I$e!inAtWp4v@!-=Z@qx$x_}%t;G=KkPMm<e+g~Wxy
zi~UeQ$Sc+_<i<Tg?<MMfthE?$550fXKo*>sqm>*<!=9s)JkalRmB1&CLGNq7FDX=4
zJIe;5`bYFHek%Pn8DMsip7edZ|A8^L*pC+-6AFL2_|%z_8o#9QC<zOmc!a0yxmsxY
zY{ZYrUMCSWm&Mj0=As(^24T}qKb0P#B-HGToaBYU0plK)%dX7Bb5&8+oPdjxpz7ix
zDB$QMldorXgwD#=k71L`T##O3uTVouu5h^(D#pQa3W)5S9F@gM6)y}DsAies3kzHU
z>nxSd>d;n{sZijTR(b?#gp{3kn=lxrEyxrKqRW^*MG&0<xr>z=W@Mj<-+KWnwF5y}
zG7YhdPs&daSigeOE#kC(0l&=IaZ}|n{!V_0-AlT!9g<HQD#^jKH!j!yYvKR(_0?fb
zeSg4rW57_54Wt{Skq`x>K_rzDMU)sV(kabGD=pn1qKI@#hm<tZjPCBv_xk<4&->^5
zJbUhTxBHx(^Escc^GpqLDqCFrh#U&oMf~i2E<9C+L|^S2$I-J3d*2OJ=@nV{`dR<&
z1WnN{yYVYu5y#m2<lrNS7WUV%z&)u?#ab#f>YK9CE|Diw);kz5{Y{ry@e?q&`vO<?
zj+BAT@(Uu_@F8~2^iNjCK%@BPaoNJ&$UsyurHIpFhc?7S|JIA|rI6gzCOLCw7ninr
z{5pLNDOtw%lZ*ULotJ!DcAb}R6I!EDn_79^bZ|%!JZ13Qa`ZUVOg?12oIhuEAow&+
zk>JoqT-QIC$+l|IK=P|=Y5Rwr4`iO}DBs1Gq_Sondgke<lPF3i_87Ii#Ml(QV<wC}
zJBhTBNyqFAiLtAh4xiQ?jVg`l&OnKc@=2kIPfn=q(<X&<RT!!kOyf|um6e81c*PGg
z5FL$eorc}Tg;~i$DuNm*5>ID=GIowP=^HX|z;(GNxm|LnfR47{+%z%T=^{7ITkH6&
zSzbZ(Iq|n(Y!;&UTH@6#d2|2iFh*7VbN)S}vQzcT_Hu23bjA9OY&DI2ja`9}V0>#_
z`=|s;O7RTSUUtixX;@LJCc;pt(M))vh0P;KkBenF-NnUa`jO|aEBzgL={f=oj}1k$
zGm#;i4Pvb_(gQHVliqfJWqZT~lb+qsb>?fh+=uk^__cFgPsswl_P=8|Uo`8x46*OF
z+=SkA$qcSsk~LC%x($iGEahG|@am2l1tsC<YAw9=u-aUuj-*26W6sZYF~1DY-k~Nx
z6-rd~>M@CuB&UIiN3O(z-todi?d6d7G(7@-A8ux|BGf+qZE%-My@i3zPk4cjli?7{
zhaO%Fk2X>V-Ftdy=x|YklhPS_G{AU=v}mmI<AA3IIDR&d@YIYo$e=v&C0Yu$>=BR7
zrGKAynAqxG;u<TVvjtxDwwm=_^DTK+8CIS}&C4@O|B2XUq0k0VU+ny@aC(vCOJ^D}
z#O(R*t%0c3s$8cJUA7QUE$vY}4uPQu06OD1JUsmQm&3xe|00p)WY9Ny&3pz4QzK-5
z*S}6?82$OiKW#!N*e8ClsLjGUh#$HAz8bfqHv99i)c8}25YYf2hRERWegA?}vRE7?
zBqObNbMUCC8Z!e$)I>1F)6H<tX82>x^sk;=C1HBA?Zxf$4HQ4_Qhoe{`Ay7ljYM+;
z`7~#I?az-#Ku-)N``J&>ZPj(__`YApOzWJ)PeDp!<>J(D>V3B{?_0}<OfRo`GW~2_
z%Bz->s6716-j0j|sYPffF0xD8uO{bk@~iigT!@+1V~2l#&EJLcw*@tSu5EfeTBfdW
zlO{5w%mVy!`Q{(V==QW8{d^y^mK<Gmv($<%6hkAc6eokvT9S1BjGHn%ZEUg)(f$<P
zQfQ<%N}}7;m|H6&#X(fFq1`Yd|Fb3ccJ^nNOJ(arM;{%8YW>A_6$^`V(Kf<4a!)qH
z-9F!$YF1Ephq+@YBfZFpRW38}q4(X*4=G@pKT{v)YBKz<=VLTc=g8mQncUyL;=Ol;
z90tVjU%!vO33geLbpSP_J4`E^WQ@`QsQdO$z?6x9>c~HF@VTZ-|5f@Y^e)%5t@2`X
z6D<h^jx~#Xp4oE95}oQbSBeR7`R4)o_f`|-LRNgGSgLMDb@#T^50Ve>xvOnJH1MhT
zV7U~dIw$MGHll|F+BKGib@mi^YU&I2hPoPAR*d_lP$BW_E*1mhHE7SGEt3|lB16G#
zaux<az>p<M8Cl}@%v$nO8sDg=R=__CU@O|NuuUi&kQo&4l=(tYuf2FZvV7M7FE}3J
z=MlA0_SadcOG|7>D6Vb9n2f0u5}4}pZY{o^sEx1G{pk?hcb7gS=xU3gC^DOh_&xD$
za9$hGJag4O<@eh)HBWuCG!ylC-4CvLp8_~lzu_Ie@*QqChTOxdavNK;{PA9~dL?v{
z-NH-7A&dS<`heDZ9+APCE=N!m>e+m*wJ#Ls`E<Gs+2jiGxPSemR=!mO`OiE_AbR%j
z`NGm9Lx^+-Pf<13eJvBAMxnnI(_T9%4E|2}st9MTAK8fDH^NrCVGK2Q0vZ(<?OSL{
zu_=`|x{f_N^>=nP0%TfCu*>+8wNpba6!OZs5k^giBPY0Xo~F0oET8JUiBldDP-1X>
zO&xZ)Gx^%7m5+>oBx6{OR#sJnVWk`v{!7(^z@PSqnkJ}b`Xi(&e%VV)K}7a`YH%Mf
z{@h!Tmjc=@o8w<afJuZjc(3~lcj_<dM_#AqtAO6Mx?jE38wqDQ!ZwsQpMF|{!FEzd
zJF>>Fp~CJ!@6ORmy7&27+~}-($Lqreg=zS#Jpd@w1*q%CH{ix4y*Z{kQghL8v9cti
zAo#4k$~UZxd=S+B=p0rP(E?8lX$edCzp?;H6H+|-r7HU&Pv9*{JW@v3cO2bc!|eEL
zXk+J>Z*nHA(C4WTPI>;)EaTq#r`6i+U?k}B7qZSmEW)3c?Y+mrw<kAuN4z^Ct$XoL
zLvTy9YgeMXCKw*Cqzsaz({F~M@Pr)?ic$I+W8ZXV3G3>g!Z~K{bF*lA14b5AMI58t
z+VGCAq(PC6$F~p}vPE&fGLPMvKtjm2S|~ZcOG6`_@aFx35t#2La@ZeL=zA)@L>VoV
zVWBYHOI|a8=aBrX6apuO`%xuGS&n$6^ZjiV^{0R?wTF?5$-$jSbGCYD+9d9dTqjwY
zZ9bp9sUV4YKm!1W8i<Ffw44IOK%+55s84)i1SHa}z)#5RuB-kv^ug_K2=!zx<0#$S
z)Npld>(#o|Jxya-p<oU6zwXQ?3X)GR-u$q(VohE}T{hfH=g!FbVgvwfmEm>t8E3Cb
zg6h?K#sMrmRF($CeQR|uJRCoWY|CT{3@7X_$cNA|!(L&0_X=l3wlPJ=qxv3OC~{8<
zoGUH9fF|_{(!Fu>ROkKVr=$n!MZ3)v?7^q5cFVMu-d)Qgjhry2qC@0XFXSaf3oC{h
z_L}UUvc@W{BKeWT0e7F%l18PTxlc7F6u)@%vr3mUN#Lgt2tbH@2HC~9BV-=K%U@=H
zw&*KYauQ;6<o|iR1^atjRFr<Kb*j!C^Z%fBsRE}xZ6h*Xp*V<A{ceLo5&ssj48jB8
z^(&jD->3B-3J`Orotk2Kpx^i-c_VQ%Y?itoY?$(dNssT;sb}}=&!kWCNpS@#lVr&2
zR}+7i$n;DZB@vky-C81X7_0a=$nJkyyrlmit^NF{Js|9jQU<jV5JV@mPsC5xSJgEe
zi^IFU<|~^xr9WbY)_^Y4!M;;CJ54rxh_V?$=)+?egx^@jBN#J`C0Zc1Q<LX?Mv?{9
z>*wkG3-4G9LiiI5bY@~jf4noI&vcz{uAm`EfXg2j+&RF1@<&$MCM)(wb}KIsq0+&$
zNT&SH1iTYrc}U>GK*qWVJNQ@|*3}%Mc)+oFv$^g-NA~1p3jol%jI9PGo-uQeKQYKz
zvy%SiPv$%IrpVRkc+-Qp<fI_(KZ+PI!ukU^I}hxdH-sxxU%M@wK5F&Ubz3@$hUjL9
z^n_Kp^DjFS+KAm;a*{E?&!Ob{V=f;UU44KX^@dL7N`&kWE8h2|08jjKP`uGa{Iep5
zD-JnpCX@fgz%@vie$TOEJdLi_`^NRWbF0F-07K4d-?=gKKOReB+~K(&A;9X|n+wOH
z(`bhT&MXvFDWrcZphO*m;x48jM24Py=&4_k^VoTh`5GV{+My5s3b~XVj&SigXdA(P
zF>xAxIU2#*90n-Tt9MTcnSU4TlQ=nC$vQm_4rKTCI+PkYVsHg8Scub)x9JBA%c(dI
zRupI6-eW>Qp7zB*xyrl86okyUpROpU$HpXfvDoqbodI`w7?Cd*F4R9{pkz;G%yzER
zA@%9gm#xxZiDY$6O;uPEAEYSK-CLtnA%O^0`n7u;$Rd=!onVvlXk)s|ygOmQs!->?
zmMmNj`|;uz0tFTIZ`8M)91~qCw`|XMsN|Vl|7jkIW5<V%1GCzz@_5BgCVvP&-P=ME
z_o>?I)vRih`ZRHF61VY~1w#!fe-p>px3R*s&QPKso%1Gwi#7f8&bXA(ah@csi?<ey
ztyJQt%B|-?zEeyIsSlSRzk%0?*D?_VX=uC@003+2EXik4?{F8cLZ7mDY%w}uL1#)N
z89Z;et%}CBw$&F@u&ndLasfVAXBf@>ihHmT0Xspt-GEVBt5n3&1e&&>{`{OJ@r)jQ
zL6%MgfED`gM;<Frk9IVu<8<u~Z-@9u*c&Ca+5o)3!lNbZAMcZZ9p&?pLjg=r?B$GA
z;>9OS*<y_@bf`n+e+EzQnwx^|_c~o!)m$~ZQ+w;K&qeeZHHdzuF~>dO;*0t3o%XM|
zavb>n!rQ%0z(^-NLZF`*tp}rY`Oz(Ne7tjSk)okAEE#?8LYZGxY&juoyVxU#&_(oW
zFC8@0Q}b6tWE&54DiEk^F?gke+R1a<Y4;X<0TKwH>pXn40{kNE#yrZ6LB)xuZLR5u
zc}|pBq_>4FX-M8l5ponp)VDTXql$0C#h_z%yJ;ZrpxsvkUMpzv3VxX3acR8I&bw9D
zlgT~QBwwZQpI)<&C&zafia!>|{31$GD;dSh!!05Y{Q!95X8536NEL0%C@;h{4;Jdo
zKwi6Dv(xmOatI=w<iO{k43?&GVrF06L4^pCCp;Z;UKVJlJ*clUEaBkAt3}wg-=}^H
zoZc4fBhsAQn+ScP;@;67N+p;yTDYPDkyN3w<RMT#_BlN4-?bS0_G-zODMGZVW+gnR
z@oMkmfaWXEg31V~Xj$3G>pWf2pn5p`;1sDq>wY8QNG~3kunqn7<GqQ&v=w@rkWn&e
zz_UN~!?f4(SEZS?BJDi#e}7-zc8%~~30`}QPTwMtS^1RQ7UMk7`)E}6VBD$NO$gIq
z6ZWnfY(uC|>^rmyP62d}j7H)~6CfODF^E4OkPM%uvIaC>fGDUy=e;U%GI`SY1DFCv
z3>b79N+AJ=vywhE{{^w&PTJJp7W0bo@|Rm&>a3i}!?z(>^n%v}UsNiUt(`YYM|WoE
zSYC7d<CVqIgU?~po!xnnDeu1vC%d}`ra#^)8;$P#oiz{%W1=X<Ub3?ys26*7=v9Xf
zA7Vmlc2YeaV0`e`Qn!!9X~1SRdCa>3PXeu+j==6bBn`|%-bszOAs;6D$9@9hjkMkF
zjy>!V`<oR`p|ftf9$V#p9`2g)zGrt88xE}gSeSMto<;8X+&HZy6G^x~o|d!}58jSY
z8_EfjE;?Jib7|p=r!c(3W;Z9MP%$l?VGP+;1*0!KRbvycx9@g{-bTU=WzmjG#D@K?
z3<c+l*<h9K#n|!jagvoeAB?$`9y$w?7O$O>$C)0xZVBR@Hth8CXj6o0>z~z(q@~T;
zqSByXU6DRj>Wo+k<Fd~;O^KO9*+!Na*VN-lde$rRu>ka-&B2R%1=}q1#akLPIW+%(
zPd_fI0w|!1Rc9KM-E9C`lZf4MfPEZG(m|(2E=DN^C;zZn+qK;{>Fv>8nGYeAS!=P*
zAO7Hv9*TpYXb@_o7Vd_tO2qf4w-<vZ-pML-4sYlFG$X8A&)wKzot5EtDceZ{%^bvm
zLG@r$Wv}Q4&&@s)F@ACU{?8dE1aqH|h@AXaPjPB0XU8OId4={#TXp40{5Ep8B<ILF
z%RDu!b%t9ya+c9*DirXC#ELiX_ZL5gY%0n5xhFC5nBV`QnA5WmwjmAwwflfnZ+9eY
z!KUl_SG)75p77V7+@}P{bX<9#pZGy=0Eh8<A&}_#wZLzBs#>oRSM+cVfh@tK=BwBc
z*l#(Y0k5+jY*V?uys)nJmV-52S7I)wuv}xmxCI{|zajSEP5LMqT74G>;J55eN%eau
z0c$HTvzXGQL5&?e_AfW|4&?Y%W*}AlJ^%jo<8@!kRXSP5K_U0>If1P9a7NhTL|5S9
z#`T_3P5q52a~j?<Q}2i>QzL5jgu$nF{}4>uIgumW9Kh}!>nuPt79ZniY>h7|%4U3q
zW<t-j%G4f`d5GIfob#`*gLs=oEOj}9(8GtW#%wveSfkA<%jRTt>WK#86WvD6kLH)K
z+X>F93JUEgkBDZ_GJbqy1iEs^Zn6rNrFdT#58%p*_@fsl(4iPyjL+fKPP%04Z3r>h
zI6H!C@PNa!joD;MP4{n%=h(g>$cm?;)K<$8Wb4txY6YrFqeLfMLzg4X{9p#1mk3z1
z5|wOSG0aLRSyOM0vuXNtf)M9`Y-1$EVA2JvC(M6vpFPFynSanEe9UI)2PxAoCx#)M
zlG0-NK7j}zjg&EI<6SxT*IXCC>3Vgi?P4Rz?lhd?%i~B-e~@`xGI5Apr#Pzla2Tma
z)Z<YIG@Dh986k~r1j0Ms?RtJX+x;LvvSR_ID>d(XOzZy^uwuq;983Ubek|Bhw+6+Q
zN`e8;jSmgj<M8`^sj4>%Y9E!V%8OLN$ho)kS?GY+Q^GsI*zr&mPV(T7!NlF$02tLO
zI7%hK-dS(<Pn}jECGhCkte0yvNLUqHN%NTSVa7XJwRF8W+(Y%Iu;p+vm6ox@%XBE6
zs&BqrB{k$ae!Zzi`Vh+YOvLjJ<_>I^nDI&Pqxu)1ic47~G!P(NFJwrJrzbab(3W4>
zo8qyNFUUUW&~^VRpRV(y>u;uF&eRVGm95@PwaXciF)vdvAJfHE;8|Wdn7LmmY28#i
z8g$p}QP?L!Je$ka1q{AB)2s%q?kPjOtKJ7@4&$W7-||VHra;H*sL`mGbZ_=u4C-HG
zf}NqNKtV5^t$$#|Wh}&TrB`A$x_m00#Ez@%Pii1oCq$D(Nr{i5%`rW?UpaNg0i)ku
z%XDibe{&rt_vwqIKRR@&f#D>g?|ItEZ~?@XvT=0dC46zWZy2ZOcP?hkqqJO85D)${
z!EODd!KZCA7$PBV6pY?<(0~$_2AxH`=VP^r^30uPkD#SCGh++_`?ZymveVjAQ!Ws6
zNVV5xn&Uk4GciXO6Tl{b$&~l=mBf=~TNa7+I-+C*N9-^M*RyZ|hRVF(uh9k^^uq3&
z%5~7Qj~|4Z&L3v-(=_@hXaBu$W!k;Lik)!$$vvO;@wpIRACQoEB1klNe-_QGPNdlK
zQy7B#EcABaoT<A%Rl*$7+-CmBh_H(KOJbtw{3qV{zLwYw2qyX6w~vhPJ>?w?@faxh
z-BZampbz+odf|i5KS=bAH(V$TOM2z5r>wsK3t)%!8v5Rpucvq&|LlD>dh+&S^9n2b
zHm^g;*Z#==oPIc6^*PC1mj2IV=6Hm*9@r&=J+BQb{+ZHRZd__PiB3klXNe-k3R-dJ
zhw4<ZC8z$w?_h*dH!~D!q?alH?MGW-H>`>&49HmJ6bIAhfxDEyEZom`X;>p8`^oEr
z*!GlpI@r7*zLT#N@**o(I^Joh>DGh{rN=F2WirT;a<2=>uxRT!|8WF(Z<;CUPcIlf
zZ#Ck9fmt6FQgyL+is>F9Lhz4W4=p_M#@fFIQCo*57k>&x&;9AUoF?V-x)`0lv+G((
z4f##%Gu3CR7U6Yetljk<GRKji#?SJX{NB(zRm<^hl3|VbdfjxTVB^+sLjS7$PVVi$
z-dqZe=8u=LKBf_P{-Sq}ASA_O@~E}5C3D_-6~XP$)#qRKy9{CPIS-?MU`|N>Kv@kl
z?;(2vwY61ALN_iv*TRm<EMcDy72dUeqNjqFx*hm<KpL@9H)nfTne}ME8%J7klO~wt
zu-{^`yrL6aLCP3k3ry~iB`%vE8ghbv92GG+9<Ll5&Qub?v>yRE#A9|Q5lnYo&>1PN
zl&r27e^ET?GNH%If|$RI+Nduh_;r|0D#QbK<XQS!TYZitTW4MEEAe4tBH(80l@{Ba
zu6vv6GLa~}`o3d2^BZA^DPmbEz2TXOY#x$7=pyaUeEs)my!*QPERhw5Lk#H$lQRMv
z5ayEi&UuF4lJ<ry6N<?*i~-3&s&+)C{wmF7pO5Ny{_E4?``G*)N3oml8wb4oa4pAo
z;(BY`se3;g=#nk6aAXS|X9HD$+?iw`7+G$ESlo^OEBjm?d*ue{o`K8(3E!J2G4qa%
zA{4V#$f^Kp$Vb5<cX_<#&B>}4IH%n!k{|khCGXD=?=nCXEAJnd4+hZqcle(_Ab2Ki
z;`Q<MH4;;nKl<Rn1Mpq8JFX)A^0<jDw@~oCelsV|(-iM#^bz%u7UM4Eo<X6Y;FN4i
zb+*B?bEHW{z?@!1{+@X{g+i>KA5bPt2Q?2#h@52JqE2+3QUb215^ezWR0RP7`l>Yr
zBT`UzT*8uE=)RF&^y<6s8abYaz0pj*BhGp!f90M}U1<c~y0&$oB(Cs?#0S2l!IjR-
z<Bo@#uQrZr3|F48okzz{+S#@&zwzz=v`e8X=G*OP`nwR9z(Gqwz|DBjc>Uc6Rj8n}
zCyncsPRAA2^(zuUPvWRCeYjze-5ABW)HR`JLjIf}q^z*2nqTWlo~eEkC3?G21fb(S
zPEavu{skE@rEaN@m3jAkB^15fR`DeLOTd`0qT{Yv=iy^XePMArtc@@whML#%`W;Ow
z(x0n18iqNBQ3)LJNJqwZo7g!{yC3XTid>0}vQ&Pj^=N2LDc`AK^pZrckS<Xqtr7^L
zey!(%IA?k=#x^1nM%BZDYmb+y^m_$)(ZO@Hb`;cOB-J1joaHE=;XT@cR~+#G9<{cI
zH<w~;ZvQOT?L#r|gMTTJ)0dw&>JL4vos?IdC(_Cpn7+}amiP6zktRl{(>Y#^hpe3H
z34Gv6oAhe{2I4BZ2mtqPLdt@A_e~0I1Br?!QCAL7NrDW1_T#ymA;%!=k6v0-;hckQ
z@1*q-6-BLJ%$>A7@1&<5;lV=>goYySI3*+b5vvRA1}8iBPht@MT;=yeu_Pw`ojd`I
zK@CMl=cI}BFBJCsT7Vl)-dbBjtmg{-sc7)M>2~O<k|1iV_Uw~^&3MV0{7A1numd)?
znl*|J>|;5O(peweAO;to^Wgpp;k`#tfozH%c=vE6Wfo*vGVyS5blKjsf&^<5qAzsB
zr&QlTdq9a<e<_(QJ1ES$p45-!AHK2Y`N$Yl-%pZASeAo7z{Vm}ci*pw*nIM&uHzzx
z@_Y>1N!wbvvYzq12R)Z$$$K+AtZjR)E}E~3H^_?eq<QyZGm8+~9-}iffL=dZtfss_
zskHxe`I7+j=%|duqs3D^yc+Bak>WF{rA>pWsy$zx+t?mY;3(~{;qT$CV`TZ3=JG61
zDdI;~{$qmTz@_P5*q8*_Yj4E!a-H0jc$^R_xPXw*vVmjIH$L!!oz97&#ytoEXg`s8
z3^hEWf<2_~WjWR`T}~D3*OO-~Nn<Nc?_;!E$WN+iYy1P!-jIRvP`fG>=1qMl)Jwz=
zlBun6IX@3~I3>RSk&&f7t*tir2abUAP+~%`%fer(=QJ_pJypan8&7rKnwHs&isy;8
z5VjPVc8#@+<ZIV|Snrr1h3NhXwH26E$nmL@0Xe_+rZEvTDaLY%5_8r<kHl-b*A_G9
z?zekBxArE~H?GcV)v?4g>hXjyblKI)WFvk$rYNI7<lWAiE0eWCB<9qey0eUAF|Ii+
zL+bf*)<B-zkm1Cz6UVktxk61B;g6Mgu#zhw=tKR%{9kYi*9R5W+Sk09p5h4~&=w5q
z*KbFwd;@=r6x$jisOf*<n^IqQR1d<--e<+i%(lIK^dW!R7y<;An$}unjO8A=os*>U
zJ!tn0lg`hZkN6w@HpQO|Ue~6wlOk7ik5@V8E!-h`w=G+Nz#*p~pj#<;Lyh?Wd5{=9
zT2`+^z(oI!r802)V(m*=jW1HLAMb6Y&DBwbb+YYv#a8EmM1&gK>IrZuVRzpKLh{Kg
zPewZW{qXmo_p3%M2nn~ymFr^~tabfMIJ-OztZgG)3Q9J<15)Hmhu_%?Hww+l{%LYj
z4u{b&%Tbcs4~*6c%CQN;LS6#rzk%B5B%2=J459V8cKc&tGm7t1u^|Z<T!-I8$v#!q
z*f?T?8tVhwxGpDgWQF~sjDe6nqnT#Ld(=cZf6M&rdAe+2^XvMd18BRagtKF1Cqo~-
zmfF8e+%~rzE{u~fpk+a6x)AFfw9fg*64V7Q+wz1l_Fbj23R1C2d<$<)xHnd!choi_
z%mtBpykm+=NvQ?A4sx)%xWQDfVHLIw>$@Nfc7AP7^dz5Oq%7zxJO?BRXOVZx#CjPs
zlpe1S8WT3xI!%Cu2K73WS)Rz|j{|rZbLqpEc+*6?+_trQ4Go7hPSbV(Kw<Lli~y;F
z#2H-mPv<;S#xMY|)K&4~Hw=(1Y`U(4-}D*uM8k|XAQ`9kPa{G`%4@X`z7$ySiCrJ9
zPQ0bd>(skl^HFU3*TR@`z_j#NUexkIupK!LIhvd;V>@q-+IYr9>bg%bpFZEY6V5+U
zE$OnD7G@FcSNX2kRyiGSR|3rd8*0cNJiR{N!X_KVq`O{eOWv8M*|?ycX?K%qod=2d
z=T}}4S8O)%y*%Aj*0jk>I)LK-vz!4>Yke7;SV^^-5Df(rQDC-Ep!%e*2;qHt&h+b?
z0{GlR#VGD4ySF%2N9+y<?UU~O%IWIB$4AuC>G=IVxTi0HaN-cJ1P;Ic2|}k{IFOTs
z5lnl|b?kVBIwicM&A3tV-+8~aLqnz0;A3T#_xRclhi`l@_Na5AUE>0}s>y=7L!a2~
zQ8*`Y;KQ5Q?z-Qp(&VHToxq!&QyE)7bVZN@xRT1BOpm%|vMwa^EV4D+vfCdPYpIFE
zR@UUZyy7&%t46=Iy#jL+VjY1?yv7wk<aQ$`p_X}u@Fh)aGNzXrlJOp|B7M*5pxm8>
zAGf~`xMD$juU=d-q{UzDJr~w1omCQNpSIiKMZ36EA#-HJ?s#vPyd@F3cStDbI^XF9
z+${gS1yYtOrN1jFVtEu@eK+f5$1;`^HUEsFc6$e+o6ON*@nrCa$OeXDf}`|3NMzan
zEShj@bG1G39YSVjYQc1Q|1iRoGTRcLg{!{2;&m0LX*ke`UNBVHehef)PdVRRv!dHG
zcXLJ#1@re-Yo^P+9vCite7sG!aAm`EC4yd9-Y~w~We%cN-;43bTZumIU330Kx9Xas
zwd=_9P9Q}S6zeYg6&mfnJqn*)-fBV@ycCH&Z}HuXo^lvVxm}GjUReMutUn&2cU2dU
zub;m)eI;Oo_loiR<y_@WO8IqHRm?U&6p-%6x4q`Wcf?lO8XO1fv+Le0N#XSudD)Fp
zfxQ%wUgwEeUJr$4`2hX-z2XMwK+kaW5@PY|TF6v-i#S2sis9bq99FIIC9Z4iLgv~_
zbudh-(pa>d&q}!9ed)AFM2TgI*mw_;AjWEPFF%>!QhT^`TL4=AMNRT(T)470&wNkB
zaqySeaJ+ZSA0F;?RTg#bpC6AbiKq8Qf=GmT4{wns>wZ1M0{WE?eQ7mD_ti)?@AW;|
z`)Su!Z{SkRGUd;6Uw6PM6}ieGH7@vSM2;aR?S2I^n2W~ak&+JX7d3?8bVjFaL=^T#
z;EFC}cLDV<*kbEt<sl-6beaPh@ys^QoQaJ^OU;<>+f=ovs49VfR}1_hiy5kwS)TEe
zpJ0F(Wvx%W=hi>D7=nTU{FJBpbNM$V*<q^ix@PX{HU5|zF-Gq%N&N+kLH=hsI25hr
z*e_C1Y*`Chjbovc`Z^Odm;?-&k|Uw40Q9A6UqWH<%*+Ed3QEju@G6U+SlejD09B^;
znCnS2+L@_QDwR|xqiz5V_v?Q+fp+`5ux>G6qerJFgM6kv+g)d{>7-CiCV=x30EIX0
zX^_ee-f4}-X!RfXkZOb)*_M=T))QGi`_8oQs6AWEUJn_WQqj*)K!WHm#cVfz{ol1T
zCqWPmM5fGP)e9=pobI;9Z$_ywRkS<qf1Ga|+-c}Y-9LtXU})+}74}=>1CDb)ZowUH
zgOi~~JnEiNZ?DZFo!X;QU)Y%K{2Ar3i=Vo!SP4Q4mqng2?YpGMpDRMo-zgsokpLi)
zxAx-CIqM;thi8iW34WuKWCDW{TZcyTbZ?>e!_mRlK^Bz-;RUmB>7B;adf|oc9f{Gk
z-aJ?fVV-(iJN%wCXbpA|M7a|KQ%mxVi>Ikk)&{kn)1k7qj_VRL?LQ2>A0<>mE_F-V
zp<!!b&GpiL{i)F~3A8Py@4Mk4My>CJf~T9w48=-P6||>Gi8-(3{sYqB;LUCvgTSC`
zB@HAxqKRRq=%}~eAJ;!EfU59xpUy62u<26ZU9WHi9^;%bB%0n7@yz{+zTEHtnJ0l{
zZ_8>uoruRwq3Bd%Y+wTBC593sa(DGA$HwzBx|a$=&5ks^erU*)IXdoof~=kYxfA$;
z9x_-1cvDN4<y0R1Y!HxWW*l1D;&_UfA;|A*Tb=H*JDyH=u71J|ZYnpu05o634*S--
zw=w@<W%Cy@=Sn2&rrdDDRe@RL7I(~x`bPnS_2&exmil=4g}4hBjyI&>86p}D=QQ-B
zK)wRbI-#tB8=;(V<~$)wr7d@_`|i!ZvFU0C=@8*~0#GG837-$G6$OyUu-*YW`3!kH
zbp~9<1APDEjxA1{Uv!OI39buDMXL_{Z9swb?BMn!y5IgoD4saM7Y|4%%Sq$3s*A}Z
zW)>W7a&R!W{a$W1NK5y~4i^jRmz674kAgi<5s~!AYlAc8xJPaF&^C>D*8){b)2Ize
z(~*gWnFF45vnPPl9r=^A&`BKOPsrEQCMGB{7J1Yiipx>cek_ErH`^r}nUg)pu6C7f
zAa<~lUMu3g@k&Oqd2Pm<u<@T!8Dy}<Sh`%UernC%0LgG{(oHU#m{aoghbAz<79`O-
zzvqUC?dH6u_aE|toG06<6TZrObW3g0WjrcsW9JGk;8T?Jw@yBzt;;;T8;rOACg3Lb
zKYI*Xbm+SBQ5#frz(>FDetC{fprPpq@n)sG8iuMI<4X|qM?fy&{tTGIh$z0Ryb-BP
zoM5U|H}XEX)5Rx!%WsGFsm?d>7WRZP-bYY?lCak2s$$Pa-1A~gkbnl^FEqZu5_AO|
z<NRV#yxwzbZ(0uPUpGNm!w&uc48+aUh<qGlBD6Pf#c0HD39tfg05FHkuMeBE;HQIa
zm`l6jbw}I1h?pZgq0>K(;3|{rZsb~g_vJS4!x>#v`u0@RNn=`ZGr{Q0fsC0a^{Nev
zde)5DIRgEYREQZK^29myx)4yLoS~OtLlit9m+Ia%q6z6OX*wG=Z#r$g#Saa<|NCt_
zpW;{~BfQxcAz1WjH3lh?wA<Fz^Tem24m!&>@8}eZ<p1IGgyO0ube_Ps$OHs0YcNu=
zw+YL40Yxh$6*{FObK<fEgH$h%!u`W;`7)M;fn<WF#2lmW1Q)!BekiF_rtYkUN+0P7
z_@GYD_cHy%z#d6u0U(PEr^EyZxLueUu<Uvbm%6Yi@(LZ$VBYge%4_e&tMo%;5HOHN
zs05%2L1yl@UuL$fX)vuH4G~d0#DSHe4d24T;l7ObWZ78r<$`@WXZIF=o~{czERuIb
zeeB_H7wT$Z&9-N>C=~5k+={<KEl$)(DK>J?@4o{wTzb@w?mV4YsZLs$%ZNMIyhR}$
z=K-#-qyPO`EIT+Vytz5D9R|Mp_AdleAz(5He|m*(G<eT=C}M;M-orB6k%Is3;S%8q
z|E4OlYJacy-xG-2NF;*6%Iai$$~jZPWqXn;F)FyRB(BafN(ATM^8gYlh(BNFbK^}-
zL-V_(MeYS{F?d2vb=jXkf5s5zRl2+)X$cAM#)H5wrS0gwPYpvuLx%{|dE*!M_UuE$
z!|~Q?T()xC0J7pYnk|Ij9P*&eMhvWn-~mG>a$a_J{E9x;XDb{2-Y68KNP1dYLCwhM
zXuQa_MUG@zN=lxXa6A4T+Qx};%MH1Ac6O1LSMzTgr~o#0cB*$mgy7Id0=E<*{w>O%
z3=cU_pjVS5*PyZcTULsnUs`@LET^xxmjoardMVtJl9<Q@yl(LJvLgca&ph)8ILJeq
z%AG{%xH-nHDz@Kbs-u9_ogL!8#B`$m_4SekM`4y%%S9IO<bb(3Lm@T@aJ)0qmz9-O
z>$svC*G{s<WF`%QHpJqxRZ3L!4g2M?c0JOY<=n7|fttFxyPJ?IN_bwMhumUoqC8l%
zVE85bu_4%gjnwUas_H=qUp?i8Tj$b)WKV$q(W95<%*xp$IF8Gi>KsJ^glC?T{I|^f
z4c*;E`3b`CnBu$kT3Y-=8ShE|+gI6m<n<B*>$QDzdXE_a@OSNDz*DJMW!#-sqo`))
z9g#Z_=^t9RrY|lokg6LS8@JD2@x~_;N(&|90*eP{)lQ?%8!$5$7xHp%en5+hEL^m(
z{P&tuQkV|^&bU?F#rgT^?%Zz>66S(QC_-6oN0j`X7D`a<__x~IiblY<Z{Jp1NY&&F
z3>coPsl5<N&Tccf3~S!0wqZP%pX@0{o;g%KtFR^o5>AR%%}d7>Nl`klUL@^YJpVFZ
z%^MWvD?pI8Vt-n_Q7tYe7NW(~YM6xzry-Jcz1tMZD3|Eh)Fd-qWQfzS+kD%I&<ODl
zbNgkepEr840BwND1FFuI1Czj$g%&BE%djuDu<5TlaZh8F10EzkOFa#GI&&bdoUM7_
zU7+~`hONNC5G2D!ovYM28trL?Mw<gOM(%i|wSzE?RW5eBYuJquz|GeDxDeVjG<3I6
zuO^O0fv`Kvf8ujT8kwBsXn|hpwjT)bKAD$8{r!=mpj|Y?ZE&XiNKi0Kcl2d|tPI&L
zTxv>6Mr=iXVj|U}hY#hlHI;PYvbkbEX)1H}9cURCq&w(o3kYTpS}eJxl~}(8e@xH;
zFSPAUT_rIi6ET&v?HX%NZy-g%isKd5XZ5zj>qR`^PnulGJzPMXrm{tTzLRJ^;Ty}_
z(o<a(*4Z(d;R_ryCw>Rd*#1r8jio8J!pQ7Ewfor{%h*3AIl8<*lg{9OOUHQx1R>ze
zerRz;?OeqCF;V2q1|oMZhDJvfKzNDwv0y7KuFj7bnHLY-M32Xee}Euc?r_%Qx_iq6
zE7a57-Jkg+Xtkhb`%4IKHM}k-<Uei_$H2_ox#S=x_dE-NgM*{R6&quX#UO-|fUXiP
zq2!kvBl%fI#W|l^zdUGOt1j3VS4%ueD&RJIrOauH6{3<Uq{8El{=<<SQ^L~)MS*Z>
z?T%Xz=}e*V%P6nV9W4jJ5ysb5!ea~tKY?@SXI*m8$=2k|)#Az(@|F3a^{H%xD^?jP
z7_A9@l9S72JN{bEHgm%Hj|AdOH6+uu<4kee8c(i7U+_yv(1rRR9vv~I_YKa@GRA0%
zI4r_}$c^eUn<-K8y%thH`pJ_g4+R8dDD+(o)xt>r<DjqQ-?Kntx|ZG<8Tra~LHq{L
zAU{Y-Ow{<dHYd^N-L&2KIXwi7qt1WqX_Hb@UxqKJTkkU>qnMe&O$IL_s;FDE?F91p
zF`Bur;vY0jmya|Rx&xi|v<t*Bnw>^w+_B1>*)dP$I#0a1>}gdvO*yl{ToGRoKdK6T
zT47%SzCxA}Oi-5DYF+_BB+?@PkF(ysmjQ&a-M-_*>+hQvT(P-=xXPR?=RocOP1%T@
z-l&}Dy>@#5=%x+5dw6nk;MgajP+|QdTayXU<+$ss0R=oK?b`>guPz2G1Erghpk(7A
zojZd7)Cw@)DoFhfeiy`Rrce-SrZ4TPigvrU#kx!to9Z-Gnd+G6kXIF`nAxCrE*3@Z
zaAtQL)Xt8j;4lJTV{R=oAvE}v2D(}Kn^okR#3A{QEc3^^Z}-0`-Ke(03ME^<u>_qP
zv6^}1W-}+_B|NI%2SLVbq7#RQcBmo*nn3RhM|_tSZ+!ViwV6H|+*vSNIaBjTQKFaz
zSM0CXJeoquOCZd_ZE7-Uo4C3IYdGF%Djy><(=k`0gpTJZB4CDzR5Uapbwc9ew6IH^
z25$-Ud84HKB6|I=RH-auat~fKv~=8xpQ@h{(s$5tAP%BZlcTfawwpc*+RdE`=IhJJ
zEmJg5y)3fEcA7kh+o-Ns5D`L-Cq8Q{vBu^Sv2$Q6Vn^N}HBkK2Sb-~b<tS5~uQ3lA
zCJYHf1?bvl{n;f)p2-Y;{bN*2Svyjo7vd?Ywp`!a+xuB33ZiBF&_C-z5~HfbJPz<5
zdaZJS^2?9IGl-Z0)N#l8JX4gj$MOd7m*dtU!&K*U<L4FDrDwz2BC2+Tz;e<^QrFUh
z=7x##k?!s$J~y*vhpO)S9m|=MM?ixX*V=o*;WCG+Nv*ar({;Z1u4M`-^!mD==8wCH
zX0Isk-ogcX_HK8db$36T(rjtFFx~#iL`&`VJCLZ!o;Th&L-g+e_y2I7EtDfXA_6BN
zq|z2ntivGcl*6mVr%tz@mX?O6PShW*ob766hF4lz%HOryta9i%W1by@mmAoaLX(vd
z^LO|MH0gE#)G}X(GdrY}*(2?{MBkt?=ayR<8*_u^kJDJ?R&%q6#5#wn&Dbu}zLTp@
zkA+5`>Cc>$ueed8i5pME)hI$*Gq&2a*`f#PcS=W+J}Qz%@!q|KsoL|BI`uR=%hYkF
z6p~+ZsIvHzDWGjumH9_ZM5JTxd0cmc8XuGWKSd`ls_hpPpP2q@R21_4n2Pk#(NPou
zYl`akHIS2W(Y|N|c^i4i$CW!E(?b24BCJcEI|GlV0_BBpOuBM5Djf6P6YR$AcX9@z
zr1JqjdL8aRlBAMN5?YjOX;>^=*yd#H+Pnk8KwWvQZ5HJ9oWU}PaEZvim1HPMz`~uw
z_x#L>(j0K#qR=5&Af*4>^W>q`RSuv?hgDq=cOtJ`m?_|+BJrI&cq=bnyx?#0iP6l9
z`68Ps9likD0jbq`pSEnI$?=~6O?-w|&g!YfaTldS9;?i%sT{Q+BEvC1aqc~2cKrLq
z>%s3Mjiw-Q>UMcC^LD#TEF<9ycEzv#7eu<dW#~)YG{ykR1CV`iXM;3UvQe$eLtdr%
z!;UR_PC2_dRyn&DhCsfJBdg`k9;>jf;0b~39v(U&N+{8w+O@fs<G-KqPu>geo}T`~
z%Q)GYv29*l*W${?$H#vj>Se7HC+V~r`AXTMtK=aPcV}>Luo~QCWljk(v7eR^z~wWx
zg$ocz^ZiboM7J#Sz08iwN3J(F^ZL5$SKR0YF9KwnDy^{x;(f`9iF=^FGkxGS1HQUe
z4$?Gpzr^Gw9JIy%`5BYz-{&<$y#DZr;ubCg+Eyp7y<1Z70>q`;CxeIBTF;r2+X1q>
zubVu#m9yuQ+i%TXe2w{AVQtw}a=5toRS1k?=qFzgOGiXRY@~{PSl1VXqQ%+%>AxZh
zrF<9y6~gW<g-GRW-=5%R8swuuHX)BY{t=k-I(Tzta1YPTGd95;^9Nb@W+W{96rV+E
zwPi#u?P@&k-XnNIdm;ivUJ*3^kD!^!ON<R&ljr3fQPWv2+=RBeFfF1dB9X+1I;AK9
z;(U0a74U@!jr7d$#{a!5)RFJP7fO1AgMTN(p20VW!JCqWg{3J4tv?fzffV?sc<U*T
zQ4T@p`geEvA&s6`Dq(6E@a`djnd}z<_#7o(V~P?#<v~C1oX+)}sgLx_O~?>|gnBBe
zp&SDcx;rWvBtZHcOOAsM1<L-fa_lRB4<H^n0Dz$%ivIpHHMK!-=x^V@`)R;{gv3OF
z-`AK9n%_WoVjxE|Ch`lsaEF+b+>_<--`KD&&YnUql<Fa-YAGJnI=hx!Jw0ieum7h3
z8{)pWoftAz^I<fM5qKF5bcZnm!|y!4<q6ydP+#__22uM7IbHzjnOqI6s1sb;v7e!I
zM2hy3oz{|x=d`@uUPQ?6td$sGzElYMPmzDAyr2jL5(Z#?0)*1T4*_Z3_}`^{l5gYr
zte7WSgR2ce=H)Eo9T*f;HR7#TKf&YUkIRs7d-@&^{eKFIJ8f?3)FT!`wkb;Z6IV4{
z!E;0a!OHE23I~GCHvCWnO9R-CC;#o+KTFB?4gm;n*@*x2TN#8q?SJt9|JhP!Hwx8U
Wl_C|xIN-s(fYMV{xdK^~@BaaOhavL-

literal 0
HcmV?d00001

diff --git a/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd b/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd
new file mode 100644
index 0000000..317300b
--- /dev/null
+++ b/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd
@@ -0,0 +1,289 @@
+---
+title: "Quick comparison: R and Python GRASS interfaces"
+author: "Veronica Andreo"
+date: 2024-04-01
+date-modified: today
+format:
+  html:
+    toc: true
+    code-tools: true
+    code-copy: true
+    code-fold: false
+categories: [GRASS GIS, Python, R, Intermediate]
+engine: knitr
+execute:
+  eval: false
+---
+
+![](images/R_Python_compare.png){.preview-image width=50%}
+
+In this short tutorial we will highlight the similarities of R and Python GRASS interfaces
+in order to streamline the use of GRASS GIS within R and Python communities.
+As you may know, there's
+an R package called [rgrass](https://github.com/rsbivand/rgrass/) that provides
+basic functionality to read and write data from and into GRASS database as well
+as to execute GRASS tools in either existing or temporary GRASS projects.
+The [GRASS Python API](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html),
+on the other hand, is composed of various packages that provide classes and
+functions for low and high level tasks, including those that can be executed
+with rgrass.
+
+
+
+There are some parallelisms between the
+**rgrass** and **grass.script**/**grass.jupyter** packages, i.e.,
+R and Python interfaces to GRASS GIS.
+Let's review them and go through some examples.
+
+
+| Task                                                       | rgrass function                | GRASS Python API function                               |
+|------------------------------------------------------------|--------------------------------|---------------------------------------------------------|
+| Load library                                               | library(rgrass)                | import grass.script as gs<br>import grass.jupyter as gj |
+| Start GRASS and set all needed <br>environmental variables | initGRASS()                    | gs.setup.init() for scripts,<br>gj.init() for notebooks |
+| Execute GRASS commands                                     | execGRASS()                    | gs.run_command(),<br>gs.read_command(),<br>gs.parse_command()  |
+| Read raster and vector data <br>from GRASS                 | read_RAST(),<br>read_VECT()    | gs.array.array(),<br>n/a                                |
+| Write raster and vector data<br>into GRASS                 | write_RAST(),<br>write_VECT()  | gs.array.write(),<br>n/a                                |
+| Get raster and vector info                                 | n/a,<br>vInfo()                | gs.raster_info(),<br>gs.vector_info()                   |
+| Close GRASS session                                        | unlink_.gislock()              | gs.setup.finish(),<br>gj.finish()                       |
+
+: R and Python GRASS interfaces compared {.striped .hover}
+
+## Comparison examples
+
+Let's see how usage examples would look like.
+
+1. **Load the library**: We need to
+load the libraries that allow us to interface with GRASS GIS
+functionality and (optionally) data. For the Python case, we first need to add
+the GRASS python package path to our system's path.
+
+::: {.panel-tabset}
+
+## R
+
+```{r}
+library(rgrass)
+```
+
+## Python
+
+```{python}
+#| python.reticulate: FALSE
+import sys
+import subprocess
+
+sys.path.append(
+    subprocess.check_output(["grass", "--config", "python_path"], text=True).strip()
+)
+
+import grass.script as gs
+import grass.jupyter as gj
+```
+
+:::
+
+2. **Start a GRASS session**: Once we loaded or imported the packages, we
+start a GRASS session. We need to pass the path to a
+temporary or existing GRASS project.
+In the case of R, `initGRASS` will automatically look for GRASS binaries, alternatively we can
+specify the path to the binaries ourselves.
+In the case of Python, it is worth noting that while grass.script and grass.jupyter init functions
+take the same arguments, `gj.init` also sets other environmental variables to
+streamline work within Jupyter Notebooks, e.g., overwrite is set to true so cells
+can be executed multiple times.
+
+::: {.panel-tabset}
+
+## R
+
+```{r}
+session <- initGRASS(gisBase = "/usr/lib/grass84", # where grass binaries live, `grass --config path`
+                     gisDbase = "/home/user/grassdata", # path to grass database or folder where your project lives
+                     location = "nc_basic_spm_grass7", # existing project name
+                     mapset = "PERMANENT" # mapset name
+                     )
+```
+
+## Python
+
+```{python}
+#| python.reticulate: FALSE
+# With grass.script for scripts
+session = gs.setup.init(path="/home/user/grassdata",
+              location="nc_basic_spm_grass7",
+              mapset="PERMANENT")
+# Optionally, the path to a mapset
+session = gs.setup.init("/home/user/grassdata/nc_basic_spm_grass7/PERMANENT")
+
+# With grass.jupyter for notebooks
+session = gj.init(path="/home/user/grassdata",
+                  location="nc_basic_spm_grass7",
+                  mapset="PERMANENT")
+# Optionally, the path to a mapset
+session = gj.init("~/grassdata/nc_basic_spm_grass7/PERMANENT")
+```
+
+:::
+
+3. **Execute GRASS commands**: Both interfaces work pretty similarly, the
+first argument is always the GRASS tool name and then we pass the parameters
+and flags. While in R we basically use `execGRASS()` for all GRASS commands, in
+the Python API, we have different wrappers to execute GRASS commands depending
+on the nature of their output.
+
+::: {.panel-tabset}
+
+## R
+
+```{r}
+# Map output
+execGRASS("r.slope.aspect",
+          elevation = "elevation",
+          slope = "slope",
+          aspect = "aspect")
+
+# Text output
+execGRASS("g.region",
+          raster = "elevation",
+          flags = "p")
+```
+
+## Python
+
+```{python}
+#| python.reticulate: FALSE
+# Map output
+gs.run_command("r.slope.aspect",
+               elevation="elevation",
+               slope="slope",
+               aspect="aspect")
+# Text output
+print(gs.read_command("g.region",
+                      raster="elevation",
+                      flags="p"))
+# Text output - dictionary
+region = gs.parse_command("g.region",
+                          raster="elevation",
+                          flags="g")
+region
+```
+
+:::
+
+4. **Read raster and vector data into other R or Python formats**:
+*rgrass* functions `read_RAST()` and `read_VECT()` convert GRASS raster and
+vector maps into terra's SpatRaster and SpatVector objects within R.
+In the case of Python, GRASS
+raster maps that can be converted into numpy arrays through
+`gs.array.array()`. Vector attribute data can be converted into
+Pandas data frames in various ways.
+
+::: {.panel-tabset}
+
+## R
+
+```{r}
+# Raster
+elevr <- read_RAST("elevation")
+
+# Vector
+schoolsr <- read_VECT("schools")
+```
+
+## Python
+
+```{python}
+#| python.reticulate: FALSE
+# Raster into numpy array
+elev = gs.array.array("elevation")
+
+# Vector attributes
+import pandas as pd
+schools = gs.parse_command("v.db.select", map="schools", format="json")
+pd.DataFrame(schools["records"])
+
+# Vector geometry and attributes to GeoJSON
+gs.run_command("v.out.ogr", input="schools", output="schools.geojson", format="GeoJSON")
+```
+
+:::
+
+
+5. **Write R or Python objects into GRASS raster and vector maps**: R terra's
+SpatRaster and SpatVector objects can be written (back) into GRASS format with
+`write_RAST()` and `write_VECT()` functions. Within the Python environment,
+numpy arrays can also be written (back) into GRASS raster maps with the
+`write()` method.
+
+::: {.panel-tabset}
+
+## R
+
+```{r}
+# Raster
+write_RAST(elevr, "elevation_r")
+
+# Vector
+write_VECT(schoolsr, "schools_r")
+```
+
+## Python
+
+```{python}
+#| python.reticulate: FALSE
+# Raster
+elev.write(mapname="elev_np", overwrite=True)
+
+# GeoJSON into GRASS vector
+gs.run_command("v.in.ogr", input="schools.geojson", output="schools2")
+```
+
+:::
+
+
+6. **Close GRASS GIS session**: In general, just closing R or Rstudio, as well
+as shutting down Jupyter notebook, will clean up and close the GRASS session
+properly. Sometimes, however, especially if the user changed mapset within the
+workflow, it is better to clean up explicitly before closing.
+
+::: {.panel-tabset}
+
+## R
+
+```{r}
+unlink_.gislock()
+```
+
+## Python
+
+```{python}
+#| python.reticulate: FALSE
+session.finish()
+```
+
+:::
+
+## Final remarks
+
+The examples and comparisons presented here are intended to facilitate the
+combination of tools and languages as well as the exchange of data and format
+conversions. We hope that's useful as a starting point for the implementation
+of different use cases and workflows that suit the needs of users.
+See R and Python tutorials for more examples:
+
+* [GRASS and Python tutorial for beginners](../get_started/fast_track_grass_and_python.qmd)
+* [GRASS and R tutorial for beginners](../get_started/fast_track_grass_and_R.qmd)
+
+## References
+
+* [GRASS Python API docs](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html)
+* [rgrass docs](https://rsbivand.github.io/rgrass/)
+
+
+***
+
+:::{.smaller}
+The development of this tutorial was funded by the US
+[National Science Foundation (NSF)](https://www.nsf.gov/),
+award [2303651](https://www.nsf.gov/awardsearch/showAward?AWD_ID=2303651).
+:::

From 2f3bd2b3c8cb835ddb2060659f7178e7ad58f70d Mon Sep 17 00:00:00 2001
From: Anna Petrasova <kratochanna@gmail.com>
Date: Fri, 27 Sep 2024 16:50:45 -0400
Subject: [PATCH 2/2] change tags

---
 .../execute-results/html.json                             | 8 +++++---
 .../quick_comparison_r_vs_python_grass_interfaces.qmd     | 2 +-
 2 files changed, 6 insertions(+), 4 deletions(-)

diff --git a/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json b/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json
index 085866a..ee88b58 100644
--- a/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json
+++ b/_freeze/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces/execute-results/html.json
@@ -1,8 +1,10 @@
 {
-  "hash": "02b5f6931312f20d479fb4fbc02fb26f",
+  "hash": "b2f5b1384f95aeb0e0a837650dceb10c",
   "result": {
-    "markdown": "---\ntitle: \"Quick comparison: R and Python GRASS interfaces\"\nauthor: \"Veronica Andreo\"\ndate: 2024-04-01\ndate-modified: today\nformat:\n  html:\n    toc: true\n    code-tools: true\n    code-copy: true\n    code-fold: false\ncategories: [GRASS GIS, Python, R, Intermediate]\nengine: knitr\nexecute:\n  eval: false\n---\n\n\n![](images/R_Python_compare.png){.preview-image width=50%}\n\nIn this short tutorial we will highlight the similarities of R and Python GRASS interfaces\nin order to streamline the use of GRASS GIS within R and Python communities.\nAs you may know, there's\nan R package called [rgrass](https://github.com/rsbivand/rgrass/) that provides\nbasic functionality to read and write data from and into GRASS database as well\nas to execute GRASS tools in either existing or temporary GRASS projects.\nThe [GRASS Python API](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html),\non the other hand, is composed of various packages that provide classes and\nfunctions for low and high level tasks, including those that can be executed\nwith rgrass.\n\n\n\nThere are some parallelisms between the\n**rgrass** and **grass.script**/**grass.jupyter** packages, i.e.,\nR and Python interfaces to GRASS GIS.\nLet's review them and go through some examples.\n\n\n| Task                                                       | rgrass function                | GRASS Python API function                               |\n|------------------------------------------------------------|--------------------------------|---------------------------------------------------------|\n| Load library                                               | library(rgrass)                | import grass.script as gs<br>import grass.jupyter as gj |\n| Start GRASS and set all needed <br>environmental variables | initGRASS()                    | gs.setup.init() for scripts,<br>gj.init() for notebooks |\n| Execute GRASS commands                                     | execGRASS()                    | gs.run_command(),<br>gs.read_command(),<br>gs.parse_command()  |\n| Read raster and vector data <br>from GRASS                 | read_RAST(),<br>read_VECT()    | gs.array.array(),<br>n/a                                |\n| Write raster and vector data<br>into GRASS                 | write_RAST(),<br>write_VECT()  | gs.array.write(),<br>n/a                                |\n| Get raster and vector info                                 | n/a,<br>vInfo()                | gs.raster_info(),<br>gs.vector_info()                   |\n| Close GRASS session                                        | unlink_.gislock()              | gs.setup.finish(),<br>gj.finish()                       |\n\n: R and Python GRASS interfaces compared {.striped .hover}\n\n## Comparison examples\n\nLet's see how usage examples would look like.\n\n1. **Load the library**: We need to\nload the libraries that allow us to interface with GRASS GIS\nfunctionality and (optionally) data. For the Python case, we first need to add\nthe GRASS python package path to our system's path.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nlibrary(rgrass)\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\nimport sys\nimport subprocess\n\nsys.path.append(\n    subprocess.check_output([\"grass\", \"--config\", \"python_path\"], text=True).strip()\n)\n\nimport grass.script as gs\nimport grass.jupyter as gj\n```\n:::\n\n\n:::\n\n2. **Start a GRASS session**: Once we loaded or imported the packages, we\nstart a GRASS session. We need to pass the path to a\ntemporary or existing GRASS project.\nIn the case of R, `initGRASS` will automatically look for GRASS binaries, alternatively we can\nspecify the path to the binaries ourselves.\nIn the case of Python, it is worth noting that while grass.script and grass.jupyter init functions\ntake the same arguments, `gj.init` also sets other environmental variables to\nstreamline work within Jupyter Notebooks, e.g., overwrite is set to true so cells\ncan be executed multiple times.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nsession <- initGRASS(gisBase = \"/usr/lib/grass84\", # where grass binaries live, `grass --config path`\n                     gisDbase = \"/home/user/grassdata\", # path to grass database or folder where your project lives\n                     location = \"nc_basic_spm_grass7\", # existing project name\n                     mapset = \"PERMANENT\" # mapset name\n                     )\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# With grass.script for scripts\nsession = gs.setup.init(path=\"/home/user/grassdata\",\n              location=\"nc_basic_spm_grass7\",\n              mapset=\"PERMANENT\")\n# Optionally, the path to a mapset\nsession = gs.setup.init(\"/home/user/grassdata/nc_basic_spm_grass7/PERMANENT\")\n\n# With grass.jupyter for notebooks\nsession = gj.init(path=\"/home/user/grassdata\",\n                  location=\"nc_basic_spm_grass7\",\n                  mapset=\"PERMANENT\")\n# Optionally, the path to a mapset\nsession = gj.init(\"~/grassdata/nc_basic_spm_grass7/PERMANENT\")\n```\n:::\n\n\n:::\n\n3. **Execute GRASS commands**: Both interfaces work pretty similarly, the\nfirst argument is always the GRASS tool name and then we pass the parameters\nand flags. While in R we basically use `execGRASS()` for all GRASS commands, in\nthe Python API, we have different wrappers to execute GRASS commands depending\non the nature of their output.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Map output\nexecGRASS(\"r.slope.aspect\",\n          elevation = \"elevation\",\n          slope = \"slope\",\n          aspect = \"aspect\")\n\n# Text output\nexecGRASS(\"g.region\",\n          raster = \"elevation\",\n          flags = \"p\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Map output\ngs.run_command(\"r.slope.aspect\",\n               elevation=\"elevation\",\n               slope=\"slope\",\n               aspect=\"aspect\")\n# Text output\nprint(gs.read_command(\"g.region\",\n                      raster=\"elevation\",\n                      flags=\"p\"))\n# Text output - dictionary\nregion = gs.parse_command(\"g.region\",\n                          raster=\"elevation\",\n                          flags=\"g\")\nregion\n```\n:::\n\n\n:::\n\n4. **Read raster and vector data into other R or Python formats**:\n*rgrass* functions `read_RAST()` and `read_VECT()` convert GRASS raster and\nvector maps into terra's SpatRaster and SpatVector objects within R.\nIn the case of Python, GRASS\nraster maps that can be converted into numpy arrays through\n`gs.array.array()`. Vector attribute data can be converted into\nPandas data frames in various ways.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Raster\nelevr <- read_RAST(\"elevation\")\n\n# Vector\nschoolsr <- read_VECT(\"schools\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Raster into numpy array\nelev = gs.array.array(\"elevation\")\n\n# Vector attributes\nimport pandas as pd\nschools = gs.parse_command(\"v.db.select\", map=\"schools\", format=\"json\")\npd.DataFrame(schools[\"records\"])\n\n# Vector geometry and attributes to GeoJSON\ngs.run_command(\"v.out.ogr\", input=\"schools\", output=\"schools.geojson\", format=\"GeoJSON\")\n```\n:::\n\n\n:::\n\n\n5. **Write R or Python objects into GRASS raster and vector maps**: R terra's\nSpatRaster and SpatVector objects can be written (back) into GRASS format with\n`write_RAST()` and `write_VECT()` functions. Within the Python environment,\nnumpy arrays can also be written (back) into GRASS raster maps with the\n`write()` method.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Raster\nwrite_RAST(elevr, \"elevation_r\")\n\n# Vector\nwrite_VECT(schoolsr, \"schools_r\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Raster\nelev.write(mapname=\"elev_np\", overwrite=True)\n\n# GeoJSON into GRASS vector\ngs.run_command(\"v.in.ogr\", input=\"schools.geojson\", output=\"schools2\")\n```\n:::\n\n\n:::\n\n\n6. **Close GRASS GIS session**: In general, just closing R or Rstudio, as well\nas shutting down Jupyter notebook, will clean up and close the GRASS session\nproperly. Sometimes, however, especially if the user changed mapset within the\nworkflow, it is better to clean up explicitly before closing.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nunlink_.gislock()\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\nsession.finish()\n```\n:::\n\n\n:::\n\n## Final remarks\n\nThe examples and comparisons presented here are intended to facilitate the\ncombination of tools and languages as well as the exchange of data and format\nconversions. We hope that's useful as a starting point for the implementation\nof different use cases and workflows that suit the needs of users.\nSee R and Python tutorials for more examples:\n\n* [GRASS and Python tutorial for beginners](../get_started/fast_track_grass_and_python.qmd)\n* [GRASS and R tutorial for beginners](../get_started/fast_track_grass_and_R.qmd)\n\n## References\n\n* [GRASS Python API docs](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html)\n* [rgrass docs](https://rsbivand.github.io/rgrass/)\n\n\n***\n\n:::{.smaller}\nThe development of this tutorial was funded by the US\n[National Science Foundation (NSF)](https://www.nsf.gov/),\naward [2303651](https://www.nsf.gov/awardsearch/showAward?AWD_ID=2303651).\n:::\n",
-    "supporting": [],
+    "markdown": "---\ntitle: \"Quick comparison: R and Python GRASS interfaces\"\nauthor: \"Veronica Andreo\"\ndate: 2024-04-01\ndate-modified: today\nformat:\n  html:\n    toc: true\n    code-tools: true\n    code-copy: true\n    code-fold: false\ncategories: [Python, R, intermediate]\nengine: knitr\nexecute:\n  eval: false\n---\n\n\n![](images/R_Python_compare.png){.preview-image width=50%}\n\nIn this short tutorial we will highlight the similarities of R and Python GRASS interfaces\nin order to streamline the use of GRASS GIS within R and Python communities.\nAs you may know, there's\nan R package called [rgrass](https://github.com/rsbivand/rgrass/) that provides\nbasic functionality to read and write data from and into GRASS database as well\nas to execute GRASS tools in either existing or temporary GRASS projects.\nThe [GRASS Python API](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html),\non the other hand, is composed of various packages that provide classes and\nfunctions for low and high level tasks, including those that can be executed\nwith rgrass.\n\n\n\nThere are some parallelisms between the\n**rgrass** and **grass.script**/**grass.jupyter** packages, i.e.,\nR and Python interfaces to GRASS GIS.\nLet's review them and go through some examples.\n\n\n| Task                                                       | rgrass function                | GRASS Python API function                               |\n|------------------------------------------------------------|--------------------------------|---------------------------------------------------------|\n| Load library                                               | library(rgrass)                | import grass.script as gs<br>import grass.jupyter as gj |\n| Start GRASS and set all needed <br>environmental variables | initGRASS()                    | gs.setup.init() for scripts,<br>gj.init() for notebooks |\n| Execute GRASS commands                                     | execGRASS()                    | gs.run_command(),<br>gs.read_command(),<br>gs.parse_command()  |\n| Read raster and vector data <br>from GRASS                 | read_RAST(),<br>read_VECT()    | gs.array.array(),<br>n/a                                |\n| Write raster and vector data<br>into GRASS                 | write_RAST(),<br>write_VECT()  | gs.array.write(),<br>n/a                                |\n| Get raster and vector info                                 | n/a,<br>vInfo()                | gs.raster_info(),<br>gs.vector_info()                   |\n| Close GRASS session                                        | unlink_.gislock()              | gs.setup.finish(),<br>gj.finish()                       |\n\n: R and Python GRASS interfaces compared {.striped .hover}\n\n## Comparison examples\n\nLet's see how usage examples would look like.\n\n1. **Load the library**: We need to\nload the libraries that allow us to interface with GRASS GIS\nfunctionality and (optionally) data. For the Python case, we first need to add\nthe GRASS python package path to our system's path.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nlibrary(rgrass)\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\nimport sys\nimport subprocess\n\nsys.path.append(\n    subprocess.check_output([\"grass\", \"--config\", \"python_path\"], text=True).strip()\n)\n\nimport grass.script as gs\nimport grass.jupyter as gj\n```\n:::\n\n\n:::\n\n2. **Start a GRASS session**: Once we loaded or imported the packages, we\nstart a GRASS session. We need to pass the path to a\ntemporary or existing GRASS project.\nIn the case of R, `initGRASS` will automatically look for GRASS binaries, alternatively we can\nspecify the path to the binaries ourselves.\nIn the case of Python, it is worth noting that while grass.script and grass.jupyter init functions\ntake the same arguments, `gj.init` also sets other environmental variables to\nstreamline work within Jupyter Notebooks, e.g., overwrite is set to true so cells\ncan be executed multiple times.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nsession <- initGRASS(gisBase = \"/usr/lib/grass84\", # where grass binaries live, `grass --config path`\n                     gisDbase = \"/home/user/grassdata\", # path to grass database or folder where your project lives\n                     location = \"nc_basic_spm_grass7\", # existing project name\n                     mapset = \"PERMANENT\" # mapset name\n                     )\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# With grass.script for scripts\nsession = gs.setup.init(path=\"/home/user/grassdata\",\n              location=\"nc_basic_spm_grass7\",\n              mapset=\"PERMANENT\")\n# Optionally, the path to a mapset\nsession = gs.setup.init(\"/home/user/grassdata/nc_basic_spm_grass7/PERMANENT\")\n\n# With grass.jupyter for notebooks\nsession = gj.init(path=\"/home/user/grassdata\",\n                  location=\"nc_basic_spm_grass7\",\n                  mapset=\"PERMANENT\")\n# Optionally, the path to a mapset\nsession = gj.init(\"~/grassdata/nc_basic_spm_grass7/PERMANENT\")\n```\n:::\n\n\n:::\n\n3. **Execute GRASS commands**: Both interfaces work pretty similarly, the\nfirst argument is always the GRASS tool name and then we pass the parameters\nand flags. While in R we basically use `execGRASS()` for all GRASS commands, in\nthe Python API, we have different wrappers to execute GRASS commands depending\non the nature of their output.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Map output\nexecGRASS(\"r.slope.aspect\",\n          elevation = \"elevation\",\n          slope = \"slope\",\n          aspect = \"aspect\")\n\n# Text output\nexecGRASS(\"g.region\",\n          raster = \"elevation\",\n          flags = \"p\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Map output\ngs.run_command(\"r.slope.aspect\",\n               elevation=\"elevation\",\n               slope=\"slope\",\n               aspect=\"aspect\")\n# Text output\nprint(gs.read_command(\"g.region\",\n                      raster=\"elevation\",\n                      flags=\"p\"))\n# Text output - dictionary\nregion = gs.parse_command(\"g.region\",\n                          raster=\"elevation\",\n                          flags=\"g\")\nregion\n```\n:::\n\n\n:::\n\n4. **Read raster and vector data into other R or Python formats**:\n*rgrass* functions `read_RAST()` and `read_VECT()` convert GRASS raster and\nvector maps into terra's SpatRaster and SpatVector objects within R.\nIn the case of Python, GRASS\nraster maps that can be converted into numpy arrays through\n`gs.array.array()`. Vector attribute data can be converted into\nPandas data frames in various ways.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Raster\nelevr <- read_RAST(\"elevation\")\n\n# Vector\nschoolsr <- read_VECT(\"schools\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Raster into numpy array\nelev = gs.array.array(\"elevation\")\n\n# Vector attributes\nimport pandas as pd\nschools = gs.parse_command(\"v.db.select\", map=\"schools\", format=\"json\")\npd.DataFrame(schools[\"records\"])\n\n# Vector geometry and attributes to GeoJSON\ngs.run_command(\"v.out.ogr\", input=\"schools\", output=\"schools.geojson\", format=\"GeoJSON\")\n```\n:::\n\n\n:::\n\n\n5. **Write R or Python objects into GRASS raster and vector maps**: R terra's\nSpatRaster and SpatVector objects can be written (back) into GRASS format with\n`write_RAST()` and `write_VECT()` functions. Within the Python environment,\nnumpy arrays can also be written (back) into GRASS raster maps with the\n`write()` method.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\n# Raster\nwrite_RAST(elevr, \"elevation_r\")\n\n# Vector\nwrite_VECT(schoolsr, \"schools_r\")\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\n# Raster\nelev.write(mapname=\"elev_np\", overwrite=True)\n\n# GeoJSON into GRASS vector\ngs.run_command(\"v.in.ogr\", input=\"schools.geojson\", output=\"schools2\")\n```\n:::\n\n\n:::\n\n\n6. **Close GRASS GIS session**: In general, just closing R or Rstudio, as well\nas shutting down Jupyter notebook, will clean up and close the GRASS session\nproperly. Sometimes, however, especially if the user changed mapset within the\nworkflow, it is better to clean up explicitly before closing.\n\n::: {.panel-tabset}\n\n## R\n\n\n::: {.cell}\n\n```{.r .cell-code}\nunlink_.gislock()\n```\n:::\n\n\n## Python\n\n\n::: {.cell python.reticulate='false'}\n\n```{.python .cell-code}\nsession.finish()\n```\n:::\n\n\n:::\n\n## Final remarks\n\nThe examples and comparisons presented here are intended to facilitate the\ncombination of tools and languages as well as the exchange of data and format\nconversions. We hope that's useful as a starting point for the implementation\nof different use cases and workflows that suit the needs of users.\nSee R and Python tutorials for more examples:\n\n* [GRASS and Python tutorial for beginners](../get_started/fast_track_grass_and_python.qmd)\n* [GRASS and R tutorial for beginners](../get_started/fast_track_grass_and_R.qmd)\n\n## References\n\n* [GRASS Python API docs](https://grass.osgeo.org/grass-stable/manuals/libpython/index.html)\n* [rgrass docs](https://rsbivand.github.io/rgrass/)\n\n\n***\n\n:::{.smaller}\nThe development of this tutorial was funded by the US\n[National Science Foundation (NSF)](https://www.nsf.gov/),\naward [2303651](https://www.nsf.gov/awardsearch/showAward?AWD_ID=2303651).\n:::\n",
+    "supporting": [
+      "quick_comparison_r_vs_python_grass_interfaces_files"
+    ],
     "filters": [
       "rmarkdown/pagebreak.lua"
     ],
diff --git a/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd b/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd
index 317300b..81b40d2 100644
--- a/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd
+++ b/content/tutorials/r_python_interfaces_comparison/quick_comparison_r_vs_python_grass_interfaces.qmd
@@ -9,7 +9,7 @@ format:
     code-tools: true
     code-copy: true
     code-fold: false
-categories: [GRASS GIS, Python, R, Intermediate]
+categories: [Python, R, intermediate]
 engine: knitr
 execute:
   eval: false