From 5af6b3cd5fc57bf2c79c17dc0298bb28bb0dcd5f Mon Sep 17 00:00:00 2001 From: Jeffrey Dick Date: Fri, 17 Jan 2025 10:26:48 +0000 Subject: [PATCH] OBIGT: Update parameters for cassiterite git-svn-id: svn://scm.r-forge.r-project.org/svnroot/chnosz/pkg/CHNOSZ@868 edb9625f-4e0d-4859-8d74-9fd3b1da38cb --- DESCRIPTION | 4 ++-- inst/NEWS.Rd | 20 +++++++++++--------- inst/extdata/OBIGT/SLOP98.csv | 2 ++ inst/extdata/OBIGT/SUPCRT92.csv | 1 + inst/extdata/OBIGT/inorganic_aq.csv | 4 ++-- inst/extdata/OBIGT/inorganic_cr.csv | 2 +- inst/extdata/misc/OBIGT_check.csv | 2 +- inst/extdata/thermo/refs.csv | 2 +- inst/extdata/thermo/stoich.csv.xz | Bin 20024 -> 20024 bytes 9 files changed, 21 insertions(+), 16 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index fdc7250..36f9f9e 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,6 +1,6 @@ -Date: 2025-01-08 +Date: 2025-01-17 Package: CHNOSZ -Version: 2.1.0-39 +Version: 2.1.0-40 Title: Thermodynamic Calculations and Diagrams for Geochemistry Authors@R: c( person("Jeffrey", "Dick", , "j3ffdick@gmail.com", role = c("aut", "cre"), diff --git a/inst/NEWS.Rd b/inst/NEWS.Rd index a40c727..9658474 100644 --- a/inst/NEWS.Rd +++ b/inst/NEWS.Rd @@ -15,7 +15,7 @@ \newcommand{\Cp}{\ifelse{latex}{\eqn{C_P}}{\ifelse{html}{\out{CP}}{Cp}}} \newcommand{\DG0}{\ifelse{latex}{\eqn{{\Delta}G^{\circ}}}{\ifelse{html}{\out{ΔG°}}{ΔG°}}} -\section{Changes in CHNOSZ version 2.1.0-39 (2025-01-08)}{ +\section{Changes in CHNOSZ version 2.1.0-40 (2025-01-17)}{ \subsection{OBIGT DEFAULT DATA}{ \itemize{ @@ -36,15 +36,16 @@ \item Add tantalum species from \href{https://doi.org/10.1016/j.gca.2024.10.019}{Hu et al. (2024)}. - \item Add Sn\S{+2}, Th\S{+4}, and Sn chloride complexes from - \href{https://doi.org/10.1016/j.gsf.2023.101624}{Liu et al. (2023)}. + \item Add Sn chloride complexes and updated parameters for Sn\S{+2}, + Th\S{+4} from \href{https://doi.org/10.1016/j.gsf.2023.101624}{Liu et al. + (2023)}. \item Move americium complexes from \file{SLOP98.csv} back to \file{inorganic_aq.csv} (entropy of the element Am has been available for checking GHS self-consistency since version 1.4.0). - \item \file{inorganic_cr.csv}: Add cerianite (CeO\s{2}) and - chromite (FeCr\s{2}O\s{4}) from + \item \file{inorganic_cr.csv}: Add cerianite (CeO\s{2}), chromite + (FeCr\s{2}O\s{4}), and cassiterite (SnO\s{2}) from \href{https://doi.org/10.3133/b2131}{Robie and Hemingway (1995)}. \item Restore steam (H\s{2}O,g) in \file{inorganic_gas.csv}. @@ -83,10 +84,11 @@ \subsection{OTHER CHANGES}{ \itemize{ - \item \code{OBIGT} now lists \samp{CGL_Ttr} in \code{model} for species - whose \code{T} value denotes a phase stability limit (as opposed to a Cp - equation limit). This replaces the (overloaded) use of signed values of - \code{T}. + \item \code{OBIGT} now uses \samp{CGL_Ttr} as the \code{model} for + species whose \code{T} value denotes a phase stability limit above which + \code{subcrt()} returns NA values (in contrast to \samp{CGL}, for which + the value of \code{T} is a Cp equation limit that only produces a + warning). This replaces the use of signed values of \code{T}. \item \code{check.EOS()} now uses values of Born coefficients \emph{X} and \emph{Q} that are consistent with either SUPCRT92 or DEW, depending diff --git a/inst/extdata/OBIGT/SLOP98.csv b/inst/extdata/OBIGT/SLOP98.csv index aa9a1a3..a7f38e6 100644 --- a/inst/extdata/OBIGT/SLOP98.csv +++ b/inst/extdata/OBIGT/SLOP98.csv @@ -148,3 +148,5 @@ PdCl2,PdCl2,PdCl2,aq,SS98a.3,NA,2015-08-16,HKF,cal,-34080,-51390,4.2,37.8,32.4,6 PdCl3-,PdCl3-,PdCl3-,aq,SS98a.3,NA,2015-08-16,HKF,cal,-68390,-98890,2.5,44.6,62.9,10.9057,18.8497,-1.6645,-3.5581,46.9452,6.0491,1.5895,-1 PdCl4-2,PdCl4-2,PdCl4-2,aq,SS98a.3,NA,2015-08-16,HKF,cal,-102130,-149470,-11.7,41.4,96.8,16.1552,31.6662,-6.698,-4.088,61.6845,5.4024,3.3914,-2 CoCl+,CoCl+,CoCl+,aq,SSH97,NA,2007-10-16,HKF,cal,-45157,-53965,-11.27,17.02,-1.5,1.8028,-3.3766,7.0702,-2.6394,22.7656,0.4323,0.7191,1 +Th+4,Th+4,Th+4,aq,SSWS97,NA,1997-11-07,HKF,cal,-168500,-183800,-101,-0.2,-53.4,-4.2886,-18.25,12.9154,-2.0244,40.2222,-3.0752,3.7093,4 +Sn+2,Sn+2,Sn+2,aq,SSWS97,NA,1997-11-11,HKF,cal,-6570,-2100,-4,-15.5,-15.6,0.0094,-7.7516,8.781,-2.4584,7.416,-6.1919,1.1216,2 diff --git a/inst/extdata/OBIGT/SUPCRT92.csv b/inst/extdata/OBIGT/SUPCRT92.csv index cbc8bfc..870f4df 100644 --- a/inst/extdata/OBIGT/SUPCRT92.csv +++ b/inst/extdata/OBIGT/SUPCRT92.csv @@ -31,6 +31,7 @@ brucite,Brc,Mg(OH)2,cr,HDNB78,NA,1978-05-05,CGL,cal,-199646,-221390,15.09,NA,24. Ca-Al-pyroxene,Ca-Al-Px,CaAl(AlSi)O6,cr,HDNB78,SPRONS92.2,1990-03-15,CGL,cal,-742067,-783793,35,NA,63.5,54.13,0.00642,-1490000,0,0,0,0,1400 Ca-phillipsite,Ca-Php,Ca(Al2Si5)O14*5H2O,cr,HDNB78,NA,1978-05-05,CGL,cal,NA,NA,166.6,NA,265,145.82,0.05267,-1913000,0,0,0,0,848 Ca-phillipsite,Ca-Php,Ca(Al2Si5)O14*5H2O,cr2,HDNB78,OBIGT.1,1978-05-05,CGL,cal,-3.895,NA,165.3767,NA,265,149.01,0.04641,-1643000,0,0,0,0,1000 +cassiterite,Cst,SnO2,cr,JH85,RH95.7,1978-05-05,CGL_Ttr,cal,-124260,-138800,12.5,NA,21.55,17.246,0.0028026,-490010,0,0,0,0,1903 calcite,Cal,CaCO3,cr,HDNB78,PB82,1990-03-09,CGL,cal,-269880,-288552,22.15,NA,36.934,24.98,0.00524,-620000,0,0,0,0,1200 celadonite,Cln,K(MgAl)Si4O10(OH)2,cr,HDNB78,NA,1978-05-05,CGL,cal,NA,NA,74.9,NA,157.1,80.25,0.0253,-1854000,0,0,0,0,1000 chabazite,Cbz,Ca(Al2Si4)O12*6H2O,cr,HDNB78,NA,1978-05-05,CGL,cal,NA,NA,152.9,NA,247.76,146,0.04447,-1643000,0,0,0,0,1000 diff --git a/inst/extdata/OBIGT/inorganic_aq.csv b/inst/extdata/OBIGT/inorganic_aq.csv index 6dc463f..6b3bd27 100644 --- a/inst/extdata/OBIGT/inorganic_aq.csv +++ b/inst/extdata/OBIGT/inorganic_aq.csv @@ -916,5 +916,5 @@ Sn+4,Sn+4,Sn+4,aq,LYX23,NA,2025-01-05,HKF,cal,600,-18083.8,-112.9,-56.12,-61.76, Th+4,Th+4,Th+4,aq,LYX23,SSWS97.6,2025-01-05,HKF,cal,-168500,-183800,-101,-53.54,-60.6,-5.2739,-20.6582,13.8696,-1.925,8.9637,-13.9401,3.7093,4 Sn+2,Sn+2,Sn+2,aq,LYX23,SSWS97.6,2025-01-05,HKF,cal,-6570,-2100,-4,-11.14,-15.5584,0.0142,-7.7466,8.7948,-2.4588,9.9465,-5.3042,1.119,2 SnCl+,SnCl+,SnCl+,aq,LYX23,NA,2025-01-05,HKF,cal,-39962,-41670,17.5442,8.7913,9.2717,3.1277,-0.1446,5.8069,-2.773,13.8796,-1.2438,0.278,1 -SnCl2,SnCl2,SnCl2,aq,LYX23,NA,2025-01-05,HKF,cal,-72232,-81852,32.2825,13.9082,37.0637,6.8239,8.8804,2.2598,-3.1461,13.9676,-0.2015,-0.038,0 -SnCl3-,SnCl3-,SnCl3-,aq,LYX23,NA,2025-01-05,HKF,cal,-103547,-121668,45.0481,4.2091,68.0292,11.3814,20.0081,-2.1139,-3.6061,17.0106,-2.1772,0.9093,-1 +SnCl2,SnCl2,SnCl2,aq,LYX23,NA,2025-01-05,HKF,cal,-72232,-81852,33.2825,13.9082,37.0637,6.8239,8.8804,2.2598,-3.1461,13.9676,-0.2015,-0.038,0 +SnCl3-,SnCl3-,SnCl3-,aq,LYX23,NA,2025-01-05,HKF,cal,-103547,-121668,47.0481,4.2091,68.0292,11.3814,20.0081,-2.1139,-3.6061,17.0106,-2.1772,0.9093,-1 diff --git a/inst/extdata/OBIGT/inorganic_cr.csv b/inst/extdata/OBIGT/inorganic_cr.csv index 7adfb4d..da6bb7b 100644 --- a/inst/extdata/OBIGT/inorganic_cr.csv +++ b/inst/extdata/OBIGT/inorganic_cr.csv @@ -64,7 +64,6 @@ sylvite,Sy,KCl,cr,HDNB78,NA,1978-05-05,CGL_Ttr,cal,-97735,-104370,19.73,NA,37.52 tenorite,Tn,CuO,cr,HDNB78,NA,1978-05-05,CGL,cal,-30568,-37200,10.18,NA,12.22,11.53,0.00188,-176000,0,0,0,0,1600 wurtzite,Wur,ZnS,cr,HDNB78,NA,1978-05-05,CGL,cal,-44810,-45850,14.064,NA,23.846,11.82,0.00116,-104000,0,0,0,0,1300 berndtite,Brn,SnS2,cr,JH85,NA,1978-05-05,CGL,cal,-34750,-36700,20.9,NA,40.96,15.51,0.0042,0,0,0,0,0,1000 -cassiterite,Cst,SnO2,cr,JH85,RH95.7,1978-05-05,CGL_Ttr,cal,-124260,-138800,12.5,NA,21.55,17.246,0.0028026,-490010,0,0,0,0,1903 romarchite,Sn-Ox,SnO,cr,JH85,NA,1978-05-05,CGL,cal,-61459,-68340,13.66,NA,20.895,9.55,0.0035,0,0,0,0,0,1237 tin,Sn,Sn,cr,JH85,NA,1985-08-00,CGL,cal,0,0,12.24,NA,16.289,4.42,0.0063,0,0,0,0,0,505.06 tin,Sn,Sn,cr2,JH85,OBIGT.1,1985-08-00,CGL,cal,679.802,1607.605,15.3519,NA,16.289,7.3,0,0,0,0,0,0,3000 @@ -160,3 +159,4 @@ beta-UO2(OH)2,Sho,UO2(OH)2,cr,NEA20,NA,2023-08-07,CGL_Ttr,J,-1398683,-1533800,13 Na2U2O7,Na2U2O7,Na2U2O7,cr,NEA20,NA,2023-08-07,CGL_Ttr,J,-3011454,-3203800,275.9,227.3,0,262.83,0.014653,-3549000,0,0,0,0,1273 chromite,Chr,FeCr2O4,cr,RH95,NA,2024-05-23,CGL,J,-1344500,-1445500,146,133.58,44.01,301.8,-0.04157,487700,-2803,1.147E-05,0,0,1800 cerianite,Cei,CeO2,cr,RH95,NA,2024-05-23,CGL,J,-1024600,-1088700,62.3,61.63,23.85,80.29,0.005699,-729400,-209.9,0,0,0,1800 +cassiterite,Cst,SnO2,cr,RH95,NA,2025-01-17,CGL_Ttr,J,-515800,-577600,49.02,53.22,21.55,76.04,0.007364,-2224000,0,0,0,0,1903 diff --git a/inst/extdata/misc/OBIGT_check.csv b/inst/extdata/misc/OBIGT_check.csv index 8f5d3bb..df6a6bb 100644 --- a/inst/extdata/misc/OBIGT_check.csv +++ b/inst/extdata/misc/OBIGT_check.csv @@ -135,7 +135,6 @@ "OBIGT","TaF7-2","aq","cal",-9.97,-4.13, "OBIGT","Sn+4","aq","cal",,-1.18, "OBIGT","Th+4","aq","cal",,-1.13, -"OBIGT","SnCl3-","aq","cal",,,596 "OBIGT","CH4","aq","cal",-2.61,, "OBIGT","hexane","aq","cal",1.43,3.21, "OBIGT","ethylene","aq","cal",6.12,-3.82, @@ -392,6 +391,7 @@ "SLOP98","MnAsO4-","aq","cal",-1.45,, "SLOP98","AmCrO4+2","aq","cal",,,-4654 "SLOP98","PdCl4-2","aq","cal",,-1.05, +"SLOP98","Th+4","aq","cal",,-1.14, "SUPCRT92","almandine","cr","cal",4.43,, "SUPCRT92","antigorite","cr","cal",,,812 "SUPCRT92","clinochlore,7A","cr","cal",,,666 diff --git a/inst/extdata/thermo/refs.csv b/inst/extdata/thermo/refs.csv index c73f596..e80d8bf 100644 --- a/inst/extdata/thermo/refs.csv +++ b/inst/extdata/thermo/refs.csv @@ -64,7 +64,7 @@ OC95,"P. A. G. O'Hare and L. A. Curtiss",1995,"J. Chem. Thermodynamics 27, 643-6 HSS95,"J. R. Haas, E. L. Shock and D. C. Sassani",1995,"Geochim. Cosmochim. Acta 59, 4329-4350","complexes of rare earth elements",https://doi.org/10.1016/0016-7037(95)00314-P PH95,"V. A. Pokrovskii and H. C. Helgeson",1995,"Am. J. Sci. 295, 1255-1342","aluminum species",https://doi.org/10.2475/ajs.295.10.1255 PK95,"V. B. Parker and I. L. Khodakovskii",1995,"J. Phys. Chem. Ref. Data 24, 1699-1745",melanterite,https://doi.org/10.1063/1.555964 -RH95,"R. A. Robie and B. S. Hemingway",1995,"U. S. Geological Survey Bull. 2131","bixbyite, cattierite, cerianite, chromite, cobalt, cobalt monoxide, guite, gypsum, hausmannite, huebnerite, linnaeite, manganese, manganosite, pyrolusite, willemite, wustite, zinc",https://doi.org/10.3133/b2131 +RH95,"R. A. Robie and B. S. Hemingway",1995,"U. S. Geological Survey Bull. 2131","bixbyite, cassiterite, cattierite, cerianite, chromite, cobalt, cobalt monoxide, guite, gypsum, hausmannite, huebnerite, linnaeite, manganese, manganosite, pyrolusite, willemite, wustite, zinc",https://doi.org/10.3133/b2131 RH95.1,"R. A. Robie and B. S. Hemingway",1995,"U. S. Geological Survey Bull. 2131","hydrogen fluoride and hydrogen chloride",https://doi.org/10.3133/b2131 RH95.2,"R. A. Robie and B. S. Hemingway",1995,"U. S. Geological Survey Bull. 2131","dawsonite: Cp coefficients corrected in @TKSS14; Cp value at 25 °C from @BPAH07, citing @FSR76",https://doi.org/10.3133/b2131 RH95.3,"R. A. Robie and B. S. Hemingway",1995,"U. S. Geological Survey Bull. 2131","almandine, dickite, glaucophane, grunerite, halloysite, pyrope: GHS and Cp at 25 °C",https://doi.org/10.3133/b2131 diff --git a/inst/extdata/thermo/stoich.csv.xz b/inst/extdata/thermo/stoich.csv.xz index 0006aa64e5f523cfce63e2dfa3c46f26d501378b..573d25f615d2f549767fc11a9b84aba6e5ef89bb 100644 GIT binary patch delta 4019 zcmV;k4@~g5oB_C;0kE1re^@F?5sZt4#%mJeq@DlW5@YbrORZb_lfEVXjbu9TeULmN znx5fq;qng~FqrdURF%Y~_^?5!Mp>i2v5&^sNKpZidKc1+Jt$1OfZpUZA#?4Q`MjEG z+{D4|goEqrqZ~U+47NZ2d|)sF8Qq=-!cWPnGHz;_QBlRPl5M2^f4K-rxUFmci>~QM z2J;}I`7NIaNFRJd`LA76q^kHA_UVW#8r;p=xD$C3LVj3M7eNImwOXZ+Ozb?4;2sX+qIz@=AitfSH* zwQ&o_m3q>@eix~|f0Z~7N36zTq|Wla{tP<1@9bCqy=+%Gz8fiiogl>5Ugoy{3OsvL>p306kFJ8`SFIW=XM&145I}Noo!V$$2qGKoy<-z zp@WB8_Q>K&^U>!xAq-!{Y8hiB)N8_&ir*(y6)+I|$C5O)%3kHuzy-kmN=eq0Eh@(WK+EjkG!s_8qIe`D zZN>k9f9i>K|BTd};j2Rk?BZib`D4qhC(hJ>b$rOE%HG}{>n@i{3#KxgXI5PdSvfuA9trR- z7$x^fK33v{a<%de!T_$DJTr%q94|i351*y6e*(sz*zTKsqeiV8&36Y{COPIY)Kv9N z6~rxB14nE1^9=RSVT4?Ks3-BcF&Za3zAi&&@n$tb+LQd&R`Z|ro&B|9oIVbAVq(~t zcx8suD%)Vl(`{gG;A|UtJ;zI+3=@3*d9u_BPJo;W{Z-t06cM839U%*x1&#?>h zu=Ws;6!%}*6auGRJ$IdG1H#bLSng4>I>@KYOY>pyKwz{OTASoqYTV6X$@Q8ykgxZs zK`_lnf`)+TT3;^Zq*%Y5kd_c*_=mZNehA) zest`q{`S+cBs>#Z3&l=RE$Z8h8zW2!dSqV@vutI59C!BYYXQpuIAuUWTU`atyqfES z()PJ^WOx$0>)AL1mG=&to#?{Q*F+urRcl8&Z}dQDcNHfLu4(EGD+X8PWj>}if7t<| z1)!ae++*4$-mNq@E?2Sgy3fa98^4o?D4XmiyXTd=dahN>A&Y+-aQQrt$kzTXNS2wm zneq0T>gul*SxGz}W;ag-vrB$Qh~=xJPfIG;bsZf40gZap~;(3|7B_ne=yUr^@NrcPC9mq|L#1J(X08Z{QsF<%PBXqrge?6 zm*0-(SUBwAGr$Xb-SCci)|EcNfxy^M0gl?;-_xc>%!Uco5i`boHn@a120Ma6%7f!& z|B~3A$C)!z+rDz5;~VVLH|fVVc99lMFQGM6Snx-e~shI)hZbwfiFIcYh?EiF%NWtL-?$xl~- zUlJb2@i?5@fRY$p*00h=e>5nSYRL*=P+4z+)bCC?2uu7c=~v~}$}`WBsWJClnb6F+ zrtYmSX3@RCPTFn{X@gRt_=*LJs%|gIbG`0G7`>WT`X~CPjMuL;5k<+8f*%uU$91#s zil^_Tl+jwn=2b-E)V_tsmDae%JFjIWa+MFW2bbRV>8NOvB7$2Ue?(2H^ zW)QF&n054HG>(QSMe~cdplXh!QO1_JvHj`!Elt?Hy zq)en6kbop=Alr+2B9SfD+oX1v+<7rusxFOf zt|}I5VWAKNi4Cp+mj1Dx`hzcRC|^n3M$K@=lw^{{cTrQFHbnF=VZ6?=*2)jO!8P{z zOso7Fe+wI;9<4eqr}=4(N6{<4)G3HKe=j1mFLJi)sYw#rwOwRGC)J3Ee6A+D zrE$05s0HnfpYJR}^)Uq%syQ6e4COAH5k=9fbW}j(P;2;3%#uaGr(($1?RN8r|05H2 zprdwk*OdH-j3Q~D^b^;R7GC%6Vcc*UtQ}(*xu%ImtkMoz>=lXLe}|R$KI{G(0->Z= ze{W_sw}R1C5f*%8YWZGCZXD#hhDHnCD_3^$4M&rIls#vi4sSF_u?To;NBl6)-QYdt z;tIL^qr;nfVttC%Y>BaGBb+xQpxMwh8oOcfbr)WLBi{xcIE8?I4fVL0D9GyZX5J_U za$Q3Z-k5K{xM5E8?O95t-2_Lo_qkSbfAHCOy8jL7&R-l}X6VJ%?+o;w1l&*%$?-IP@_`K5L zj4v3)?dGTh*6csJOwZFpj$`sZ6{UA0*{1ffJ_=tSIu=iKv4B{7#O17{w>+zC6I-)CeW5vc!=O;eC`2C{ z<4^gXp+j+<#9~r_I|%UFtZdZZE=9}j(J)(nD_C&}JP$r~0KAxZ5$9Gzwb3{C&X(+ zdQGT2d-O-KX$%L+cl*Qe$CMT{PJ67K_F>QIO|9jEBx}fTrHr~_+=7KPsx``1Sd&>A z>G`x=wL?$V8NxBZ50}~ff1I9L*6JWkB8a0jRzgqDCwoa!t3Wh_Z*%ldSAmwK)w^Y( zX1%CS4vhV<$spn~LIr+H|+KZkd`uzeBqyX(O$@) zJ}-ch9P*fT8ib;|h(P!p`R-bt@`H1e=H*62IeAjWQgDD z5rSCrDeO>gHF0lcq~(`18yl<#|2YnTJx*sjqsfkos<2qReHS|vg$z-driG-Y2PXY`+gXS-Q*{9LW*(LR#N&lrO9=WT{CxvH+`fU-L=-3fy6E9( zP5!+^Bl|)<2Kc*T-EM4hiGy6JYUeSwnBDd)n*=JoOhGx9f08-1<)R5M;G&M}syLP2 zs4(!XHP$Uy%V8ZEjDz*>r`8!PVX}gdbX(C(85?UO#F`$u`nHO3KHrX+^~2O!Vz#jD zv5lTu&3lhSMy*96ZF?ur4<-Rwmf@TOpV$GBH{dP%D Ze&4V;0MFpe004O~h?o4T{fUzxP8bYK{c!*Q delta 4019 zcmV;k4@~g5oB_C;0kE1re_W2l!E#KW&lIki24GVdDH*!U(@pf?nJD#SxX3KloSu4T zUR7}ZsdFO-9qQikIz9^fEWWb!AI}XUu(;c6`<6rZ>8L6+pIKPLLjG({1N`( zrd(aa&nlqtiql87u20GMg=F7ptf^rI^a4cCmn3Jaz)axun!))~VBXDr`Y2h5c;tsR z%^SaS$qZKw^4MR?XFAKHScg}k{D5>igpmf%(|`jg@Hs?qeU3f`{@X^5>!}pcoxrt4hoA>h!*;0eogKfKK zEmhn>PT1=Vf1s&ghz?Nva?pq6bl{n1L`Yo*DujB3EbJ}tLF9ff-X`@ zLQ`E;c}#|}(T;p0+kXN1I>!B&AP+t+Vy^k(q8Mo26x-HPJ3%8wMDR~qrvF>GRj?_W zbQi>Pd?^km5}6?(LegScfYfcu%7`i}|IH(8)^IgFmR!eY}TeXOv! z46WOoj7M^OFrKtEp<n_0O= ze}I8gK+!whmHxsIXxAh%xFyrr(KFnAp1Wz+W{Kwed$gRl&dznq;3{xO6%K z(ge>lNNy10Jozatg0A)iq+ZJ}wlt}4SjKL0$2Mxb0qFXL`ToLQ2=9^AA`=x;e`T~t zqj1F6tO6V-AwRjf@rGF@uxr}w-?TLPOdrtG51&9jh1Huj6U{#Zfz2?)QY?m5p);N7 zaaUl%zGgwfm2E59hkS*aLIrj2DX{#nufmzDLvc&q)`i1$1(bm$t$H`aM+JP{T%RpS zU_3rBh-`>~;l?6?one2lPbIAk)X`S4#E~KzjaX|dHDhbwTDKvEB!RilC%$!1VJDiil7fK*mqi1W_8$#`ujT9e`w2!f^-_L z`|D6y!F|d&Hp!rRz1i}p2D3=6=1-0|&=4!pGJRvz%4MZm$T@g_jzib9+1o|r66L)% z1-iyb54VaIjrmBAurE2jo__h-B@Bb)qN=P)yF<8)S&R=wQ<$zNjp^r!N+ZX;=i;8l zNV3Z|rn#OXKc;Wp3648Ae~$RIFW2gh5eN^=aF#IZkwVk23LFznwA3rR~rop@CRuBQscNRaO}uNmZdAxwL+gj z-mltl%)YK&3I6V|J!wMmx#In7Rx9bt1Kk^C$8yaCV_H@p(z(={f1FJ_>sUhsM>@)S zA7j5z%BgFR;hcyN8lkMlZuS(;<^cLpk4OVuk)^FN1$u-6`t52FnwA&~*0+QT zM`HJ5Vn3B>V8tn*Q;?lN?FfuW)4HqxEC?HR4ov-g{?2)m%v^_@lys>mY85NaXi-U@ zFc?0HZ4`P0pb~E(Z%6#3IA)DGRp;ixcGL=(=#Rf8muqAGe-|9z*NXvUjwfQCc=bSY z7D@aPahTa%b%}L&7Kv%MW3f^tH3m$5;_wCYMxiH6!CFoCjKiwhX2KV=paD`-2({q)R0a=%$+EoPxu108r#FLcvSNuJSYrY^#*7Tj_fLk2_ zvi?GiA1_Kxe?}!{ER*z8bo7)8?JXQV9r!0n*YS|ZcZdEYRT}#8Yw@FGl_ZNDwd$Ug zX5F7i$NcTkdvVoVh^&5TWj1|e@?j9SUdo;;5rq;`f&l&?bOyeYOaHJA`EaUn3LPS&MQ( zo(+yjgjTx05KhoT1SqwA57JmlL0dRQ79pAmDXd9h^w9(vW%qf|P7LWki61pyIb!#y z3#BdwK@*hJa2ep>VmNBxn*{o=DE%q|ltZH>e~dw_?95xbN7GHQ9;M;Tv4z(}1W7J1 z0U4O$jWL#m7RK!R=PX`W!w4m*$7UR@Wg|vj*SE`!tXmr{k=HKJ*QX@$5x5JdnOleL ztCsrU{MamWmP7Z-G+RDqmh~*2AY(v4ron0ogsms^4P9~?7p6)-LV*p7gZd2OWm~V!Xx|+3qoi|d(~AzZXYBP}uge8e<1c)l61}_0BlLA0 z0AI@Iha|wVplbGEl4MGUQWvRfLwQ{1e|c!qR zbLQa~%}#KRT46$j*3_`aPAGdzc*uV?DL3mxWM0-gP%%hk@X^rb4mjuXCdlG(E{~L) z5R2JSM&vO&8Me0)WhH>?`b8+2+t^p=*X&316t^kOD>>PCMVf8*(6_D2oT{3nf3fhQ zIQj9wdm?If;jYk;vyR9};@FUlODW_w1)Wg9UrGFlNGF04YGrha@FX5iC#f9B;R&hpEl{dmAdPyP#y6@&@^wB%q%ZcWQI^=8J@UXb1) zw!@g~0DWaWe#?!J6!)j41tbiFf6ux*@A8TD$O?C)NB(<1d=|5qmgZK0&J*}|dU$oL z;^rIRy0FPcfR6AA>h2MsXZOrejYLhe{4Nf6DT5qqAFo{-DPeI%#pi%h`J$w6Uf)sImHTgR!POpuxLd=YqtTLMWEwfkkCo}uFd)5>as`h703=GPMM^_D zf>)?%ktCwIS_}K-Lca4cr%Ggc2j^y021q;#00T_LR2LoUoyc#QN+luWNvhN#-A6zF zJ!Ug}V}ylTdd=TTkP`wEe_Mq}b#P0+0paiRY;|P&@Kp3`+Qr>jiY(Qa_XT;j3Pf~k zpXP!G7vcr;0x&$BvanmK<)LkX$FrmU`jzml=MAbz$`t;QUMeTf?B0Gs>;yYVWu->Y z8=&-lSSY%k;me6}6wjE<^f81CVVPd1gPG()< zHFTl=cGUERgnzBpRp=2gsTCb7HS)GXZ=>_Q5VygJfIS_JAYS?$eJk!76AUu3?`tT+ z;TX6^v|0&)+pqn_3sAz+S*e>moExO^ENqvq^g5TUGj zs>p~EPF1p7>2$~;J_nT-X7^Lf-mjHlJ+9|dMLv0cAszc9naLS6J+`ZoT{Qv$l7-84 zt)MQlq>$|dYAubH3F0EmcR3dpWS3B><(3k}1O@M}W3bUHGk}9W=87%#C=@_Hj?gVK z9ir*xKwm;`f2GYthVFq+Xcu`UoXuvadR2j2hG;_7!P3lI%BhFP@Vsm8QlNz3RxE{E zzr*PQv!RXNIZ*6o4y0HlzrSCZkLaGInMjPx%dYCF`F$qMdK0Tz>}TL}R|jl07Zd ZilQO;S3%aJ005QhgQn^6r+SkhP8ef?+7SQ%