From 0e2b19806c8c0bdbc090054a8eb61bdeaaeaee3d Mon Sep 17 00:00:00 2001 From: Felix Garmirian Date: Mon, 23 Oct 2023 18:08:27 +0200 Subject: [PATCH 1/5] naive fix test small=slave --- src/interfaces/interfaces.f90 | 27 ++++++++++++++++++++++++--- 1 file changed, 24 insertions(+), 3 deletions(-) diff --git a/src/interfaces/interfaces.f90 b/src/interfaces/interfaces.f90 index c09ca59c2..55e7ed37c 100644 --- a/src/interfaces/interfaces.f90 +++ b/src/interfaces/interfaces.f90 @@ -128,15 +128,36 @@ SUBROUTINE InitInterfaces ! b) - dielectric <-> dielectric : RIEMANN_DIELECTRIC = 2 ! a1) - dielectric -> vacuum : RIEMANN_DIELECTRIC2VAC = 3 or 5 (when using non-conservative fluxes) ! a2) - vacuum -> dielectric : RIEMANN_VAC2DIELECTRIC = 4 or 6 (when using non-conservative fluxes) + ! am1) - vacuum -> dielectric mortar : RIEMANN_VAC2DIELECTRIC = 4 or 6 (when using non-conservative fluxes) + ! am2) - dielectric -> vacuum mortar : RIEMANN_DIELECTRIC2VAC = 3 or 5 (when using non-conservative fluxes) IF(DoDielectric) THEN IF (isDielectricFace(SideID))THEN ! 1.) RiemannDielectric IF(isDielectricInterFace(SideID))THEN ! a) physical <-> dielectric region: for Riemann solver, select A+ and A- as functions of f(Eps0,Mu0) or f(EpsR,MuR) ElemID = SideToElem(S2E_ELEM_ID,SideID) ! get master element ID for checking if it is in a physical or dielectric region - IF(MortarType(1,SideID).GE.0) CALL abort(__STAMP__,'Mortars not fully implemented for dielectric <-> vacuum interfaces') + !IF(MortarType(1,SideID).GE.0) CALL abort(__STAMP__,'Mortars not fully implemented for dielectric <-> vacuum interfaces') IF(ElemID.EQ.-1) THEN - InterfaceRiemann(SideID)=-1 - CYCLE ! skip + IF(MortarType(1,SideID).EQ.0) THEN + ! small mortar slave sides have no corresponding master element + IF(isDielectricElem(SideToElem(S2E_NB_ELEM_ID,SideID)))THEN + ! am1) big elem is PHYSICAL and small slave DIELECTRIC + IF(DielectricFluxNonConserving)THEN + InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC_NC ! use two different Riemann solvers + ELSE + InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC ! A+(EpsR,MuR) and A-(Eps0,Mu0) + END IF + ELSE + ! am2) big elem is DIELECTRIC and small slave PHYSICAL + IF(DielectricFluxNonConserving)THEN ! use one flux (conserving) or two fluxes (non-conserving) at the interface + InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC_NC ! use two different Riemann solvers + ELSE + InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC ! A+(Eps0,Mu0) and A-(EpsR,MuR) + END IF + END IF + ELSE + InterfaceRiemann(SideID)=-1 + CYCLE ! skip + END IF END IF IF(isDielectricElem(ElemID))THEN ! a1) master is DIELECTRIC and slave PHYSICAL From 48e129477f637f352dbe1f4673c83809376670ab Mon Sep 17 00:00:00 2001 From: Felix Garmirian Date: Tue, 24 Oct 2023 11:17:42 +0200 Subject: [PATCH 2/5] reactivate reggie --- .../excludeBuild.ini | 5 ----- 1 file changed, 5 deletions(-) diff --git a/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/excludeBuild.ini b/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/excludeBuild.ini index 67ed5c0ba..2169008ef 100644 --- a/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/excludeBuild.ini +++ b/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/excludeBuild.ini @@ -1,7 +1,2 @@ -! Deactivate reggie for now -! Fix by implementing dielectric <-> vacuum mortar interfaces correctly -! Delete the following line afterwards -PICLAS_EQNSYSNAME = poisson - ! run only with particles PICLAS_PARTICLES=OFF From 04f7768a5fb00a61eb12c8434b97a04f39d0a0a8 Mon Sep 17 00:00:00 2001 From: Felix Garmirian Date: Tue, 24 Oct 2023 16:30:43 +0200 Subject: [PATCH 3/5] really reativate dielectric mortar reggie --- regressioncheck/NIG_convtest_poisson/builds.ini | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/regressioncheck/NIG_convtest_poisson/builds.ini b/regressioncheck/NIG_convtest_poisson/builds.ini index ebf13abd5..a0e9e2cc1 100644 --- a/regressioncheck/NIG_convtest_poisson/builds.ini +++ b/regressioncheck/NIG_convtest_poisson/builds.ini @@ -9,4 +9,4 @@ PICLAS_EQNSYSNAME = poisson PICLAS_TIMEDISCMETHOD = RK3 LIBS_USE_MPI = ON PICLAS_NODETYPE = GAUSS -PICLAS_PARTICLES = OFF ! activate also PICLAS_PARTICLES=ON as soon as the mortar dielectric reggie is fixed +PICLAS_PARTICLES = ON,OFF From 2c7ae490eb188e9fc3fdb1102a87d3e75090c2cd Mon Sep 17 00:00:00 2001 From: Felix Garmirian Date: Wed, 25 Oct 2023 12:10:47 +0200 Subject: [PATCH 4/5] fix particle bc for reggie with dummy particle --- .../NIG_convtest_poisson/Laplace_h_N1_mortar/parameter.ini | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/parameter.ini b/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/parameter.ini index b2949c3e5..66eea7c03 100644 --- a/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/parameter.ini +++ b/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/parameter.ini @@ -106,7 +106,7 @@ PIC-Deposition-Type = shape_function PIC-shapefunction-radius = 0.3 PIC-Depo-Periodic = FALSE -Part-nBounds = 6 +Part-nBounds = 7 Part-Boundary1-SourceName = BC_z- Part-Boundary1-Condition = open Part-Boundary2-SourceName = BC_y- @@ -118,6 +118,8 @@ Part-Boundary4-Condition = open Part-Boundary5-SourceName = BC_x- Part-Boundary5-Condition = open Part-Boundary6-SourceName = BC_z+ +Part-Boundary6-Condition = open +Part-Boundary7-SourceName = inner ! =============================================================================== ! ! DSMC From 500174b6e10425c3b465b3e7200cf2a86fc09039 Mon Sep 17 00:00:00 2001 From: Paul Nizenkov Date: Tue, 14 Nov 2023 19:07:26 +0100 Subject: [PATCH 5/5] Skipping ElemID.EQ.-1 --- .../parameter.ini | 146 +---------------- ...RT_HEX_PERIODIC_MORTAR_FLIPPED_002_mesh.h5 | Bin 13074 -> 0 bytes ...RT_HEX_PERIODIC_MORTAR_FLIPPED_004_mesh.h5 | Bin 47634 -> 0 bytes ...RT_HEX_PERIODIC_MORTAR_FLIPPED_008_mesh.h5 | Bin 320274 -> 0 bytes .../Laplace_h_N1_mortar/parameter.ini | 148 +----------------- .../NIG_convtest_poisson/builds.ini | 2 +- src/interfaces/interfaces.f90 | 59 +++---- 7 files changed, 41 insertions(+), 314 deletions(-) delete mode 100644 regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_002_mesh.h5 delete mode 100644 regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_004_mesh.h5 delete mode 100644 regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_008_mesh.h5 diff --git a/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/parameter.ini b/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/parameter.ini index b6d217335..bc2815846 100644 --- a/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/parameter.ini +++ b/regressioncheck/NIG_convtest_poisson/Dielectric_sphere_in_sphere_curved_mortar/parameter.ini @@ -6,16 +6,11 @@ maxIterCG = 12000 crossProductMetrics = F ! bad for hyperbolic equations (free stream preservation), maybe good for HDG?? - MeshFile = MESH3_sphere_center_004_002_NGeo2_mesh.h5, MESH3_sphere_center_008_004_NGeo2_mesh.h5 ! --> NGeo=2 - - useCurveds = T -!MeshFile = sphere_004_NGeo4_mesh.h5 DielectricRadiusValue = 1.0 - ! =============================================================================== ! ! EQUATION (linearscalaradvection) ! =============================================================================== ! @@ -26,56 +21,24 @@ TrackingMethod = refmapping ! OUTPUT / VISUALIZATION ! =============================================================================== ! ProjectName = sphere -!NVisu = 4 ! Number of visualization points -!NodeType_visu = VISU ! VISU... Equidistant points -!VarName = ElectricFieldX -!VarName = ElectricFieldY -!VarName = ElectricFieldZ -!Visu3D = T Logging = F WriteErrorFiles= F printRandomSeeds = F DoCalcErrorNorms = T ! =============================================================================== ! -! ANALYZE -! =============================================================================== ! -CalcKineticEnergy = T -CalcPotentialEnergy = T -CalcNumSpec = T -CalcInternalEnergy = T - -CalcTemp = T -CalcTransTemp = T ! noch nicht im trunk - -CalcPartBalance = T! T -> defekt - -CalcVelos = F! T -VelocityDirections = (/1,1,1,1/) ! x,y,z,abs - -CalcCollRates = T -CalcReacRates = T -Particles-DSMC-CalcQualityFactors = T - -Part-WriteMacroValues = T -Part-IterationForMacroVal = 100 -Part-WriteFieldsToVTK = T -! =============================================================================== ! ! HDF5-Output ! =============================================================================== ! PIC-OutputSource = T - ! =============================================================================== ! ! CALCULATION ! =============================================================================== ! -tend = 0.1!100.0E-12 -Analyze_dt = 0.1!100.0E-12 +tend = 0.1 +Analyze_dt = 0.1 Part-DelayTime = 1 ! delay time for particle push -!ManualTimeStep=1.12E-13 -ManualTimestep=1!e-13 ! 1E21: 1.12E-13 -! 1E23: 1E-14 -Particles-HaloEpsVelo=1!300E6 +ManualTimestep=1 +Particles-HaloEpsVelo=1 CFLscale = 0.1 ! Scaling of theoretical CFL number ! =============================================================================== ! @@ -86,7 +49,6 @@ DielectricMuR = 1 xyzDielectricMinMax = (/-1.0,1.0,-1.0,1.0,-1.0,1.0/) ! DielectricCheckRadius = T Dielectric_E_0 = -1. - ! =============================================================================== ! ! Field Boundaries ! =============================================================================== ! @@ -97,124 +59,32 @@ BoundaryType=(/2,200/) ! =============================================================================== ! ! PARTICLES ! =============================================================================== ! -Part-maxParticleNumber=5!7000000 -Part-nSpecies=1!7 -PIC-externalField=(/0.,0.,0.,0.,0.,0./) +Part-maxParticleNumber=5 +Part-nSpecies=1 Part-FIBGMdeltas=(/0.1,0.1,0.1/) PIC-Deposition-Type = shape_function PIC-shapefunction-radius = 0.3 -PIC-Depo-Periodic = FALSE Part-nBounds = 1 Part-Boundary1-SourceName = BC_outer Part-Boundary1-Condition = open - ! =============================================================================== ! ! DSMC ! =============================================================================== ! UseDSMC=F -Particles-DSMCElectronicDatabase=DSMCSpecies_electronic_state_full_Data.h5 ! when supplied: doQK=true -Particles-DSMCReservoirSim=false -Particles-NumberForDSMCOutputs=0 -Part-TimeFracForSampling=0.0 -Particles-DSMC-CollisMode=3 !(1:elast coll, 2: elast + rela, 3:chem) -Part-NumberOfRandomSeeds =2 -Particles-RandomSeed1= 1 -Particles-RandomSeed2= 2 - -! for chem reac + HDG -Particles-DSMC-UseOctree=T -Particles-DSMC-UseNearestNeighbour = T -Particles-OctreePartNumNode=80 -Particles-OctreePartNumNodeMin=50 - -EpsMergeElectronicState = 1.e-2 ! merge QK levels when difference falls below eps -Particles-DSMC-BackwardReacRate = T ! automatic reaction backward rate calculation (when foward rate is defined) -Particles-DSMC-PartitionMaxTemp = 1e9! maximale temperatur für allocation von kb=f(T), sonst abort -Particles-DSMC-PartitionInterval = 50 ! dT for Intervall generation of kb - ! =============================================================================== ! ! Species1 | Al ! =============================================================================== ! Part-Species1-ChargeIC = 1 -Part-Species1-MassIC = 100000!4.47908080018E-26 +Part-Species1-MassIC = 100000 Part-Species1-MacroParticleFactor = 88 Part-Species1-nInits = 1 Part-Species1-Init1-SpaceIC = point -Part-Species1-Init1-velocityDistribution = maxwell Part-Species1-Init1-MWTemperatureIC = 400 Part-Species1-Init1-ParticleNumber = 1 -Part-Species1-Init1-RadiusIC = 0.!20E-6 Part-Species1-Init1-BasePointIC = (/0.0,0.0,0.0/) -Part-Species1-Init1-BaseVector1IC = (/0.0,0.0,0.0/) -Part-Species1-Init1-BaseVector2IC = (/0.0,0.0,0.0/) -Part-Species1-Init1-CylinderHeightIC = 0.!0E-6 Part-Species1-Init1-VeloIC = 0 -Part-Species1-Init1-VeloVecIC = (/0.,0.,1./) -Part-Species1-Init1-Tempelec = 400 -! =============================================================================== ! -! Species2, e- -! =============================================================================== ! -Part-Species2-ChargeIC = -1.60217653E-19 -Part-Species2-MassIC = 9.1093826E-31 -Part-Species2-MacroParticleFactor = 88 -Part-Species2-nInits = 1 - -Part-Species2-Init1-SpaceIC = cylinder -Part-Species2-Init1-velocityDistribution = maxwell -Part-Species2-Init1-MWTemperatureIC = 1.160e7 -Part-Species2-Init1-ParticleNumber = 2000000 -Part-Species2-Init1-RadiusIC = 20E-6 -Part-Species2-Init1-BasePointIC = (/0.0,0.0,0.0/) -Part-Species2-Init1-BaseVector1IC = (/1.0,0.0,0.0/) -Part-Species2-Init1-BaseVector2IC = (/0.0,1.0,0.0/) -Part-Species2-Init1-CylinderHeightIC = 70E-6 -Part-Species2-Init1-VeloIC = 0 -Part-Species2-Init1-VeloVecIC = (/0.,0.,1./) -! =============================================================================== ! -! Species3, Al+ -! =============================================================================== ! -Part-Species3-ChargeIC=1.60217653E-19 -Part-Species3-MassIC=4.4789897064E-26 -Part-Species3-MacroParticleFactor=88 -! =============================================================================== ! -! Species4, Al+2 -! =============================================================================== ! -Part-Species4-ChargeIC = 3.20435306000000E-19 -Part-Species4-MassIC = 4.4788986125E-26 -Part-Species4-MacroParticleFactor = 88 -Part-Species4-nInits = 1 - -Part-Species4-Init1-SpaceIC = cylinder -Part-Species4-Init1-velocityDistribution = maxwell -Part-Species4-Init1-MWTemperatureIC = 1.160e5 -Part-Species4-Init1-ParticleNumber = 1000000 -Part-Species4-Init1-RadiusIC = 20E-6 -Part-Species4-Init1-BasePointIC = (/0.0,0.0,0.0/) -Part-Species4-Init1-BaseVector1IC = (/1.0,0.0,0.0/) -Part-Species4-Init1-BaseVector2IC = (/0.0,1.0,0.0/) -Part-Species4-Init1-CylinderHeightIC = 70E-6 -Part-Species4-Init1-VeloIC = 0 -Part-Species4-Init1-VeloVecIC = (/0.,0.,1./) -Part-Species4-Init1-Tempelec = 1.160e5 -! =============================================================================== ! -! Species5, Al+3 -! =============================================================================== ! -Part-Species5-ChargeIC=4.80652959000000E-19 -Part-Species5-MassIC=4.4788075187E-26 -Part-Species5-MacroParticleFactor=88 -! =============================================================================== ! -! Species6, Al+4 -! =============================================================================== ! -Part-Species6-ChargeIC= 6.40870612000000E-19 -Part-Species6-MassIC=4.4787164249E-26 -Part-Species6-MacroParticleFactor=88 -! =============================================================================== ! -! Species7, Al+5 -! =============================================================================== ! -Part-Species7-ChargeIC= 8.010882650000E-19 -Part-Species7-MassIC=4.4786253311E-26 -Part-Species7-MacroParticleFactor=88 +Part-Species1-Init1-VeloVecIC = (/0.,0.,1./) \ No newline at end of file diff --git a/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_002_mesh.h5 b/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_002_mesh.h5 deleted file mode 100644 index 78e3744d328a244e2d53a1542b51a445546b61b7..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 13074 zcmeI2$#Ydj6vppMK%xX?a6p_M$`BdUQWlg|8caarCQ{-+x$&V!5-TJKY6Vtotiqj3 zmn>bnc9VZVDax{(hdAS$(sCZ~{BED`()8uNmrzV)(p7y=pRd0@r%!j^zWq9Rp>Olu z3ocrAQAX|YlOvMb)K{C|!>-7;SJGyBAHXO$?i_$lWN)A ztkdlE?;IFKUueeWP{``}v&xR_9k(+tw==S3cwpDwsoUwyy=2yQQ}CC$ooz!;3{1n$ z_9{F5qtmj}DC}sw7lgPT92t80*#Ym>)Z@FQxjpaQ)a|Wyd$sI>e4h&s=CLcUskMcf z(YX3HW-){$z1{XjnCX1awbcaP))owP*xEpEuQkFl%j&^FFqbX%)dzTZ@_Czv-ClEc z#kT(YHnWz)tYz7YOPc%7Uv8spy{tJ;>#A{ibKaO&Kee_w&j&H=ru z@Jrmp!F;EFzc7sYINsl#7yG@|IQWSDerh~o!!MHle_|N(5&NAmjP^`4=`-rWEp$$K1BcN7gB2yTdeyyj+;%FU&T>3jIkmqjvK7Z{E`&M8$N#VcY$xB z?7?LWR$5kDVb@zhKESZ~rrV3Z3-p()I9i>H`3{nJFG#6A{d4gFAMOc;u}J63xwzGM zG`z=JCo5JC}u41m zySmBshX;ms4n7mwQ{OxKER{K+#nTn4E1%{03oUdHpm}5eB&GWF*Y|tgrABT~kMg|k zR@GURt;U1zAW8dOs-OCPcRQbQOLa|Mu9e4SFvn*fXL1B4$Ak4$%4dp8PwS;iikC{& zORBy~7v4(BP9eKPi$eD}Q#(by^HndrQcvuQS0yJC`{GTsz|(fAE<3^{g{SrvmvHfg zOIQ08SG{oYs-*a;mtX0U!c)HLB~!oR3zwagFTBJb7xP3p;kfuG)>dmNUv|`@=x1VI zxOkE&F8jhI#g`N=zVNg^DW39WM?C3;^z|z}wV&d$FI=*)pV%w-j+g8wc*;-h7kXH;dQIK&a&*tpKyMBF1IXP{T05#>Ta+sT=AykvC=rJEY~g5CLPb!hQXoV zb;~`ykM3VzpRvip-#4?vA6q+PC%?|@&+PEWX7=LWKZ9qVFRv0`yIAL|8|(ImjQlT* z;{TiIaAW@QznJKRT#|kwNb#e>CC7u$Bhn>>yVOZg_Kr>!cY8-tzHpa?eb-C z-$gs=lHKK-zY}lvb<-z#>NkCl@J}0-8u?R?@=f2Es?Lclk4IPD58VuN>Lal>NvS^l zo8eb|T($>;{w$Bnxw#n*{~*{l&h;ZM F{{v?e!RG(~ diff --git a/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_004_mesh.h5 b/regressioncheck/NIG_convtest_poisson/Laplace_h_N1_mortar/CART_HEX_PERIODIC_MORTAR_FLIPPED_004_mesh.h5 deleted file mode 100644 index 8b7bca95e6a6a46cefa943a1cd64d01853298014..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 47634 zcmeI52bi5zwT4d;N(e2568ZoEf}sSAArzTOgA@pXP(sb11Oh@5X$Dk`C>RwH0Todp zDj;?&AYw`Eh=5~56iZapsHh`y6;a`S-#PpJS*$rHC-Ex!T=<`7vHtb$Z?CoX{>z!7 zlh;m}bkK%_wj5ORF>qjQt=i!FDgSw*97FGD5Msu=%SpfDmlS_<|Hetb+WIw*4KI9f z$=@-c0Y(2LI2%~N4{m9kkQd#v{{QC=Oqo1>rdL)jdyV^S z>8AhQ^>V*1UAIxEyOi6YHl=>nuQsr@iSF*B=bnAe=?mu9hIDj__~@#QS-&>siE^L! z>(_Wfcuw$jX|1<&1L=MAtF2wzK>22zHg|#J!AiE4Pih_Y|F+6Qz9;jnThB9p@|knb zI%oCsbl0oYG`Fi!uU7KRI{lQnYmjGgd9M59nXzDv^3=-yJ-Oc-Hr932{L^22-rUmF z>aTBR*{BxzN_VTD?^9*J`qkE{zvmkN@pT;8*uVabg55<|USqW!572wj%6z@o*4<>@ zzkegqK;0XT2lVT1+-9}fdJ3y!Xu12Wd7$@_*|pljYV7GenKk3ENxqk-l%iK`vg+#f zhg*`Fo2^=>1=Wx4%~!2kREF>R$f|X0;*_;lt@}@n4k@3POG*QLu0Hu(HQpEVCrs4S z+4ufw^}gUf{CBVGOcn6G=!3nshx>eR8MBA`+@I}uxLHRTvxob9qzn%?8(qfi;XWT% z#_Zuf?~>u+W>aK%xY^7yW)JuI>@sE#_j$Jr4>y}9!^6!MmN9#{&li<3d$`XR%kbuY z6`NgJ{1SP%*<~`bD*`vWvao9cH*;)TDkG2M>k7L$aI+g^X14@x_Wr_d4?MEvfqU+4 zGP65n)HVA|VLgGHeL-e+f8b`GZ?-~49RJa9ATHtPu7>~Wdd$iU6kDr|J%W`ksA;{rExO|!1R%{D7+ zO5kQ&%gkm5Zni^VvjaDCe`ei*o9$89yui)I$jlZ7Znj@xivl-u-)4&gH#?}XC4rkA zA~U-paI+%{yC!fm_inZ{aI=#NyE$;PQ)Fhh1a5XlVYdfv=5?Aa58Ukh!tM;*>}4{u zp1{p6F6{om&6dc_Rs?Q#d0{I9H@jM9)*HClwS_$zxS99G%r749f3uqk8yL9R2W4i1 z12_Aa%xrkzX1B}CIs!NQoXpJMQ&QLLOER<3ft%eWGaDDUnfKAGD{!;>WM)$WH~XQ? zY-ZqQKa-iw4&2OpY1SRM+3#g$^8z<}L}s=yaI?o`W{Uzhs};7ma3iyUGP5Otn+-1P zionf=%gn9`+^nOprGcA`l$qTexY_8!ZVB9MoXqU@z|Fb}TOPRC6q(tbft$@NtS4}@ z*)p^H12^j~Y(?N^^JHc#12vu$sPiY4?4BwLbXrGspkL_XS%t zGS>DkzV7~8zAtFEAKo(`#Z@Q71IqV5t-jA*#rI3o8=0$bN$O$WeVJEePmk-nAH(Z~ z`MSN@h-zQf`1{1G%1SnfE9$%M#c}vQcHVA3)+wFg2bAxRhv_+bLSw7A{?|9AKE|U~ z5Bu)Jl2YNn%lXGRQ|F&r|Cx|G^AA2^>>pfAT(IEmQ|jLZ)@lczx!?sSpJ_UElF@{T zGfqBBKfswV@yH7n&aKs)Z^FrEU+}`Y^Uu|fuIlkIb5B2Q-nor+dwk-8^E_Aw_0JHU zjcg3J{BPgu*O<|-eWR=9m61+U*^N-e-WD_AN%uj>0X~B^I_ee zrSp%$V_oNA9XhS|?3z0{D2Lavp=@V>yz@P*eCDygKZN3eC+QJ zn(Mw}e?t}L{)UO|CgX5_{_*4}fV^w`R{1uPkNy2d>?E<+pRa}YD)#3ePg2*tyM{m6 zcbEBaelvfuwuj7zneCUV;p-LG(JSVCiP!fRV(ugM_n_E`VzIw%6uZaR-?oZ#f6lX~ z%m@8-=zO$+4*L6vY79~w`}?t&_dfQwWtro=K0N|Tb*yr;3h@2SS|VzIyb#Euhle>}gVbnZPR zr}KSBI^UPEKVR?MpRc!jD`MXJ*xwk%zD9BWm!U^6zd2;_guNuxD`&%yNdpGvyefM6){+3C1j9Bb%f?}^X_BT;+?r)OV)-oUTH(BRf z)yb3nIrg~w^ZoDhJ5<;2y4YW*V((S#?^DtpEf)K`UF;|^_vib;`KHQpf78VH=9K=X zi}^j_!_4OU%z0vepHK~dR>c1HRqX2)`}>@9?k)EBaj_%CVth|@D z{moLG`}11aC!N-N_VtOsDtCXrj=m>jfA3R&&Kdjjv+H%l{<;+V-iZCZzpTyu`EU)t zmvVo8m-_RT!|x%#bB6)sJ$EOW-*vIS8&t#l75jUSnD;tf-~AN(-iZBuL^|(t?9ZR& zxxYDL{%qiIf88>FEy3abc9Z#Djs4xM8s5{`-^q%--nf4+Q0#oMzjtb`bH)DNA?Evv z>v@ghcgcJ`d@uNLFLPzS2VIwE@Kl|9Ph;KF6z?L7b)B6DCRwl>%CUY`xp0lzT!h=vF=%ly$7-G zn>2SDu~_#CF+U%%?oCbVUN7c-iFIGB*zc2AccEghBi6lKbNzmfb>AT7nz8N&#Qctk z>%C0O_d%?Cp5p1USoa*oUPr9^I?eSvKh}MT;!!f!bsgUWA8b;0iJ0GYvF^EwT_e^# zMX~oK)_slU`h6ekURv_S_1+@p=Q6JMrHZG?T$lTKk-)%B~ruccpULodQSg-j7GWQ(c8y72{ zB8zojrP$Xk)_tX7?`f=iYm>UKmS%*E^}6oMWWJWM?j?#3mc_c4Dn3XS>%K;@_aN5& zwCZ}Eb7M<%zPl{eJ-@`U?#s)*c)kzMf33`4uf)2qQ@o8V*8NPAx)+Ffui`$Rt=RWz ztousoynnIo>lIIu#k!X%_PZ(8{etEW5@W6Q=Zkq?V%-Ib{p`oO=Sb%`*1cTuL|Lr+ zM#ZCLvF;a})IC?s`x5J(sd$(y)?FxFhghup2F0Em>%K{Ghb-2;Q*+l7i|bu1=KYIx zU#Qs6eyn??Vt?K0x;(#E=-hi6>t3aJds(dewI+2Jh&g|(>(6C>2E@8oD;_V4b>FPm zdm8J$Mez=@Soa>)7$_FkdzM&-SgiY2#r_P9b+1v}C5v_6rg%?Ttb48E9c8iZy-n(# zA?AIFb>FUdq%78bhvH6Itou&Id&pwl>lBZW#kzl24XzDqK2^;766;>C*k6mqx=R%w zD2sL9rP$AYtot(M^ZU94nP-4}>?Ut-;NEA~4()_sp+=Z|%7P`q2=a`o$X zqrB&`Ua#-uQX_sY94F>IrLOyXMAlzyfXs*IzE^gDEY`hQ@hDkb@2fP|?^xFBxhIK@ zkdJkb7xSLRx_@d?_kGgrFN<~GuXtBktos4QUPoN->B{GODAs+On4g(g_giA?ipBMQ zQ1O1USocGUoj=ySNwL=v>&{g^-$SwPQDS~(V%={xsrzAR_LarDA5lC`7VF-k*fnF_ z!<5g@d93?pv0-Ad?)Syk5sT|RTFm{qt`A?gkID9t#kwC??0YTNJw$W;oX5H|#e9uo z-TRx=oh9b$7VF-sc&seeeXVrP<+^^3+}|`YKX0+_R5AA#>pmprwZ`?%7W4BF>)x*T zd9qmdlZr>l;(8C(T)$sp-9yE!iFJS1r0xk~zF%V9I~0$R#kz~6^XED1^>evW=5#^Hp?m50U=14a}EY|&uVy_|A{j_55DeHCJU&&%kpHCF?HIH>arPzBP>%LIB z9mQhZ&nezp7VCaivDXpn{#0|h*3=&_=6gKW{Y04;>z*c^pTSsnx#GQKvFLzQ<$T7fI)LPOSTR#a*&k_X~=>2eIz+HTMKD)@px%*hKkQ_sfd? zK8SU@r1QSSx?fS;DT{T#s@T^&)*aWR?tx<7t62AIiv2!_btg#Y{fl+)QtZ#hSoiCS zeUHbwV>I_fvAEuS#k?=E?mdeAT*kVSr1Sj}>)xx_HDleqq&rB=b$u`SejY8$b)P5Z zeTj9ysd%g`);&l%KOeE~w-o!cG1mRIV$Ws0UgL914cGBuwwsvul)CP3lySte8G}av{mg`Oz8zo(=`+dd!jEZ%?quBZ4 zdiT*>U&~l`XEEQqvF=V{$BV`F&Jc61SoepD{h1u={y?$!FV@{(bA9i|x;u)^5Q}wp zXj1oZF+bC>?h#`CjB;InZhB9DBJ(wmbvvZ%6pMAY6LZh8?zUpbiN*EqBIb3(x<`un zGdb4%sbcSathxpd_7{_W5oQ~80&5+HdV}Zc_(inW=*U+ zT+I6t>uxSKM=Y-QIbvhQV%=lKb`pzqS1R^B9@o2>bbi0Yy4#C+Ut--&#d6&}#C)H| zx+jW_5Q}v;5j$9n_4>ISChHdS^~Z*axfkl1|4Qbbx0CrW+gmz6gRyS6*p6bc?(a*? zdR=#jtXrI#_JhTIuf@8*QS7~sb;nBQeTsEY5%XtbtlO*D>xgwXmhM=wSa&0_F0oj5 zg<|i0th=9de)ePCd1C%-jCCJY>~+Mt8#bxCf!KKIV%?uB_Oli1c1q`GKh`}{tV1l; zeME7WEY@9Ly4hlJy@SLih{d{pQS4_c)*UF_M6pAHbuHv z_iu`4%3|HMq?;rb>prU3HC)$+*YhQr_a)XHpt;ACan{_w;K8rp<=P_d78V8SgbohI`37i>;Jo<_a)Z- zyJ{RQ##-(FCg!gXV%>E#_YkpI_a&;kwOFjXu5^A!#JX#Xc^$FtqfP35LmIyiTz7!X zy{xCX)8%8`3srY3u~_$|Vq1#2E}ut(G}n6&>;8}G9wipn`%sw|>u#XAez(WEe=YkG z>t3X~Ta@|KeN3^hd93@VCUyTT=Dmvhypifo6N`2KA)WUx*1cGDhl|DaZd~S4*K2&F z)QF!8zZCbonDzR7+AH0WV%~Qjp8E$e?|rO$iL{%G#kxaP*Uw^JE?_=FRig{02 zuj~G;6(Tyf3lt&zsbJi{=iK zF4ldk*ibRouy{2 zAl6-5IzQ8~?(davnpmv6qDkFfiFsdQ-R)Jw*CX!p4r2bSjdg#dxqdEV-4W9H=cU5+ zW1jxoA?7{cFxyVd*OJ4}^iE>O$;bE0iHZly+#h>?g3dRV#r}5ITz@vk{=TpNW{JiA zMoO3a+e*xPz~TP35%cxnaDTgs9WNjI`$^f;xPR_rBU$Y4IhuR4SnO|<*hymU&-a3B z>?W4`+e~b_#yH&Ha4|ne9G<(2n6F#x@7v`X#r}4e&ifMk+e6HsJ+Z$%#pZ~`{zi-C z{sxQr-O1tphKWrB$h*HC#eCn#{=DxS%Hs9iOFHjW>~C)|>h{-|YiuYsMc#+=o2@VA z&n+Kjw&#j<+soMJZ z_chh<^>BYato?zQ@3q+9zMAW2Gxqlt>Aa`0zx|~1I%0qOH|g(xq!};sL4SYMx!-j@ z%xrg;J&67NSj_iY?C${8@b^ctzvn4lPZqE5^QH4TVt)rV>FP77W*5mc%UryH(olwhhl%7VqRN(ejg}%5bI6|>yW(u%=|*$ zpJg0=j(#TF7a;GoOw`;t^07bHTUQqQ+f}jO)v>=xn(IA{{e4gU`5ub>O>WZPkEQWv zJcs*xK(-G+-u)e-dbAuTi|dQ&i(bvDn`irSo-*{cWw-^ph5d4+-ngQ2qYS(^bRwLF{j4sLTELyQFd5KgbT0`EX6!a@F-c z6xZ>wvM=%a-XnIDbg@5`^((*XPp3M1mo@h*T5$AvoEVQTd&cpMw~FDd{OZm}hi|O< z$@P;zbM%b2^XSQw^_ip3x>hmzoCkg8=v$3vj-EXDtb?P+qr+#OV|1-NdOUT}XHLv` z*0|51AM;>(e%+ce~;WH=BKI_|!qsM1GG2{5mbBwN)M~`o{A8>T!gA?Q7 z=(9(cIdL8*#-qbi508!*j~)+4pLweoJ@e7E^62w8oILmyosAJ2A+arWe?e7Y|1mh)Qd8~Qw-J~Kz3J^IXv$;17CK6CW=tk0a7 zJeB_!?N3jOeS_z95clbee(=j03i6<5y>RluiL)M#M~@EA{#0?Be$eNB;N-#Qy15_n zWIa6V(Y5mE87B`h9-j4#*261b?oVR$ zcsL%uvfNL^czD+LiBnz1>7P90&AgqSdE`M4?`_dPJnPBB{C542C+p$lCl7Jf!z*9* zZB6@OKKXfm;dpq?pLshy>mm<&IQu|9nP)wD(6{R+=SL4GKb*Mo>He(drTx!(`po)D zSMIx97d@Oj%*(u;p8X^bdUzF=eloAdTdcQTKh^xSUU=nON1{r)}yO@**9YJ^q=!*p7lMfGxZNgPoLDqXP))sK~Fz$ z^zceYUv1W#>l5S2PfT6vWS;emqo+?edU(#C`Ttd4^|6}!ll$TMh0_l_=g&Or*@tRg zxn6Mec^|Sq_k*51=;;SeemHru9^OZv?l0>06pkKFUGit1_2fZc^_lX+bN-*F%>w%*urY`w2 z&w9qu(+?ayJm=56ot`@6K@V@|vz|Q6qi;BRI68dhS)a%AJoIqt664`npYsz}zPx{k z$wS`E$q(;M_doOC=;1s^#CUktlLuX{iyls%tcSPLGcW7m)FH;h$qy$FF*-PUI6mv) zS*3@lrY`zi7d@PLSr5sdGY%IE8w)>Y{eCx10w@(|MpbL26OJBEo~(zr(leiV=!q*| z)|>V0Po+!O5j{TVhm(gG4<|pIeIRBWPK>8+*2B@^$%7sbM+YYlF&>`vtP36E@G36* z+0$bG;aN`}^5fy?;p8F4!`tZ@&w4od@R<|i$%7sbM+a}`vz~om9?u6HJ)F8(56^n? zpl2Lj#sB2>r7rn%pX7(5!)KoLIX~lYVmx{AnG@rC)BQ|8@Jd(qhv$g==;2vU9`x0C zsY{F}Pu9cH;o<1p^}}_~dN_G=9`uY;m$>q2KQpi5vY+e^oVxHTPWzwr^po}EL63){ zhm(gG4{xVuJnP}D_8E>2PagWf!_mRXQ^jT9;O+F)c)5jI_Jm~Robnu>(Kl65a z#>qpBCm%j@VtmgU_CwxQp1h2+55#zQ)~{(l8sYFpi!)aD3Lo zvz|QY8HW?&$&>YPblIbCrGry9>*3@@PmG?pl_yWu!}C1&>gwq$>*4M4GmpCH;p8F4 z!?V6s49}jre?}o*u4-Q5=FMXiToR~cHgFbWg_^i*Ixbmyn5B2kUGq3WpPc8D(CwlV0 ziScmsc*gOW=NMfpkG}FP`ez@|kq=Ichok2@GLFw2JsurCb7JyP7k%dFE8k|l%*T_L z`I!?lp7oid$J0N$%&WNEZ&}Ygc-F6OKj`q}CC2Cc)J0$Uvd`$5pY@re&v~ji?F0Ip z2Yu$~8E^I6!O@ck&iu^Lf3XC-ZVYnG=&I>!0NQTR{H%L;kx${u@vJ zJ5c_+MFVC2J74fYI``i;g7Y`>a_qnD=f8#Mzv1Y=*_iVU5gRH)=UT&LuI0b~>9I{^ z&f&jh1oz)l+B|UQ9xijwTgcqoma?q(+P0FR%k@0wzrW=i+sNF@wlep!U6a_ob;#V? z_A+aCkXgHW z8*5mX_r^7{eh;xdWuEK5J?Fo5=D&evt^dB1^Nx|Ba~}WgE_ALlRyIzyk1Xrq?qy#Y zd8h~X-`3hs2KV2}^55Wsd#(OESa9!6r!4E~d%VWFWLXdQUQdvD%@bv=HA&_*PnLDb z4wAuL>tI=z%-nOQ$lUW(nR}Zi%X+VQx(r>e=RG?_hR)aGP#HS!?F<<@UyH+J=gAsqh#og*7-3qbhCAStPI_8IzL{9?gX8mC_{IW&gaO` zb?f|O8M+sgvGwJh%Fjyy|m=W z`}Hz$=RlYDn|Uvn#yzsnTrc=V8uJ?AT<=$iEtXv@d!;PrgTGp1uaf0_mxx^|TOxan zEcXki4(p(f=f1YgeVu&XFXuF0EPK6-ePqAk?)@?uoa+U5?{AR7x!&+MYV2|u{7pK) zLI%G|=U2+$Z`S$MvaIJm+ec%EG||I(-;w97(!NCo-%saLW$?G@{2CejK%GyMT`POL z?47cA$g-Zk(dByAX^zLPmpR8$nRC2L=9=%8IS0D;=-fG$%B;hq-(So{Sm(!dZq3JKp7RNrHMh#Ff%|#&-opJndyh7fx%S2~?@`td6?2d1 za=jrMvu4AhxlP`hPs*GNK3M0j0rz{xJ>~dzG52_h%su9OpAz#L(M^~*=fd4Pp0S6P zf59_*`HMz9V=sP(ul_T8`HP;}H~)pt;H%$Hp8^~MbQ0tKDb_#mQU4c_(d!ROKPC3u zIG=xcpjP95A>7LMENkS<<8b`yj-y{)PhRxI_@2%d{Uu(_UrkPa{ue3LeDC-balf3hB)>%w!NaP%uX8~X%DpE-3ir_Z*$ z)ajF_@IG}5e`wiLuOs_}o*0jw7*AbdJpB{nv!3;3Jv`Ti=RV=+iSg)(@$^HCrypYc z^0I%#_&(23(Pe$L4=E4)!4^EfAJAt`KRFK^J@et{GtYX)vmVYq5aZbgVtnTGlk>ot z$8*(+qpQZ3mFF?FeM=eopr^i`bt%$&OL%IhEK zS*P(eQTf8CU;XC;`Z}p!19UewxP~8~agPk?=&b*i&A)3Hvwm&N6aL~%$J(ME@8945 z?xuee*U%5>cZa{p_feFMU}(u=@dLe=%x?T_u`%|a_;26+7+G(lUt>Adzv$wxJ$?8m zH@s}*z|rX=dI-Zzvusd&OqJfoAhvGql@8H z_E~e!|Fz?-?APWwrhDoOUAwBf)~{x&s;*jy-QBC-kc0O-aB$W9H5)^G)K+z#yQ=f= z#Xg^Drs-Em%MQl13hlf$)@l08V{@f-FYhK`{O)w zw4P`1=7ahjaPWWqJXJ-%e>%^O`}OVj-#O2Mt?kl7ZR-~^RTceoKk&bIp5g-kW4`B^ zKCYbx?|0B4{fe&s>*H&gv(i^({r=^AMR))8^KCMHaagz6Tkmtz|KmD#oIdWEi|0X71(bvs_g*e8>OT-+FPA(qpIpIFA9fKI(h@$9av#ACG^uSfGEC zK4TJR`!@;yhw=Qy>*cVbkIL)nAJ^6N`^Dg{n{svjdH?*o`$c8n{jVP0AO!A3ACsO*DDe}cO)uYwpzhD1(&T&$hriVd2G}YmGgMK zQNikh%QhoqJ%Y=+6|8q~nJ-zgdZKcE*=_|J5M0)mkTnFC4J=q=@MKNFtv!^G4I?UE z*--@>9$a=JAsZ1~c1FRPgUik*WTS}6`DK?CY;`Ow{ zKDcZuA?pxa_6s4a2`=*$Le??3Y_5WJ3NBlKkaY)&Wg8Z3KyVo^BvrD8;Igd?))-u7AF`(4vR(xn7F@O)AsZfC*0*3Ig3AUG zvgY8jp#>WiTy_*88y#GBe8I*Bm)VP~CAjRYf{hC{~V|xNKCx zCI^?=?0?y~f>j5XO(10LgUenaWF3Ob-Xvr-!Da6eGQX#+ zTz4{imvstW8IP<>@XB~(wZSXnk<|sS;_KN9)+4yA9UVwOcBxD1E%a$!z zLvWd6k~Idej8oPWyc1tLCfTsyvW*HhJh*H#LN+3}tXsjFgUfmnvQfcheF`=@xU8O# zjSVgvP_UNZvIasnF1W0*U=xDNnh4p%;Id%_n-pAj79pD)Ty}oJrUaK=Ldd2Cmo*ox ziue4LeJHz*kW~kl-Au^Z2bbMJ$T|dfo|72wD5! zvU3Sphv2digsdjG>_$S?CAe%fA*&588$-zIg3InFWIck*9wB7C zgUeh$S$%NXvxICwaM??QtRcAUbwbt{T;>|enu5zdB4opY%RVP$!-LDdA!H+h%Uo|+ zb8y*jgltrB*{lT{9b7gSAsZW9*1ljZ!DWjOvT?y>OBQTGaM`kiY+`U($AV1?E?a|; zO%5*WT(Bv@W$P2NX~AW+1*_s`3f}+S|FSIzS#@w(H$v7vxNHYP)*-m8HzBJDF582U zbqp@+OUOC}mklIjU4qLR2w81#*`b82F1YL{Le?X=?07=fJGg8ZA*&BAJByGF2rfIH zkTnFCT|&qjgUgx;SyOP?b%bnKaM{g-Yym)`yTa2A9V2Haxg&2q7C0T-Hd)nuE)ZA!MV1%T6R@ql3#%BV=QP%Z3xOmf*4r3E8;d zvdakBgy6D~gluAP*(gFbDY)!bLN+@Gq!CAjQ9LN+b9tc8$O@jIi+{a^M3A*&89 zdzz574=#IwkaY+yn@GrNg3I0}WF3RcJ|JYBg3CT7WL<*GyryNf!DUkkSzU11FNCZ| zaM{evfvk6MSv4W64=$URkPQefTZoV~1eYyN$QpyoY6w|VaM_B4Y*=vFYJ_ZfaM{{~ zY(#Kb7edw?T(%J*8x>r(86g`TT(&hK8yj5KgOIfZm-QlKqf}h2bb+Y$T|d<^(JIB!DV|8vW~%J zeF<5o;Ie^)tV?iN10kyoE<2Qv)diOwMaX&tmmN>YdIy&cBV_f#WoHqx0l{VG6S9Wj zvP%eAV{lnBA!`aQyN-|z3og5vkPQzmyMvI82re5-$eM%89wcO=g3BHwWTS)2#uKu! z!DSN&Sxa!)D}-!ZaM_!LY(j9^dxUIaaM@%+HYvF5OF}j|xa>PZHYK?1CqgzYxU9-+ z0#-HGANyZ6yVnG)I=HMIA!{F8wg4gP5M0*5fYk(-EltQe2A3^Q$T|g=twPAU1aD<< zS#9u6uyqMpU2xfkgsew!Syw{VJGiWlkktp5bthy4g3EdmvWDQYK7_0>xU8O#H3gUT zCuGBd%LWm$;lX7?2-%3>vPMGI99(t`AsZE3b|N7g9b9%AAsZW9Hk^>P1eaY%$i@Yi zT}H?z1ec8@WD|qSMiH_}!DY7+vdO_^cM-BF!DaUmvT4C(Erg7}`Z(kMFMEQJRR@32w9Kd zvYELD$a)8tRTHxM;IerM*?{1(g$P+gaM|L7tTDK(hLAM{m#s+1h6R_cM#zQ-m#t06 zMg*62A!NqO zEB8X#=-{%s2-(=+vi5|mCAe%6LN+eAY)L{kA-HT=LN+nDtRo?t6kN6jA)6dr)|rq^ z2`*cokWC9Nt0iPr^Zc>@Wm^!k>fo|&gsgpV*$#xPLvUGdLRJ%8wg(~W7+ltukaY?! z8%W5y1eY}svfALXLkU@3aM@9WtVeL!@r0~*aM>_IRv%n;79krDTy{PoYX~m8gpf4` zmo*czrr@&c2-&dUvYQFn@ZhpL2-%3>vay7$Ik@aWLN+S6>@h+%I=E~+AsZW9Hi3|} z1ed)+$i@Yiy-CO>1ed)>$R-AtO(tZMg3G=nWRru-z9VE)g3Ep)WYdDns(ej=Rn7ay z{+G?}YXYn~xU3x^Yad*;03quTT-L#W)dZIMncvcTy_j08x>r3A|V?cTy`2E8yj3UoRGBymt9E6#s!yMM#v@v zmyIN36NAe}5wc0aWw#Qt$-!lJ5wa=4W%m)XX~AVJgsf`5KlZ=u2|`vKT=p~}Yad+p z0wL=VTsD!A)dZKlO~^V1mwiCUIt7<~O31ndm-(JfRvTP4m5|j1m;FM>dIXov%xj>m zcW_xXA*&BAo0pIc2rgTQkTnFCEl$W9gUf0NSyOP?iiB)faM@~vYH6a@tT-Jk-wFH;-B4p!&%XTAV6N1b3CS((X%LWj#Nx@}< z3EAY}vY~`*N^sc`glt-HSrZ}iFYgcHsrB#wz)mJ))xl+F5VH2cW#V z*%gGWV{q9ugsf9=*^Pv(OK{m}LRK4GHinSZ1()4V$a(~qJwnKO2bYZ_Wc9&i&l0i$ z!DTNIvWDQY*9lo;aGCF6Wlh0l9}%))!DXKlvf;sH-w?79!DZ73S#xmNZ-i`AaM`T9 z7m$q(E}M&xjSVhqPsmz=%N8MITyWW&gls}^*?WX+VsP1HLN+P5>`Ov6Ik@aQLN+D1>?cAtEx4@8_e3!MG2)E- zzif8j6Tzy3%i0mL_Q7Qf5V8)zWgQGyO>o)Lgsfw5+46*}Q*hZTgse;ORtA^V2JZx0 zmyp#3mu*PMdIXnsC1kyW%jyVOeQ;TKLN*|{tS2FB2rlbG$Qpyo>Iqp>a9MvsHY~Vo z5Fr~LTsDM|jR-DlBxKFOWycV*QNd*=60*_3Wv3CcvB71-30X^U*@c8`TyWWCgls}^ z*+@b*F}Q3LA)6Fjb}J#999(u6A)69hb{`>|7F^as$g21^(1UpL{g3PkLRKAI_B0`D zA6)hVA?pxaHj$9k1ed){$T|j>eL%=M1($tF$hriV`8|!SHn?mmA*%~6`-PD82riqM z_kgnA!DZEitUkDGUP3k?xNIRp)(~8_I3a5cE~_DAO~GX=60%{zWvdag;lX8V6S5J( zWnBnab8y*4gltrB*=B@nba2_$glueZSr03kWCCO8$ie= z1(yvbWRru-h7z(V!DUAfvT4C(O@xdOMrhst{PU|B-*?HXgAXHjI3a5vTy`NL>kwRa z86m3)E*nY6ItG`GB4nL{%SIEjF2QAE30ZA$SqmYn3od(tko5>Idzz5-4la9%kktp5 zy-COh1ed)>$QpvnCKIy8;Ib)%tSPu`8X+4NT=p9w8y;LXEAL@tBZAA?5whmsvV{oQ zsNk~23EAl2vKm4*Hn^-KA!`XP>qN-L1($UpWD|nRY6;oI;IcYGHYvEQJ0Y7KT-KA2 zO$jdRL&&BDm(>%ps)hg9|FZsstU9=C5Fu+HTsDM|bqFqNBxE(gWle;vV{q9pLe?p` zY&aq75?nTdkktm4H50PB;IiuoS&!hdn+aL(;Ig|2S$%NX{e)~laM>e-tRc8;93g8A zE}KBenu5zF60%{zWp5L*;lX7e5V8@$WuFtW=HRmL2-&FMvY!aq=-{#{-XqJ#2A5S6 zvXr2Rb1eXmWWW9sS zh7z*+;Ig9#*?{1(6A4*EaM@{utTDK3I3a5aF1wJB4GS*2jF1fvF1v=1jR-EgnUFOH zm)%9kMg^DMPsm0Gmpwwr#s-)9y{xPyxa?U%HZHjAB|*xp zl}!#_`My^+C3xlgUfHzZmG66Het^$^`ngzspTL8EE;9X5I{u&dT&&*3nf!>28IS*4 zJ{Rjhoe$>B2mfq#`lEC_OY!+%|L$|<{)5j8d#RX@T7MoKgMvY!Tk@c^yLTt<*^^!x9Ncg4(WUF z^o%VB9k};igQQz+ChEFr_q`6_0}Q)v+Vk*+epOYTuj^hz4&SHW;GukUW9#wG{r1~; zz|iUS{`sZ{9^xOig37(5hL;=A&|>n{mG@zNj4k&)@IvE9RxDF30-kc^Qy@{QtiG|9A%6n|_wB z>-=}V_9~G18uY8Y_}W!)%@-{$G;?G=cL?V3&(!(+SHg#1{{*Y%*|*(3WIw|C@cb9T zo1SbMOgDc$rS&1RPKEvXknz`6|J0vF#`BQ(qt*crnda=I+jbNVVa~ytdVfObK%zIJx$FJ^N*UA2T zSo$;o7!@L^8J(J#d$aDVEy9~$nnCs?` zMeHxVZ&SMgEWK~R_z~>CjCXmC-9zbpgS@3+>0Oay`%ACCmUZ8!_syc`^u7jjzvVnn zgw@0JR<5I0Va{QCUxEEtnrG+4HA?R!@;py^Uxsl2-V*^7JB zxyt??!@;lm{-wW1Vc(bfTaz5emi;|Wp5w~?yf2me+t|A3#r|rEm4n&eL!574o-1o0 zdl2T{&-v@hvE#_=V^fY@-|X*xY8_klw>ih<{?rp8bucc5PVpHwNaM=Dy#aW9L4-x02`US$fyv*gcfql{of(klx#zMD^xgn-ebYM%HV~HH>tViTrgsO9 z9Y=b7KeH5(-fOA#`bzILumxf1?ZvTU%Xz+70hQNBm zvcF?su1C(_${hRJll`5^v2&X1b`Z5o!P0vZ%)OT02AKD;^qv57{?fYw$G+C4_Z*HL zTY3jmyF4tt*TKBD(t9??&VgPZ&eZ|r)$p8NYi1$Rdp^hI^|%0bK+L`CaUn4{n2!B$ z_`V_gI}qkvWq)UI?3`wQ1IY6@`y0WrbCCUA%yGHDOJEJri~U_n92m_0`r=rY=j?AE zSPd-uJDFqmM)tQidGD3>^LgYrw(RfnqAu5MFKU;8r}s3N_nh<|%dvYSz4he1Tk>8- zj{6|JBRO_#>D`^$jmQN_j-=Y>oE#;P|UsSaRV_V znCsCChwt~Yzdc~fz_P!~IJUR!ZzuBJDfM?VIU5n#-z^-M`@0o(SoC6lw-JrO?5_uo z6?xA7wu3qM+22JRyHB&fKIFL`xu5UgcmpE)yPad_G}mogYP~09x!2NrAIBRL>HP=C<@LB9 zc6iLa>+t|_L@?K*D~_#r&i=NBIak@=XpWuJ>~9nD+-uq2LmWE?*e;dM_m+bFGj@=vC-v;EpUfR$1lH=I2zsERs9Jy}mQ@af;z3ahz zJxlM^94|wpw+nf%mAp@q@UpBy`m^sY16vK2-lvPY^!}ARKZ88i<5}|D zjxD_}a$MOj zeA)V5>Ls4HEA~c(z04Oqg1KJ~BJ4N&TLs4}rT$(eXLTa`dx~S{Ap7e`?IE!2?=_Cg z{k;zB9=*7q-ypUPW`8db&T00yJdT%3{k=s_XCnK1n`8G-_O}eRjj-%*636BK-i2)% zz1ZJ-g#V9_%ESKNCtSDeZ)F^=QSPe`IM$i{eaNwUE&E%FTIV_Y`-tOme;>oTMlbd^ znb;(l{p~?4NMwI&;&4v0zfU+`o5=p|;&^!?`&)!s=QR6U0_Gg#zTcJO1&H*nO`hjU z?`IsZMWlBu$KG?&yAZWrOX*z}=3J%sQ;wa}^!|-J_i1{+;J6Es-Y+@!H7dRHQ+pIF zz3pM{we)_?vGbSSg~@ZS()$g^e;lKX6>$f74(a$6j;){YcaXbN_8a%u8f{OX6_9WPd+#ybh85 zP2t!%$o^)d_HbDC_Y246{(gmR6usErZ^Yk&+22&cb7g1_J=AL7Jp2Awc za1UjFzyJ38?@h4m&tKD*`~B_%I|tJ-8^^9&_BT7dxJKFEPk5gy^*0B# z?zQaiXY!oq?5~= zmOSS)`}+~*94tck(CeSiMWlCLyzb}penah(CGUKA*CQ{zQ^-3SmR^6!>|CXH9%>zP zdKV+fKltMq!?rq_pa z2C#&z5Xjq{yL70dCvYmf;lhQ-$ESELS%nESvS8I%Kld3{I#&` zuOqD7pTCdu{^X<5pTD>B-d17qm+7UB=j`uOnAb=4w<2}URra?9{rNpu_V*IC50}>0 z-))uq^LJpyU&ETw-~8|++l1%r?@O3_GW%PGTIVGDdyPEjKKpwP)>7)v z`*pd$*4G%gkIMY{8uR;fG^4+#IA2#-_BS5387%wr_qVmM?C(kPocru=9PFV|f9ti; z-}>rlH1Ec<(eye(nb-@~v6Oa1xYz&XhK zMhkg+z|#8=>=;;jyHe|1rS}2ywt%Jge%J#g@47|L>3xtq$CloIz>Y3?{d+j>q4eHI zp7-|jZc+4?*V~pb*DbwcskP_y-UB-dmgDvBRC#TscMN%3!_w>bP>v(L+rV6p^xj47 z?y&UU2|E&&-a2X>TYB#xZ#!6eeedU5rq}O>%h!kB4}VLZkIMe@_fEby_8~KF#jzXD z+27qT=PLX2caqL&_IERRzE98o{zl#=uoU825zdx`sEc+Wp-d3>e&))+&j_hw2Sh>GGFyE*9sPwlh?5pW$Mt|27bCCUA z4RgO_e-Co(y*T?DNnSTt_O~?-$Cmx=4s#s2udbrjy^-E#*kQ2r-ox<D?P1alO=iulG%l?MLoYU;j-xoUf+27gZ-Bs$(pZmKHvcCqHbCBzHCbc`l(t8H% zU|4$hq1L%d@9E^-S@QZbF7FZPJrw3#rT0{7cYvjL7;FeEy?au-H7vbnk>|YRe(~=? zct1$5Ki@8|hd4d$F?f4foZHJJTfK%VoQ{T)r7 z_tWg}7+ASKf95@jJRg<*nqcotM>G067DsPb_IC`-InDn32RzPw_IEUS?zQaiMDlir zWq-%RoP%7qBdOgPmfj;^2g1_p?|Pi8^d3&$j0AS@@_5d(-X*ZZ0S7%R$h-YVXu?tqp}`n!Csq=W~|3SICg?%e~mE5oBbU~ zt$QQ;Yas8IQh(===lW)U=fcYUod=spo{vg@=fhr|j%M_C3df5R_BRV*?*4e4)9mkT zYMuM+ZyE-0@0!!}|u=09b344J&AC>iJhCM$W z%~+4KIQDzn?5`hQ=QR5}j#~F=_IDt8&U5y66?uJN+22T5xxcGn6Ug&X>F*lYbJNj` z{^sSFssES$cEwu{%l@vT*85uaw-0%)NA|Z1?0QaUf0VR(M)5_tzZ+oB@aH}%{oM$A zb~>8T-y9r!o}9m(aO?@o{%$78@1wH6J;`$|v%ei-{(pqCzgwv-_jeoYDQbOG`Wp>< zdODiXUr)~Gdf1;2bGL)}y;b(NA+=p$*% zC#Rzs{q@1&9Atl6!*+*de>+g?>uUCQ0(tJq?C&1(99#A`2Ie^Ox~!wN9+uuMVSQog z?MAKdBhq^md0wmO-2!%XY29w8*11aW=H%@OOYdf|ePHRWrFK(Tdbc9ac}ef4u#qM2 zP1L&P>0OQEYC>=2zT1`By*Wt`s z!O}aHTE~{&{^WTdO7Fwu*>ieZVD+%{?#OX_BE1`Nz5%fGK1!`?p56_~8wgA92Cyqj z^ZX#S&Ov%RaXcrH-u0>77na_CQtO(hcVqJQf~B_$?23~2erjE#^sYyq^P<;>*WBM= z{drFBcxsP3yEO%6*7GYW0G{n{0?fI}{+6V_{b1SO z2jqE8O-}+_UHSa`HAdriDJBYU4BNb*I;@VB5!|KdOs)6IY{paST9(5y~aBb>0PiG zM|!`cb})JAou9l5*}Lielswlcylb-vy=qz9-Lh%l@Xq>?ix1ri&E>p`WG*b@_Nh+ zyNbNZ!}XX2HnIYEc0Fc=9mI3?H#ZLFFZ-K~JYQ$Bzu94qE&Kb0Mx2A}Zw~Ux{Z+%R zB(L(Yzd2#e6~MFo%>`@ZIs5ybleveozjow#4QGFi%-?pf?5~QkIj7m*Pn^#=$aR~C z^Bo3D?>FQf3`_64%L0wB$#uW-p@Ji zM5On9@|@=!?=pA~ho!fsn1l3ArPh6w-et)H{|;__IQPrJ&MA(4k}VHAH<;`34zW6s z{e8m8oWJaE1!|jM+24vV=PLVqjau(R*~B?SkA-D_tHHeIWPhu1y?Gzc{#GZi+}|3ov&gGF>~BrjnH9jZ{jCIZ-Lk)* zIhkve{dFSG?~k&-ZP{1;3?TdaEB*C_Wq&VmKIb&|)mog-xliwR`I9Hs_1Y_&U8$ljppo_bFHt zEWPci?E*{hv*bAk=^YO{q2%3!+D*wz@4C#h-&>{kNotRUrMGi2f9ah-o^zGnaj@e{ z-oJ6|*wVWKUcaA7?|S4>{yVt!;a*!Gc0zINlWYUniNV}|8^W4+&i-EJWbT*j?@`!s zD>(G zIMRDB=Q{z3o~ z&iRfhdAH-(dsuqwV8_GKJDT&I49oG}4RcTCcpq=WyFKSSn#g(Hf#Wku-W_2&(t9)K z`?=)33HFchx|ZFjwU^viyD`&*)5**6_J+BZ>AfDuFR&c%EimW4GG4tssCA9f z+Xr?ek=|W7K9#)m?gn%3ruS+bKf=;G3g*0|cS~xIhoyH9*bzi}_vH9A^3q!mTc>m# zU5VqTlJ^>z^OD|8sr6b&@7}P(iS+Kn@p&b0Uzok*crV5A11!gT3CuN0@4nP}U8eWi zqCdSpT#o_dozHWQcR$#MrSV=!?KJYzdl}3%O7H&EdLK&fRphyD=^aGg6|nU7hk3pn z?|IaI4@>X4Ft3I59!Ra%WqNNQ&wY^I2J$Y4rFSsQ`AhHF)J}z^_bk|bu=E~Gt@oky z-b$YPKD|T9Ylfxw5SVL}-gBt6m-L3yXQ??W)J>GXQP9#Qi8+F>s_-s7k}50>60 zSTiiWZ^0U2Io?xX_LpAYllVPij@S25o0i6VG__vm={+9i8m0HmHoPanT(|U|h~t8i z_avB(9Pbg-egRAG;V`fLir4#)ukZGf*U{zVxi``~jJykBIo?xYzBf+qq11j!UV4v& zod-*=?{!_H^o}6Uy^&tO&-I>_-ZMD1m-HS??Psv`4uM?;OYhm#x?j?JRMB60hm+?$ zHNEF>?C(F*JD6I(S5EIxm}`_?e;(+5NpB;0u3LKjdF8;8*PpR=FL@83);UP;K$v5z zc)eeJ2s;Fp`}tUybD!ReaQOO@4gUXH`p>hxa9vExYZ-qcs+O`H+kM)(t9m=*5-I`fSH%xUevx1OYcswYhdZU znOf&Gy~D}#H8j1W$a7B9dkXA4SbBG)_5)aYcZNBC>AjWOW>|X9C+~DvdYj2}J<@wy zQJdbLg(JN^U?Zumc)fn7!487uem)iE`;zothQswpuirCzy`^^tYMuM^ZVz)#)B8gk z-ZNmnKS=NGI7Y$Jdk1VJEWO)Odo?V*-C(0(>HQ7X0L$@i1G@*7-Z4152Gctl<~Y*3 z6}9f6^md0`1xxR*ZFuWoW64YJy*RFgrT0FVj`VI$?OU+)ZU(y*R`Jft*LPBDFS)Na zCGTEXdLJauYd^g=!7hiTcN1#ghNX84m~)@rd#D`=OK&ZC_rcQpFnNCem)=J>wwLt& zo!U2G>D>@^Gc3K2Q|o?7@0LY>>3xDczZXsKH5_{{N$+~pzFG2a1apnjJC0iSOM1JK z=eni$N%H!XyjODUHJ#or)H(+_-gRM)t>X3ior9m79RSPyyft|vVd;Ge$1WxBc#iEk zy&F*L*wWh>b^|QE)opmUBhS5--e++5bA_Der#W^U>0O&z*D}3p!EP*h=fyh+mgDV2 zp4WMLpU2_PIMVwZ$BrYtYfyU=EWIz1*PHi4j#rQKJYO5$)u_FVKTq$=)cW(5^uEHe zV@vPK)V>T$??m$YgxBl$1(Z-J%vb@Fy1()$?4u4Q^xp!OA5dRKut2bJ~J z`z*E2ReF~v&-+VyUn9@j^u7f%FTKlB`ywp8?~u1^X}ptQ&S`r8LhVDa^u9@+bDG}! zVcv(*yEL^g!P2`N%=ydlzE7=tD7{OP*8)rLOXRs8>HV;%O>a%%Nbg7F`R|-^UoFPI z+K=<)eqM~)+e_XF)Vdz&y$f~~EWJxn>)fYz37B(Q8L#`fLmS=> z9Pdw@uUpCcGsnL7Nw2+~0L$^tisK)!9B(zuc}eezZFqm>eA|?~cXRADliuH8_LAP2 zaQGf6y$ixzqx7zUw;wFWJ7>{f?yIplMv|A_M%Y=f^j1;pd!+Qv40DaryLubm*~lA1 zUV0zEaWyQxGtsP$^#1nC@4x>Omfl}sEwGB$*M;9XpS|RL@)z>l8|iH!&-Xs*9RfQ8 zmfoMK9S2MA?=a^+y-u)el-{4nb8n<~_F`=5^>y(i^3wYwwNJp(I}PSI()$aw&S`q* z#IXfjZ{_-(6=u)rt%liOdZ*%e9G2eiVGqF4`y=OfkEgdC=i9vGod?!JUV7(-Iga#x zgX1w+dcTHw4@>X&oZo9Yz3n;QW+m_Zu!qP??|d-Fk=`$HJPJ$ig5<3f_euBF*PQ=O zp3}P!dH$Xy=lNZZeeF)~!sNBU()%gihhgbml)M#7-cLE7_xtpIRm@B7s}IQYbv3=S za6Z>Oy_2bZq~!e!=3M1Gf64j0m!$Vw@|?f)E{5YpSb7(S=}7N~)IJBx@h(Z;5~cAj z0<)L&zEAC=u=FlX-b=9b{srbeBE6HSod8SkM=<9;$J>G07Fc@UA3t3FGqCi&4SN)p zW}z3ZXk@pTPyFr6LKbOjVwNp`B@jAxQg(L5iy~yid8t>ceAK$N~_hV{Z^Yrda?JFg3 zZ`fO~^xi_P`yjnH!JN|^?|XQCzn0$j$aBrpy9Be4XB#i~iEPI|ZwjynDdjpfVW+`@&u;c@JpAdnM;{ z?sL3*Qu|8DdmDM~jr8^-uVcyEANG35I|Oe}SdRBH&gZf-iM|4G-@Zo(t9S%`AhFZ zc&~w__i6Hc&z0T*INTfQJq$JpR`GiMHnFd^hvk00mGgOhq<1K__LANw$@6tJy@SYe z-O}3#dl#18U>?>G$Pj16|GUsz%(t9kmAHdSvMBcZs^d3i^*Gzhk zFKR1Z$9PiV$bEG>%sEZ(jqIz-Vd*`Zys5DCo`~ZkSb9%_eF01Fanw5Z={*MaJ}kYr z;JpHt-XqBS9+uuya4bn~dWRLgr1vmt-3RGC8|F2b-m5v+r6upNMStl%4aX8C@9D75 zs7>!sYTXCvZG^c->AjwFT~_iQMV@P(-ZOD5UhQ0FTLlr;T^*Hocr{iL9P2Gy(f_OUdek7c^yjLb77yB zyd&^-gXMS|IG^*9-t(w+jnX@aJg>_f?*-&7TJl~9YvTP&j`!j=yn{KP^OD|+sQtX; z9YEfclJ{cr7AbizDQYWT$9O>D$bEGP%sEZ(nVjo_lD8juUY9xE%Wy1Q@?H+}^T*P= zAGOYXdi%qi)AXLpxh^bu_a^VllJ`;^*5>tlCCt3^)>G>mrT0jf*GGC!#(N$t*LP3y zet_k8ufXAYq<19j$kKRsr`9z}?_MzHFTJPXJs+0d-N>5;OYiwO+#Bh=7N(;zUa#M4 z*;iY^az7u>`QC%2x0zadNpBzWzJaB81bMDodas8a5#!Z+T^rtAIiGW%-m9r~jndnj zyst~%QRF$c^xgoQqvX92Zyk9#-d@Ex(t8uN&S`pklIIU%bG+A)XV2-q5oUim&o{K; z-I4QsOI~_!rPe)^-tEcz9hTlE@?4|z-U4$y(mNXOmarUe56xk~T0)HeF`=XmfnN#o&ig59eF>) z(t8ySuebEx3v*AVcL!=+^Ym^5v%mEA!+R<$y<3p?BP_kw;P84&?>}I(#dYL7??kP8 zJiS}OT#xh)#CsYny_=Eu6D+;=;jlKn55UaJ@ph&5cUXEifw^w!-5u{Ku=KWYzG}E$ zhSPdY+z)e3)B7-NmeP1P!eLQ*H-$NW>D`O-4}+!mQO-9<$vYM{wd8#a<~VY^8{qgA zmfja&Pr%aKlk=YhOYak$Z}yV+9+-P0$NLb>ain)W9A5kB-4NzJP46z8|72Kt$8o;d zO5VF*ULWaw5~ic#^>yKK?i*cUdB1uD?^Ce!K2EK@q<1`Q){^%Q*tg`R_bJ%#u=GCB zhWBrHo%{6ulUmm(y-&ktDS2;$d0nRW8JK%0y-(xa1f1htrVCh{MhjX9a6=BY4df%aT zU08Z2!Tda*9Phg@-?Qd;--FppdY8iCK1lE8Fz>PH{eW8kjc0mSD*Dsw!}&AJojIEKO^t(lJ|4ioF(t8ZFm>qe9nD(Kc?0-O7E)Vxi`}L1$lm+ zPmcFX*j%t2?;ChGDvfs`&gZHVI(LtyEhTGUp&uIK!PBlp$fFy}PAKTx|~$y-gH*JXNtB=69Y z_b1rWu=LJDt#hB=xnRy|dVi*NeOQinR`TYBrFR;6*5-JBg_)P$nW=S+(z`6o>m$9t zQ~Nhqdc9!3gX!g|_4@sVJl7+=Gm-ZfSbBf^dHNcqcNVsbw)9TH z>)uH3Y~;BQDqgSON$e|g_@f!I3SO^|^!~>A>?OUkQ{-pFrT07XT(|VjLEikZ^uF7M z_h;t7xlixR)VfCLt%mt|b?KcG<~5VvxnK*x()$r(2G1C)U;kLleR}8S4_u@4wuAZE zdFh=8wou7CFRVk!`*9oIsd$~2^v*|Z4J^IwVPC=0J3nkNEXTV5%wBT5U*Kg){$;%L zvc9fSdKaYDYcRb(6#eP-K}yv^{O5p@cVXDVCGVGQc)#R)u2Fgyp?0Z~_gnH@xAZPb z-XK_xw*$<-r;y|Q9xrc^{$;$MaX#lIy^B%1BrM1K5qYj#dY2$?e^`3&X8xYTp?5aI z>u72l-pQQLc}ee*oX@}Om)`fu^Y3}2cd24*>Ae^4vn4NomDsmt|1#e9igBd(FPzW6 zbC=#4*r$cJa*oAe_MG0OVfL5vOe=r#euLL_OYbt&E)GlYvanBJj+dv_-j{~ z$)hn|QmWo9<|^-#Z^Cq>cLkDHf~9vw*ko9GSAw|@(z`s&xvGqp!>V_RIZf|GnAd80 zS0=e5EWNA1K8B@tRak#mdOu*i50}RK0S@6>Kc&`dHNDTmyf>xy zuhcFNOYhU0KM+d z>D_?ZWywqLhOl>G>HRy*c}ed^MQwWjSvb=BFwFNf>8+)9&60N$*rbxTE6jb6-c4a% zgXw)7hjX9aM_|rrdN-rC6D-HOIqV%+j&}>#9q9XNbgwK`mpqFMeQ1} z^gdejm-}jK@^&XLz1zUN2GjctweExT-V1Y$(%X&N)nVyGa+Wb6(QB9kr{$(mR?w*Dbx6 zmfl;)^R*(qI~HS0?-zID`&yRbV;Z-muqTj+dua?=CRsCB5_V=HkxMcyGdCFX_DwrlaEZ z{ZI?hj#}S)`r!0cyYioj$RHR%V6G{(t8Kqc}wHHlswlX zz5Q|Q2up8YnBz$AJ=A({O7F$6T3CATY{NT(JYOr)I{=40r*}V?{iU~sT3-Xwdm(Hs zSbG17*YoH4?$7yNAWzOy>wf81j4i#p!(5N_-jBm;HN6)V<4ErlZFmQ9zUNEc17OZo zdI!Rsm-L>C!`IsMo&#GOmfop+f4**MybYXh0(m*!gJ2uO(t9Ax`AhFvIQ)4@de4M) zhNbt1HoQYQ-*e=p_Yl}du=E}Z^L*(&6~{WT^ghUSejH;<@2_~*DUG+0^F2##dJl*F z9hTmMVa{KAPsY10EWM|}{t8R)FKu{__9QI5r@*#^rFR(2UeY@!c5rMHpuxkl+dgIe!H={=G>*DbwglD7pcy{E!FUwS*>ogJ3mA)L=^ zA-!i)>vfslob;Yc z-e$1$o(FS{(t8lK_LAN~FdY@I@2_SiI@3$O9-Ts-*H(Hjz_BSTz30Q6tMnd7t#g{* z6B+OD7_Z)0+VGxEt@o1jUWlV>$$JsZy^-ESsddfMdjQP-a=g`e&B^heO`g|udPm^c z1eW977v_6`^bV!gYb(9`!Ca5@&e4YVJo3Dz(|a+F9(DggJld-HlrBFX`PA=CzRHUAhhL zVDj95>Ae<5H&~AMI+%|1?n3R6u=F0rc+ZdVQdTvR|JY03HyWt*8cgpf^0tAc_Xe12 znciO19tBHpADDAr8L!@J_>XIp-h;{Wo|E32$lD5*-kV|W-SqBA?IE!A?hJEY(t8E9 zu2FgqC(mm?y|x(t8MbUi<02oxHKI^xjeQ zlHP5pJ(#@oc7wS_>Ajm;ugmoARP>kLG34D-^4FA=Ve;;RrS}n-Yn0wiskN8%{vD>HGG1Q?S0GlS zmwY|iojk9t^gf2;PFQ*$g*jL0t)G@fB)yO0xTEBK0_NUG z?>5xB=IPxCW`8-})$p!b8t*>jc}=JHpEz!Z<#@-zwuhy6OKQE&)4M**^+@lkZFu)3 z&ucopPvRI2OYeBt4zTpDL#@|Vdb_}!m-L>A*W+B@&gA)ep5AA0+y+bUv#_0E>0OK3 zdRThbg*h+jJ);e8C-Qodm)@swSexGGVdkZG4Qjo2rnf)L`Acsz-UUnJy@Wj9H>Gz1 z4(BwzFT!?$rFS)IcZa2SO_=kS-YeViUQXWXD{yDYVRO5PP*UTM?Gt_sO&8^iGD^U(WNQ)b@j=cM+J^bb4Q<);*ryJIPxCmfp|EbB)sb zDa>)CcOh!~!qU4S%zIdRU!vA)I=y4aTOO9)FUj+MpWZKEjw8MEQ@al=z1zcl-;~}L zsP!I}-uuW~4wl|2E%i zOM2%Z&oxT#Pvl(*OYhGx?-A*pmD&badgp>U_v!tPTGuGOvykVWOz)56S)1P9VCJQF zCTa)6(pv>P6qeqvsC6yV`xJRg!_xaJdCqBiXW~W3Yc;*U{lw=D!_qr5=X3tj`#I-x zEz|oPc}v04I|~lyG`+LJyf>xyXS{=8>HQt%yrlONyk7h1eF5*1u=LK3!~1S}Z{%ll zFP83;KjPS*y!8G8b6(Q>F<$R=>3tdR60r1E_=XDXXAXXQF`CS z>ou6(UyJ_IJ2#GtOWtc3@5`m}en;)VlJ^IgYn0v(@%lQG-k->G-O@V`-Vv}I@4PS_ zmGS!h$C1S1^ulRp#Mhk9Yb(8TQEM;hokE^_BfayJcM&YTS1{gJW4wBgYQy^t=X37U z+m7?OU()*}dG3w$E>Mguy|?3izT`a=?-DUq$Nfbyj`YsK`JB`Ae#kypq~vW+o;{~` zL74sJJRj1A_j9~n)9Ia!TK9N*Kf${&EXTVrd9G1<7lOGS>D?di(xvf!ir0HsdS|BA zYdXE3<6Q`r-bKmtexKe&V2&fb1KaRU#(Nlf>HVDr@E(@lukbDiOYdSh&V{9SahP+E z-hOyJSMKKz@g788dVe6n`+a&pDdr`;OX4`E77KL>yh4N$U7UB-lbvAUwU_L z!}~VpJC?lkE=TRbu=Ku3UL!2ME0A{{u{YBDF15}Epy{l2{9He(aymP{Gylau?wUFLUFxN7@o8eutG~Q?N9zkAu|4OZEp5BG=R>RWU znY<=gde??Ij`VKYhIbdv=VuP2cO7bvEO{3tZ;q08UGk1DdH)8p=ZcrDQuS0Zw%k{H za6bQyc6zTTdGnI@FXYWu@^-=DwV&f%9p)OP_nBfX)B7aMv8DG$yj#F>o|nKoJ1pmU zJsc;((z|}qOL`y2;d-QZAI|4zXr%WcyxmIP@kM{>-4MqKCGScw??dT*lv?*+dY^!~ zM(OQ!9#KWA=IY{r{sde4b`y6@h z)AZJocN{Fe8^au1dN*prJCJ$)v^38fQ0p9|_i4OdZ|U8Hyi-cvt}y2yy$?|9{{tk) zdq2$2z{&A$PObYTy)EQ9FX`QaykRBpmM~wF(>t~e?}ISsD!p4#>wZb^Bjh=M>D`*V zlVCZ|+rWHHu6P~eJ=8h}xv%bmIkxm}OYIm~dLJXtHBWDM@=k`Ow;Rl!)B6u<9b0;D zhxyq%>Fq)7-Q=bBpX9mb>D`{ZpJD0U0d@u~y|+;7dZhPem_4U=Cu;uzOYa!+97lS4 zk@pKMy***heR^-8*8P{>QLqI|Fq<_Pq6gv3Uh1~uh;K2)cTp0 zxu36w&0q5FOs#W}-iOF@jncaZc|VrCd%{jFd7IntUJG*$a-Mgj*11aWedM_()4LaW z(_lH?yu06rc&~yvFX`<|t$QQAmy+joncn{7{QyhvzOeIP=^fFA_X?Qv zlHUEOJsp)6tJE6lY_??Kew2TSi+X3iJ(*hn|4zA|?}d4Nr1vmtozwKbL!Pha=^aR(`z5_c!p?=I_rx~5r@)+p z^bV!gxk~T*Vs=fKi?bQ|82V6I1cPo&m8nck1dI}n!MlgRrLmfn+L&Ov$)r}q8QcrS*zM(G_! zt=E2fk0H-_N$;uTeF01FX|M}m={=+k?-4NPD!pe=dpRtaRN$&*| zTw3xDB=4e<_d@a@bI^6WXi7sKo? z=eb`S-s54eTY4{}_7Ygm^S9*f2g`ZBoIK|xy?c}AXMR`4>owaSW-oc4+zX~7z0I8O zdRThDCvN~Oy;qUvK1lDCFz2c=UcG&(bxzZ}2h4j!datH-6fC{d$lDi|-fPJF2$tSm z$y=i|-u+?jgY@=+*>ifY=X}?}(z_3Njw8LJ$omkM-s@n_eU7(*TK7SEcZPYbruRl_ zuZ5*|Px4%k^xj0?2e9xy7HY47rFRGNT#xkLM&A3d^xh0} z9O-SO)_YTW&xiTGB)y}ly&9I@?a1?5yqbdMS91R=h)J_4tc#w z-g=n#*z|UWIga!`!})H2rFTj8RUcS-pCHek)BE)Q;O#BI?W(TtZ73A?;%>z$?#11u zSn;C4p}4!d2d7XRin}DZ!^Is!2_E3$4*xO9e$O|4qvt$_e6-x>p)2d%W6m+w+H3D~ z6H2aF<@0@BUmo=?ao^>{>b|+}SuKN|9^5S`K}$lPs_K?f%`rwR(s{^on8HEPu+JWv3f_i?~HAHr;+bk z()oIylgD#n?mN9|SBv$|BUb0YeWw$v%>QLLU-_no?p@62L#Kiv0Kc|0TL>wQ_Q*5SSjskWX~_nlI#-c9a1g;@P;Jiot_ zi`D&b-#6s(w3z$8CRXck-$hkh?0J|5vzT1-(^%= z&y4#{)bw-T3B~H?4dlrrOHlp9_o0cLg!`{a(7;#N77>vAUP;d!0NE5_8|-#A@$Z zZ~1G1=DtIVRX+E9U$u26-FG?pE+^)` zLk&0h&lBDEC)L(7<5)5G9Y(C~hx>k{+BzHVyRv+j6?5OA<@u}#;rF&4!eSarbd+)x-srG0w_x(=3bvE30 zIMvqszH-1lPD9xCR(U&y!q-gNgJMYZ)Fa^K-qTj$<=N0shgG4~x!toF)%ziRf> zeLoQ^4}9zUt6s5Z;(axeeCyeA-!bKJPxy{5R(s{XpQ~T(sr!B;w!N*ld_QdCJGy-9 zUE;oD%j53w9Y?Irr2BrWes#~?_fxUT=j(l6?CDtVALLulmivw;kGsTtpZ~4rs6Lb3 z_XqW>Gw!}0h}HdY-*?;ijwj!Gm$>f)^6e3G-wDO)dms0GOa1COci(|xwO8)@rr6`L z-uJ|6|J-*{c|0NJz7scl;J$CDww`nM9UxX`(tTfV<2y*K_Qidtl*i-YJC#`VbKh4~ zThE*OzARQh2f%$_5_>Gx`>|M^d-t7I9*>3Z;2-suZq>P;J(kd@%>D! z&Y%0vAdg4Ie4l3&D-ZX5Mz!CIx$o0ryV-i{b8u?aRu+G6d?{T$gYNq$=^hbt-&w@! zUb=6OYQGb6-)F^Y?^$p8PN&+sNACNLboKso-`S*lNX&g_7kf?2eIHZp7h>-Fs95E5 z--%RP&!zkRBwam&?mL%s4~x0)oMQFNxbMTN{Zh<*pA@V6;l7ipw$6t8elK0UbKG}N zeJ0NwVz9Pe@nmaNoVj0>s$clt zBDTM+w|sAIR$SKZxnkZ*86Xz9Uxe68Ak>zO#zC?~3w$Ud(;3Rqg9y z?t6_`<#XSw#lDR7UM2R7bnbhKd}k4J-<9S2q?r5uOSSbbao-!n>VCNIm2G^l5UcMa z+;{b6Ps_Ky>fOGESmkry%T!zMZ1;UztoF}+FBSVN);pQL&#CX5+;`PxZSK3a*t6k# zvHHClzL$&DI^6f7HojBI_sZ~HR~~gQ-FH2)=fr%y7pUJWV!qx>#A;vM_k6L>W4+VJ z_X_EJy&KBo1u^$MPqnM4wtV%i<TJ00W%8X? ztnM%Se3EpF$9hi`dsjO5Jx#3k-hH=d?veYRAYGje_ual(o1d%WrCTg~j}v=4e2*2Y zb-3@})UWo`edp7;x+r|Nk*+ek@888LpYQWA@_0+k*L$>BJ=5;HzG~|Xx$lDVy->{8 zyNh&nkKA`>v08`w9x0E3V(xo{SiQsC_phq0XWD%ik?%jn+;?~B>b>i}yNT60-1jhf zyea0skBim!Ywo+IYU>^5zDvmW0x|c!P3LM#vGP^P&(H7H`bF)%`yMKfH^kic5V6Ya zzALG=o@w_zO1fHy`}ULX6k_hXmuhS8-S=SC{#(p_9}=t2WcOW7we{KPzK2U!_rrbf zmG9(Y?z@j_>sfW*16BL6nEM{kWO3g;o6PRJzjSpr-1mO@P8PoVskZjiefLvsZ}{#j z_L-Rb?ycI|d-vT(x;h)~dysUse}1kG7OQ;jyO(M|67%)$DfWe!`|hFI+Ee#kRp;s) zG56hHx>}q29wt`Ls{8IPj}OJ%ckgB$?z^jM>kPT?8uC3`%zclLuI{D#9w}DuCimS% z9)rZ(cMq}J7x&#swe{?~@7nS`OU!+bmag8r?z@xvEfVj^zl(h!o%`-8R{P?)dxe`JO4}zR&4gO(<5rD*5^OJ+8^(>)k;f?@8yr2Z_}^a^I~~ThG4x?$+dU--yOy39=Y#!@~zKM_uW~#x^M3LntaC>^Yxyj+VWt%_4|))R9pWN ztUhD;XKS%~rrq~=)m9ew-AcMT8}2(mzT=6x@0QZdVe2j5K5cxr6|22>-xK9q=f!=u zkgm>#`<~vc&Ck{5(#;;en~BvA6k)ADe{UgH>u}#=)vxx{eYexOIyHPxm98?o?^$A% z&-ZzgHohB+)idqBN2#{XxclxT-&4eVz2{0-_sD(E5vz5$??&?YG1j}GSiQsC_i)wL zGwr^+$@gS2_dQ>_dhfdLKg4Ps?z?^)-}S`4lFoe(QEk18-FGkfo+ReJy*gK8iIuNP ze!hRF=ohv3?z^r$eu(w1BUYK+cVE@kGwr^Aldjg`zMso?OfmPpNVT>1?)%p^zH5tp zBAxplpxXMJao-K3tNY=;U(0t4G55VhwY3iST}!q1#(LKj`&2sjy9 zeZQ6O=wj}BnQCht?)#TEzN?A-B%S+SrP|NM+;sehx zwUx!+&nt_Shx^W}+K0s4cO~iSbJ~4}(tjg1s+jv;quSal_gz)BwWsd8qSy~&?z@0$ z>(5`@cLnLb6m#Ee8;|m>uUgxoV)ZPz@A9gxd*r^$HT~Ro4%I#ozNhFpIx2jxm#+58 zeQy-2eR1DqMSTb~O{%j2Fn z&r3B~-1mkiv-_SQ-y_5KYU%2}x$muFwWsd8WE6?5Oo_1+j(tbA4S^YeS1 zeo@bs`_3beItT7Mw^(KN^)95^N5tHB0qJTT?mM-7hY@q%hgDm9@4j=Xw$6e3&M8(u zx7mFcRqdnUJD+rQKiqdZ`Tj0^|E1b`w(zazcMjFwC6Bt-{4;x##eE-YGW&WLm#)r+ z`_3fap~Lra)z+T6?`&;+XBDfT>+8OcsJ8areHWIl&W8KWD&L{Pw@0KU8PAN8b_`a^%`WfNwJGpeDhVQ)ct)F}C>wQDDwO9Do+Wx58%Hq$_q+;da z>wQzTdx^R4#L|r-=IfnbzV&my-S^+Bt-W&JDO6i~>b?_+Egx51G54K7x{<}) zcR~5q&uDhvM^#(Tg8R;_+B$3QJHA-$i~GK(+WL2Q?mMn@^*(jqMdVvQ=i7bXRc$>B z?mN9|>#Vu&WMZ{1?)!mi_Y`yAv85YH%zfW)Jj%De>Rvu5R_`_U{ex=j8Fb&Vntty4 zv1;qz`Md8J($#;D?Y@JgtG#mHkHl(U+;?=<)-&k7qlwjdao^8WySteCjw0RIV($Bn zbaf`(_Y1M==e{GUww?v|9Z_sPG57siwY!M9@9@%1D(1cerK@}7zMqQK{czvmR9nx! z`;IJD=hS__Q|+!|?mMh>lZd(R0O{)aaNi%q$^+l}Jo=q#Zx!?R)zD(|h`H~Vs;w;U zJCt;ji@EPB($)LGeTN=iKc8L9eSiG=KYfQ$zuJ5E{Z{?z9=Y!is+~;CeP58So)7mO zR=#z1-S<1y-XiAf{Yk9$#eF|fzq&{6`>k|iiMj7H($)LGeeV&g&rA3Hrj76SVzn>s z`=M&tD{*_nq$hm1=KRt*`e>v08`wj;wz5th(>#(v26s z!!`HPeMb_jeD3>M8{bdG>b|+}=&G$})qQ)V8$Wz|n|*QLwYA>*z7yYC;}@!}EdCsQ zELI-wJC=O64c`x?8z+1}mv5Z|_Z?HUwO8)@iE3+4-FJ{!J(un~j%v3P^Yy+j-MC_Y zfB(>UxbMnZ@4>O&_f%W=$bH{!`nm4}>bG_HZmQ3o`g?ZVcU+P%e)tmCTK{`Lr|8BCl?|4mS z_uX2)?}zUg($#%)-$})4Pg!sIzR|{afY=1mx$lIkt$XRd+spU8@SQ@sdKTPwaFP7teP zy>j0d%CE*TkboyYU|l@-zVi!=fHh?#45A< z&ZgQ;#N79J>1rMR{C*?f`kZ#(Syfwm@4k<#w$6e3J|;GmnETGD+D*mW_ZjKxez@;2 z($#0E`;I17&lbLY^&CB-+H1u8JU=W}S=@J)CbRp#E?u1s_Z>>Q`kujk=N79yb>D~D z_&z8$MflFH+S+^feObCX8}2)wboD)w`_3y?>u}%uReQCVulGK&DaCxfe^zbnz5CuH zT|IB^yO4DAin;FsVzmzUy{nCHf3aDlbKgZ(JGGem-YMO5V(z=Rbn}V1?;K*aSNPT% z?^SJO@pE;%Sb4bblB!)_%zbZ_uAVpdU0S;N#oTuZvDz#5?Wfw>Q}?}DY-Ta{T~@Ul zh`H~L(oG}gzDqS8?mM$soqPAaLA7;{-1qvXpZl(;+I7SCc)d4XiTB2G($!wM?+RkI zFYbG-Jn9~~?=@mGin;G9s$Eaa*L#wDUk={|rK|hqzSD};xyQG!p5LqFaiy4_=YKU> z+;{mVv-_SV-_cVP!beg^& zD&P9=GAq-%s+}o(&y#QcebRoO*ORWE1@~Q9to91u>UT-Ar+%(35GxN~?*^(}BYZEA z?{ne1p>%%=-}S|6uiW>-W>4Ms9If~CSno!vT{C&de4mYZlT)M zrE}kZsJ8Bh`}UD;HZk{IM68}IeEaJ8Jx#TjiurjywaMb^-MGo@zL!f^=hS_-mTva& z-Ab(X)O}BG<9m|UTmOH7^5r~luG-pr_q|xP^}M<7w$jZRzT1h_XQTU`pxR5se7(nu zRc2rBj;gKm;=X;Qt97{V-=&*F%=dW*vD$n0J*JKC8T#I={(ld?-knuj=f!)A!MwFmBdq-rk~^Yxylg8KipxbJSNT_}8yRc-By``#ws$HI4a z)y@^ZhqduNO04$E_jym%E*QQiNmu*lzIV#^(eT|%wR4B>!K%GT%-4IkSnZ4Z?xWg2 zi}^kuC|#Ww_w6s=N5c0t?ccD{`92@e#`h4h+86iTU%m^3?|#zN_m1v+K(jVKS3k*j z)bQO`wHJ!{diN2lb-3@r>Q~>d`FbCau0G@4cVFo$v-=(-R{4Crd$sXBQLOHp`yQ&= z`hLxQ@0PCKL+*Q+bhR(;`<2!^vKYR##=TWrS^Pb@yI6U+?-A;^w3z$eAzgi+=)Olv zSLeWe4;QPwa^F2wTYKuhyNK00!hMfc?J{ESdy91SIpe-PI#>6{xq3_M9Zov;-CwnJ z#@%-(vC8be$I5ppG56g=x>|?(J|o}z!uJ@}*514CL8`5D;J!PG)wAlp$E$WpG56h7 zy1F0k`+|J$4d3HbTkCM&LseVPs{3v)R_`YFJyEqwh`H~!($)QN-u&ylX~hx?u>R(s{XeND>44s;zgJ`>rci`{KTr zs&;uX_gz!E`mAu@%bK3z7wO!0H`UfV)_qqOt97{Vm8xA<%zamrZhbNLy;!=+ z?7mlsRX+FKS+(^pao^3w>b|+}RjOT1%zampZUZs*y+*p)7x%qZY$q}I{flZVi~Fu9 zRv!5F)$c1-5Iak(-kbGRpWVxgZ7AlxH_D^Vf%{%3R(s{XtE#s4)P0u~t9OL^F4M+$ zX|au@bKha~vj%Tert;O}_;bEn<-2eAuBF;KwKCB)7Y^YyMQR_k!z;pBTm z_}(Vp+IwH`I;yR6;J%BC)wAlpi?#7xPOR>S`;H{v>%;esX7Am1J=NB;>b{GL)w{`k z7ZE!{%-6fHSltiz9aX;9g>OIk);ipG1J%~M$$dB1jd{V|lYMor7HZ?Wpjf@L-S@7> z!_U=CVzn>syQymHnQ`AgH~oCQ3y7UAov(L(u~ozO9(mMWx$nJVwJ+}bH`Ug=-F-LG zdY_B+&ez6wUa{4}_x{GieIF95bML-$tF}Hv-FL30pRadLvA)vzdgl;ZJA5CMNA11) zJ|b58;=Z%1wmw7M_j&z(|MFPxY;AmJ6Ao|mw*E}Z z&-2rbhxIcZQ~)ulIS?*1dG!Y4ogKp-k@kjC6J1-1j-L+81B%dGe?~ zgYxx0qV--J>wQ7BbuZm_M)_VIzCF^_zPRs;VwKr_r)$>czCH52)O`Euy)lh+r>NG? z^DOecEPNl6N1a{weOc_U;XA$h)gHL-v|_b3_PKngZsR+Ld@l{(SLIR9s{6hswodp? zp?-DG-S>}TwJ+{Fx!B2KzTUaydrA1dA&>e@cHe)CZ6W6SJc;_PE9Smaiq*cj@5F6< z=aKKl;ro_6>U$seeOqkv@V!p`{usU!h}9n8ThGy(s;#~9b2YzwFACpxrK@*@`@SPq zYxDJ9Estr$e7!S@)jfCL*Hv44;JypX_rmb~mvr@gll#6WwyBu=-X@Q!#N2mWvC8be zAE>tOoBK{AU9H2P->apo@0;BBeX-hm_q|EAQ;WIpEMoOubKgO#t#j(W<4aff!+md% zuDXU%>8AXfY0zQeci9Zsy?2k!fwJa!dx-|xk?6La4u zR6B{7`wk;kS@5mr_p>ImpR31Ls6?5P5#OkcM@7JoW zy?5WCq^tLk`~KK?xbM&-4F3Nzx$m>8ohW>N{AzGN_Z>#HbuZm_9lbZsiucA)($!wM z@9)HFU)=X8`A!_Z6N=UQz@8|md{i+I+nO-JdhNAmqc z_zse9Wp>|j)UWQF`}Rs#>+tu+xYF$|=DuT#)!zGhKalU(Vt$^76|47}`;MpDI;ZaY zQFA}scYNt~l+Jy}(t7WV^?sn*%HqE7h?NJvef9jlEp~*MpXblS>iKZriRAltG54K7 ztoF*+`<`lRPu+Jkt+$?i_Z`^A_baiw=k7a+JoXoJ-#==-_r!Ytt=c*V?)!#VWp>{I zVuy?QdfyeRb-3>o^4(9&eJ2yEz4!ILrrJ6O?)$1(eb3;&ue9-fORVmP`~FcL`--{m zG-B1yeP2{R0!} zeWz4yoelSWU%Gns-S-^5vk%w)l&_Zg^Pk5qreD;5Yw5n@$#?kheN6o-v-{2|-+CXo zZ;y1f4);A@zK4bH9ICCoci(a3JA#<6_Yw80v*EsTtG3<;?)$iObwAwqBKaN~zH_Oz z*5SUB$oKbR?)#AX)wAHf^QyL3Ez2CTkCM&3FSMSnES4+bGnq6 z`_8A@dJnnpGxDu-@4i>b_u%lIU$wPY?)#u>D~tQyCsrQ#_SNUnyB|-`nJSK=>}L=jicxjviHQopJZQTdeMf`}S|+`=(gEOWb!+ zt*t&U{k^f6*ui4H-aAxV_sxBmkWR6GQtY3<`sjZb7OO0L6 z1vTGa*D4?HlBJ}p&o|GvwAevnp6_PmJ5$W_Eh}B-TTradBwu~>i$9CiedDY0EiP8~ z&GXHx`F^_Q`IeLJP;t+ZSK2@Y9A7F-#IkjCwy0x?nrU>T}|wEG55Vv9(L3srlTnES3LU7o+`#Oj^HS0DXiMzMPT@m1$weX-gr&v&al>KXKW z8%TGAnCDwvtk&lFE@<|^^KB?y=9^lq-c5Y<(Jv+!tM>+9m2V@l$HY8eKY7%-@_c`j zu6{r5`8F1-GvxWsR&DLG=i5}e%r~)EeNW9-AN^uNvHJd=ugbTX*fV0D?|gaGo_fB` zr8`Q@^KBt^lbGl0tJ*q$p6?v7+Iv4&TdMX+G55Vny1F0kJE`WoDyREyrP>?C-S-sL z9xLX)CyUjdy6?8CeL>89&ycS6*?lM0e3$UuMzysE?t7wYj}dd-aH?HP8OYi%#qiSmpJm0aZJz322?Id018%}Hs>GMd!m^89x7JraNjuNxvUoG=FvotNXEuYU^xxzKi5h?{Uw!k92hoJm0=zwFjQ> z0M+&t^L+ca$=9na8%f7kAN}HEu?>w?zP(gid+Pc2lSh39c)kOrt8?J_4idXc%=7K7 z+S9~5-+^Mar+(k|P;KqK`|c%Oedf7u{e2hPDX04$s@i`^=e~QW_Ea(V-CeBq)O~kV zZJj^&-B-HWXZNkY$7kE{JzTZ52kyI@YEKb!-y@~V^D;ntU;oaAuRi+4K(Td<)&2Op zYU`YOzJugZ`{()U?=9L!df%6$R9k!C`HmJFBv$$Ad-w7T7U79`h`*xUgwfFA3y>xZf z-1h|OY7gA^M6t`n-1j80UNQIGrdfyk9xPp*4fowzy7R=`_gLxb9JuevVwKN*PZj%6 z%zaN2%lq+!_I@?#`0AryJSw)TvAQ4oskY9B=j)?8nGT#HrvWj$k_0ccx6Iw>0cek>vBpqLU^o#ytD;leO|4?oHe#`T{qO-e&u2nvM|8~A~wGPjBfmr=s$n$O3 ztgZ6ZS9$zXy3BW{vaDcP`bc-X*z(4zb{W;y_sX8H{y$oq>zd~~K)N~yp6^1jdY^i} zO`Egf`7V+!^WCg0%UPB_^1EJaS!0!NUDeiSndjR;9<`@_e-Dwao+Hn9vDiPwJm1=? ztuyKQ))lM0_vdnL)z)Xc`>rKj-4FM@RJz&&_q|N)0x|bpL$!6*+;@GkT8I0tpxXMq zmHXD;r?^Snk1M3Bb-3@9V&{vw@4v*p7As%g&l}6Po)7oET)JAD`_`YuZXCX+NLTM1 z_q|H&A7bu%wb&lt+4OQfs)bKljQz4Cp&M!I^(y6?4Omx#IVbzRuJ+G;Z;-Cm=Ds(IT`cClH;H{C=DvSv*5ST4OIK&ZeecjdZy5W0 zi*%LGeQy=JD14V`GM8_C)iblESlvtay-l@uhVQNN-9Y+UFYlq-rK|6`-1iQ#3&VE> zdDOnR?_b60Ub^p{s;xhJ@bxY&UG0JIb3f^>7IWYJV!StWX{bK>=Rmz9=BbIkYOTc% z5u4Xo@rC74XT$R?E>?T$`RULFsCrJ>Naj)gE}hd&M%}e#$bpboE8Py~XA- zR(t_@)EV-8ONrIKc)t7P(NE0tU8(!Fu6p==dqBEcujhNPO};&pWlrhpi+sC^&0(zg zL#nM8M&;wZv7GYN^WpjKldks4^VR<+avjxqzK5l|R65W1NSl1yE6Z$_rH}gTBsROT z%J-;he-iV2D=A-{NzZqebamf6-+ay9`~7`Px+|pfe2hSx%Yf?OIP>O^F1Y9eQ)dedcvY zzSY&&_vPg_`BqhyX;oWaGrShCE;W zJ z-(+I7Ue7mB9@mI@zD1?0caGNC&ty)RweU(ffBSgpil`WKZwIOdcLr9ogw%AQo73QzF&z|KKC6*wRO+k_jIu%#eBVAtM(N!U+?hJ z-7e<7pG#Nw&3(TW%kwg#_I^z1>Wlj^s@Pb@ijOaky62v6M6p`0=lf0`FN=A;VWqo8 z%=3LOUG0JA`=L#~;gw|!>FSGozZVPo0CzE6<K*I(evogyUp(KB@~yq{eBVn~ z-&1?OLDJQE@qEK5U+saPtHtD7pFQq-so0@n?)!Vy)?T^qxALvugSqcF@~wOBzCSgY z-S@j@ZSMOo`5qF!!^`)9@cmk~cZBZQLMjo?)znv*?q5(@4?|avT7d;-zD^p zSVH%z*30iNhm)@Ek^7D!R`Sqen7y0fMtG~Y(R(vSk-?|^3?>*(KXUp@AFOPe~Jm1IC)t*+q`l@~tNcWnU=bNxi zzFU-~eili6k?%IK`Z*!6;_u3%?uX}lOYD9z&o{9=eiZY3Z%S8t>iH&V*5UakZIka( zWvQRrLY6-I#g$_9b01)pZ&0%bp6?B@`y$`uYlsrOY*3*;l3}5Jt5}4)2a45 zG5394y4n}NAJa?sikSP3Csy~|eZNs{ea^V=b7DPW?mLZYYyaH$S?Ox8+;>LlY8~$T z2eCTO?)!ym>%Hr~&xk!F=Dt&^w$|ppV``u0i+!G1y2|XnGl^9`_q|%r<-ua^`;>ei z6?5NNR9pAlee35*%^SY6N>}&DeWwzu%G<4I+n@gobU8hPPsoR_dM;NG>oHc(<=m?MGV%>4&D^@?&)YoG zJr?=?Bv$v=^WCQVae$cTo3~BAWt64}Ce2)}!-&@qL&YJr!quSa(_q|=Z+AH^6R=VfK z+;=&#`b=@(n^jx;;=VVC)fsZ%6;%7NnEUpVuJ*uvSCsBKG51|wtTMarwW_VPx$nx- z<@;q^owa^;*6RF?FV^2!oxfF7`;q#2zI)_R=f(4_D&4bUo^NTf+5^vbrtZfPVxDic zHu?UbEcM@h@m2Z85UYQ;##iNANwu}7p6>zWt8?Y~>hCF@N&k4hzo@p(kmoyI`Hqf! ztGCJbMf2Z0lBM!}Ep~n6TSK*V-#p({%{}sbXK$hW3y>#TXc%T-(Z?EAM?n|!^> za;;^le1pWUF;@GxwrW3!eAmgN_SEzBNLx94|Ng4l0g>+_)z+SRzRSdF5B%J(tJ?R) zd|z&ouFkXj*5BPZefX}Y+83qseK|?@t^N$keXmr%+AH^6N42%5?t7bbbxz&)L3vIW zzU!;D{w&RXN0#r6V($A-^{YK_-wjk-d*#0WkgncW?t7m+rw!kYRLk=#`On`UETemQ zU~@0)KGywRT&#ZHM12*je@0cl`g1PNcY%Cs|2*H{J?RzKebR{ZzMcXQ-BtJxRNx0yWZ_gkLt1L^AfXV15xbhS3mx2agI z!|&S&^1VsS_vKdgt8?JKo2$0o?e6=$boKe>zFSCF_t$;Eh``##4=iYs{ zRP8`9_kCHq`aE~vt)#2_=Du5t<$39&y|160$==s_*+i^h|z4Cm0#cB^c z-!}4iQ_SS@{&$n%xe0wWP{X4k&;=JrBcAT-|=gOm=LC<%hSna*% z+g=|37V~_kN>}&8^IfQ(Q|KShx1DNhPyN1~ui83S?t8IVy-VD82i3kI=Dw#(SNrF_ z=gD*O@ZGV=?7qjVw$|pp$BNav(|vbR?YrT7v~+cb-1i)LRyJSnE~?G*vY+<;1l#*M zF9(U8XsqtX3G%4D@_fgM)gE}hUFGpk z>H4Z{_0P%jsC(}D4il?=@qByAwKgH&5N%nQ@qCB2$=6p|>Sqts7x~T<>uap|k@BebhUeQutoGjXO{;wM@Ap05 zKGN0w@O%fT=LGu4^Bvagso%F_R9olDeRmgoFMMZ^@8#ipfONHg?z^`<#}D7bin%E_KdMQe;2B@ z_R90EBabh|Jm0C()xGq5=Ze+ZJm0#_9(cZs+vM9oS)R2lm2YFQ`fmyNs(hR4zMbCW zv+@ajsbTa$aWEuhsJZtoIs$dY8EGwI$O3dD-vBit?yEaNkwLYHjYj ztbET8-|90$_#P;adUv|-^s?v8+CP7O_mxL| zR=DqtV)gqE_kBUN^=AO?yP;T}4fnl4wY7hKf0vc6_R4*4lCIX_zBh~gAm+ZWsP;ZF z_gzh_?z#J3r`lSZ`%a=gu6yn0`4;JFZSH$(`8I#AhA#d5eb|M1X6nBUtgm`r>N&bt ztk+oa0m^c(iag&k&3ZlGZStryi?GWke_Zpt>&NlhZRF)4cOXWLTY>=_ax41m& ze)#^aAXaC{^UbOAe2G})wcGKl`=8WC?rmoW*-f8T%PniQf{qlNp z^yu<>a`fmJQwK-S^_(Z?dZ+s2=;_mKL{Fd8Cr>?d>d?a(d*br}eRA~lLC2UJJ$=wI zcH;D*o;f2t`@qrBAC8{uIZw`Y`Y=bAoa=CWQx8W!_*47;`}{_q96kM+qfgFt`Y=bI zoLux_PF-@Y(}y|wzMa>4T0jIr{8_J~{gAgFZQW`mm4alcT3k>N|1z(3iU8 zT&E9n^vStSA9Re#(bESVV{-Id&v|mLmw&zuoGec}~SU z2aY~DdUTBFnKSnO%KrFge|YwXqi4Fmb(WB!!`sC=*F`{SAh@R{8 zW6lUi&zv!JaP(ZqmpLOGJ@s&O$$R7d&vp7_9eQ+(=$SM2Mt}04r#~D$^>Fmb(dTvM zjMSmW7mhADdah@l`sC;tdz$^CZ*t~bPknOqjF~4nbNZw{IeJFchb}pDuBSdZ zdd8mr(f*__b4GHKgE?dB;OM!Y`IDokKJ}eAeYnn;Jm=^-jp!N4gFZQPuBSdZ`dl}< zfNQGEZ{!^htek^z>m) zp5$Do4|DX%xlSK+jLFf{2OVSb)H6>#ywixj)AjU4pFHc)F>;+bBYn8eeM^oW9V2?? zjOgjZ9DQ=GXFdAl=zBYSpO_pyzVxGia`g0Jjy^fp>BAg-a;|3|^vTg@AN0x5(+3@6 za`f~;$Cy0z%u^3%WZ#*iPtNu9MV}l!eb6x`Pd)S0!#j=WJ6%s-^vSax9V6G7Q_skC z#^gCi*J(u0h%fr&%(xL|&UMBiy^fxdJ~?uok?Wo2 z@J=K3jJf{g%;}f<QY|4#Gd=;^~8UGlt+uG5H~KB-HNK1b@3GtcYH z8PVm)bw;jtn!~$|^kHP3%o(!}96i_Jo#yCMhdxL2$(hq9KTndQXH0!^=Ja7C4|7KJ zjBxbHxz0$Q)WOkl9nPE)jy^|p$(iSM=8WhV;rKFVgrjFa;OLU0=X%bQbDciS(RJcn zCl9*h%(Fm$fzv-ZdipS;XU>S8K5+EO(W7HT&zuo`_CcQT^V&oH>1Zi&T;pmg2r#~Zl=8Qef^Fbdt`sC=*aUFef^ynDTGiQX;hp}hj!B?0w zqNg5?J~?`HTt}ZAJvv78%o*YIVI0sQKYHp@pBz0pu4jKZdghGinKPyypVY(Ab1vcN zlcUFv5j}H8^z?zFPoDRSJoKR+j-Ebn^vTiFhY>wC+mU=dm?!`k?0=wED&K3rEkK zG4sID)4$~({o&}Thi88{_ly42(I+{2bd2bkGoq&tzUY&qr#|(`(Py9kRsWt2_Z?0i za>B`zJoUeHz33R}&zv#)|I$3vr#|^_-5)>J%YMQc$;ljja<2DsUm4^1Vjs}s3#Wf_ z^z>(rJ~`Lv!yJ7lP9O5{yeH2-si!|2J$=ZN{o&~814o}c_4uY9j-LAnN1q&h`l3&c zo<5mBIePkVJ^RDaGjEN}d27w%9MV76n>_ocp8jz3y^SB&$)6lOzNt@+o<8_ye>i&j zz|nW&^e4}|@%+NG4;(%H;ppj49`wo4(1q(1p? z-Jj=?^|CHFdB_PzpB#PONA$_jPE`oPgAM^At9 zpihpTKB-Th`d^xd`Qx8X<9DPgI6RAt#(X$y3ku)WgxY`o?}HN8hRseRA~lVI&Xp))@DtbsswL?4Nr2!qJnn zHAa6pdiuaKe{$yd=JNv2K5+E-!qL;85k2$P*sO~_aP-N4>;CwYhkS7IFlR*1oH6zE zOFbNYPrt#R>~QqS(c{O6o;f3W`oPgAN8i#lxzHy^PoLhy zANs@5(+7?|IePk#2Yqt%^htek^x5Zs)xYHv&ncWd=r z)F;pWy*!uA{;*y+diucGpXBK2lj}{6o<6DX#OXsGp5Nr1a-yd{ebCc~Jm{07r%&pW zqt8D7tNxsa><=doIpOG&qbGmrlcUc*zqCL7&@*Su{&4j8!qJnP5j}H8^wh)ACr6Kt z>*$lCN5_bsIV1Y)gFZQW>KV~9XGBjQ_78n>^wg(5dFr`Nf98zn=?h1n96dU&qfd?= z9V2??jJ&Vt!$=?IjOeL{qfd?=9oNw(M~{vXJ#$9%+(-JMPmUfPBYNhH=;?zm`sC=T zPknOq^kK~YaP;(nqfd?=9oNw(M~{vXJ#)s?(=YXK^y~{9eRA~pF`{SAh@L)h^vQeT z`vmgPhk7`A`oPgAM^7I{^vqjhlZQTV^vUVp%F%rPgFZR>R^QMkM^B&D^(JR>^ynDL z!@Mh-Jku-{_xBXM^A2e>bahJIQmw7?0<6f)N?)i!_hNmM9(}&`nSgBzO?2| z&W?4Bp6k?goJW6joyJbNl zlcR6xnq26Uqo+@6?6N=TTKz%~r#~Z{Jjt2UzcogGIC@5S_J^mA-|wJD2S?8cN1vQI zeR93Y(KDvL6Q>V(c;1tD%88!-^g&M_@}N(So<6BhjvgJp^iPf+-_$2ZPoLJC?I%=MA3gArF1f!O_znjy^eh`Y@tr&WN5qaP-O1qw9@*NRFO9sZWlcKDplH z=;@RCPMkjE$@2rxeSo863bRCp|oKvLA5fsfRP7@8*l1{LJZ}oa=COjQBD~56?WwQ%@hp^hFQP{>f9{pYJQ0 zkvwqpaB?svPd$CmlLw9-j*c;T>N#iV$pc3Z$2axx)YAt&dEn^b^hrIuo1V|})Whk^ zn4JD_`k>E!Ko3X9m^}6D1A2Vn=;4_sdFtsyfAYZ5!?S<#)YF%dJgs^2T)^}BN=|<` zeHh8Z96dbyC-0`mmp%v=8b2Ja-iH8z|+=)2{iKJ{=$`k-sw&&C(dNPhC99-eyU z=+l?$aQdVk-c6r5(Zg9k{gbC2U-X$1Jv{p-PksO9eAAyiaP)BIsfRP7rw@Aaz|q5* zrykCT9vvgT%+bT~<$Cheqhm~8^ziJTyqmtKIsf!Q4`)uEBlG7hfAI9`W56}L|JL%#0(q};AZ|OSl)YFH$%!3||jxl-a*%$QWhogsQ zp5&>g5B>fzn=)TJIyJ#$7lBb+|WGbegD^VGvz zV{;#wqi2LO!t?uq>xe(2$uCwc1W(;CAUJvp;~@?WcO)&Ii%$vk`> z!RZHQ&U%w`9geOw{+IWs)wfx1`eq*1MICy0Cq3sN?-zaG)TbVvdVJAm9`tbfq#mAn z=IF@@Z;j1;$vo)csb@}qa>5zm^hrHD^#hyxik>`h^zcsqKO%VQ=|g{V!qLO&lX`gS zxi9F+2}cjdH}&w;(+52{;ppM?NjRa=8|G-;gvp?CN{P1pi?px~Nxj)J252sJ=4;(%Fko}XVzBP~M7|z%!Px90= z&-`#kaC(lP`MmWdG!;=l;^4JaF{z%#%Fz^r1g_;OODmKl!iK59$77o&nA06r4Qp z?4LaK?8Edj^Pq>*C-v~uvk&OW2}cjl zJjqi}ANrFMjvk)$y3ihpeGL;J$y*_hd%Tt51jt+ z?4SJC>U(4V=@0MLzqiBw!|B6%lhYs0KAfzn=1DbtKJv{d( zIsM`EAy4iDdU*Cvo_h9y{`kVt!!u9v)YFIlCnvl$b~#t%Og(e@lM~Jer%&qPsb`*gIOBh-KhG(B8R6t4Px5a1%!wY( zeWZW#)DLN%>`#6;{geOFJcF9gi`2vEpX)-OzUblfNj1 zq#jQHPU{`eJg4}ghtr4t$x}}s^tmqd@a&&F^_(mElLw9-o_Ug|o<8&^4;(!_`zPx9IWxwR-k} z_0k{St$#ju=;8FCfAZ8bM^7GjYyAJYKh)EoJa9%h`IXLW6+QXk=;6#$4{wdl zdg+6noN)AT=BbA>qDRL_&eptH9~|G-^~gzo^ynDV7d<@tC-0^o&>;_;dh#T1U5`HK zSr<9c!zNZhyrqlv!c*T`XQL+%96dbqBu_nk@Fyo6Jv{p-Pd$BptNuKX_>u>n`IC3k zcUy03oxkw@4r!k3PkwmjhquN!=jd72AaXM{dU)zveH(qZJoKRsJ)Hi@(ZAaz5A~^s zGiD$3>5CrDJoWI_7&+1NJmz}g)TbWKh@L*^b6x1+t$B05lDDom`CEM(J-%@C@XV7u z_4H|t;ftP}**|$VeNWS$KIq}h$&;MxaQdL{wqEK}4`-wgI!1h%qlagn|?oaa6(}(=s)|>s&!|4xaPJi;i8R3~HdFq+dKYh6l&;H3f>EW4&b;Fsb9?poK zJm@nIdN}jc!x_<|W5ky^dU)nZo_hK)rZ0MU_D`Prfz5MEfAYZ5!_hG&Pd$CmlLw9- zj*c;T>N%I_$pc3Z$2axx)YAt&dEn^b^hrIulb(Lva?&^TaK`L|zFSV}Qx9jP4?0GC znWKlZ4_r^4dUTBGiyof+lmA*h^_d6Gn7-%-G~Zv56Fs~&Z=TEKt?NynR^LXCFC0A_ zU#=%lJ-XHyzUaxB{gZdn!!u7lFX$N2GfzEvxSn~?!_hG&Pd)pBo;+~$@XV7u_4J`X zdEn^b**|&e`R@YgPaZgWIKHWer=C9O$pc3Zr%&qPsqbmdJ$iD&(ZlggJv{aFK~GLN zdN_Sj5AURBecf`RV?@s!UpT&8@0PQrYxWgAoc?h3f$R9f(Ze%O^3>CZ{^^Szp8b<| z)AM;mAN256j%I(7r=C94bz3hwM)b_-56}G}4;(!_^CV9_eeh3T^ziJTyptZDd3t~0 zK9eu?fx#HiT4lr;ppL+Cwc1WlX>9i$)Ej`r@nuO{BZQ-#FsqDQ%@iC z$DjVp8PPMxH}&w;Q-_{BaP)Bcq#oW$KOoLIdivnk)96wUPd$Cmcgsni z)WhjdA9A5*&WN6Q=1HDoNyREmUIZwG>IQ_E^`fhpXlX^J)=>tcPFC0BQKi`w5o<8(X zU-a4Tn}aP)9=jLB1<_lrDm^l*IXpFH*SK~Eky zdU*Cv-boM7JiX2P2~LjG!&6Tl^qB`eoIa_Cr@n`LjLq`}M^FCDlRWkG$vklMg4|?*z(ZjQU^3>;hfIM*YaD2&=JoWTJPaZgWc=k{JANpp$ z>EFuJaP)9|Qx8vl|LEV!6aC@zp?`Av!+WAXIpOHxnJ0Pb>63Zj=zE(S**|$F{eXoB zuRQbM3uk{)4^KUL&}Sa>aQdVkp8B3<-RQ{)M-R_D$x}}s`jZom9-jS^che8(a3A69 z19_6uAKt&i{=m_*57|F?>e+|P14rK}Px92$hyLV&qldGu)WbXJ;hBf$19i$(j9=r~ZZJeA1shaP)9|Qx8u)ebAE!jvh{* z)WcIxF7(V9(KE+4_3+eFhn}2p^l63c+uhp{;*&p7me{ZudnFmfE`X{G9ytPjY4}Np84{-Ew zeCeM&_4GlXzUblEKY1rTJoDu9(i7(jJv{Z~p?~H<52sJ+;i+eT(31y_9-euUr=C9a zCnp>|Jo_h4{lE_U2}e&(e94nM_4Gkc9yoe<_D`OA{ORACH|s(V$2axx)KiC^JaF`I z`lKG-P0#vM4`-x5b4L0zrw@8^!_mX3OFcaG)S*WQM-Rt0_3+fw2YqY)U%J18;{1~n z-YpM$a;2U=nIDdx`qsGc|9zjsQ_s1hKY8Hj;rON=o_hMAZ{_H)|MW>coc`G-_4uZq zKBu=sPS%?|^{sWsdg17q<4gbK zspmR+^1#uV{;#w(+54d;q-@7mwI^WsY8zrjvkJ0>fx!U4|?W|=$X?e z_3+dW>Nog(j-H%w^l*Gq4^KUP(4&K+hcl)go_fv~dUC?i!|_c$JoWTJPYyVGIDJwN zPd$B74^KUPQV&l(eb6&!M9-W)sfTyd=RTl^Gt!?qW9sompE=RP>63bR>Zu$26bv;K zIyimcAU5jA9>Kj=?@>$JoKSIb4L0zXGD)L96fws^Ldeac;9a3e8~eRC;5}7o^_*V&WN5lzNv?& zo;vj8fuo1hC-v~u=X&W2M-S)UFVH`E>gj`?oN)B;?4P`wzFVH2cz)5t=?@>$JoL%+ z!s*}JIr$c@?dU9r- zXTZXPo#Et2Jv{ZS7k%bI52sJ+;i>Ow){UN=aP;uZlRWkGp+7m{=;7HvdFpw-=uaLv zdN{tRho_!C=*a^|52sJ+;i(_cVgKRi$&WAnlc%0O=*bC356}L|yXommAM|kcfjr4m zPapK%)|>s&!|4xi`8Urg^VS$V^Pq#%AI_Nh;poXjpVY%s&$`gl7mgm@Kh9h7)YFIl zrw{#EKO8+g`zKF5edy1e5k2$lpS+vCTOO{Hhk0vko_FR$ znuj|2lOLY>`#0|y^yuK|;piBXr=C9O$pc3ZN5_~v_5b0|dy_Ht6JPpM52rtz5j}a} z=;8E9J-m}1PG9tVUcl*_dU)#TgML7>FL3m5`lKG7dUBy}%^M%|aC}n_Pd#<$$q7df zr%&qPspt7ZPhU8Ccx#Mv1y4PF=ub{KdN_Sj4^MroF7m+9Q;#qGlc%0O=$W_1#+P~a zPfq{j^yi$>2R)p;sfVYYKIqZG(ZgF~>?b_+tb_jKgrkS!n|gTa>4UzNqu=0;=;8E9 zJv{Z)p(iICJ-mP8n|gTa>4UyC?~tE9sfW{_K5+Dz6FrbvYSI>yWi&;HC?efpw@(s&!|4xaAL!2rM-R_D$x}}s zbm@y8p8b<|(!=qg&w%Fp6L{`J^3>CZ{+S0oJo_h4eQ$IA=}#UwdU)nZo_hMwpFD8% z@a&&F_4N6z`tv;EOCEUUPu@-6ZM{Q!e}^oBQ6HH}AdFvk#dEjz0Gx^_d4febAE!PJcLk zQV&o4pk`g@@r9#@gj`?oN)AT`lKG7`rZ!t;poYaFa49Jo<8Ww2}cjl{>i)P zc^>J59?m|HCwc1WgTC8(vp;$`{o(8b{mBDI56?WwQ%@iIr!RVV_D|kP56?Wk&3gcz z`;a{KBO^TW~iG&!<=^3>Dkx9ZRHh%b5InLl|q zeYf=v>HQtjJlUW8@XQa-{fDDxT?6C1B~Lwl=-(|*_D2t=e{%YF%R^o2;q=cw=+Sk` z!}Zj|>Cc=#=y@LD=;73*9-eyopwD%ohtntZ@YLu1B_|v`Jo6+^J$*6{96kB7fAZ9` zF8Y%PjvmgxM@T(9_4Gkc9yoe9eNqqarq7({;rP%$dFt^+PdywxoH6zA)YE4`e9pq@ zlQ~mQ|J2h5JvumgIAiMJspmeTCnp>|9N*N#Q%@iC&G~87cleybm;UJC^oM8v);#*d=|lhI^oKLT)3?*}l6rhopZ%FLqGwJY z`X}$ECqI4A!x`z%oH6zELEkM8eNqpnKYifnnKPni-aqbp^3-#k{^^Szp8b<|)8j`U z^l;uomx-Spk^aGgBN8R^e_Nb^uffAYiW5C4z; z&G!-TmcHTiY4s1DditQJ9*!Q)n0k2XSqJ)_CJ!7v9N*N#Q%@iCiyq_4Gm4Ef0NC52t^pKGfmc+dL2O?2|n8^y%b_o;=w< zdFt~%k_V0+o_Ug|o<5l$j=rbar|h3R_5bnvi>418Jvs3uPx92$2R(V<=;7HvdFluB z8+^y;PaZgWIKHWer=C9O$pc3Zr%&qPsizP5T4RTG(Wj-0`-N}odeetFIq^jg&;H4~ z>G7ctdN}u&IU_vv^l9ljte2ehNj;qY=-G$V!&5(`dFVrb-Y;5AM|kc zfjr4mPapK%@??MXaQegP4<`@2x4GZ&?2|n8^r3(HqK9YyNBWZojvk(QlBb?N^iN;(@a&&F^@HL( zlLw9-o_Ug|o<5l$j=rb4FWEnN>U%rnhodJazT`=sditOz4;(!_`zKF5eSWL{JdgO2 z2cG$pchh%U?~vZ#A;vb3{P4Ulu;p_wb zlhYqgAM}|MJv{p-Pd)oUe?~ZZc;-o-ditPaU2yd9?4P`w9v}Lkhv)tzPd$Cmb<30e z(ZlJVoc@fdhyTC!?gRR^t1AC~Lhpp$A@mw*=p{4(K{^6T3%&Oa(tGbxMIZSqoj;wAy4S_i`7}@8#GBVV-S_*dUcPte{EhQbcb|AVpXTY!|EkZY{WpCcd%t+P z{#)<;HvfNnK79E1y4{a+Zl2z}^U#KVCT<M-hEtg%er*DytDdgT`lwB zFVk`Ko#o$p^T!r%oUWFw%eQs4%m**iarE|Yy?Ni`IOBA+Y~AMRIGGP#rsL>*)H1z! zJ~)})JRh7)Z=Mg%=ai18^J$*Gi8rr#`rh7OKIV6y-B0`D>CHO_oOS4GSwGFwaWWsg zOvlmrsAc*l-gxtL=fJmhwaf=E({c3nZ@qcvfOCK7c)E3)r#H_BZ=9}{`7}?*$#i`8 ziKF-VwBCIC8>g$~_HSJ+H@~T$oy+${@0ZTEeelMutCsmRPsho0ymjd~y6^Gk>CN-O z8>g#fKF!l{vU6yjj+3q1JRK+V!OL_UosU|kZ{m$NPj?P{TUX0`@G>1oZ~xXe@y7Y6 zW&eDEZ|iEA4_>C@=8pA^J$Kb zH(zd^-u$|ArsMfpS1r?<=Y!J+9Z%=eJiYm&>pq%C$J6_B*Lw4OaQdX<>FwY8Cf+)H z@O0<4sy$!JIW*4)zp0=0$J6=KodbX4^i{p|+5H%&t5^9~e1EOm{nJ-{)Hyfb=VKlp zyiCW@`K;<%cN4!FU(QV}Z`FCWu9o+&^4Q_qupGPHw$< z`@-voj;D7&tvAnyzkcX=di%G&j;D7&=F|1jJiU4S;EmJOGN0z@IJwV3pL87Ey8K&j zo)6wQT`jkN>*`gm_n3Z;UUVGY=TSecH_r!eoW81;_g?!Or>kW;|4sckx8~`5+Xrvl zx@wtE^K{%Q*ZFTR=g{Xvw=VzI`O^8|?Te15w}0!Kc-_0_@ZBe#-u|sO@4lF4UAlGY`e~lt zJRiJ%=y*Dx=INXG?h{YvE4SXfb@AONo{p1SZ{Ft%uOB*|?%bNEH_r#JA3C1Sr+Iqw zM=zg8yguo8x^+)~Cc|Lf3((&~6Z@u~LHT~1^`n0ZoT5p~YUO#j^z5QF?#BVRpCm%e$ z&!_d~`RISse%l{U=TCPI{Pjb})4QM6o9Dy7*TvJ@zx8!IeY?)z{iF9ewBEda_;){e zdi%HDyz}9&A3C1i{j}aZAO8BGjKUH;a$u3GN?TUX1?Z`$wP-rx4`b=7kF zx2~3(*FSy5EuV{>%lo$P7rl8t_|^F0qn3A;f9qm#X}x(qc;j^SDzE9^ zx-0IGU2n8LpY77|^quA3dh_jXzjQpk{afF}Z|cW6>{#}Rr}L-pZ9jbY>z~e_-sfOl zI*zU%KCL&;2XCCds+aRLiuDGR7KKhX9csd`o zOmDuw$Bom~^090Br{iQk`jF{3Iv=%6-^BaeHBWaAd|T&F=Yuy+SFiG#^Re!VTkfxO z;BQ?zo_=ikx86J-{;ToDM=c*){;jKJyt?~Q%Y6>5tL65$4jo5t|JIv#J~-p_RlVFV z>o!kso)2E8&nx{A4=VP3% zmaWUbb+ya~zsk#fTyeFp6<7Ku-nc&1vUA|yx?1Lgm+3fq`?tP{H_k^b`{xLJTUX0` z@G>1oZ~xXe@y7Y6- zbev4bcb_=A^J$*meES=xt7Sg?TUX1@p?Nw^wr=xuoXiKm%Jm%5`KV?3%6EBQ=s4qa zwY<0I!-v02=Sy$@)|+=ebeZnCq4&D2H_r#BA3C1i{;hA~t-}XTcRp&F-aH@NrhfP| zPv_5vj+g0p`c^%+tvAmH*X!cx?ce$)-a35n^#124tvAmHx2d1@$J6=KosZf&bhWIX z=IJ<@4_>C@=zP>NeG_lIdAf7p+qzojgO}+zdi%HDymP?0KXg3Zy3Ny@=Yuy+SIc~w zr{iQgzWc<{`+Qn&zWt5U)pGl{u9o{8^g~z6y>9DjnUD2kI*#7{t#9Iu^HIyrhi~g@ znGas3N1>8Ho(=9_2T z=IPDz!5gQmWj@W*aWWm>ed6dix%KAT-#A?@^Oajy%iaTc{m^mrtu^uQ+~UZ&&dn{xBLU*mMOY+dWLu9o@WWjc=D{;jX$ z=`x?~`o7|alj+U#X^xIJUv8e>{Eq6+JUX6kUA0Vao)1nRbUd9;^YrGo*7Q%u^Rup6 zrZ>+Ar%yVb&Zl{L^Zgz(kB+BXS1r?<=Y!KH9Z%=eJiU27&C{Fb(>%R-KKLEW{ifsT ze43|k;*VbZo2Ofcf9w3|eDLPe@pQR)dh^b~dNQ34-MVU--aH?iKInKlpXTYCc=MX4 zTSqO^`O~+S`>Dp!@pQR)dh^bqd3y8KZJypdAG}(|tNAof-^ANz^K_ZNTIR3jgYQ1^ zbn}|0H*X#tN8euVKV3h~)0^jmSIc-cpXTYC`0f)=m-(yZ?G5|wKJj!u{9E6|n=j+l zzQ_5u-n?=ArheKVPv_q{f9o|*KYBUe=IQ+LeDJ+4p3bLvdh^yXu9o?z_0v4PdGql4 zq~qy)nx}8#`8H4IE4R*{&IjLp;^{cK_2#XEm+5$V_tSdwd~j+RuWtX=H}TftgQsuG z&GRvDQ$IKvujWse={Wj!egEsHd3y7F@M;;a=F>cV6R-c~=`w${%wNq1-+kig<~2{> z#N%YVdS|^qTW{X@p7zuJ_C@F4I)CdmPshov^QZH{_qupGPHw$<>#WOsjO#}&)0;Og z>qm{J^D*A~Cf+(qKc23i*7@tddAj+{ z(`7s#e6NeAtDC3GcpTlj^zHgRte@uT&GW(cx_COD=INXG9ZSE>)A`D+^QZH{cb|AV zPHue@Z(Z~8^keJ&-+J?W&@fjUKdYq|JFD0oBA_b$gtSlX2^Eb9|boH{ZXnjMLR}`?ju@`Sd(GuE*Ea&3lgU z#_4L=y3NyZG9SE5$IPv69w*F1ew=8xl}?mqE! zoZNcz#^rVM==y1%j+1fga&vr|r*GnoH&2)OtJh^dc$sbE_|pc=fu>r+IqMqvNeh z$J6EJ>CN-O8>g#fKF!l{oty7*wRM}P<2vKzb>n<+GTpq^)yA8rJBQ}!IGGP#rsL>* z)H1z!^^P@v-@$yGb(^O*ZyvtK)iR&v={T8=w;mlw-(L3FJiU27c;j@n%%^!ePUeHh z(Q$OSd3y7F@Xmpbr}Jr^-n@Bu934-$Zu9i!`QX(uUd^X@dh@%R- zK6o4*PnVmgH?PKTFa6VHx^z{^YrHV;MFo-&8K<#CcgW`(`Ei@xq0j2yH7lw zPxJKV&BJfk=NC`!ep+w7{q;e|)7!uGP5h>Q`1OA2{M!eQ(F6&1v?^vEcH6Q$@e#~p0&YzFmJiYn7?S~J4&pn+#z5QF? z#Pj8Yr*ALktDn}J=Y!wWPy6HP{OQ(je|qzK+YcZ9{OD@A{aaVd&D$3pN8hgVZ@qaw z-9H^Kw}0zuxz9mAbhX^;wyu`>SWl+o=BZ=Mg%b416}`7}@8#GBVVeS0s@hmZN)XZO?oczW~van_-$Wj@W*aWWsg zOvlmrsAYQd?jz1PT`gO;c{)z!gO}+zIv=%6Z{F`HoN>BZwr=xuoXiI=({Xe@YMI`A zKTpQ#YT3GcTUX0`@G>1oZ~xYtZy)1yna*D=)0^*FKXjSShi~hfc=Pz+>3tttZ=Mfs zQ$OvGr}L+)`K#CE?gvL#cjhmz8|QE^YrHr_nlIW$kl$$aoK9Y^P*mg&u_`+m^P zS6jDvdh_Ptdt5E^X`YUg>3HkWadhXZ|+-PYCex^?KywLgx&i8tOnUFNG^ zm-*mjx_R`jTVFSh>)gC`o1-`1+`8w2&d0pg`L~Xj`QYhYx8CD?aJ?>`u5O+#h6!)pGNj`Z3S?c)Cn)f8%s@Xa22sjg$R3L5-(d*LdsA<7E3% zPnYR?+mHGD^-q`SzL)XVrQ_-5HBWDz4<1Ly)8*#r&8zYHq|0>c zHcxL}jqg5;)A=+{U&qtCPrh_@^K==%XtTUX1?Z|X-MGG5J>-u}kvYMGDRx>|M)`0Zu?bR2!F z&ZqU}`S7&ql)7AQE-8h{OUO#j^UEMrg#^dDe+7F&?-R9}dgJL-#Y&2 z>XU!#cgK0Kds~WG~c@Y@X-&S*71CrZynD^KkZM)^P%Hg$MbjY_}1}!nr|J? zr}x`Bo=@}ZbUyl7<>k3j+ZP|Y+&Z4We(G^mL^&z*;cW?L0 zr~7HW{hPNg9nVMI{nKSWtGezNUFOeMZXM4bx5~?XR5!n>H}r#(+m9~qZ9l7d`dIQh zoxgGSr}cGx;?0-weAM{XjrV+f>v;FkIKFkfb(?P;&u3LH=h?b(K5G56z8XJ#S0}6b zSLCmcRbJLrTVJ1ax%KAN&C_M)uEw`+yw}CIj^~4u*XexBUwv<^xi5U0Z{525)$LE0 z`Ka-&8}H{3-#T7D-GA$NKI-gJL z-#Wg1x_>_S_QAJ~=kIyKw~pu2eCy5YlRuu0ckXn2>v%qU+7F-RTmN;f;=GW#`&o6tsCcqlUv90!O5-T`Pe_cbv&QuTgUTJw?AFxqsF&x{CoB1Abjh1{dE7W z+;u!Y~5A8{BG1IU2fg^@mcL>@o62;r}@_LeDtxJx7>F+o(~=0 zI-dWkuKww(ysWEkf4a;cx5|r;8m~{f+&X@HO+T&U^|Pus?3a%|m)cDp{i#UxxqkTY zr{nq1@vYq51-~+|8@JX)~nAKeU+EzLG6C&pDy=( z#J7&``-pEH&qs}KeKo$OfBrc4hmYEL>&E%0@vR%@gOgjw^TEljH?MAEzZ%bnzRGo9@cexq@vYN0tr}=d{ zAN!SCw=N%?+&Z3*e(M@z(Kt_Ou^9&9~n7OCP;19nXi3Zyn$L z<6Fn`(LcU*eEZ;A$F~o@bv&Q#<$Dz0I-XDS>vTT$>s(v6E}!OG$Mexo`_u7!==j$0 z{W-$7j@M83-#VU8_tQF_PxGz2Uwrh#r*%A^=3B?}(NFu+@qFm`$EI9p5^hzvl+uI-XDSt>gLhep|=$X}BnWp7*)?=ku}XD{j&Ev>)@DZ(VI)d-^<> z*L>@0{{P$O)AvO`KEKwbH5(p9p8HMd%9oq znr~g*>u&1DyzakszOsIBt*ebU-#T8_53Y5!@#b4^eoy;xe$BVu{(IUFpXS%;_N5uTf8w~m+fgKJ%Fy!qDga=%YnSMzDUb-djBZC%Z$`PTdX>W7a^$Md1%TE`pLk6NbV z`OtB#<5zk4-o&@AUX2&uI$rL6T37RFew}Wf^R&;_&6k^ReN`{NAN14ybi7Q*x31Rb zj^#Y@t>ag@_}0~YdcUpX<>p&o&0o$(|0}L^yz%B+$FK6TUtH_z)p+r(H@~O-xIfLe z-u}*0A2J=!hmKpP8`lr6b@gg|IUlvm2QTYq#Vz-jjx$bI%XGYKTyDK<+`5ck<=SuS zYJSc^Ez{kvxb(VeM&2PSSyu2D;d|Ovrzxmeja`)5vs;+(U zUFBuJYV-6Xw~m)_a_ee7`p379U*)AweCuj|e!3P8&!_uo9nYuv*75D5pVsm1gKr(r$9oUoI-XDStvA1?`|bOvpVr%dPy6B1 zeCula(toc@$IEnl>uUb`!MBcIv*~O*46D}pRMENz4iBbY=1gmnT~H=E$?YR{F-l_ue_)I z@M*qv=gVI|YMGAbL&vp_UyU#CbG58bJRkasTi)|@yxje?uIAT#>v-9^xYpIin{OR2 z`#j=WR~v7>b-djBZC%Z$`PQ4?(|+9F=38(7J?)22^R4rz?`c0fmglhf*4uwi`{C34 zI=%Z>%bWVd@o#=LzI=Yw@}@pl+;VQ~cfG0SLNUDL+jG@qffei*6Dog_kZjDT9?0m+MkZ^KJl&Njd%a8 zg>v%r;T-D2SOULt}<6GZ0-{;_vT|4Bnnz!h9KKjSEj_0EveCv2V&9{zkpI?=~ ze;(2Pbp7a)j$fyD|Et>lsr$mm`QTf(u5taSWjdY@9oKsEd)iNbj`+9U{(IUFpXOWN zUe4FP^xx~!@qFm`*75xHgKr(rr}@_L?enYh-(K?^()FWHx_;K_-M{Cd_4aW-c>Z)e zAN}B4$Mb2vbv!@os%1Lfx^!IYc;nqq>v%rRx4!NDvSj`6;ZMi&q2pV}^U)8!bv&Qu zTgUV3ep<&{xB1rbe7c|3@qC(Z9nYuVC#~c8G~YU&Pw%&NJfG&*>3sCF;yhpV-Z0*L z>wH&v@xisOUX2&uI)0UlZ(ZF!YPrwNJkP_5JABt0g^n{$SIcy~`)1u0xA@cXtGtK% z}f#UwGrq^Y7=Ob@S!sTUYaO4*F>wzseiV$9!2o=C`gk-hAtLS^v1!SM?s= zcf98n2K`LSPl`;h(f!8u;y1c~D(PMA1R;v|U^|C0Ko(N2~)MdIX%VNaPD z_Ed=(J9T2l4oHmmX%gdo+C=m8(`W2-2|j1c^%;Ii&G*b1J5xfRXUX-<(brjXJ&b*w zE!VR~^Yn9M?Cc3X=gRdt6Z)F50~7k1v2!Qr=gswb67=)u`g{q#7s&N6d@q>mVYm~`@c;7UEJ9VyamN+1B^F;IX`K;a| z!N+wf3V=UlDR`@L7j?wz=A;y#Jy=N&xExf8?mdh~hM-7ldJ{qmi$`zQ2azYmD^ zz{G5n1lEBScf0&rHUYeM(mnFvI-;&02qw zn6=)Jn6Wn|#^X&1`mFV*yXoWo=EU^;mc(lkf0k%|_V(uqT(9@mjJ+-K7m2qgntw;M zcP9Qa@vcPk^Vt|?J{!Z#vpCFrHpa8-K5yK7R^FZH^X5CBt@k9F*Vl}_H(?*+^}b(H zJBQa~Y|ab!fn2{o@tVZwA4+^M!T0c7e>l;+zP2(pjD370*W=l>=jfvun>BDB&-KR= z^iSsc6N#ImORYY6T|#fVi=z@oX-`m&$IbHPd?T=GHVR;jl`V! zHxt8rD>01s3eI~7*Y%4tXL>v;F+Jk=(Z8Ls+2idO<}S&@9GV#Bu8Co8otQOnx6Adc zacE+=TPCK*yC60Xnd-kCFu zb-Le&X6$K+2PTeAG(ThW8Jn@KM6XBx+lsqYwL?d`x0z;_ni}eJAnv ziSH)ped%*SmK2heMI5D5S zlO*H#-6i@-qvP6d=AJAub5EXdPR?cKo+5$k^E_q7PL((yaq8VTey53c+61os<~=dY z=@zE@rGGoztYIJaOFu)#PM@HkDc5IA(C1xs=G63gH=ZT+%sFdf#?F?QIcHDI*f|pR zg|jc5jPK{AdFP7je)@Uo{kX5VbEe00CBB&b9GHCW#CZ~*NPXU9`^0$;aWnUPiJ5!; z#LT@wV&+~jVP80(1KfqeTsU#j#6=P~e$Exwe(uY~GIyBnmpfJ19ZF zWUeofpno~@FO`};zYQ*(dgfdvF=GcOX3k|3Gj_Rzec|j2C*v=lIaf$rDRIR_^Uf95 z{aiV7X6&=+Z#;TG?ki5lUnO&{n)poSUoE-$k3~NuI!?x4J!98Md^+4Ulbiok^lL@O zT|3vmwHtSxTwgbV^Pa;Uk@-JL;5>)=#96nW*Y@)qI@jstdWqTNZzuYB_1?yHJ-yBM z*6iW>2^>HAoWTte^f~k2NzmsEZ*k614BsMwn=|>{1Wv|J&$moWZ?{S`KYPA)0yk&lC z0}?oX4~q8S1g`xqyNu1dZah3Eo~QFi8xNfOgqwHgLlOrk9=aRH?;_D2mcX^&%zb!b z<~|~!Py3mGN!<#{+kc zT#pA%#(PekQ}g{i>qEbu(_3VXd7iu9ex9Eg=1IG8PtNsI61b=4`e_N=({ufd1n%t1 zm_GMg#`Jme)Z>Ah-`L}UJ8Q0=nQ%@xey7gZvl6)Wo4LAXh zfqQGPUD3(g*KGJgK-__~bGpD_p6j_tsp$ErEM`uKyx|d&e^7oZh*NIj47}9uM4K=6XDE@6Ppm67Ca@-+MFm zz67rQX72kFo@0I5&&>Tm0@wTbV8%X_zTWX8Ucz_s5uGxn_n?xni1cXY1*Hi7H< zn2h~h;ya0*yK(&J_VL{?`ecm>&NzF+H{){f8Nw zJ+>eHM;V(lvd(s{=Pc=-Bf96J>wn6exi_AZ_R}}+pTmyFza+-7uYS2%uj`?z4x@z_87uop-@9>+^PEdBf`$Af;})Z=l{IcBkPJq zN{yQ}=iP{#p64Bko1W+0i1QrcJlD9>gr9fwX_McZjN|9I#HfYa z{g`FH_e>s-(CvuFlbUpCNOmH$Hv8JU6($-v?*?GiLso6K6^^Prp;Rvt*o) z_0AgYY>B&M{G{1`_iG-`JlxDZdt&CEBGKnHYn~&4>+_zT&zZovPq^9Jxe~a(p9f~_ z+==rh&a)fG?|jkDpTM=Bc^3?GfrNSbrJp8yJ8gn)zjw~LjmL!&!`^8*XXiz~YxMD; zJMRl;&h*Bo@7Gl_Z}x!e`+bp&T{MBaSgsFB;I5eY7f+45Z0bv-#>x18Up3$FyGzFF zQVHDgb3N}-+~u;~rK92ceS4XV9h}f7j^C9s|FRj!wcpIWTw?Zk`NZ^cg~aUfiiw$f zr37wzxpD&6b$?&s`1v05eMg@&zenYl({&g7a2+?>g^5;z&(zgL>?-#c^G*G}O2_tKp8ZzXX3du#69brSl-@tb>h z-2|@vX72S8v&Y{~OfT0@%pPx$n7O}`z)deVOyIhn--5gUUd8eIX2!pjI5P2#1bxo@ zMhW`YGd91g>2oh`oS<)IZ2nv__wFW%xpy~BG;d$HUeCOnXYS1s*7@CB-y%W3RjzND zp#NU3Z=INT=P>i`9A@6d!^}H-7|#{XbAprc{+#sQ?CpJt=KT)Fb-&K#gPHrGME6S{ z?++yG<2LE__Y?H@WbFM3`ipb@k_7#wxqf-#?3pvne0L2qzX68v=UMy0*%xlUGj5xh z?}is9aQvJruKnDX@wi>0_d_4=7bSW>x6hnAB<`3vG|~K>qTM-x>-FxEv0?7IFx@Zx zZW*67?8AQV9_=29dnN9fXg~UWGIsAo`^}vDE@Ssg)~Ef9_x%&N-p>Ov_P_-0LAic# z0{4(yKQw`RSgs$Qz&#?@!#pxE%%c*+JUTH9-9G8|*?fOqnm>0ucmKTUL%;ob?SA|7 z{s-}UOyY5g$0n|snD^k-67yc6^P%&3e3&OBoKy2p%-E9>xL)tc8GA|s_tacJErEM_ zuAh;>Ju}zOO5nQg_ce~+bu#yQ3Hy9@uE*m!iDCQug8tl$jmPB@<8k@KcswsLJw87% zJ+>eH1sR(?wjce48Jjb*&WmzAXG!;*(LFcK>+7nSJNK%2y65c0VaDSn3w!b8@$h?r z{)f?KP5SwAJsvMjOph;1OpndKJlgCL*Xz9^V{=BhSLS-o2=}U7&l%(XDA#kwxUMgo zu`4AGPRzZ+@uS=4tHX@PA1`eGUZlS!H->p5fk8*@GPq=N?b7E&nxEP%)`yxyC!Dt-4ZkR0ttU!nl*xmG(=K6&R`ulP{&pF+5uJ5^b|Cne$y63$8X3sxd=Kd&o z_Pm{#Uj8XDo}E&O3Nj*J3KJkdeLlW&bdwya9H)s8?yK!?iPfp6 zo3nj-0yp>W-xKzMv+mrx|4884Z|44IV$R^liRtCP5_1MWNzB~;PT;1OpC)i!PjAmi z%pRVy8^@16XYiy1ea`x43HqGj6B6_}>z^m+bG8pn=zH$be#bwqAARoK{uysS`tdV% zyo7a5nClZH=qJwgi4yda=K3UwBjP(ApH2+xxtMqHu=II%jtBkYsmG&v&&?w<{(waD z^t)&59*OSPxy;-LCc0nx%)NiYK2DxqPL`nGGh_En&`+7MQzYp3&e;7D_VLk$)z`gJ z4@>`8uE&FZpVZ^wT;{vuR2gqS`WG{P>WsIa`!aJ6Nc4W_Gxs!!-VgnB89QylKF^Tr z(%Jmr&^t0sp%nAD0a(&hW{T#VIdxCzhT%R*RKX@g$k=$$&z7*sn{`@BJ}ho%2c zuE&FZ!(5L?^UmeMIh#i%nx|hRV-HJozwXP-eRQJxrO(_)ChX(J>E%WV`a?4I@C5y) z8M{e>{?Lp)B4Hml%h<5`x=hA~rQbYb<3T?-W8>jmem7&cNVFgQRvEixqW#>LnS1L* z?}t8fe=pJd`Tc0ON!%`R+eG`(?~t+EC))4OXm?E9IdP{%^LL4M*95NDyIaQYp1|EB z*Y`}|&Xw!=-oTwZ*9RtWUH9{h-dy_gXyrd8Xex$xY*_w&?zKTi)@JRZDwKO}iP9-0`I{;*un zn)HX~dY;KgB&NqlCZ@;c9~Eu(i0kzpov}G1+#lq6&ItFIT+bQf9-HeqV_esd%h=q* z$0z1q;rP+*^9f?*l-(ll1-&J(`q}yjdukF{* zYy0){+Wq$PY9IFdiuC--#2+PIm3VdHj}xy=ye2`n4&6F4_n^egy?A2gUO6#yub!aW zC*3}qKOko|zoAZ_nBPjx)17nobB1NE`_%jC|Gnqy(&J$X-03p^lnLD1@>}+;37o&L zjhnIg4LcrJO8EQM<8i74-9G8|`OeIp-_Wnm+>0jKkM5k?@7EAf|!nKzy@ z_U^=t&6+dzp2UnzPw!2zK!&h_Dmk0iDd&Cl3Jm$8o} z_j)sD^m(qHl+cHM`Ob4aj6Ur5m|o@`JKnEM;HH;PCvaW&dExl=dCgj%Nz7XFoiJzl*~FaX zQxmh+=MuQN7oShyrnkcqvxhGwny1g+<~x;-^=8hFWo*8)rMuw8O5!Pr=OnuRYP7Erk4)%u#=eozXP-y&-wgY$1n#I@e>;Ke`mZzg zHwo(;o$F!r`M0?qMxV#zde*Sd=Kn5ZI}^C?55xRZ zV!Y@+XLO&luAQUvZofY7f6f}i{7YiI{xvboza_?t?ma{Io_W#IyZ2P{FUk1q5!dV8 zCS$|gc47WKYYy`t%i47BS^6h4K0S5)pV6kr4=1L__WN%-Aqp|18?i6CX;vKhb_OHqZWy&1bCp+3y7TXJOHC zy`TLvK8$sam+N7UpFF*sAbFS*CWbjt@|@9$lV>l@pCtM)Cru1UX8q z*BMg}qpvgNdKmh9^ZVgL3HsM^{iOuo8Jpi6d}nNaW1Tr3AIUoJ%Juw4Y5pu3ALb1i z8>ZKDF1YU3xmc(7d)9c*XX$Lo^V_8PvqvB19Eo9iJ^CBd+glR)&@bN^o8LV8u-`l8 zaXd6}=fwO*YW~X6hPg^&m|pK{(Qcf;b-yzAKNo4cz&n&lGSX1%`-}W zO`i4nOz@q5uW%T?ug=&oU0*VumrA@gV}~W$Z^q^`He>S{>wcUouJ_|yt<&|T<9(UL zWfKP{+Hc0@9W!I|F0zh(==1O2%sbIK_QUt+j1ANE<>Gny#LkR;J<)zMHt*OOn|EyY zbA^nZCV}hyTrp$ASm#Q)zH$P0m0Vvnf%BX>r{?>49H|3 z>-mnv_2;Si{ybbgbBEDq*Vo9{tYM$myw)8ab8+-E)oV)*< zz5Pf0W3n6pTDGp5|0cc1)iOWp()n-FjkizheChnx z>FtM?+lS6i&4;d*TQ_f=Ze8E_FP)N;>P%j-D$o@Rcs zdFpj}9oM@3wNLBj(be)g-MR9|w{Dz{Tc=yEdAe-AT5g>WT`jlHN1wQLy78?;e);?V z$PGH*_SvKx=dYIOeCb=&&pO??;Bo78=Y_W}z58jM|2o~e?L*&L`)Zy4I^DYMLqE3s zTj#${w{H8a)A=>;T)WS8Iv@P1K77~n+BzTW;@9b`@#TE%XPwT+y7+ZEAN-cS}fpF```rSot7?apV(_#@WnN3GGnx<<$Amp{IB>$a}X)|+>p z&C`2bdiP1kZ!P-uTE`6CbC^_(|P$2$i)zV+ta|K?lQC!J60eLna;pVswB z=hHeLy8Sj!$Lo`hZ~eRV+~SW}qvN|DeCy5IZ}Y9|lg_90A5@?CA8yd~iRaTgAG-5t zo{rZi9p8HMI}hLW{5MZ;{@67^i-=z1tc)D|F-8pR1_0v9d zK0BA^g^q8%`R)@>cRsB&S8_TpZ1}*Pk&zc)A7!Y z&cF5Mx9anPr{k?l=ihptKfceOj^DXFZ*+X?&D(GDbi6+4_|}hF`e@zvMeBHd^2fK{ zy!|#$$6J?%aBpk3IbV{QZWf_qz1%la6<8bbRYa)boJXfAeo&`fdG)H9Efcg>Sui=ifXXuOB+T z^{-Z+_@mb7czrhCdh^bcKc0@)CmrAVG1VvjyKD63_1`@G2jxTm;U4PxX&*YDJ_mg3 z&FiOmx_$AXJBLlWKHG=Rr+-J`TW?-J&C~6R58XLz()HOsbUs^$?0WCe@vS$npXTZI z#fR=3HtG6lANoc|De|+oKZC#(OH@|b)PxJI%m)?ES@mtIN zrQ=)IPwV>Gq#wET-#+wSm)?ES@y?BoZ{51B>%aBpx7M6Lz1O98pLD!)qvKoGPwV<` zy+23%`}0C?-g!1p@6R>8-vjjKt=l{u@4nFSoAh26Pq(kuox>(wKkY;3vvc`;(ebS} z-+kig&Zl+f(|Yq;%k$hk9e-@~kKd$spLn`+Xx%w%()HOs^rM#j@uB-1t(-_7*{Cn)%74X|u21hX{;l(&`@A;a zx^?M%T0hQxS^9K;o2Pd_d|Gebew%OIx^zCRALnzi^xvN^KCK_Q@a;qI_XEA(59@U6 zwhx_;bED&1Z{E7i)9s57-8r=0yndRe>Si`>-uRQI$kZ~)pGOJ zZJusl=F^=|>&+W)o_^HfyZ+pwmVLg|viF8sK4OgzUf<2PuAkQR(|Ys9o2TnjEq9-E zyuRu9)~(z6QFT6iTK9QsA9}A#?>^~x>(TM8>!)@7wBG!&%loi-daq0GKIwSpNyoQt z-PZNldVgN{_vdAu?mXLv-k)oFzX#~eTeo>S-hH9tH|f1Do^D^QJD=8@*H80wy!Gh# zO}g>+p*tV7?0nR6^ZIF?j#tZgwcL8Ii=*3D>&~b3=8ZQ`KWcd|sb!xdwT#yv9p8HM z#+#?(`O@*NA5s0|y&u-`I2o^&`K$RgPsgieyjp(y^1f8dN37Amx<)^0jgIgB@vZBp zb^Wy7yz%Dgy)M1`q~o0n9pAcjTi53%{mA7U+lS6aA9Q@{&HMgoo{rZi9p8GN1HR9V zj@Kt0-+J@*+dLhwPddK!W9qrZe|L?J*JtysH*dfE@pQaC>G;-vP<`TmxJJkKy7<NUM)9o-R9}`Wj@{cwBEe&=IMC% ziH_f->!*F__Qi+p95(6tX&*Wt&pjRAdh_~eo^D@!=+3A0=JnG&y?N)+JRPr=@oKqw z{WjmaTIQpc``qY#&aLZ@&Zl)gbh&vto-ZBWy6=nDeP6VW=PTpYa`XI}Z(S|(QOmxU z)w1t#wM^F^olon{8*iSD=S#=8-n`F6^K|FNhu-Hy$LoWRZ(TpF>!jNXOr&y`Qux+E}j1-oiCpMIvuyl%f8gB zd24k3&C|`7^{-y-Z}D%veeiVa()rUpFYEO7!PE6a=d(`NNAqgqeAekWSwCu-&PU#) zn~&eQoc}s~#T~NibAn%|^O27&|JM1e^1`p`W!?5!jW2xP7k#eN`#IrH?><}SzfQMq z`_N_Us%8FanQmP=AG*9wZy!9JFP+ai-9C-uTjx*bze#t$`Qux+E}j2zbkCW*JZD?# zCF^tbx!Ryx7k`{S&6f{Z{yC6ZUZeBZr*-IRxpnjCYI&W$bLodKzIEet+$P<9#jWz1 zb@BY=b-L$(FTVBFc%A<`-Tmdm*Ymbc?>_m^t+!6+-#p#A?L+4)x86QDx^?L?9k)(5 z-aLJ4Ie#^u*6ojuTc^7}c>bGoee%b*Ze2S6P5RDdzx?s7TbIs%lkRiPAK$ul>HOE} z{F0{m%T0miFI&5#~QS zG=DAm_!*y((;Rrnekb|Gz}5FV_~iQ?JpU!z`RjyHXZ(1_JKir0PO|XF-~aue$$!vp z{@S&`uIv3yySNPe{q3jRI{*JF%gFyX>A2@{p defekt - -CalcVelos = F! T -VelocityDirections = (/1,1,1,1/) ! x,y,z,abs - -CalcCollRates = T -CalcReacRates = T -Particles-DSMC-CalcQualityFactors = T - -Part-WriteMacroValues = T -Part-IterationForMacroVal = 100 -Part-WriteFieldsToVTK = T -! =============================================================================== ! ! HDF5-Output ! =============================================================================== ! PIC-OutputSource = T - ! =============================================================================== ! ! CALCULATION ! =============================================================================== ! @@ -82,29 +52,16 @@ ManualTimestep = 1 Particles-HaloEpsVelo = 1 CFLscale = 0.1 ! Scaling of theoretical CFL number -! =============================================================================== ! -! Dielectric -! =============================================================================== ! -DoDielectric = F -!DielectricTestCase = FishEyeLens ! Greenwood 1999: A field picture of wave propagation (Maxwell 1860) -DielectricMuR = 1 -xyzDielectricMinMax = (/-1.0,1.0,-1.0,1.0,-1.0,1.0/) ! -DielectricCheckRadius = T -Dielectric_E_0 = -1. - ! =============================================================================== ! ! PARTICLES ! =============================================================================== ! Part-maxParticleNumber=5 Part-nSpecies=1 -PIC-externalField=(/0.,0.,0.,0.,0.,0./) Part-FIBGMdeltas=(/0.1,0.1,0.1/) - PIC-Deposition-Type = shape_function PIC-shapefunction-radius = 0.3 -PIC-Depo-Periodic = FALSE Part-nBounds = 7 Part-Boundary1-SourceName = BC_z- @@ -121,46 +78,12 @@ Part-Boundary6-SourceName = BC_z+ Part-Boundary6-Condition = open Part-Boundary7-SourceName = inner -! =============================================================================== ! -! DSMC -! =============================================================================== ! UseDSMC=F -Particles-DSMCElectronicDatabase=DSMCSpecies_electronic_state_full_Data.h5 ! when supplied: doQK=true -Particles-DSMCReservoirSim=false -Particles-NumberForDSMCOutputs=0 -Part-TimeFracForSampling=0.0 -Particles-DSMC-CollisMode=3 !(1:elast coll, 2: elast + rela, 3:chem) -Part-NumberOfRandomSeeds =2 -Particles-RandomSeed1= 1 -Particles-RandomSeed2= 2 - - - -! HDG - -! for chem reac + HDG -Particles-DSMC-UseOctree=T -Particles-DSMC-UseNearestNeighbour = T -Particles-OctreePartNumNode=80 -Particles-OctreePartNumNodeMin=50 - - -!xyzPhysicalMinMax = (/-5.00E-01,5.00E-01,-5.00E-01,5.00E-01,-5.00E-01,5.00E-01,/) ! - - -EpsMergeElectronicState = 1.e-2 ! merge QK levels when difference falls below eps -Particles-DSMC-BackwardReacRate = T ! automatic reaction backward rate calculation (when foward rate is defined) -Particles-DSMC-PartitionMaxTemp = 1e9! maximale temperatur für allocation von kb=f(T), sonst abort -Particles-DSMC-PartitionInterval = 50 ! dT for Intervall generation of kb - - - - ! =============================================================================== ! ! Species1 | Al ! =============================================================================== ! Part-Species1-ChargeIC = 1 -Part-Species1-MassIC = 100000!4.47908080018E-26 +Part-Species1-MassIC = 100000 Part-Species1-MacroParticleFactor = 88 Part-Species1-nInits = 1 @@ -168,75 +91,6 @@ Part-Species1-Init1-SpaceIC = point Part-Species1-Init1-velocityDistribution = maxwell Part-Species1-Init1-MWTemperatureIC = 400 Part-Species1-Init1-ParticleNumber = 1 -Part-Species1-Init1-RadiusIC = 0.!20E-6 Part-Species1-Init1-BasePointIC = (/0.0,0.0,0.0/) -Part-Species1-Init1-BaseVector1IC = (/0.0,0.0,0.0/) -Part-Species1-Init1-BaseVector2IC = (/0.0,0.0,0.0/) -Part-Species1-Init1-CylinderHeightIC = 0.!0E-6 Part-Species1-Init1-VeloIC = 0 Part-Species1-Init1-VeloVecIC = (/0.,0.,1./) -Part-Species1-Init1-Tempelec = 400 -! =============================================================================== ! -! Species2, e- -! =============================================================================== ! -Part-Species2-VeloVecIC = (/0.,0.,1./) -Part-Species2-ChargeIC = -1.60217653E-19 -Part-Species2-MassIC = 9.1093826E-31 -Part-Species2-nInits = 1 - -Part-Species2-Init1-SpaceIC = cylinder -Part-Species2-Init1-velocityDistribution = maxwell -Part-Species2-Init1-MWTemperatureIC = 1.160e7 -Part-Species2-Init1-ParticleNumber = 2000000 -Part-Species2-Init1-RadiusIC = 20E-6 -Part-Species2-Init1-BasePointIC = (/0.0,0.0,0.0/) -Part-Species2-Init1-BaseVector1IC = (/1.0,0.0,0.0/) -Part-Species2-Init1-BaseVector2IC = (/0.0,1.0,0.0/) -Part-Species2-Init1-CylinderHeightIC = 70E-6 -Part-Species2-Init1-VeloIC = 0 -Part-Species2-Init1-MacroParticleFactor = 88 -! =============================================================================== ! -! Species3, Al+ -! =============================================================================== ! -Part-Species3-ChargeIC=1.60217653E-19 -Part-Species3-MassIC=4.4789897064E-26 -Part-Species3-MacroParticleFactor=88 -Part-Species3-Tempelec=1.160e5 -! =============================================================================== ! -! Species4, Al+2 -! =============================================================================== ! -Part-Species4-ChargeIC = 3.20435306000000E-19 -Part-Species4-MassIC = 4.4788986125E-26 -Part-Species4-MacroParticleFactor = 88 -Part-Species4-nInits = 1 - -Part-Species4-Init1-SpaceIC = cylinder -Part-Species4-Init1-velocityDistribution = maxwell -Part-Species4-Init1-MWTemperatureIC = 1.160e5 -Part-Species4-Init1-ParticleNumber = 1000000 -Part-Species4-Init1-RadiusIC = 20E-6 -Part-Species4-Init1-BasePointIC = (/0.0,0.0,0.0/) -Part-Species4-Init1-BaseVector1IC = (/1.0,0.0,0.0/) -Part-Species4-Init1-BaseVector2IC = (/0.0,1.0,0.0/) -Part-Species4-Init1-CylinderHeightIC = 70E-6 -Part-Species4-Init1-VeloIC = 0 -Part-Species4-Init1-VeloVecIC = (/0.,0.,1./) -Part-Species4-Init1-Tempelec = 1.160e5 -! =============================================================================== ! -! Species5, Al+3 -! =============================================================================== ! -Part-Species5-ChargeIC=4.80652959000000E-19 -Part-Species5-MassIC=4.4788075187E-26 -Part-Species5-MacroParticleFactor=88 -! =============================================================================== ! -! Species6, Al+4 -! =============================================================================== ! -Part-Species6-ChargeIC= 6.40870612000000E-19 -Part-Species6-MassIC=4.4787164249E-26 -Part-Species6-MacroParticleFactor=88 -! =============================================================================== ! -! Species7, Al+5 -! =============================================================================== ! -Part-Species7-ChargeIC= 8.010882650000E-19 -Part-Species7-MassIC=4.4786253311E-26 -Part-Species7-MacroParticleFactor=88 diff --git a/regressioncheck/NIG_convtest_poisson/builds.ini b/regressioncheck/NIG_convtest_poisson/builds.ini index a0e9e2cc1..7b30fd44e 100644 --- a/regressioncheck/NIG_convtest_poisson/builds.ini +++ b/regressioncheck/NIG_convtest_poisson/builds.ini @@ -2,7 +2,7 @@ binary=./bin/piclas ! fixed compiler flags -CMAKE_BUILD_TYPE = Release +CMAKE_BUILD_TYPE = Debug,Release LIBS_BUILD_HDF5 = OFF PICLAS_POLYNOMIAL_DEGREE = N PICLAS_EQNSYSNAME = poisson diff --git a/src/interfaces/interfaces.f90 b/src/interfaces/interfaces.f90 index 55e7ed37c..589462382 100644 --- a/src/interfaces/interfaces.f90 +++ b/src/interfaces/interfaces.f90 @@ -64,14 +64,14 @@ SUBROUTINE InitInterfaces !=================================================================================================================================== !> Check every face and set the correct identifier for selecting the corresponding Riemann solver !> possible connections are (Master <-> Slave direction is important): -!> - vaccuum <-> vacuum : RIEMANN_VACUUM = 0 +!> - vacuum <-> vacuum : RIEMANN_VACUUM = 0 !> - PML <-> vacuum : RIEMANN_PML = 1 !> - PML <-> PML : RIEMANN_PML = 1 !> - dielectric <-> dielectric : RIEMANN_DIELECTRIC = 2 !> - dielectric -> vacuum : RIEMANN_DIELECTRIC2VAC = 3 ! for conservative fluxes (one flux) -!> - vacuum -> dielectri : RIEMANN_VAC2DIELECTRIC = 4 ! for conservative fluxes (one flux) +!> - vacuum -> dielectric : RIEMANN_VAC2DIELECTRIC = 4 ! for conservative fluxes (one flux) !> - dielectric -> vacuum : RIEMANN_DIELECTRIC2VAC_NC = 5 ! for non-conservative fluxes (two fluxes) -!> - vacuum -> dielectri : RIEMANN_VAC2DIELECTRIC_NC = 6 ! for non-conservative fluxes (two fluxes) +!> - vacuum -> dielectric : RIEMANN_VAC2DIELECTRIC_NC = 6 ! for non-conservative fluxes (two fluxes) !=================================================================================================================================== ! MODULES USE MOD_globals @@ -135,43 +135,46 @@ SUBROUTINE InitInterfaces IF(isDielectricInterFace(SideID))THEN ! a) physical <-> dielectric region: for Riemann solver, select A+ and A- as functions of f(Eps0,Mu0) or f(EpsR,MuR) ElemID = SideToElem(S2E_ELEM_ID,SideID) ! get master element ID for checking if it is in a physical or dielectric region - !IF(MortarType(1,SideID).GE.0) CALL abort(__STAMP__,'Mortars not fully implemented for dielectric <-> vacuum interfaces') IF(ElemID.EQ.-1) THEN IF(MortarType(1,SideID).EQ.0) THEN ! small mortar slave sides have no corresponding master element - IF(isDielectricElem(SideToElem(S2E_NB_ELEM_ID,SideID)))THEN - ! am1) big elem is PHYSICAL and small slave DIELECTRIC - IF(DielectricFluxNonConserving)THEN - InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC_NC ! use two different Riemann solvers + IF(SideToElem(S2E_NB_ELEM_ID,SideID).GT.0) THEN + IF(isDielectricElem(SideToElem(S2E_NB_ELEM_ID,SideID)))THEN + ! am1) big elem is PHYSICAL and small slave DIELECTRIC + IF(DielectricFluxNonConserving)THEN + InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC_NC ! use two different Riemann solvers + ELSE + InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC ! A+(EpsR,MuR) and A-(Eps0,Mu0) + END IF ELSE - InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC ! A+(EpsR,MuR) and A-(Eps0,Mu0) + ! am2) big elem is DIELECTRIC and small slave PHYSICAL + IF(DielectricFluxNonConserving)THEN ! use one flux (conserving) or two fluxes (non-conserving) at the interface + InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC_NC ! use two different Riemann solvers + ELSE + InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC ! A+(Eps0,Mu0) and A-(EpsR,MuR) + END IF END IF ELSE - ! am2) big elem is DIELECTRIC and small slave PHYSICAL - IF(DielectricFluxNonConserving)THEN ! use one flux (conserving) or two fluxes (non-conserving) at the interface - InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC_NC ! use two different Riemann solvers - ELSE - InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC ! A+(Eps0,Mu0) and A-(EpsR,MuR) - END IF + InterfaceRiemann(SideID)=-1 END IF ELSE InterfaceRiemann(SideID)=-1 - CYCLE ! skip - END IF - END IF - IF(isDielectricElem(ElemID))THEN - ! a1) master is DIELECTRIC and slave PHYSICAL - IF(DielectricFluxNonConserving)THEN ! use one flux (conserving) or two fluxes (non-conserving) at the interface - InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC_NC ! use two different Riemann solvers - ELSE - InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC ! A+(Eps0,Mu0) and A-(EpsR,MuR) END IF ELSE - ! a2) master is PHYSICAL and slave DIELECTRIC - IF(DielectricFluxNonConserving)THEN - InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC_NC ! use two different Riemann solvers + IF(isDielectricElem(ElemID))THEN + ! a1) master is DIELECTRIC and slave PHYSICAL + IF(DielectricFluxNonConserving)THEN ! use one flux (conserving) or two fluxes (non-conserving) at the interface + InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC_NC ! use two different Riemann solvers + ELSE + InterfaceRiemann(SideID)=RIEMANN_DIELECTRIC2VAC ! A+(Eps0,Mu0) and A-(EpsR,MuR) + END IF ELSE - InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC ! A+(EpsR,MuR) and A-(Eps0,Mu0) + ! a2) master is PHYSICAL and slave DIELECTRIC + IF(DielectricFluxNonConserving)THEN + InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC_NC ! use two different Riemann solvers + ELSE + InterfaceRiemann(SideID)=RIEMANN_VAC2DIELECTRIC ! A+(EpsR,MuR) and A-(Eps0,Mu0) + END IF END IF END IF ELSE