From 84af7878e522c0f59d69a746f904064d250eb567 Mon Sep 17 00:00:00 2001
From: jain651 <1Amitjn@iitk>
Date: Thu, 2 Dec 2021 23:08:57 -0700
Subject: [PATCH] improving coverage for concrete ASR related objects (ref
 #272)

---
 .../ConcreteExpansionEigenstrainBase.h        |   2 +-
 src/materials/ConcreteASREigenstrain.C        |   6 +-
 .../ConcreteExpansionEigenstrainBase.C        |  25 +++--
 .../ConcreteExpansionMicrocrackingDamage.C    |  11 ++-
 .../asr_confined_kelvin.cmp                   |  36 +++++++
 .../gold/asr_confined_strip_kelvin_out.e      | Bin 0 -> 147200 bytes
 test/tests/concrete_ASR_swelling/tests        |  29 +++++-
 .../concrete_anisotropic_eigenstrain.i        |   8 --
 .../gold/concrete_expansion_load_10_out.csv   |  13 ---
 .../gold/concrete_expansion_load_50_out.csv   |  13 +++
 .../concrete_anisotropic_eigenstrain/tests    |  93 ++++++++++++++----
 .../concrete_expansion_microcracking.i        |   1 -
 ...cking_excluding_confinement_effect_out.csv |  13 +++
 ...crocracking_no_isotropic_expansion_out.csv |  13 +++
 .../concrete_expansion_microcracking/tests    |  60 ++++++++---
 15 files changed, 241 insertions(+), 82 deletions(-)
 create mode 100644 test/tests/concrete_ASR_swelling/asr_confined_kelvin.cmp
 create mode 100644 test/tests/concrete_ASR_swelling/gold/asr_confined_strip_kelvin_out.e
 delete mode 100644 test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_10_out.csv
 create mode 100644 test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_50_out.csv
 create mode 100644 test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_excluding_confinement_effect_out.csv
 create mode 100644 test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_no_isotropic_expansion_out.csv

diff --git a/include/materials/ConcreteExpansionEigenstrainBase.h b/include/materials/ConcreteExpansionEigenstrainBase.h
index 848ec46d..e069a205 100644
--- a/include/materials/ConcreteExpansionEigenstrainBase.h
+++ b/include/materials/ConcreteExpansionEigenstrainBase.h
@@ -112,7 +112,7 @@ class ConcreteExpansionEigenstrainBase : public ComputeEigenstrainBase
    * @param sig  Stress
    * @return index
    */
-  int findNeighborIndex(Real sig);
+  unsigned int findNeighborIndex(Real sig);
 
   /**
    * Determine the term of the A or B vector in Equation 17 of Saouma and Perotti
diff --git a/src/materials/ConcreteASREigenstrain.C b/src/materials/ConcreteASREigenstrain.C
index 914cbec5..01cca11b 100644
--- a/src/materials/ConcreteASREigenstrain.C
+++ b/src/materials/ConcreteASREigenstrain.C
@@ -266,8 +266,7 @@ ConcreteASREigenstrain::computeResidual(unsigned qp, Real scalar)
     else // hydrostatic compression: retarding ASR rection
     {
       f = 1.0 + _alpha * I_sigma / (3.0 * -_f_compress);
-      if (f < 1.0)
-        mooseError("Wrong retardation for ASR latency time calculation!");
+      mooseAssert("f >= 1.0", "Wrong retardation for ASR latency time calculation!");
     }
   }
   else
@@ -301,8 +300,7 @@ ConcreteASREigenstrain::computeDerivative(unsigned qp, Real scalar)
     else // hydrostatic compression: retarding ASR rection
     {
       f = 1.0 + _alpha * I_sigma / (3.0 * -_f_compress);
-      if (f < 1.0)
-        mooseError("Wrong retardation for ASR latency time calculation!");
+      mooseAssert("f >= 1.0", "Wrong retardation for ASR latency time calculation!");
     }
   }
   else
diff --git a/src/materials/ConcreteExpansionEigenstrainBase.C b/src/materials/ConcreteExpansionEigenstrainBase.C
index b0ba7cc2..18de8564 100644
--- a/src/materials/ConcreteExpansionEigenstrainBase.C
+++ b/src/materials/ConcreteExpansionEigenstrainBase.C
@@ -51,12 +51,14 @@ ConcreteExpansionEigenstrainBase::ConcreteExpansionEigenstrainBase(
   if (_expansion_type == ExpansionType::Anisotropic)
   {
     if (!parameters.isParamSetByUser("compressive_strength"))
-      paramError("compressive_strength", "parameter is required for expansion_type = Anisotropic");
+      paramError("compressive_strength",
+                 "compressive_strength is required for expansion_type = Anisotropic");
     if (!parameters.isParamSetByUser("expansion_stress_limit"))
       paramError("expansion_stress_limit",
-                 "parameter is required for expansion_type = Anisotropic");
+                 "expansion_stress_limit is required for expansion_type = Anisotropic");
     if (!parameters.isParamSetByUser("tensile_strength"))
-      paramError("tensile_strength", "parameter is required for expansion_type = Anisotropic");
+      paramError("tensile_strength",
+                 "tensile_strength is required for expansion_type = Anisotropic");
   }
 
   // Initialize triaxial weight table
@@ -165,17 +167,15 @@ ConcreteExpansionEigenstrainBase::weight(Real sig_l, Real sig_m, Real sig_k)
   return computeW(N1, N2, N3, N4, N5, N6, a, b, sig_l, sig_m, sig_k);
 }
 
-int
+unsigned int
 ConcreteExpansionEigenstrainBase::findNeighborIndex(Real sig)
 {
   if (sig <= -_sigma_u)
     return 2;
   else if (sig > -_sigma_u && sig <= 0)
     return 1;
-  else if (sig > 0)
-    return 0;
   else
-    mooseError("Invalid value for sig.");
+    return 0;
 }
 
 Real
@@ -184,16 +184,13 @@ ConcreteExpansionEigenstrainBase::computeAB(const Real ab1,
                                             const Real ab3,
                                             const unsigned int pbound)
 {
-  Real ab;
+  mooseAssert("pbound <= 2", "pbound outside allowed range");
   if (pbound == 0)
-    ab = ab1;
+    return ab1;
   else if (pbound == 1)
-    ab = ab2;
-  else if (pbound == 2)
-    ab = ab3;
+    return ab2;
   else
-    mooseError("Invalid value for pbound. Got ", pbound, ", expected 0, 1, or 2.");
-  return ab;
+    return ab3;
 }
 
 Real
diff --git a/src/materials/ConcreteExpansionMicrocrackingDamage.C b/src/materials/ConcreteExpansionMicrocrackingDamage.C
index 66f0fcf3..b9a1500f 100644
--- a/src/materials/ConcreteExpansionMicrocrackingDamage.C
+++ b/src/materials/ConcreteExpansionMicrocrackingDamage.C
@@ -165,8 +165,11 @@ ConcreteExpansionMicrocrackingDamage::computeLinearExpansion(const RankTwoTensor
     return std::max(0.0, strain(0, 0) * _eigenstrain_factor);
 
   // otherwise we use the principal expansion directions
-  strain.symmetricEigenvalues(_eigenvalues);
-  return std::max(std::max(0.0, _eigenvalues[0] * _eigenstrain_factor),
-                  std::max(std::max(0.0, _eigenvalues[1] * _eigenstrain_factor),
-                           std::max(0.0, _eigenvalues[2] * _eigenstrain_factor)));
+  else
+  {
+    strain.symmetricEigenvalues(_eigenvalues);
+    return std::max(std::max(0.0, _eigenvalues[0] * _eigenstrain_factor),
+                    std::max(std::max(0.0, _eigenvalues[1] * _eigenstrain_factor),
+                             std::max(0.0, _eigenvalues[2] * _eigenstrain_factor)));
+  }
 }
diff --git a/test/tests/concrete_ASR_swelling/asr_confined_kelvin.cmp b/test/tests/concrete_ASR_swelling/asr_confined_kelvin.cmp
new file mode 100644
index 00000000..9ee6a292
--- /dev/null
+++ b/test/tests/concrete_ASR_swelling/asr_confined_kelvin.cmp
@@ -0,0 +1,36 @@
+COORDINATES absolute 1.e-6    # min separation not calculated
+
+TIME STEPS relative 1.e-6 floor 0.0     # min:               0 @ t1 max:             500 @ t6
+
+
+# No GLOBAL VARIABLES
+
+NODAL VARIABLES relative 1.e-5 floor 1.e-10
+	contact_pressure     relative 1.e-4 # min:               0 @ t1,n1	max:       865610.43 @ t6,n4
+	disp_x                # min:               0 @ t1,n1	max:   6.5003353e-06 @ t6,n48
+	disp_y                # min:               0 @ t1,n1	max:   1.5779584e-06 @ t6,n1
+	nodal_area_leftright  # min:               0 @ t1,n2	max:   0.00061261057 @ t1,n1
+	penetration           # min:               0 @ t1,n1	max:   1.2681453e-11 @ t6,n4
+
+ELEMENT VARIABLES relative 1.e-5 floor 1.e-10
+	stress_xx        # min:               0 @ t1,b1,e1	max:       858917.53 @ t6,b1,e1
+	stress_yy        floor 10.0 # min:               0 @ t1,b1,e1	max:        19997480 @ t2,b1,e1
+	stress_zz        # min:               0 @ t1,b1,e1	max:        18601279 @ t6,b2,e23
+	stress_xy        relative 1.e-2 floor 1.0 # min:               0 @ t1,b1,e1	max:       3462.2872 @ t6,b2,e23
+	stress_yz        # min:               0 @ t1,b1,e1	max:               0 @ t1,b1,e1
+	stress_zx        # min:               0 @ t1,b1,e1	max:               0 @ t1,b1,e1
+	ASR_ex           # min:               0 @ t1,b1,e1	max:    7.947504e-06 @ t6,b1,e2
+	ASR_vstrain      # min:               0 @ t1,b1,e1	max:   1.6310732e-08 @ t6,b1,e2
+	ASR_strain_xx    # min:               0 @ t1,b1,e1	max:   6.9259599e-09 @ t6,b1,e2
+	ASR_strain_yy    # min:               0 @ t1,b1,e1	max:   2.4583968e-09 @ t6,b1,e1
+	ASR_strain_zz    # min:               0 @ t1,b1,e1	max:   6.9263808e-09 @ t6,b1,e2
+	ASR_strain_xy    # min:               0 @ t1,b1,e1	max:   7.1325043e-13 @ t6,b1,e1
+	ASR_strain_yz    # min:               0 @ t1,b1,e1	max:               0 @ t1,b1,e1
+	ASR_strain_zx    # min:               0 @ t1,b1,e1	max:               0 @ t1,b1,e1
+	total_strain_zz  # min:               0 @ t1,b1,e1	max:   9.9985036e-05 @ t6,b1,e3
+	total_strain_xx  # min:               0 @ t1,b1,e1	max:   9.9995609e-05 @ t6,b1,e1
+	total_strain_yy  # min:               0 @ t1,b1,e1	max:   0.00052599114 @ t6,b1,e1
+
+# No NODESET VARIABLES
+
+# No SIDESET VARIABLES
diff --git a/test/tests/concrete_ASR_swelling/gold/asr_confined_strip_kelvin_out.e b/test/tests/concrete_ASR_swelling/gold/asr_confined_strip_kelvin_out.e
new file mode 100644
index 0000000000000000000000000000000000000000..c0488fccde160cfedeb108420ddcf158bf2cf915
GIT binary patch
literal 147200
zcmeHQ2YeGZ7Y}>yJql)(RWcw@mM;X7FvEsG8OE_MiNUdh?Sv$>Y$&U=tU^oKwCvHc
zGYNa|J<HzELdz=Sdr!LacI>;Dex>#Oe)j)#lHN&APkPeRxK=?S*}c5Ha>Fk_{3;=;
zS&d9@(5W<$;PUdS0>8ZYOs&$eaNh-fIU}^OIz;&8GpM3jh``2xgS$w0jXWB|WBfdt
z*l3wj6^-8bzI^D`XqBuUpLluIhhJWN#;REiR0e+esaqDV9t18LuMi(kre_UYTueK8
zu2(4`9Tv|Eb2L;MB;!Zul!pHz935E<tAlqE56p{YcubX-S2p-rg(n^uz97<}@6{s_
zOt<+n42bVDf7SrpNYnh8zny0d?L2G5JuAScUqfC?cxv{Hq!0OQ_6%R^RnGDmzLuof
zEF8YpnjXH^@>!f*htH7rBp+CAGWxI>JR@bP(nM&zymD2<F!0L*zXI|Yl`M|c=~Y?{
zIK6VgFI$VH`K5q)96T=@q1MU`AdRvZtx98%4bke9!1F|uEi41%WnTE@i%_Yd{HceT
zAK^7bY~&9@DmTQU*zob-)+?7mWl$qAxR7TOxn3tzXf+Wa?(ib02r;rjtU3;4NE>VL
zCLv2i%i~qiD8nG*GN>GphA52B)0Mz|0r<6nUw%?QhD2$(mtcBwz>mL>hnRlIN!mB?
zpN|bUB@~~+kNU$6T=)*oU3@a$L--h-HU{O74~zW9`<z@I#C+vw8wCHB@n?WLpGI%<
zzI-Gk{loW_rMJa@0L%;F&EJP%oxP8bzX8JL)7!oe+aCThf8RlPm>0sEzYoK@`o0?w
zHn&0x5))6zBTtau`~~;fkmqm}e+_>$Yyx?ZN!XH;Aq4Ur*^Bao{KYcCwDT7}3!;wF
zSO<`A=HZcV*5UVZ;Rj%N$hUcT%(sK^Wom(6e$1ya22B&zFY7p5n#%kJeg&}&mSJNO
zEyw2DI?ouFOH)~FEb>CV(69;v${OYmWgN?bFOPhr9WW%snANA;U3vrtk(P#}CEEZ9
z4_8^_H%|{6ETmTmeud#DfSy-T2oG2JzVKI&yx%fE*zWKb<_E*URbCB$c}aLZX+$ya
zMbL-4%)=rN@|v~tnS`xS9m8^TF<%&#OC#LT^5{kGyef<U<moyKdltiTba@?uZ7qLc
zdhrKWMPK;i={gHb%9W#w_n6+|@Z&FYx{cu>PuE%4<`|Zv%j>cdJW)V%x-Bp)$6se*
zgD@;dm)D0S;TeCK(`^F}dAiQRcEGT-J{aYp6i*Zo>E@;Gik<PE3(J>VY3fJ*&BJ0@
z^I`e2z<vUMVOS+~R2%@vhvoTOmim!@49m%H#hqqhd7hzu=PwM)rS~HV-gD{Y`B0t@
z2xuO*7Pt$=$<@&MXqBkbv9iO3^%M1${q{3k+HCN%h=+`U%SpU#Ht$0wIggjI1;SXo
z4;coRllSdt17R%UA@ksJ5^v{O2xAcs83~t@c)Oz@OhIhOd46(nDgYwi`3vdnrH+rV
z{uUx}h@f-5E&8sC;Aa&tJD2s3UVvXQQg(biKQ`Q|{IJdDFD%P-#F1@}ZQS~<#@)LN
zXjs0tfva)%^|FcU?`qurd2Ql0bUE%8o4Ac!je7vg)H1)Hx*GRjS(~_xU5<<LVfo%B
zuEsqS4{?iI%7!l=nU`S)vT?YUZCs;ZIg5MNHm*^yoWwmc-ZrjLu$;s_+SoR(QLvoF
zy<;2KC|FM79-D0&*C<#{;vR2f8`mgU&f>m+xR$a3a%I&9WB-c3u<qw2j%+8E0b0fd
zxpFZswjZ`~MI-ZRcJ$sYwsAqOT#LKSHZI7OYjOA5#s#@@E$&g<xFA=q#XV~q7v##d
zxWC%Q1-Wu9?j75>AXl!%{SD&Mel@5enL!>7ZY)=n9n?kqg*s_Gbx1w|EQsTOnT}P!
z=t~c7jPnuvuz$&4*}byOy+|FB4S?v!|1dm%Zx7KQ`y=>X9Os}s+Y5)!%)_A@<K_-n
zxA;tcY0Q#<wL1sB4oTm|;ZEc8HG#En6ra~JfMQ_4s?^&DE`Oe8#Fb5dT)(1EvjB=s
z1D0=A%!zBVg~QDO%lxa1Zdumfh(E2nGN{nAp96UOvEJOS=T3PC@VIzIAE3B0fZ~dP
zMK#M+9XzcThcyEhaOvf79+y6ffdO+!dYW_CG+;LOejaCX@8fYAN1tNNfGLDt5Ql*Q
z<4OF+9Qp^0;=)s`84%C;Q}hprBJLI(a(xouG<OC39Hm<)+4+N!zdZl>=TXxQDFcVS
z3gD2B$D{RgkDKx<tNgL%BTvtR9u7F*0N;k1x2xFxqOAq>0@}9FUd7voc)JeTSBk<9
z?H0x1hxUMy@WXL7j-PQ{To!&fj>YjQj*(Eupx(f_iOTSkz^@AYs=}`t{HnvR2K+vT
zUrqScf?sX;)q!7K_@OSU2S0E4`M}Q?emEpT?Pr}2)C;^W!Mx!;pC>+leBStc@_FR*
z%jdNq#Kn6)&uG`edp_?xA9x<{{NQ=P^M&UL&mW#QJfC<T@%-X>#q*8l8P7kScM?bo
z-t#=<`N{K==PS=sp1(YAc|P+z=K0O@8u{*w-#o8*9`n5AdCK#W=ONEKo@YF-kVn>K
z=}y`x3-9aqzS50!-`5ZOTquuNUy%3cMnpGX{#f={$B=L6kM#%Z4Av8@8(1H(4xrql
ztRp>?X_RA>U6fY|{7^1Y7E%6C=1|U1wy-{1*JHlEq8y<N@%0bo24x221=ERj2<r{j
z6|5gs;D_}9Wgq1oWgO)i^MX9Y^pkd@JKSV?fq~p0R~TT%g4OF|b*xQbrAi+oi??}f
zdO!jsnA|pZV3w1sWpW)Whnc4cn4gM_GFSwPVKuBl2Xkh0lFh=8a|CwlDr4ilBJgUV
zZ)uAd0Or2rFzXd>H=KEZ1VYG)j#+>N62Qto!aP7CP_XuIVHO|}(qQASHxH2LW$WJ(
z1H{p6*o)sX0{ozyXcq0nk1z&+9N6{e13(rXT-+-X5*_$^aRDHU4lX``Lkb410VD!~
zam8gr$0C5eq*??>boe57v?n%+Gkf3?{c3hD`(<NCrtHajHAyR|Oa@jd2~UuOcIwnE
zSW?YiI7fBG=DD2&R<l#bY~ndgFA$Wp(nd$iHA+bbn3t9W>LO#KS&cz20Wa6-7-d-E
z(?bu_@IC`UV&pzBy*dzP_4GdBYPn)iI4jqg?&@h{43c^(i&=4Rl_W4WzNcKLl838V
zy|=fwmrr*el}2SyfgC~;sDwMMM%uu)vA17CNht6Zf~xgCSRdK=7`X-)0QeZ#Xjp#%
zQ}Q}i77LN3?a4H|Boqlch@)d8SRJcTurjmQwl6Y7vASqeoj`(+EQO9Wus&RMiqNsb
zh`<sTR?Wpmf-*-GI>j1dVu71vK|;7(&q_00QXr{UPoqVG^&*ww^}?Y9pdu?+Ec2dp
zsf460m{PE$&8mtkw;?*1c7XK~g9rg&2$o7yS7JH{Ln{ZDs-_at6{#{v6d<4g{&L^&
zhK(Ds5uZx58cBWM`hNBN>eutdu7K^;K@}d#>ZACKnn>CPb?r&#;O^ZNw&Wy|1|W&v
z*k2HMb?Xq^*{uXy#p*zS4U+D#Y|4~hU}cC-B?)2`kW`6ZgQmW|P5qk?_Bl&%a6GGs
z1-qP@B{X4?OLPn&Y(~)6uc=@CrVSgoM9|p_|Ns72H`;2Tl1V~f(T`cp2nZ_vv|jnM
zDgXB*_UY&5BuFICQY#FqIM!6_EbgR|S{kiP$wtUy)duXqI=PS_ZLlMB+GrVRtZ=Ck
zuH++xrIHST-GjS1?eg;Pu}2V0#->C{C2c$Q4D7(IL*k=3Mu|kOR%?gIG;AcaTyd;y
z@DLU|yLz6aR3fZGOC*X|y#X>CuGT6B8OiYEBog5~<OqU|8=7hpA^uJ{L7cdS=0Qd%
z@?=aXl?3+a-pMh~tbHUB2ZVT%!U>Xf8jh(7Cy48X_~h6@sibvak8a(pbLP<P_NwRb
zO=3OzK<}wNtJAP*y}O$Uu}QLuCz15==6Y-HNg<YNrVW_zngRs%ox-u>8VKzo^ffFm
z;f;etlBpz0jUCyptXvI6XvH&YCnkq9G?wWyxx}O}TpO!V%5@1QDrR@W8(3%&qgcft
z8BZ4G>!5qAQs7_+E-@)Y3y4muhMB@h+XfLD2czg)^UIza6Txao-v)1OFEJ@(Nr#D*
z*@z61mrA;~4F#<i+*uafG00w;ET3VV2vE#gGew8bjF*@cYQgZp8jK_?Z=@15MY(=q
zY+OJqi)0Nj5R>U+BO+AsW;DfAJtWBL436ZG4)hG5$XrP$rbW+Ekuyt|=~hgZWA!X-
zInc+zkQUZ)L%$}<L=Q??I1iyOL&m#iWK6dp*wm!vN`4^>Hn844?rB?He5Jr(MipSR
zD`hdY$gB@u1$32b26Z2zb>8|DZex}`Dge1)y8fLiK&4T#@iLVuo91_7S{q7mpg9@S
zy%>Tqa)TlY_Ry#kG}>si)0ubS#q=E9u2XbSw`7JPS}^Tt^w#=RN-U5@<Dt{W44LJ2
zouY@)8TO52gXs|2R3dOEhKK(TC1q^SUEtvdr=$`*N=hp6C@CkRjrI0FI8OPym6TBH
z9E?*uO3F&72p^~XtCW<(X?0@_|L{gq2_7W{oh6TwGK#SM8;_E*XbOo#j4{p0NXnz6
z?2DiADKg$OlCqRVYVPAkQuY>98M#MqK0qav=ov|Ql$4bi+d8pYct%pzje~HV;zUW=
zUp-}89DHlz8A(|d30mj=+aoE*>yM36^x-L~#CVUA@+c{z06V(zC@JgKP+WWHQBszC
z^vp>)s&U5W$S@ZcPH-bB$9rJR2!3dDQi-0Cl*im>lwXH8SW`SBDdd9bl1hXRd^}3Z
z#6P_AC@F_kj`u(|2~biFw@32yygs0j6t_J#i{mI<(D*M{-N)}#^epf;in8Smmhyj4
zNqH7{W5Js)I{gu+Rl_D^+%CZFJ-3t<GZ)5m`xMguW9=M&n1ET{9Nc~&D!^jne|*-=
zX;{jh1>VR7(<Nq6D8CuUM8)h*3_;HxA2SBRt(ZMNbXeh8;LUSV%%&kH3%t$SA_uLi
z*<(zI+jYtZzjn@?hj^?_=GTX)q_Vhr4#mc!q>u}yON>tOC@B;F@J`IiG`DX=b)ux)
zzAQmF!4H0cw|N7bdE<({C>b6lg<LRQQi(YaWlmSmeTYgbi*r86ERT{hi45L(l$6uv
z2HEt_KP|zDlJd;>I4Nf=WzR_Jy`Se`p`;uiiOCg~+m+OZaf(q-{tU=0&mTDvdQsk3
zC~iN+C_)~sQrkcG1eH{1r;eRM+Xi-UEx}lq-LL4a58}kP@Rt}h>h^P%%#R-66J|uI
zq*bSm-MR;Mbn)aljHp-0;q$bSMj943VzTV)Bq*&`&+@|G5{~%Hw74-fnI5(${HQ&9
z1YruB&$pN?yZxXNI1n^Ot0AY<^64_6nDMgnnMb}z^qH0i@x-LSae}P%xnag6idix4
z_TE=SoH0GNH)2u<->}8QX>Ez&6ER{#rr~1cYP%GLr|?5t-Iv8VGioC6=oCxo5z{z$
zbP9a_+Ux*Vk4|x_15oHZI>ox`3EwB>(J4+vEz<zTaSDGllRFNG@aPn4!4^JS;n67`
zog(;*agR>1tav!A@aU9G6~KG&wALwETsGkur?~5kagR>%=oF7man|CY3h~&!oN?)I
z_P5a~?%3Jt8K+owmpryF&wPq)y+#@F%%^0%rg(IUt!({YpIGY@w;bu?(J5K4DV}kP
zXPn|0r#NdHv7UK!iZd?#&Hgqz#WPM3wY0;t?#o(@hI-b0WnPz&VV-e{wb%%sdGzR%
z%;oHF_ciJicOIDI8K-2urg+9FAKo}c=y$+9I>oY{;FP*Yr#R68*kXSGpIGhg#;r3b
zB~Elohrrfd13R`3cEUW)8wvj1`%@bDHum;w=vw-?bLZUs9dIW)C8TZFZZ2+4fovcT
zJ^NEYNxk3wDIT5Tq{^TidiJLv7fhG%m1!S@PO;d#n#CO|$Sfy1#k1$pg{tR7r+D@}
zy1+x*O_Uy;67Sjbn3?2Y^X1u}k{O4b`w`ohg-*%hYGRL0QE9FB9g7*!dvuD&_9em)
z<k|Cx4Up*)J__*Y6z5ta^ANT=#SLesdB!Owp~gFpPWf-FCHC09GH-0WurhdTU&sa1
zC0y&d+4l7T9C6mNmEMihGfsq5j5krE@_3b8EmLdd%B&|A<-^I@&i;ZDiUgf>IE6P~
zN|ioFEmyG7_D+|QN@}9NEM7t&p{A4cVhY{fII!^som>_zhtH$L%M9{JixZc`oU=p;
z#<2>6R_7Lis^}QCN})1f_H8aPJqPzs6uNs}&p~m}+)E|;1igWcb|NiCZ;a*kvz?us
zJJ*qojh1UPEn68=k=YHV@bCUg4y|&TrOMn7Q|LgX(F=R98m$ITH8PKt$vwtVs^pO}
zy*!SUsjQ7GNKkld6Qk3{$Rb&Tfz`?Mu@Mofczg|aiTMzQqr?%CVp13<SI4q^imXv8
zsaKz^>-%XYY_stCKo(BK%TU&`uJJJM$TBBdPxyKpkJj^OJxMK%R;FYl<gsc)ZHroF
z*MV~A_Eq;Eyw(d+!8hij)U3fJy@x#fzi2&YU)wR3IHvGJ(0Z(%wck4wXIc8YwVu#J
zhul28H5RwKZpc_K;qTRY4Fn|~1>2+bJX+7A^_*23tSL^$dLM|^OK_G|-XF8<(R#=Q
z)Ab(4djDsw_W>AnTY*JEt}wV~<j706JCP`QH5&o8!^kLiq!R@WEW{DejW)o_6TH>x
z6mVR$R?WhvWE8AWvj|5}rBUdJ4HhWs)vA~*6O3Us@Ih%gCen-+zD-O`QOGn}U9`ia
z>HWYc(#Yj<!$c^+@T!4N1{lEJj&6%fjLB3Xk1@pR%nQPdr0@wE?0xBUU?P^W%1BnG
z1iMv?*=I4t5JVe`TP^J%dC{?IQp9F?<!+^tK`IRk16Gy7rGmj|yj05+S{7_`3Kf)(
zk+PU9yWftj0n1^G7BbCp5@|LYAVzGW<uNg=QWm3AMay*wvPgN1OdAowjn>4}lxWM(
zDH+*EZbT8(Dp}hq!qJr)%S_GcUHnu>v{q-3>tv8|omQ@hvQFU;_>@PKfmOOBHC`B-
zHVra$w9L-t95IE$5!6R(wT39DbvRIi<^md5xWw?#hF~;nh|(&NTviCG5?L9vvvRe;
z;wwgCvTTP8O&v`kG_ZcIf6fi#U~46VkH)~CPHTXM%zT}sn7YhW5mTzn@1zpfn|~B6
zjLX+KJqrU}CDy)ya)mZrCC8}-F+4O9Lh9C+RO0eRI3oEN-$>A^T_Ci@$e~H+*{LmJ
zvTR1s)C1wO4yEkv>+9=szA>JG)#$Oypkvi6e2j@FMNy2{IQgiJ&?zK0Tb4e6K6Y?C
ztBA$g%j(?H&xZbMtU(K{0;zT?q-@r)fr1t;hncE*Bq-E$pQ+(qObRnK5Fo+s(0SKm
zf5#=Y2RsZ1Eh>YGqSvZndNe|=ve41OYl>W<U}FfyXnDL0CPuMWqE%@OG8J@rXeUnq
zRY;))!5Ec7u2YJm#kGe-Y~;!~UJ+Zi4ZuU8br}hg4ywI=0~U%<f_Pj`Ec6Wd%}qL1
zuTsJXwPfKe^p!zZ!voeJ3s>lYxbx9HjsYk^u3rP?P5U5R)`iy;H0i>CSq7ijV`X|)
zu2V$G<Z88c2$Z&39v!Zf%b>Dj@?6sP)JAMgPz|cUm_p$^)EIPHwavO1BSE2cSrMy;
zVwTCn^*9c-5*kvM*?J9~7Y1K1A<L9tOpFaJ@<4>xP=cWH9H-E?F`odZh|4C($T21c
zHfq~hP2h^`nZE@I>gD0ELBqURFbX$zZe@xnRx!vWsRqg>f?TtZ1zIGbK+rT$v{x)p
zQ=li}>vAz8L1+i<k)qFlAX~SA%k(hIhG;vBSd4Fg8pqc2bA3UclZ9?#loYE`lFyDA
z(`0!gJj<+32P@n_N!ic4@aD^0{_Ip}Orfy4Y(Y>Rsnw|rQPHMsOJNv;dds{vWOffG
z3^T<^Epu2sj79?CRnMFtYG#$?7az06Ly$=K613BXkwOvFe5oBa)}&;Fm$DjGZ`yFq
zyp#pB|D^e{-%uMV(k#0^B9%zY^s$kO*^N}fEj-Cg_{{t)2+CsQib3*7i?I|4Hfzwk
zW|>2}5(#oJK7x&w#l|SX`el|lax0bi*KgwAq|vAKo47RI1A-9|(A*@%*pQ;6rIH~r
zW;t?t3uYiv%CN`~nF$6av_PN%f=enCitFt^$ue}V!BAz0jR9n1jCGjLy70QJp^~js
z#W_^;bp?nbWSJ7g<k_PrI)(N+N|<fof(<EjlwF=+T|Xge6bO=b6~_`bMS>}Gofyg#
z;UUU)M63q7<+c|Q8v%l2BxpBzKpzP@m1eU5gdn!7ZY^ay4Ye`vq?84j`>1^r&5<x1
zCB_0M^@s>!BA123LKg#O+bmO=-$IsSqh*mgmC{~`P#$JZ!OE8sv|Ed5u7z9Yz)BB!
zlt<zmoE}DeGBP|h8lZ*CVWvkLpx6i<Ns$&_rj0=hiyb~@ksxGQ9&OQf!)A?oF)*)&
zO*-C0CBS%>a}I({9Fft`%?HCCHmcs<-d+YtJ-r5=sP$&S!?oIIFgHZQvTqm$g7w%C
z1u>$uYNfNiA6S<yvTQr}ffSNN!2)6KEt2}aelVa5S1G_&`Kk3Q(27mW7IylYL~ck;
z5!7OH1RaebMHC86p`E@amqI~}4f2o@bkt)KQjf6~8VNdTaR?6u^$&rPJf6-hM}X}?
zt0Orz{s~VZAqW-_ID$rn<6NuSBHVrmA`5mNhhtLVjRO!=LHhJKMZsIEpkD_Sg=AY4
zLG0VWh#rRZ;dU(TLd`6^i7Ha3g$|PqCxr{P9!v_B4UDmb@kT1)cRkf(V`1%sRf})q
zgxQoqu&jZNM>8kRFi?wyiKcKKV(S$h0WBKNG1%jw>l%;N6l%+L*e8}sI=2hy*xn{<
zjvhdkVPBZRU};`A<yeGH(d#MS?DaNDp}+uKP+bOsnb}1|6bcb*+-GE#QPZs$f~;Jb
zfI|VW>?+x;Ph~<k#6*ycFf6PtSxlk5X*i*Ug4~0NWt~c^lz~ymsyP%PHs%CLjx9+F
zZyd}BT9TqAEsh}H!R8iHnS=&g8==-YSo@^W%lP_YPK(?NFJ(dy)O<9Vn9&<xT3rEa
zF0?wz0O3l?oS;p6jae3|0;pwxS_$L^BGiPX;cO&YK80r)+f`U$1Z%Xwlx?d`g%gZp
zb>Uj;S!Gk%2qpLd%-*(EYs2O0&e&6M%l>y_x_7*c3`<$)mhCtgGQxToSkVX@-xaWG
zEE3j#k+pB2rjej9HTnS*N+_rj=xM6mQX-DDe^#E8WvhYRy2_ODXs{eYdEoL;11ov~
z0M@ZgXX%BvHt%73ox`3NCSNLn)y9Miz?e{#0-9eP7*;`Cg9wwX`HgX~Qz)|@C!+^m
z!j%*d#2P}Q+x`o!DM-P|R#~Ao2QncH8XlGn*c^)oZ<9(ocIxQjf#t|gOd(ELJDd=f
zN{n^Mn*Vm)NRX3hl}5?N+Y}C1+NsDBZ@V-s9uk7ALM;bVhs}m^^b{l5Za>0;A1Yyj
zY%Vck1J(gmG}!2j%;Hv$6XKvoy@j;d7>QcHaIGNpI}Iw!A-I^DlBqe3oJHX3K?G3V
zm_kgNTaOMPV|rX2qBnr86vmQzSk@2r3|{LA97#zeMuG|2SWTo}7Ohprs$<QlV!^mH
zJHut!?Pp@ReKRi-1lu4dhgeKzm(cFq2Xa~h>SZ=Q2G@&eHKH{#ZN-pV|0ca0tb#oq
zU^$knG#1UYnFF&-x;8AVc~-)TI~iGX3T7~f$7&)CQP5@!ZN9)jj1wRBw|@XRwqAYL
z?XZ}_|6<P}3Pvdk>={FYsE}bx7)PzVu$4#`nE=Jo(BIXWIOL(ilAxtc2v&?Y_HJS(
z7fx;oL0FCljjUb;D=v+~f*yEDB@MiBACjwwvtS%t-)Bw`Q~=a$SZa&=UWByPn0THb
zh#}tEUc$thCkWN)pWU5vd+*EChS+oup-|*{ovHbCA~ckTZtu<+TYn;jAT(-KwDI&W
zsid*s-Eby?c7{MwNuWlh*Fq;8riPqmm}c4S9ZuYFrBNi<X%#@urIIEtHw=7D;R#x{
z3S6aun&Kx^1vr{8BP5}5aJ~KK>M{~EvyB+bh$D!5E9y6LX_+oZkn3X`*K@-ZZiJT0
zF}5}?l&vWU#<NQh*ZYx$e*#IG5&UnB0wP%GcHkCXHX6o0a@<8~Diq<L(89?~(44F+
zDmGfBR9SBg6Kd4OVFgFfc4`5up3un!vvG=b$ICSSab5Qzt8`IjYLn+3k75Md?E@d|
zEnqv%wPKdj)-bu?-2g6A<;wgHb7n_a_zYvAxc*Hue3cbJS**sWO(;>Rq&=&SbKIx4
zPA?@00~F|0W7k+mr(8|hcD;NlX>)>Th=hIH(9D?VVO!(!(AJb?nnIhd9KH>^b9mOC
ztE)IMj+=+KZr*EVT}WAW6Du)q8t&dEi$hUJ4)X9{GaBt4#o~ztp;{QuK$Gr(+nD7*
z{Bv}Pv3<34L=2ycAsFPB2}P(WuyeuHjtDJfd1R!Hjl>!CXt~~c23Jf9O$1|br<G+l
zHfsbeX%R+bjUaUDOqz}|Rt!NWa{)E|B=xNqv^fw$uIjmeE`XdUVy}or${{@_+NyMU
z(G<f&Zi5|9!V*P}m&qIvY@<Ulv|NTFE`|T}0k{Gde5qkkLNwRXk)^+}U?hRQC$FMW
zJ^T8)Fr>l~P{1x6ScIs?P2e_NS)}x-t43^2u$68%5omCkpry1JU1G$>?c-&;g#i-+
z<R}j@a%d#T6EWr&Zlscy@<gx_Ic{GpeK4-=fzD<WfrqYkya;8C`k$Oaf<@?fVuIx;
zk1`U574-$IP4>cAI@U-K&Y>~R`||{aPX{AGeGIEm!ErvaC{}K<O%qC4(4HJVg%$)Y
z(kW1~d4dMoSHt`&W6|!+GZgwZmXhnll^H=%38e&Wn?Za|E`d4r!bol(-nxG@$`^<n
zHm{)`F>iImWZ7N`8;NFxpn1Ybf#4sk73tKmTlc_@E{?-7j++OM63Gx;WdsX#Wn8a`
z9O8g)!b^CTKZwyNbQo~*2qlIen2H{jBAfmNqfv5#L|fi|A8ce$chDPJt7Y?;<^SN!
zV|c835WiXiSzvL&Hy2#!x0>_<Mt&c|ivQd+oKaHA1U(!i08^#B#ZP$i<>b7Ce?+I9
zp1@)ueI|z(Is5>+#4u>o!~R5>g?8d{jbpv~t{axTFq+$!X}P_#q=#;>&QWtCrj*^@
zbFf@L#uWw+Xc+&r0h|Zu;(Xx;vZ~6g=Kz})SheVG184$dm8a*x3k*duEw&R2r+}7$
zD-^d|I6XZF>rxi69^38KIZw~QqEy7_lz&0bL5SrB`Q_<3SRNe>y)V<K#)UiNp_872
zpLY}NIk-5U`apUP+}e&T+Z&L@$ish2&w*}>>PA++WFkJ8<^Q0~;lF4ipgXJVT1>$p
zX=lHo6S4UZ&g*>udfDCkq})6J{~+wxoKDH|Nc7*PQ$p})WfpV_Z1CYs+4dC7|3jzz
z14g6&MV;dI!Q(%EX@^GzVD&a#!pA8egbLtitU(n(W0$q*Ok~!-Wmv(js}1INz&hCS
ziCOjxE07DOOUxt}3~T>$Gmn{X0OTvCaYI%nqZB%#cl!oFw=Wm;=oC{e!8_qP1y)tj
zJw6dS*3MooF~@+sw^da}$px(fAKIf?hsFU;oN#?*Eyi;5%;Sfsm$NiPLmqncGIGIm
z{W~WLJwr273FDoZ_7FY{$sg`vZ)3hkFK5~jp<NNG=&g5Lab@$5I2Osx=01;JmRWrg
zR(RjW&3gGmIH3`yLbI|0$hs!JhhsA{?Kg;NU^DH^Awd_XZMnQ-&qkq+RkC=hCv4^d
z8%AW57FK8TUt|$3hukQ|l8i~O2q#Du7|B%Jhj>ktjP0S&4KI|SqZ1#2hk^;hQEueh
zwDuF31%mbyDV!iN5a?jnj0}#bhr)*K&ak@`zPP7>-DbooXf!2}km7~?gxq2+tFr(_
z@X(eZo_vCPYRDOt!n5r5&x6tr{mQgDC2RFD10spS1)H3851aj1a!!Jdg--<-$fs6i
z@(A1&0xPNEBS|h73@A}NJl{ygScmy4BRqw6c2)Dl5T%k3I6NX0nqljv)V>?Wq8b&R
zj?HJyoERl2bTte4E+&wJZ_$Lqp=5ls#S<3Z9)`!N;HW|%%heS4>Lz>=)u^3?^H2fj
zVCb~)VF+WvK+)2VKcZ+07i?2qHdbC*myJ~$h`O!@a3e$}hqsDRQc9KC=lS@eFrkpt
zp&&Wg6(C0I;TJVb39`A9N_=|gS)JZzAXxR~KCtL)AbcEM?-Q<u!-T_OWw?(4zK}!y
z^!Sw=A7CNfQI6{)^g~!ReBLb52c-t*1jyH1!SbYv@pg2?gjI#%8dx9#Vq<{szv|(;
zP3FlIE;e{l6@0E6c6wToqGX{^1l2v*8w>{r=v59zDOf0i=y}1~J~&!Xt;P?m+Wp|1
zR1y>%64;}Ii`(@mLE50&ZLlGQLO<ANzbuU90^ipYlxNH`7A|)tKTaYUSU5w`c)?hY
zh2!w8gu}6uQvp<Q$c99*8X0yN$rnZlQK41uKR62TWZAKN{oj20gyJ#D>_7G_epFI@
zynX0)%}^E&HPxsT|FG0z0|03Zh0V9D&K=`s(xOP@%1pk=C8#+B9^#eBgs^a}hXi2-
zA{^nWg##GyJTmx-MU<8wOu`{-@&vV3XPl^ZBH>65%>;7z^qm}z@3Z+>w8#{ypox~#
zE#OuPRrsemFoo7uSTmyHQV2#}JoXR{@{Tpc#9AE?A#}WqDU2lFJi)V3p^JqNOqaM6
zW@2g)W5*_z5gj=jNpAB&dT}X?#=ZXVMG=mNaO_UBN=>!BdkC5hO~ny}!o>qHKo&Ld
z87=aGdOG~Z?p&gix?$-}i_?N2IYrrcaHbnl$P*;rg(ZdpAw)AtA)JOqi%L*0)($~S
z@fMOg<R&GEr2%_Uaoikhw1-eHF`w$dJv4C4C7yi6Z`7v5@LyaC?Kg#lsz`!hi-as&
zNDs$Kkp%5DGhD8vl5U}a9XeP_Zze7pEt6?~aVf-pI)0rNBo7YZ!4LD|iL>xkEl#cl
zXi}O&Bnb8+Scwa#4-umae3nuhu9Bn8K>#(Q9A@%eE|CN?{ic-Ap-Lu2Cc=yI5V}kJ
z$UdtbCQjNMcbA0b*Nh<P*nnQKR?@;G(GEer8MM+^F4A>@Ak|nd((^JBG}aT_evTkf
z<Ow?La=Mx_Uc!t79en~yY$+=ez@nst&#yr-!)JZ;vKXyOV}KJs#ir0c!TRC}I_zbO
z>tkoq%NCnu%MJ}+RDj|MI_S2G>$GRmZ3lujD<>>Ja$tFiD4#ht9u%J}4vd#@G>btV
zX&eWL&9Ze9!@mf4SQ9ifo}V2b&#$seNZ{`h6ZyM%ct$=fPQoS_pC`s!=76tJE+d~3
zhXV_-S#~(E5I3yIWMCmSg$|l&@q%ql&^WLVn?j48oOwMIH>@D%NSa0}e1{fAz%41{
z35pq5@N#I8Wg|hJvye_#`(Snnnhz|*X4&DuLfo(-lYxcU6gnJOh$iT+fdyYE+!xh(
zrwK);Fa}OlBx)Wf4q>X0D0-kPnxLbh56s@S?bx=v%U^or^F~FFd`sPOg`e0IIvh)h
zCFr)X6rW|cjiq>k4#!fiHw#oT8VR~<EJX>rc`=#GO5Bxzd<uQCw5r?eO*{`BTC&A<
z+JT^FXomaZa2b<cfXaiy;S1&O9GVH=9gL~8Ka%nc&CF_;!&|<T&AWvBGlF7OfYQ;7
z%6k}^*^JB_W*0s9?0JS}JSW8WZ9GFW({e<yHN|3_VqOnDL$gflv3Xs#;Gt(|*0E<`
zhqggZxkm%=tpJXs+%o55PS9pSD0+&mDSFS)jOVRr9(smmxTaH~)#<<{OV|Ipyi}Jm
z-;F~vB<MzRbP<%dPgbYO%o>glU@=}NvLMENYrP~Ac?^E^T&9FI2*&zmb|c(GU`CKE
zBe49yi-@|ch+-9kc*B%g(zsjUDO6}-5tc0>{)zAbI858Z$KLJo4~dFNAz37EoIv4|
z#-Rv;`Zy~Z7Ea*=mD-4Cd5i@SbBAz(Y&@%g@07#0g3aj|?_nf`_QSS;9VFe`h6Z=*
z9^6?L+%d@EG!9H}#ltF2M;ydLak(iZO@u^(%lL^Q(wHcV8~Pn^5Cz;vqr|mlus>!H
z`A90xXq&Y&d{-ng+RcR$se$jVX)V^vnI}<rp+v#{EprN1_rm)&4#!ZGvhYbRn3raC
zakh)*xReSn<$>_+0jq49JB9Pm?K>Ahm|#()Ocg<`;$+K*N@J9C;e`U<n2wM|>f|v|
z<~cRq3+JJUAnf*wvyq@I6NF{8(dJ*c2ZFw?PH<u6gUKUfmPrc9XBo{XN+qtJdcqNe
z?O+C=pusOWK#_pWJrY_09V;jDBqGb^1Y!DGrHqxUjfE_t`HCR~=`JByF^;Lr+SlGG
zL5)`9DA&R=jWc@4LsB$qZKP5D_#5F>PamZn!pF9Q!fOg>`B<{9I|}x+YISCtX9#%_
zDNN8Pi13+8);yuY+d~7ZaIXr8GYm3Mk@dasEXOP5)<n$RVhCo^nG@6dijTHggPoNW
zQg!dff3mZp7437nW%z@+lWD_o+aYcol}GDkde$HdPmsy^aig*G<5E75V5X{GxL^Z8
zeJ7H{K$<3ehv^d5)=<N^{hfatP559R<*-W`rUcEV;B33I!dn{{oWaiEI5{lihi}Wm
z?s;9bc{3rrrZ^exWs#umXirp@>39@0o&k2Kn+llRiAf>$k96>TZSw_bMzT@~H=-6$
z%5Fcq569myy=;FpC5DGc(CS0e#vF?nfa6$-4$k2CfM{4K0?R^qsRV%(Q}uX)Dvj0t
z9!zJJ30e;#Me&f!GHiRd%x+dwXm@yFZ+P1>xOLl(k`8S<1_ySPbnMiTy9>buje%`D
z1Unx5x+u)%!7x{Rp}_7n1MEUq!1_w9PR}vZUJ5fGfSU-yX%ku<jMI2p_%bMEp(A<|
zLC$<{CWM5b&=I{gL5uWdF@*_GqoYX=TL;VV=*Fh32)5paAQlWtHnEDBm!;!-<U3!=
z+&W#$!WLThoD!sobiJgJps}xk4(m&}nFR?N$d`SsX%V+n65O$8+pe8Dh6cOzabAp+
zsa-@&Ym=!NLyXv{RpDf)XKx@Sw2-|p%iJK3ABoD82G}<gr-Bv&JRya?LPp^B2wEsD
zF+8+RA?ZNyde997&HG-MWg!*^+CRioF?8A?mV}TOcmWRu867Z=h8=<g;gim4*f6V6
zSeROXAY0GR^<_qcpm8*Bub>4Mvjc*rf)-OKCVPW1g=PgUkfQw*nw5{3QU?2+S_Nxq
zVFRT>6@$0<8twwE=j+{AQs1HB+N$x!Xgon)-C2wTpiHcyWbWh$q6N|ZxiliP?ARnJ
zW~|5AIBAvUSZ9PRq|himMq+rA1%g)UlOvmjLeay~KCF!eS`5J$9FAI_QN`yQPLIi4
zPON+}9<j{9r)upUNhYx7c^L`X)llForcij>A@sX!rv*e2wCoDY;j6@aUshxaGqb{q
zNuge$!kxPju#eYv$R~<mW=2ynrHocs+8Dyru;n6taXbVoES!enKyE~(5KRq^hT!4}
z+Li|jJqWh&BP{f2CkEYYZE@xZ3h#;#f|x)vWk-n76xPycaU-`JR9J1RmcWh^<Cu+4
zjC71Qp!FiatjG@0>IO-$Un-Hu#!F0%F&&e@Ac<3p7!eZ(Gb!rWXq(m$c_>`#5rVMJ
z!LIpPOidvKVY8J&rRNW>Fy>fPDHDQr`*J`RvJBlh*IR8&3NC=+7bo+l&=`{_Cd+it
zXsPsM7OkQHs+baodsEO;04*YwZeC#I<@Euq2<i;-r`6D{AE52sG9O^4ii>|lW-Gdd
zc6M9JP%=LF$pv7D0G48wo!GIh<iH7WN)T3yz$|d2P8%Cztp>y>DYqNOBVyH#XX3=D
zfI)ibfWR~gX5Z$L!hx~an)9`ZRc6jaBoG}D0j74lCrD!RH6$i;5|+M5P#X~eEjFC^
zWl!^oDHJ0?uoK`5EG{twV-)1DS_?WxhwxJ_aadu{L+2(2jilCMC`yf55QL!sbUcg+
zbOS+HZ4{eDg0u(gKyJhoia{3(hipeDfY?B#wi@+$)Tm`Wv*YU1S*L|L)@U~>sZ3QA
zHV#7j*qP3Mlm4Ndj5&S;)<cVmVX7kB9dz15J9X+7Ea7<L%wz>!b_WfikzhzLN6T@t
z%zV)39Ilh=608S}d{Tw;kS7R7;o4j5CY&I*YJ%=L<<n(C;dNQZM#9=6FdtINdUD!!
zRtS<05SR#9-w97)rX3M6g00ZPksl7VDAbg$pP5O8(_1`)lZ;%gzAj<Rx$yQd)1Hl(
zs>kYda7GRs5Xv8u2Hg@(BCFHlia|Ag9)sJj;hL!cXZ7&x9a!UlrwPK54p!rO;iYUM
z$S-=c(9KyPNC%M?<dBE1WtA^}nAOwj6fE$|-a-*5YvH9F5hn{*=!xwxQl~P2GR2(-
zW@01M=Aep%&L%6vaz(Pn!sbW{PvIa|r-38O6<QcJ%Ok-&Vv#_|vd|RrC&V})nFWGM
z7Wb8~8ifTLg%czTYBewrfP=;1&|2sm!wM9Wy+t?=9n&JkBAlS{B>8Z;PDdtwVBv&3
zlJA6xAgE*GRN7d*OaqzbSDfnD!Le{gDp+xZ6Xfa$neef*gF$%J(<OqyV55zv!x<}x
z@kT1?%~w;GkUL4<Z226YV$NJAJj>+kB*;IyY{IIdh~OdVvKw<~b0b`^p}p4V!IB0G
z`|L@#@GPUT&ty`~athy@?W={pWdvRx{AcjNr7G$P<l%sa10D|i12_;6RQzeZ@@G@5
zlj7*PbvTdP!vPNmJRI<Fz{3F#2Rt0`aKOU>4+lIP@NmGxfqw%BTD%xYvf;6>{xj^W
zd^rTxnO3FgPa`W?%UEEY*N2l5^FFnH?s30g4loxkZ8%$F%3<dGY<5Q6m?waJ4^)|5
zjJa@O;>m{JRb<X}uU+N(nLL0kJLf4hhdEo+E9PokcIIr6*MDAl`h+=sVDGKa*hz@J
zx-8npoceOwp4g(XfYWn4EPjAF*{WZuu&@3?EZ&fPz#M;9v*oshD-m1FZkLlemh^kf
zpy_i6-uT3yIU1Kdf5?F#f?o{QFo)R!Kg?|OhT!;JSD6DB=fA2KYar--{T#EeyMNfh
zGid}5eR7D|ePHj#L1}p@E?Ca&Y}Kzs{s9vy$}TV&@V@0G9#br}pGgDyUr$}hVQpqB
z(tGe5#m7H08!*01MzQ}MW(|g)E2B8;V`d5Z59&x!8_LW;x3)jUfvc=bkGU}9`^z=D
zR|CFptW<3MDd6)IW$lqQm<zGFnP#QcX-NYbb}rJuFD<F>s^_K47qpI_mejvUCFP=6
z{&`wbzjsHjr))5DKe;mb)^IcTpJQ)aI&SV>-HV}(BRquoe@@sGbk5xWt1tQvGIvjj
z@6+DG{c*}|3wN#d#;H!EH>K{%8y5a`OK!7t=jv|hE_&9|U9g9xyV&oR?$X~{x(jW$
za+kAo7cOk+&Y#=p9&&wd>sY@Z)8ZC9UVCl-SpGSWj_j$9p7ZfM`1f)kb=B|%US7H1
z{k?DXye%NI_O+I==Yu$GDH2<8@pufha`Vr5`n=m2^5Koh?A%<qX|Z!7I2`T-9fcca
zMTd9u$kVlUJA2Oi8}UPvFe|<o^WlKc9{zr41uQVVADTQcSBlRJE5Ex1hXD?d_vg`B
zI7j~c^FXWjGFEOr9FM$zKt`SY-48r$SOM_Uuttp#yZXFbbPupo(Ehex0CqIQRa=9Y
z)P?y9acjwb-4HjPo>dfallE|Jg2Q(VCm82dp5PG8c#4~5Q=GM$VC;eA1P#wy5{z55
znfiy6rqA;aC73Y$IDOu}2F1HmTP#zoQin8Dw;21`v+`-mHrt}F#Y^998{DGey2{C~
zE*Ec6^^0Av+6QXW|4QxeyZY7W^o-u$fAVI2m4Gcf(@VTz%BPPSdd$1Tiyi5aBimlw
z|KP`f_Pfdlv7?)%o#>nW@we;hrRUJ4oh@_Y7;`@ClFWbjG-l(y(u+4O3rPE=)h7zo
zglEixh}yp69=&6}?LVixZ{<^I%X8iFIXX3!`8@62r^|dNqz&+{dh_LsztU!XbMr>N
zBGohW#eE7quCpUUU#$Q0P4WtWm#)9ewHxq6<eX!T0grrndf5ZOqcdhI+aP}1Y1u)<
ziCx1BAx<tgz7xUE@;xHh?Am35(o{7;W+y{&0*N=788Vl;4J*0lgDH-?MKL;<;?(05
zC!VG_xw*jV=?N~70~vJ=zTQ#z&#;U-dz$`H`^jRE>*X0Gn~nhtPa4!H1hDmngNaR0
z-tXUE@dojA%a*$km-laz8*yXwn|_F!{ad#|+{8*~{nB;!BJTKO%_!y@NYGT@0*QOb
zHw_8Ko$N)>kUx=N!ti6%Kjk{byI%*C-@oePjS&+A%I_;3_EGuv;Ld+&<m7qaes-?J
zy+FXnvZsR~+`ixCSNje{WDouV;r9LB@z*wgB6e*)trfv7E8__cH~da;)Zu;v$D|}v
z{Okn9iR~z^+{(ELy{=<=2T}L_jui8^p*WBW->zMHVBJ@}4)v=E*y>U0k+Fd5gSY;g
z1h}Hum`c|Ge_DC`#R<Trv!-9IiuhfLUX2l_e0{D7;y0CR?IRerCXJw!X-SZ|7(y^$
z$5x782U8qdn8SP&r#9gHOHds7Dfisp(toNyr2oO`bHAMdd>HxNK^E{yogPK{V7_GW
z{!bA@^6fj1*zLhA%Ey~8=c2o7cx*3%V=ItwUB_M>PuyetTT+xhrx@LkVp9kFZO8F1
zA()1<<}dHxt{pG#8RZL4^&T~2$F`Mo``sCsI6i%!yw%xajqC9KJO=8jOssV1hlA-}
zOWVvIQam~BZM&s|nm#QWQ0|@Y>^9F=2j=zIS^tUacmwO4fBHlCQ*XW;c?_%@64N(t
z>07+pgMSYPn2R4@k#3!AV9sA|GUZUK`+$BWef#BME>>S#C4KyU=G?c-FEzdq2Dqc^
zgUaKWvz@%lz0pqvtTVC3T`6-W&*Ii)F4Y2Dv~^hX+RUkK<?e3nv>C8`&5-DN%*ps~
zUOtR|h4}32)0LSMetGUZ-7yjI)}HMnnPYc{o$X$aC3y7v$IQ`%f9@;4KY(E6-_J0I
zH=NVAU-TP6=IPhW!5VwIcb}9&uxGJ9nSI}GY^dC(CCHY^#q7zmxU^yLY>GWvGCSkH
zc|Q8mc8b6DW-{P?(<81@jC#tX0sY`7MJV3R$!tY>*NRfSQJ>j>@vpU^_;@z62EvB~
zbf8%MDYFFp$5j1}qSqZ}2D-~Nq1gGpbty6z(zg~F)8Gj5eT=GXd*JggDV3Ie##~7K
zXwr-1VrfZ=J~cW?N2Ddm_8t6W(a+ZL(~{&q&xVh>#6JhO<m0=gzB6-|D3|BErDpC@
zhr*haF?T=x%E!!3DKB+I+3JwF|B0S+wwb$6cHVD6@5GxtCCvR#ynPyB?mn44nze9e
z+gr-gz1z%gDlc_1+vTq<{SQ^Kbf+X*x-U1ebe|n&>0V>?zBN@X{db$}%<)zGuU?#f
zG%fC8r)ISaw&la|7!Xu+Z;b^#jCO4E=T&R0T*13N`1f!i#lfKssjH5^_VUUT_xGg2
z^A?^10fWm(R^}P_HDHggCMlKzHthFc+z7zNx$1;W1Z?8<@@zMBr<{D_PjG-V5HWqw
zs4&E|7uB*6yuIl;!N;;T1Rve_f#Ra;6xR{DX^)a<xJPL{=<{E$Qe3}~;`ZMtE^k5c
z=Mai(S2BmkT^iJ;>pJG}myHhZSkVXE@7i{pw-DSpo+MQr1(<7b|C9X?KiT&5M}j^V
zyb=8pwv9*hUHxzb!KtY%!TGhH5}emzEy1~`dQ<#kAbnnFKE*OY^!d9>oPRU=yn27?
zzEF<hGmg%~5z>3rjvVS+y@Yi4?Q9#bOz)7oE|;pskBf3i`)`lwc{y^GRB>We^5y*T
zsnWE->DTo4QnP=rY}0sJ!RD`~tglh|!nBkf$EG&bXA4PL-r%ZYOR2)mzuG-{?VAY`
zQ%Yw1;T`FBEhXrDclGa&lA2d!9v%36^WkRChX;3wS+p<t>5NlVBD-8@c6Da|4l%VR
zH9Pa#zxS2ncara&oa%MF>iK4A$A{%8o#CInuy==#_dKqid@g8So(i|GGds6!t$X=&
zXTU9enkK6dPp&^Y3vvI$@+%SdRG<4j;@*#&|3t9XS+KW)`=jb33069Nm0;zve$@SH
zDaF_KD8AXqVJi1Lj^bZ;x#vL~Hs_vqq4@S3_k0|O2Y5O5;QQeK)UThu-1*ZJkZXgi
zqhchuUo_kOYZSP%%f~$Y0Wi1s=rWy9?rUuic!&7$?me3k>nENMMfBbNb0LDG8_XrR
zVC+eP^ETcmI5%wu#rrMkbL;x{a1!-@qon@U_3hy}>JRn$hW_i;Nl#w`p1szseiy)1
zrFGrz0<PV2Apd&6pFa*<)D>{uo`HQ!Anx9C`5fYboZ8cf`xTe^5)92AN-%824+ML@
z`<3F7D2ki=P~7#9;(jh%(n{)1-c0dWPKv8rQQV@TxT!%#-J;7bm-S5toK*2y_4R;1
zUfeUQ2jHS7-)uO6Sgq@eiHNG89g`6w?<rd&D%W-1iO9Aaw}xQfvbP9!E*U|v)2NyR
zyD&E?ZdFnLr*)|R;PVu3MpOR+J*j`?xzwG|nPAs8G=8V(hAo!O9Bo*AqHT+@(M1YY
zy|-zb#&>d=QDv65s5E*;qwLYUT2!5_8a5<Dmj1f;_B;8<Uryh)yHnmZy2$}s%Z=JG
z{hK}MqZ3|SUQz0o^vFKr7He<p4CwGh-HG>y$kMLnDVy}^?Vr+f9{H)6-;(*vrKIui
z3e@&uwmnGOQt{MpX<JWSzLvgZSK7}Drk$0XZOH6ydNsCe@rr5dwnhvwJS&=(cJRfo
zr|axxR*x<FL(ieZm{q<TtCif68|<w0s(*hZ1?-r86l<z30=zKn+=)NYuDQSJ7u^64
z-i*tp1U&R?LzPs-=excNK)k<f<ClmJhquchn6@@A!TDDr2`)VUkl><+uPJV-OL6s8
ziu2Y`Tt0^4&)oBwYpHwvK8h<gQ(Qdf{nTO)Rh|PGbxXAAyC<%GM%_ZQ+BQnc1@ir-
z;jA|=0XG%i*Z%<E9JOp{A(Z=h>lc#xwW3AMz3A>$eo7RgGVHVrF``JBH;61Vv=qVi
zVP6pJ+VeQYIU5OfE^>}w=*p(ly?YbEE(LoN>~Q2F!LA=M^!YMc4?9IS3aC)b>r9>0
z9|9^A8FX)K_{e|?MF%vfP=f{RkgMHicL9UG%{jFQV8|TR^y7&0a^K&G_+{Cv2E>UI
z>uf@tls&!%aqr8QMDS3KYXlFkmr!i6hQqZK!z)o7v5sOWqyy+3mi$67s6Rz7F8qkK
z>4DxkdzFbBlpa{GL-4rYmjI^hEbwSI;QBtdM(qV$`+K&;Ujf(s`gGG3#H&}@c1AoS
z?T{bw?D^?;5ziH`b%VH{SDirc#kE5OUmhPoaqe=8b2JpEr&65t8|UAf;v}%w0^OHC
zETcF*nc_?}q`zgcO`8G{7q7{+8nCsZLi}OO*DP79eTc)lXM2D+nTeQ%ILUkDM8t_J
zUt_$GiASoZQk*w~;Ng4)DK<!?sNneKwTik&{7ms0$6ucB9Dk1#kyHQDpHi$hhl{t?
zD97cyzKDri+hyCzV?FLqsPR|&zTdV!%ekfv@6ThPu1ob!d3yV$d!5KOyK{d-+S~7a
z^7Qv#7f`;{;PHCqm%!qkMqpTe9i<2V8S7dD>sC3MF6n8zp0)0j5|z$RSi+aD2mc-p
zFc)Lqb$-a`nG4llWcwj@Cg9^UC4(z77d89x7OnW2Id?wat>rbY09ISFL*c`m{o3bt
z!Y`cwKWUjT)sHz76tZyVy{U*LXC7Y0oO;>+!jypBfU@!pzHh;tTr*_!_D`B2j;ML%
z5_3ZFP5xr7S|H}Ry0$QLyzH`XPB3Q>r*@m#lR0`QqnwxWDZyO1_b^9ttq;5s(SzXe
zVH)P3_C}YV>(3$h%ZZ)LzVoHz6Tb+fc=Q3YCn#jF!as^)$=l4%HA6-=n{$}rf^y7u
zcu&bXWhwfVVbX#Axj9;j<0mp(k>1>56c-I(HemckQz))WV%B2#H#I4qyU#2^|96uq
zUOmLjK=*@H1e=#uSeGJmG4DQA#)b08_u1i#zXCoVTs7yj(ago%AK#Z2YLu1~G5J}x
z8k5tK!k-M^^{3IE$kRznQua+cVAyEp?hv_G@`;%{pn1)iJI&mIZx>G~ZRT!SZCYKk
z^`^?ew^xJA)~704W`DP7oO$>hFCH1!3BvnZ=2-M9z{3B>p5YemWjD<1oDhDQ&mK$v
z+Al5LT}N2D+t;;p*E3s>i}ZS!*@e;F+3Nk}I#|Xp5oPJFd(`NTmp?hX?L>#PA!XwV
z)jMXiKlA=zXWi#tuYY^q%_BbdU3RAv@Alx|!-15TZ0~-}@os?i+fr9o?+bQTx%G38
z`~7ktAa2#ffQmI3z&}o|t@9b+`Uf?b^MG4J&)@6_nA~PY&E<eeX+LePgP3~jk9vqZ
z3p5Kr+)+8r7jgT$eA9`$@a&-k3(dPsvE*Zl_3~3>E>djIVdF~F?SGhJ*VPnjMp5)V
zMX^a4=47*{6;8Iw$((E+y13enRe+@wr!M^gSiV}{1Mz?r&ZbJ|ArAQXs|AR%Ey)KF
z*^&nrAu3xvIz_PNhJged#(YoE|J!#Ii&1}n?@iR*d<w;OV=1<{M$vaI#SRjR)#_0E
zq&~$aRi&@g@pGp6pO&5}dtgaMwy4x?jixRtQ*@kE`>f)|hzsFTL(b$g+26HG?by2V
zfSgB0r4}D}uutB+`=$Bsm)dr1;^dT*!vZEW{3}n&*7jYJbIZ;)pVxoX$=Tj%DW6=M
zJ#72_XDK~WLPuY|U!b{n>pG?GEL_yQaP_&ZUUv^j$#Zg)vPO#u$v2V{TZS~5-n>Aq
zY4`WM%bEOcR&0X0qbm8`<WYqNt^A?cqjxFohkYK{>{_-WM<VO>WOnY~vS{Qdp@5t7
z?ESqW;<LqnfqixV=Ifm==nyYl@ALri(v@Fp61*|+bApEk#1lN4k(c1Hz<LyOeM>PP
zp?B$6fv(hD;0gD93dKB06!W*@a4q*-Lb2dciUmk|E*;D7`+jSOhdR`^lm<a5#Xz1X
zwLF@-8nAGO`Jsye%TBLeW(8nHrS4cyzzTg{7D-0@tW(=rh!JHEuR&Dy`jpnYa#aoy
zccaq-2>L&-MzFzOd8u1+lVb2miv5aH^!}1!!^RXfTPRjc<ovmM*Q6?QS2OnLC;R6y
zceUmHAADs1clxR8Gd2P4>0WYKG~oVH-&}|T+}GpF-CGd5fALWl#D4i}>_+Ufz11*+
zGm;|+u29|}xIAAkg3H<lQ5?9R;^>(aXC_fh45k=6pW=Lfii%tmW6Dz;3ii>kyE%0^
zz>XTW|5>(@c>!<!@U!<qz-#MA|NI-^)tj<}9O$mg4&H-U`$6YF5kFpE?;T>zPN#hd
z%6<U55V*VFLb`w_R}<`UrwPHH4-yE5-7H6+rxLm~!*Va8?wVl#fp|TVr8L~P_X+m=
zt{TB^``WZvw&%jF`6awtj2#|aXYQO!+cd*j#oI$?TS$CL{j_`w+oI|(k7oZ-?PmJx
z=YMqRJEU~_HlN(S)1<!zY}JGtZhFs4AM^Q)6Sr~=O^@_1P~_v>ivvRakNkA8Te-Bm
zoog|9e=C%p^LfLEZB}n!E^e)OtMHm3%&uu)ru8j;D{bDBq)jW7Mba{Ee|0};Pa9@$
z(y{p`l%J(77<OZPWuGQ#`)`bDxoAv5=F5TwW(_QKo!R=)@T0-M&I3DZ-A2l4USP)z
z+v$_91K_!08<$K4JY4O_!zjT0O<M$B1U%5WK;ymywXYr^1~r(7>s=3ol!Nuch+)sh
z5^Ovf?X3qwnr<`@Y_^|-4`~MWaeUrvW;(?&xc(K~-!7y0)k2DK!zq4OmEzb@+;d|;
z+0*j89}Z;H{cQN3VfD9U)a`y{-=UcN8FhO#-@2<g%qxWb89S-v48RMY>bI;zc~}0l
z{CmVfY}c2F<+C#j5I<JWmLS$DQHq4C+45*Rbl2?sMLmK&-wY(!`=eNbVX@N)cAL<g
zx|v=CdrZ$xu*-K(=<}Y%3HBVcl_1o&fC{|=m-_3M2UO^}dikc=?*b}>ElfS9i3j|u
z=bB&t037@Kf}qa;zkD@VxgPPalu1nxZyZhryK2~%ce?L>L~!5Uc?7>6)r;Wt)g=6#
znV&x;?wR>V5S+2~2F3foQGDP-{Ts@tJNXR7!@d+ReM|k@jY<#nEnTrj&1&g^K7ISV
zJ3R#OaIXEI!n#(Uol#f1J_p?PVt0dvfEf>-UzQ<0Z2we>cx{q$IO2^)SBd*Z^k9lz
z&K@VvyJg4;cHdZ(;s~&#Lip|@z&?wZ7)^21Lhd=NZ^h>^K@>+Hrf2|rEqPuM(mh(|
z7tkJXOfH%9Bf#+)y;%+B=Z}o1Um{+t(6J-py+TEEBHnqLewm=Idp?5S7Q^^o-ubpu
zA%Zi0Zbk9IC4w{Sfjt-eXAGiznsM$reLjif?}aexE)z}h_Ej4G*ZdUkjik@pjpA!X
zM)}#nUDhPN-L@*%sY%7YNlD*7rqLMR-zV|@JO=7U_U~3^n=;+&U8r*Jt95B_?`5or
z-107<{P<amo&<g!Sk|+i^&h#OwXSN8-(S;<@#V;4VBOO{HYhs6Xb<6^d+_~ofVr4b
zCSqWfrOXBO;*USkNC4-AmK%`){gm?;+HO0_oUhcr?|^D~09i%q=a-qYNAp}%`40um
zS>-_=mN_$W!{v4tk0W-OvSTW9x`9G<So#KW`PbRN&U&-Ba)vYq;(+8ArI-`nW;9P;
z`Wa%=@!#KNj<<;^r5!wv;L}Z0n4_;MJ)YjLKfzgs@0lZEbEZA&dYRzW`pcPvNiRc|
zt8)-sS^Ya^f2E$U)ZV{R99fRpGjhWk@7}d3HvEg(d9%3uw~uF2^v=(0hxh$bxgW*q
zhnaMse`nxpihX^VB&63okz(KtW&_5r{Fb8cd}b|%KbJ~z-4kXB`Vae&V$X5R40KE2
z+yKNeQ>{ynx!Cn^lkBx;Am7(4D`x;cU!GXKLUZO~=juzJd+$$6QtiL@u3$y5x7KX(
z;?OSZ_-RRzGpBESGmU=^?lDQ%DjV&C;2tu#!*iXPe`3ozznQH!jZ9?smYrbk-(pH1
z<0?UjpBQv5?GJPR;A7{_)~7}$1}7vvu<-wUqK}1pMEO1z?%|*1v~+iAWa%EU#nPQ<
z{7xgLCwO?TZ!G=43AS{Ponh$?47GHJ{AKA*F!o2&;urPV(`oUUv>`nr3KZG+D<6(W
zu&3@1sdDo$@8%Jox12u8m{0z>2j4FTnj0n6@;#}m<Dg9QoU?rDalLO21bp5je9VS>
zih$31=9nK{GYH%Z^<Qq-2)KN4{SLJNSIvLS-UVEFcxuz15fA;`VKw6E*WrPPr=~x7
zKu}RjO)&B3ZGs8yWd!3B&QY9|LUAR&XXUB*U5_dLl#k+y{uEE-p*VdQ#U)aTTjnt5
zf~7~!?db-%r)T!oQGkOl?!7P?FyW&uFK+@4-SEf5;ed&6O3YIrzW%Z7=lDG5{sRvX
zb4=Jdil8hpouEO_5Y$&XLQvPW3&qb@Q5^P^;;d#A$G@XE=1=bVHxwg6C?=Mr_??Ey
znRoWGmMIUVSK~V`ZTzl$>fRq#?fGQWENT4B{2#}uRnlQy!u{tR)TD-ebT4_v!vd)l
zy6>r}Z`N8`GHzu4$unQ2T-p7+&(KpnQ+DK=JURNp-sVf2%5qO{)HkJh__W^ZHx^47
zu%u0olxlxA|Fq!3`K3>nZC?J5N^>`_O-?CTCcf{*Z@x^)d#6&#$-m@iUZFw9N(;Y!
zl9IFM=8^#?6H-dAujJdm@~UP}R&~t5*1ga?dymEW_I3h$=&2$T8?I{&Sm={!b;|)3
zo*ecV?5L*-Ym>|Vil{t%QjHjycKijwo%6y7E?Ux&;L_DM2ri3)b*~U^S<!YBO9fLb
zn?SKV*jLfN%tMMLb8`Mt?s*Kw;xNyO;VNXKSgt+AQbs-ZzU#?<l|FMW=zerklT*yO
zkZR|ymuduZ9a(qW(87RmyI(v>1N{8UKF6`1Bt9Eot0Uqc<2DQd9GbgF<pYR0mg)u&
zZ2S2Vf^oBT1of+K5!Cq>q8P^2x5@B6h^NaaP0;}RvCuuXD8-wXDTW7A{Ngc1sBe!F
zdyNcJ|HwQV+W2}In3s8!xUQyjX)VCyDXmi`0;X^7Tgn$OZF=hL(umW(EL<3Ie*C=0
zh;v7cyGxLnd5)lL!aaih27N-XZ|hKsF-Is)`jz6(bt#TXqB!vq#T|nv4$4DuR2ap1
z4>Ia5C^PirL0Hc^|7GP`+CJbev1(f2cYsBACbat<u*ma43kRcH`sK*!h|PYByMx#`
zsVvw@7Zh#MyVxj#Q8&R31nz#Hgc0ojMLa>-raBZy7Nl7COM?C4z@7@>`}<v|?si3}
ze>Jcd!Seyz_fhOa;|<6NX|e3p*3r#J-D)v*>YPtn4_d!Xb0qJIO$mEiRH^(nq4)1G
zEvg>sQ7+<%B>l~-m&1FOotvIsv&*Ex@)ZGFcLen;+HOz!*wJ6TsC2YddStF5fv-0V
z4+wqseD2G~&(ofyzL@tNqfO6MJR|qO3D21eOVVF>U!1_~+vnxm^8TQ-`2`nEeC}5|
z?ZD~czkfGsK6A24yFL4Mf1S2rV2dw&D&0;yx+I}`%bInV6^s0e%$hle*_FHPM`g3Y
zI@98{1A}h;0Cvpo<%+Id0eJd{#J3gEu4$MuXFK4Yah20`fO{8=ZH;!-y$hPm-j6u_
z(vm5NGf$k!j`(fS%|nR$x0PuGpBDU>;2#xRQuj3uFQOfF@64wexPBGzl8WL<T;B@r
zZMP{teoFE3F^bm~{8!tVzboGvbr=1(UlLR~qwd_&zkHeo^9KtCP8>Q3)@v`QR`;vG
z&j9mvz3?@xGha|N@%s^15I-Nk5au%$6s`M+S%%nn?935}O)Inn`zh3uwvS5_{LB~T
zCBQAa3G)_!;gx9pl#Q-P{n>2PzXRBJ!M}fv{RDfD-b_$F>k!3JGYJmJXc<spLCneQ
zcd7?en0M*eeQ#I?Ise#--+I70%lWDK<057QZe07ORz1MYDc(~i6Ri677~+qO#!N$8
zb}X(k;_{X&VVxvkp4Vpy7HoTjV4;^qDDH-Jm*8J$lNZHh1u5PEJ0d(U>=jGhpMqT$
z!%YmOIF5@K4R+A_KIz?_WTXfB89K-xRRp|UBV_$0z)M|I%7dNN@7yKLj5xsarydnO
zMQ~~1JBa)D+=un1e&-L4yzmFX&ArYNoa5b};M`W-2+nIzg}Pt0r}*+OiVwhU3-RZ^
zo=)-k1&a5;t_#oS<=RHw&p3Q_3(~*7;V?beIp?q6cP*$W;FeLXhVI3DZ8#OW8u5pj
zeJ3EU`0BU{arvNZFA<mBOa7X;3;%YKV3F*k9F`ZFw1c<{jVIxj7rIO0FE5k_%MbWA
z6YS;4&+#08XTiE-e7=1&eSUlub%)Y$g&K^f?r7sWl#KG16W;kB9k6Xx(_astOG`^X
zP;Rp9!jVnZ^(avHOOv&w<}Xao_DQ|)PU9Y?y=yMNm~tX6p#0$<OExLfKCsf55g5m_
zj`g3hj+HNmz`DQm@XcPsc0KDSl3FJXXN~LNu<Sh7`{w|2DR1p#(S2qy7dA@$zurF*
zu;H4p2L|R+u3L`==WE2ASH$#DjOY(Ian;qS(C;m_`PJ@RaGpT14RwY{M=@v8A77vO
zb{b+>kF~d%)0)}?dVSj)ap9+bcr&LeDnc7x9ER8yOt#F4{c}ri$!kC?IdTrHV~s!k
z>#rN?BVLVqm4i7}v*@Z}!{Iyuz)^WtF-KO+c)Wg1ae_bfXvG}LlXqPI`ZEa*$di}Z
zulW3<qBDIdUj2^Qlm2-9_BKZ-ew>5ZRZ-DSKkF&Qg>{+j@SYkKN>HqJh)D<ftH1l4
z;>3PT64INtkm4LaW+TR5_K@O+dCXb}zow@b#WF1OBlsI?B@@&?Ey2vh@GCzfXsBv)
zJ_d7f?!G;*?lc0v*C`)2?HAzly`^b6CNUT1lq<M6@Aqj*(WUD7bo&GBtx?sYZuPR3
zfV8B6=U3!Wnw<kMa9^80{fzSy;GcZ8(UuH89=OxCtu8Ry+<onl=?mQGpIRriqPc&@
z;yZiI-P^xEXtw@zVA}Q;F~;#J#^2s5i4P9vy@I=SSF`h0!5#9sncWxNLq6vtknka?
z#{KuyoxIG_f94oV_X(+`d#l+w14#cf%k+jQjn+9zf5ETD`0>~0#5bL{Hf_lCDrIul
z+0H-b5$vf4`iyzKhIjLb&%+)}FrGiaKlkAOhXbjr_d^}Xv*Z8s_W$l%14ib3TtyQ1
zRlvx6eLlN9Ari1$Kf}NxfR%1l_YVh@+%A@M53%=&KYI~W-I#$WTb_IYvHyl|*Atu(
zaFgH?rUk*pb*~Uyv>+G7^$8TWad=uv-Dz_v?od*^*p=ei&J?#kr?|5|b6FN2azmDb
zxjew{)h~_Sf;(;8@^{Ap_YNrE)Bt$ka<M$c0r%fHbK8&Ll+Tk8AAd6`0P*4ZJf9L=
zxbZx}G`~QCDWOXVCV%=d#l!;~9;68NA&j3qaW=)t`zTKMfnpEx-iOICe{$g}F$KI1
zt}C>oo%Hq>*DI92P&f7HoT689HJL6Q*|g!CrRTnvj)`K*e^J6eRW{)Jw%h$OQmf_|
zq6~d_Kw4qx;JsJtk59RI<4A`j&FqwYaYOrVZSZ^ZwX>#;dpmJdN}HvJ)jik&DG}o;
zk6C|cOmk+(^J|}{4sBk2)pt7s>YYj{+DN4;eY;yqg||VQ-b}61yw>UQza5d+O(~Mn
z?OoXKds3<|$#FNR<Id(@ySnt-Q!{__lD8*34{JJ<*?DE#o%};!o$8g;BSD9H0jAF>
zSA7Cv?6BsK5hI(=UW_=X%wm}T1gw;^I^w|JC#Dm8(P%5dm$d?@dp)d^1^*XGaGn5S
z8k|pnxV;4D?niM;EzTW6F=aQ!%?!ouaNYpM+Xi-7#8mBnx90vmIR`n;KXA{hI#AzM
zNBef02kv!#=k7@XcLn9_^%da3=j$U#y||G2<w$VvKbNfi3-R|Qr8gly^7`fq;=?fk
z5`r1kcM(j@cZOhc@mmy|(C5jc4^j7Ze~JS?qBwdT!IYBR^ETD#^QdTw6Du*lADF#Q
z&H6FFAM}}iwRsnCpIA7ia!J4oVVl+c054bVwI&wuQor$?D<k%O*6bHVc4hCAh{_(T
zE)lG?s5L?6{zHP&=w<|)uaHsP6+!Wgj^dw>DV}^r@xo?`Z~IbAUqbP49L0;N8FlYA
z8ovD8J-{+2PUwRH>$g4;4EEH!J_9bqjtBIfapxhdf4l2_p~8t0i071ThajHzUR#0S
zvA**OUOo%w1i<qXl>-T$PJc`B_a+2Sjtijfn;R%93lltb1?)%=@AS+h>fUvP;Mt;>
z9(cY%L9uGS7R$XWPAah?ti{+RMb{_%5VTEOIj&-rwm-F~T5`%~MVl6EQT6tB@mqKA
zN`KqrSe}R$+0s+*1n2zjn|%REWgApHeB#gaaZ|dsefIEadc>)-t#+=<7tnFj;za{4
zTuOWXc-^w-F)Pz^ojSj<?$UwGh5p~yD1Y@fb26-H#?+#<)0Q5p_;SU%WogGNwr)Dx
zX9{yRcJl1#;_|c&^Kxz|^irC3O4ahmKWZ*#QvaCHqI<hJ%)Y|ayUqNw9N1ZF2G?#i
z3GA4i!<$~H19-A>>29mhu9?#2XF1@`z3=+ty3<{E3biSQc)L*jAjI=i*XKaIIJ-QC
z+ja58bS=SO?)^k?<K#C4H%+7AHihP+?w%<W1A9^I{gz@_ev01JD2C6X7}A7duxFn6
zJvVC^b??Iz<Htoa>fU0c$vd$=jK5yCq6{!l-86jwV9i&#`#uKrsqlIU%y-=NDOY-A
zbwKZZZQElw@6&~b!@AVF-e<n)40c*@r~7sx2>Ue&p4kiQ(c$^IU9^6kZT6J<4>>{b
z^z+^n!@#Z!;SOsH5IlRm7R7dO?gu<SU9~_!g}dK;H#h{=quzcyzva{?0Tu2%Y1ET}
z^^SKg_5EH4cELNR3zbUT1$egN_fG;5EBG%Rh?uiIvlTJVt2?a`^F~drLa<)%sRaF&
zl_%IBK7``N-W01R6Rdw4?2QnwL0BG&k|y-|PirZ@%uAn_;^I{VyQ+VezV-hAyJv$L
ze)WC-1b1F(kwQO$`>lVTFUJ7B+B)`B7r?(Z$ydf8ZX9}I65@(sxfda>3Rzu-;83qK
z1p5bpJr$l0_-Zmid2BqzMkM|!IRoznf4M){X%U0KzKYl~JH<M~D2C*r=)++HZ%F@%
zS>IRt6Yylc2KiTDey*PK0ekeFvn6K83jv<_D<x?by7MJ#Y5u3YEQ9X6MQ(uo8qj}w
zH;PMPeJb*+))xfp-wvkkCIQsl2-eAhfBl_s&H#pc#_>0|FLnFTaQ+!lG+YXt?||VK
z1yTQs#&Lc|g&MP$oPP7mwpB6ld546KN<XkB@w+Fc@2v3EC{TBO^R22^eooKUd(j_y
YX}`30Ul*@YYb&f<F7-p6K9BtV2b?&VdH?_b

literal 0
HcmV?d00001

diff --git a/test/tests/concrete_ASR_swelling/tests b/test/tests/concrete_ASR_swelling/tests
index 1c1543af..23021cfb 100644
--- a/test/tests/concrete_ASR_swelling/tests
+++ b/test/tests/concrete_ASR_swelling/tests
@@ -14,7 +14,18 @@
       custom_cmp = 'asr_confined.cmp'
       detail = 'with anisotropic expansion strains, '
     []
-
+    [test_strip_kelvin]
+      type = 'Exodiff'
+      input = 'asr_confined.i'
+      exodiff = 'asr_confined_strip_out.e'
+      cli_args = 'AuxVariables/T/initial_condition=308.15 '
+                 'Materials/ASR_expansion/temperature_unit=Kelvin '
+                 'Materials/ASR_expansion/reference_temperature=308.15 '
+                 'Materials/thermal_strain_concrete/stress_free_temperature=308.15 '
+                 'Materials/thermal_strain_steel/stress_free_temperature=308.15 '
+      custom_cmp = 'asr_confined_kelvin.cmp'
+      detail = 'with anisotropic expansion strains and temperature in Kelvin unit, '
+    []
     [test_strip_isotropic]
       type = 'Exodiff'
       input = 'asr_confined.i'
@@ -24,7 +35,13 @@
       custom_cmp = 'asr_confined.cmp'
       detail = 'with isotropic expansion strains, '
     []
-
+    [test_strip_in_tension]
+      type = 'RunApp'
+      input = 'asr_confined.i'
+      cli_args = 'BCs/axial_load/value=10e6 '
+                 'Outputs/file_base=asr_confined_strip_in_tension_out'
+      detail = 'with uniaxial tenile loading, '
+    []
     [test_full]
       type = 'Exodiff'
       input = 'asr_confined.i'
@@ -36,4 +53,12 @@
       detail = 'and with a mesh that captures a larger axial region of the cylinder.'
     []
   []
+  [max_it_error]
+    type = 'RunException'
+    input = 'asr_confined.i'
+    max_time = 5
+    cli_args = 'Materials/ASR_expansion/max_its=2 Materials/ASR_expansion/output_iteration_info_on_error=true Materials/ASR_expansion/output_iteration_info=true'
+    expect_err = "Max material iteration hit during nonlinear constitutive model solve!"
+    requirement = 'Blackbear shall produce an error message when number of iterations exceeds max_its during the volumeteric strain residual calculation '
+  []
 []
diff --git a/test/tests/concrete_anisotropic_eigenstrain/concrete_anisotropic_eigenstrain.i b/test/tests/concrete_anisotropic_eigenstrain/concrete_anisotropic_eigenstrain.i
index 4fb5cc59..a0c8f19d 100644
--- a/test/tests/concrete_anisotropic_eigenstrain/concrete_anisotropic_eigenstrain.i
+++ b/test/tests/concrete_anisotropic_eigenstrain/concrete_anisotropic_eigenstrain.i
@@ -174,7 +174,6 @@
     boundary = back
     value = 0
   [../]
-
   [./load_x]
     type = NeumannBC
     variable = disp_x
@@ -199,28 +198,21 @@
   [stress]
     type = ComputeLinearElasticStress
   []
-
   [./concrete]
     type = ComputeIsotropicElasticityTensor
     youngs_modulus = 30e9
     poissons_ratio = 0.2
   [../]
-
   [./strain]
     type = ComputeSmallStrain
     displacements = 'disp_x disp_y disp_z'
     eigenstrain_names = 'concrete_expansion'
   [../]
-
   [concrete_expansion]
     type = TestConcreteExpansionEigenstrain
     expansion_type = Anisotropic
     function = expansion_vs_time
     eigenstrain_name = concrete_expansion
-
-    compressive_strength = 30e6
-    tensile_strength = 3e6
-    expansion_stress_limit = 10e6
   []
 []
 
diff --git a/test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_10_out.csv b/test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_10_out.csv
deleted file mode 100644
index 8c9c6bc2..00000000
--- a/test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_10_out.csv
+++ /dev/null
@@ -1,13 +0,0 @@
-time,expansion_strain_xx,expansion_strain_yy,expansion_strain_zz,expansion_vol_strain,strain_xx,strain_yy,strain_zz,stress_xx,stress_yy,stress_zz
-0,0,0,0,0,0,0,0,0,0,0
-0.0001,0,0,0,0,-0.0003332888992563,6.6657779851259e-05,6.6657779851259e-05,-9998666.9776889,-3.2014213502407e-10,-4.6566128730774e-10
-0.0002,4.4434077036047e-09,4.9997778296148e-05,4.9997778296148e-05,0.0001,-0.00033325113760079,0.00011664889449788,0.00011664889449788,-9997667.4302543,1.1935917427763e-06,1.1933734640479e-06
-0.0003,1.22186401892e-08,9.9993890679905e-05,9.9993890679905e-05,0.0002,-0.00033321005022565,0.0001666383444531,0.0001666383444531,-9996668.0659748,1.194013748318e-06,1.1927913874388e-06
-0.0004,2.3325086939729e-08,0.00014998833745653,0.00014998833745653,0.0003,-0.00033316563773969,0.00021662613002188,0.00021662613002188,-9995668.8847983,1.1937809176743e-06,1.194013748318e-06
-0.0005,3.776213761212e-08,0.00019998111893119,0.00019998111893119,0.0004,-0.00033311790075152,0.00026661225150905,0.00026661225150905,-9994669.8866736,1.1953088687733e-06,1.1942465789616e-06
-0.0006,5.5529182033402e-08,0.00024997223540898,0.00024997223540898,0.0005,-0.00033306683986965,0.00031659670921935,0.00031659670921935,-9993671.0715499,1.1921511031687e-06,1.1922675184906e-06
-0.0007,7.6625610200449e-08,0.0002999616871949,0.0002999616871949,0.0006,-0.00033301245570235,0.00036657950345744,0.00036657950345744,-9992672.4393761,1.1945958249271e-06,1.1952943168581e-06
-0.0008,1.010508122801e-07,0.00034994947459386,0.00034994947459386,0.0007,-0.00033295474885778,0.0004165606345279,0.0004165606345279,-9991673.9901013,1.191277988255e-06,1.1924421414733e-06
-0.0009,1.2880417860901e-07,0.0003999355979107,0.0003999355979107,0.0008,-0.00033289371994389,0.00046654010273522,0.00046654010273522,-9990675.7236746,1.1927913874388e-06,1.1950032785535e-06
-0.001,1.5988509969359e-07,0.00044992005745015,0.00044992005745015,0.0009,-0.00033282936956849,0.00051651790838382,0.00051651790838382,-9989677.6400451,1.1936208466068e-06,1.1909287422895e-06
-0.0011,1.942929662099e-07,0.0004999028535169,0.0004999028535169,0.001,-0.0003327616983392,0.00056649405177801,0.00056649405177801,-9988679.7391619,1.1938391253352e-06,1.1937227100134e-06
diff --git a/test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_50_out.csv b/test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_50_out.csv
new file mode 100644
index 00000000..d4fe60c0
--- /dev/null
+++ b/test/tests/concrete_anisotropic_eigenstrain/gold/concrete_expansion_load_50_out.csv
@@ -0,0 +1,13 @@
+time,expansion_strain_xx,expansion_strain_yy,expansion_strain_zz,expansion_vol_strain,strain_xx,strain_yy,strain_zz,stress_xx,stress_yy,stress_zz
+0,0,0,0,0,0,0,0,0,0,0
+0.0001,0,0,0,0,-0.0016655568500029,0.00033311137000059,0.00033311137000059,-49966705.500088,1.1641532182693e-09,1.3969838619232e-09
+0.0002,4.0407154020135e-37,5e-05,5e-05,0.0001,-0.0016653904730541,0.00038307809461086,0.00038307809461086,-49961714.191624,1.1996598914266e-06,1.2002419680357e-06
+0.0003,8.8894966257367e-37,0.0001,0.0001,0.0002,-0.0016652241210173,0.00043304482420349,0.00043304482420349,-49956723.630519,1.1994270607829e-06,1.1990778148174e-06
+0.0004,3.5254263760279e-36,0.00015,0.00015,0.0003,-0.0016650577938874,0.00048301155877751,0.00048301155877751,-49951733.816622,1.2001255527139e-06,1.200707629323e-06
+0.0005,4.3078694856466e-36,0.0002,0.0002,0.0004,-0.0016648914916595,0.00053297829833194,0.00053297829833194,-49946744.749786,1.1991360224783e-06,1.1972151696682e-06
+0.0006,5.0447289533329e-36,0.00025,0.00025,0.0005,-0.0016647252143287,0.00058294504286576,0.00058294504286576,-49941756.42986,1.2006494216621e-06,1.2002419680357e-06
+0.0007,6.0356005287523e-36,0.0003,0.0003,0.0006,-0.0016645589618899,0.000632911792378,0.000632911792378,-49936768.856696,1.2008240446448e-06,1.1993106454611e-06
+0.0008,7.1415772009698e-36,0.00035,0.00035,0.0007,-0.0016643927343382,0.00068287854686766,0.00068287854686766,-49931782.030144,1.2018135748804e-06,1.2016389518976e-06
+0.0009,7.9626511844921e-36,0.0004,0.0004,0.0008,-0.0016642265316686,0.00073284530633375,0.00073284530633375,-49926795.950057,1.1936645023525e-06,1.1934898793697e-06
+0.001,8.7016754749121e-36,0.00045,0.00045,0.0009,-0.0016640603538762,0.00078281207077527,0.00078281207077527,-49921810.616285,1.1992524378002e-06,1.1960510164499e-06
+0.0011,9.8374790946429e-36,0.0005,0.0005,0.001,-0.001663894200956,0.00083277884019123,0.00083277884019123,-49916826.028679,1.2004747986794e-06,1.1988449841738e-06
diff --git a/test/tests/concrete_anisotropic_eigenstrain/tests b/test/tests/concrete_anisotropic_eigenstrain/tests
index ca6b002f..622b0a3a 100644
--- a/test/tests/concrete_anisotropic_eigenstrain/tests
+++ b/test/tests/concrete_anisotropic_eigenstrain/tests
@@ -12,73 +12,126 @@
       input = 'concrete_anisotropic_eigenstrain.i'
       csvdiff = 'concrete_free_expansion_out.csv'
       abs_zero = 1e-5
-      cli_args = 'Outputs/file_base=concrete_free_expansion_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  Outputs/file_base=concrete_free_expansion_out'
       allow_test_objects = True
       detail = 'under free expansion conditions, '
     []
-
     [concrete_expansion_load_5]
       type = 'CSVDiff'
       input = 'concrete_anisotropic_eigenstrain.i'
       csvdiff = 'concrete_expansion_load_5_out.csv'
       abs_zero = 1e-5
-      cli_args = 'BCs/load_x/value=-5e6 Outputs/file_base=concrete_expansion_load_5_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  BCs/load_x/value=-5e6
+                  Outputs/file_base=concrete_expansion_load_5_out'
       allow_test_objects = True
       detail = 'with a 5 MPa compressive pressure in the x direction, '
     []
-
-    [concrete_expansion_load_10]
+    [concrete_expansion_load_50]
       type = 'CSVDiff'
       input = 'concrete_anisotropic_eigenstrain.i'
-      csvdiff = 'concrete_expansion_load_10_out.csv'
+      csvdiff = 'concrete_expansion_load_50_out.csv'
       abs_zero = 1e-5
-      cli_args = 'BCs/load_x/value=-10e6 Outputs/file_base=concrete_expansion_load_10_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  BCs/load_x/value=-50e6
+                  Outputs/file_base=concrete_expansion_load_50_out'
       allow_test_objects = True
-      detail = 'with a 10 MPa compressive pressure in the x direction, '
+      detail = 'with a 50 MPa compressive pressure in the x direction, '
     []
-
     [concrete_expansion_load_5_5]
       type = 'CSVDiff'
       input = 'concrete_anisotropic_eigenstrain.i'
       csvdiff = 'concrete_expansion_load_5_5_out.csv'
       abs_zero = 1e-5
-      cli_args = 'BCs/load_x/value=-5e6 BCs/load_y/value=-5e6 '
-                 'Outputs/file_base=concrete_expansion_load_5_5_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  BCs/load_x/value=-5e6
+                  BCs/load_y/value=-5e6
+                  Outputs/file_base=concrete_expansion_load_5_5_out'
       allow_test_objects = True
       detail = 'with 5 MPa compressive pressure in the x and y directions, '
     []
-
     [concrete_expansion_load_10_10]
       type = 'CSVDiff'
       input = 'concrete_anisotropic_eigenstrain.i'
       csvdiff = 'concrete_expansion_load_10_10_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-10e6 BCs/load_y/value=-10e6 '
-                 'Outputs/file_base=concrete_expansion_load_10_10_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  BCs/load_x/value=-10e6
+                  BCs/load_y/value=-10e6
+                  Outputs/file_base=concrete_expansion_load_10_10_out'
       allow_test_objects = True
       detail = 'with 10 MPa compressive pressure in the x and y directions, '
     []
-
     [concrete_expansion_load_5_5_5]
       type = 'CSVDiff'
       input = 'concrete_anisotropic_eigenstrain.i'
       csvdiff = 'concrete_expansion_load_5_5_5_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-5e6 BCs/load_y/value=-5e6 BCs/load_z/value=-5e6 '
-                 'Outputs/file_base=concrete_expansion_load_5_5_5_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  BCs/load_x/value=-5e6
+                  BCs/load_y/value=-5e6
+                  BCs/load_z/value=-5e6
+                  Outputs/file_base=concrete_expansion_load_5_5_5_out'
       allow_test_objects = True
       detail = 'with 5 MPa compressive pressure in the x, y, and z directions, '
     []
-
     [concrete_expansion_load_10_10_10]
       type = 'CSVDiff'
       input = 'concrete_anisotropic_eigenstrain.i'
       csvdiff = 'concrete_expansion_load_10_10_10_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-10e6 BCs/load_y/value=-10e6 BCs/load_z/value=-10e6 '
-                 'Outputs/file_base=concrete_expansion_load_10_10_10_out'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6
+                  BCs/load_x/value=-10e6
+                  BCs/load_y/value=-10e6
+                  BCs/load_z/value=-10e6
+                  Outputs/file_base=concrete_expansion_load_10_10_10_out'
       allow_test_objects = True
       detail = 'and with 10 MPa compressive pressure in the x, y, and z directions.'
     []
   []
+  [error_message_check]
+    requirement = 'Blackbear shall produce an error message when parameter(s) required for the anisotropic expansion_type are not set by user:'
+    [compressive_strength]
+      type = 'RunException'
+      input = 'concrete_anisotropic_eigenstrain.i'
+      cli_args = 'Materials/concrete_expansion/tensile_strength=3e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6'
+      expect_err = "compressive_strength is required for expansion_type = Anisotropic"
+      allow_test_objects = True
+      detail = "Anisotropic expansion_type shall generate error if compressive_strength is not set by user."
+    []
+    [expansion_stress_limit]
+      type = 'RunException'
+      input = 'concrete_anisotropic_eigenstrain.i'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/tensile_strength=3e6'
+      expect_err = "expansion_stress_limit is required for expansion_type = Anisotropic"
+      allow_test_objects = True
+      detail = "Anisotropic expansion_type shall generate error if expansion_stress_limit is not set by user."
+    []
+    [tensile_strength]
+      type = 'RunException'
+      input = 'concrete_anisotropic_eigenstrain.i'
+      cli_args = 'Materials/concrete_expansion/compressive_strength=30e6
+                  Materials/concrete_expansion/expansion_stress_limit=10e6'
+      expect_err = "tensile_strength is required for expansion_type = Anisotropic"
+      allow_test_objects = True
+      detail = "Anisotropic expansion_type shall generate error if tensile_strength is not set by user."
+    []
+  []
 []
diff --git a/test/tests/concrete_expansion_microcracking/concrete_expansion_microcracking.i b/test/tests/concrete_expansion_microcracking/concrete_expansion_microcracking.i
index feaef91d..a3625e00 100644
--- a/test/tests/concrete_expansion_microcracking/concrete_expansion_microcracking.i
+++ b/test/tests/concrete_expansion_microcracking/concrete_expansion_microcracking.i
@@ -215,7 +215,6 @@
     microcracking_eigenstrain_name = concrete_expansion
     microcracking_initiation_strain = 0.0001
     microcracking_strain_branch = 0.0002
-    expansion_stress_limit = 10e6
   []
 
   [stress]
diff --git a/test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_excluding_confinement_effect_out.csv b/test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_excluding_confinement_effect_out.csv
new file mode 100644
index 00000000..c5d67efb
--- /dev/null
+++ b/test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_excluding_confinement_effect_out.csv
@@ -0,0 +1,13 @@
+time,damage_index,expansion_strain_xx,expansion_strain_yy,expansion_strain_zz,expansion_vol_strain,strain_xx,strain_yy,strain_zz,stress_xx,stress_yy,stress_zz
+0,0,0,0,0,0,0,0,0,0,0,0
+0.0001,0,0,0,0,0,0,0,0,0,0,0
+0.0002,0,3.3331666777766e-05,3.3331666777766e-05,3.3331666777766e-05,0.0001,3.3331666777795e-05,3.3331666777795e-05,3.3331666777795e-05,1.4499509991099e-06,1.4499086474625e-06,1.4499933507573e-06
+0.0003,0,6.6660000888748e-05,6.6660000888748e-05,6.6660000888748e-05,0.0002,6.6660000888777e-05,6.6660000888777e-05,6.6660000888777e-05,1.4491886694574e-06,1.4494427793416e-06,1.4491886694574e-06
+0.0004,0,9.9985002999314e-05,9.9985002999314e-05,9.9985002999314e-05,0.0003,9.9985002999343e-05,9.9985002999343e-05,9.9985002999343e-05,1.4486239808259e-06,1.4486239808259e-06,1.4487086841206e-06
+0.0005,0.1428142986695,0.00013330667377563,0.00013330667377563,0.00013330667377563,0.0004,0.00013330667377566,0.00013330667377566,0.00013330667377566,1.2415219436091e-06,1.2412315177967e-06,1.2416671565153e-06
+0.0006,0.24992502311761,0.00016662501388366,0.00016662501388366,0.00016662501388366,0.0005,0.00016662501388369,0.00016662501388369,0.00016662501388369,1.0859200827703e-06,1.0854118121956e-06,1.085665947483e-06
+0.0007,0.33323336498923,0.00019994002398918,0.00019994002398918,0.00019994002398918,0.0006,0.00019994002398921,0.00019994002398921,0.00019994002398921,9.6538584111399e-07,9.651599317908e-07,9.651599317908e-07
+0.0008,0.39988003918566,0.00023325170475776,0.00023325170475776,0.00023325170475776,0.0007,0.00023325170475778,0.00023325170475778,0.00023325170475778,8.6878901309091e-07,8.6878901309091e-07,8.6868734881509e-07
+0.0009,0.45440913696524,0.00026656005685475,0.00026656005685475,0.00026656005685475,0.0008,0.00026656005685478,0.00026656005685478,0.00026656005685478,7.8932390990299e-07,7.8932390990299e-07,7.8932390990299e-07
+0.001,0.49985005247748,0.00029986508094533,0.00029986508094533,0.00029986508094533,0.0009,0.00029986508094536,0.00029986508094536,0.00029986508094536,7.2347009927962e-07,7.2313118449233e-07,7.2347009927962e-07
+0.0011,0.53830005858223,0.00033316677769447,0.00033316677769447,0.00033316677769447,0.001,0.0003331667776945,0.0003331667776945,0.0003331667776945,6.6759120272028e-07,6.6790406276998e-07,6.6806049279483e-07
diff --git a/test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_no_isotropic_expansion_out.csv b/test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_no_isotropic_expansion_out.csv
new file mode 100644
index 00000000..c5d67efb
--- /dev/null
+++ b/test/tests/concrete_expansion_microcracking/gold/concrete_expansion_microcracking_no_isotropic_expansion_out.csv
@@ -0,0 +1,13 @@
+time,damage_index,expansion_strain_xx,expansion_strain_yy,expansion_strain_zz,expansion_vol_strain,strain_xx,strain_yy,strain_zz,stress_xx,stress_yy,stress_zz
+0,0,0,0,0,0,0,0,0,0,0,0
+0.0001,0,0,0,0,0,0,0,0,0,0,0
+0.0002,0,3.3331666777766e-05,3.3331666777766e-05,3.3331666777766e-05,0.0001,3.3331666777795e-05,3.3331666777795e-05,3.3331666777795e-05,1.4499509991099e-06,1.4499086474625e-06,1.4499933507573e-06
+0.0003,0,6.6660000888748e-05,6.6660000888748e-05,6.6660000888748e-05,0.0002,6.6660000888777e-05,6.6660000888777e-05,6.6660000888777e-05,1.4491886694574e-06,1.4494427793416e-06,1.4491886694574e-06
+0.0004,0,9.9985002999314e-05,9.9985002999314e-05,9.9985002999314e-05,0.0003,9.9985002999343e-05,9.9985002999343e-05,9.9985002999343e-05,1.4486239808259e-06,1.4486239808259e-06,1.4487086841206e-06
+0.0005,0.1428142986695,0.00013330667377563,0.00013330667377563,0.00013330667377563,0.0004,0.00013330667377566,0.00013330667377566,0.00013330667377566,1.2415219436091e-06,1.2412315177967e-06,1.2416671565153e-06
+0.0006,0.24992502311761,0.00016662501388366,0.00016662501388366,0.00016662501388366,0.0005,0.00016662501388369,0.00016662501388369,0.00016662501388369,1.0859200827703e-06,1.0854118121956e-06,1.085665947483e-06
+0.0007,0.33323336498923,0.00019994002398918,0.00019994002398918,0.00019994002398918,0.0006,0.00019994002398921,0.00019994002398921,0.00019994002398921,9.6538584111399e-07,9.651599317908e-07,9.651599317908e-07
+0.0008,0.39988003918566,0.00023325170475776,0.00023325170475776,0.00023325170475776,0.0007,0.00023325170475778,0.00023325170475778,0.00023325170475778,8.6878901309091e-07,8.6878901309091e-07,8.6868734881509e-07
+0.0009,0.45440913696524,0.00026656005685475,0.00026656005685475,0.00026656005685475,0.0008,0.00026656005685478,0.00026656005685478,0.00026656005685478,7.8932390990299e-07,7.8932390990299e-07,7.8932390990299e-07
+0.001,0.49985005247748,0.00029986508094533,0.00029986508094533,0.00029986508094533,0.0009,0.00029986508094536,0.00029986508094536,0.00029986508094536,7.2347009927962e-07,7.2313118449233e-07,7.2347009927962e-07
+0.0011,0.53830005858223,0.00033316677769447,0.00033316677769447,0.00033316677769447,0.001,0.0003331667776945,0.0003331667776945,0.0003331667776945,6.6759120272028e-07,6.6790406276998e-07,6.6806049279483e-07
diff --git a/test/tests/concrete_expansion_microcracking/tests b/test/tests/concrete_expansion_microcracking/tests
index e03d5a60..54ab47bd 100644
--- a/test/tests/concrete_expansion_microcracking/tests
+++ b/test/tests/concrete_expansion_microcracking/tests
@@ -12,82 +12,94 @@
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_free_out.csv'
       abs_zero = 1e-5
-      cli_args = 'Outputs/file_base=concrete_expansion_microcracking_free_out'
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 Outputs/file_base=concrete_expansion_microcracking_free_out'
       allow_test_objects = True
       detail = 'under free expansion conditions, '
     []
-
     [concrete_free_shrinkage]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_shrinkage_out.csv'
       abs_zero = 1e-5
       cli_args = 'Functions/expansion_vs_time/value="-max(0, t - 0.0001)" '
+                 'Materials/microcracking/expansion_stress_limit=10e6 '
                  'Materials/microcracking/eigenstrain_factor=-1.0 '
                  'Outputs/file_base=concrete_expansion_microcracking_shrinkage_out'
       allow_test_objects = True
       detail = 'under free shrinkage conditions, '
     []
-
+    [concrete_free_expansion_excluding_confinement_effect]
+      type = 'CSVDiff'
+      input = 'concrete_expansion_microcracking.i'
+      csvdiff = 'concrete_expansion_microcracking_excluding_confinement_effect_out.csv'
+      abs_zero = 1e-5
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 Materials/microcracking/include_confinement_effects=false Outputs/file_base=concrete_expansion_microcracking_excluding_confinement_effect_out'
+      allow_test_objects = True
+      detail = 'under free expansion conditions excluding confinement effect, '
+    []
+    [concrete_free_expansion_no_isotropic_expansion]
+      type = 'CSVDiff'
+      input = 'concrete_expansion_microcracking.i'
+      csvdiff = 'concrete_expansion_microcracking_no_isotropic_expansion_out.csv'
+      abs_zero = 1e-5
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 Materials/microcracking/assume_isotropic_expansion=false Outputs/file_base=concrete_expansion_microcracking_no_isotropic_expansion_out'
+      allow_test_objects = True
+      detail = 'under free expansion conditions without the assumption of isotropic expansion, '
+    []
     [concrete_expansion_load_5]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_5_out.csv'
       abs_zero = 1e-5
-      cli_args = 'BCs/load_x/value=-5e6 Outputs/file_base=concrete_expansion_microcracking_5_out'
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 BCs/load_x/value=-5e6 Outputs/file_base=concrete_expansion_microcracking_5_out'
       allow_test_objects = True
       detail = 'with a 5 MPa compressive pressure in the x direction, '
     []
-
     [concrete_expansion_load_10]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_10_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-10e6 Outputs/file_base=concrete_expansion_microcracking_10_out'
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 BCs/load_x/value=-10e6 Outputs/file_base=concrete_expansion_microcracking_10_out'
       allow_test_objects = True
       detail = 'with a 10 MPa compressive pressure in the x direction, '
     []
-
     [concrete_expansion_load_5_5]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_5_5_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-5e6 BCs/load_y/value=-5e6 '
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 BCs/load_x/value=-5e6 BCs/load_y/value=-5e6 '
                  'Outputs/file_base=concrete_expansion_microcracking_5_5_out'
       allow_test_objects = True
       detail = 'with a 5 MPa compressive pressure in the x and y directions, '
     []
-
     [concrete_expansion_load_10_10]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_10_10_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-10e6 BCs/load_y/value=-10e6 '
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 BCs/load_x/value=-10e6 BCs/load_y/value=-10e6 '
                  'Outputs/file_base=concrete_expansion_microcracking_10_10_out'
       allow_test_objects = True
       detail = 'with a 10 MPa compressive pressure in the x and y directions, '
     []
-
     [concrete_expansion_load_5_5_5]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_5_5_5_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-5e6 BCs/load_y/value=-5e6 BCs/load_z/value=-5e6 '
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 BCs/load_x/value=-5e6 BCs/load_y/value=-5e6 BCs/load_z/value=-5e6 '
                  'Outputs/file_base=concrete_expansion_microcracking_5_5_5_out'
       allow_test_objects = True
       detail = 'with a 5 MPa compressive pressure in the x, y, and z directions, '
     []
-
     [concrete_expansion_load_10_10_10]
       type = 'CSVDiff'
       input = 'concrete_expansion_microcracking.i'
       csvdiff = 'concrete_expansion_microcracking_10_10_10_out.csv'
       abs_zero = 1e-6
-      cli_args = 'BCs/load_x/value=-10e6 BCs/load_y/value=-10e6 BCs/load_z/value=-10e6 '
+      cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 BCs/load_x/value=-10e6 BCs/load_y/value=-10e6 BCs/load_z/value=-10e6 '
                  'Outputs/file_base=concrete_expansion_microcracking_10_10_10_out'
       allow_test_objects = True
       detail = 'and with a 10 MPa compressive pressure in the x, y, and z directions.'
@@ -97,8 +109,26 @@
     type = 'RunException'
     input = 'concrete_expansion_microcracking.i'
     expect_err = 'ConcreteExpansionMicrocrackingDamage requires that the elasticity tensor be guaranteed isotropic'
-    cli_args = 'Materials/inactive=concrete Materials/aniso/type=ComputeElasticityTensor Materials/aniso/fill_methode=orthotropic Materials/aniso/C_ijkl="1 1 1 1 1 1 1 1 1"'
+    cli_args = 'Materials/microcracking/expansion_stress_limit=10e6 Materials/inactive=concrete Materials/aniso/type=ComputeElasticityTensor Materials/aniso/fill_methode=orthotropic Materials/aniso/C_ijkl="1 1 1 1 1 1 1 1 1"'
     allow_test_objects = True
     requirement = 'The ConcreteExpansionMicrocrackingDamage model shall generate an error if used with a material that is not guaranteed isotropic'
   []
+  [error_message]
+    requirement = 'Blackbear shall produce an error message when '
+    [expansion_stress_limit_parameter_defintion_error]
+      type = 'RunException'
+      input = 'concrete_expansion_microcracking.i'
+      expect_err = 'is a required parameter for include_confinement_effects = true'
+      allow_test_objects = True
+      detail = 'considering confinement effect without specifying expansion_stress_limit.'
+    []
+    [expansion_stress_limit_parameter_negative_error]
+      type = 'RunException'
+      input = 'concrete_expansion_microcracking.i'
+      expect_err = 'needs to be strictly > 0 for include_confinement_effects = true'
+      cli_args = 'Materials/microcracking/expansion_stress_limit=-10e6'
+      allow_test_objects = True
+      detail = 'considering confinement effect with negative expansion_stress_limit.'
+    []
+  []
 []