~fluidity-core/fluidity/darcy_weak_bcs

« back to all changes in this revision

Viewing changes to tests/darcy_impes_p1_2phase_quadraticrelperm_velBCinlet_strongpressoutlet_p1satdiag/reference_solution/darcy_impes_p1_2phase_quadraticrelperm_velBCinlet_strongpressoutlet_p1satdiag_modrelpermupwind_satfesweby_2d_A_1.vtu

  • Committer: Brendan Tollit
  • Date: 2012-07-24 11:24:05 UTC
  • Revision ID: brendan.tollit05@imperial.ac.uk-20120724112405-fp4yicflsgtd2rt0
Correct some test case regression answers, mainly for 3d as some have changed slightly.

Add a consistent subcycling test case to stress this code path, although I still wouldnt 
recommend using it.

Show diffs side-by-side

added added

removed removed

Lines of Context:
9
9
        <DataArray type="Float64" Name="Phase1::Density" format="appended" RangeMin="1"                    RangeMax="1"                    offset="1012"                />
10
10
        <DataArray type="Float64" Name="Phase1::DarcyVelocityOverPorosityCFL" format="appended" RangeMin="0.0011856467684"      RangeMax="0.32043101139"        offset="1060"                />
11
11
        <DataArray type="Float64" Name="Phase1::DivergenceTotalDarcyVelocity" format="appended" RangeMin="-0.055243936325"      RangeMax="0.62011832284"        offset="1452"                />
12
 
        <DataArray type="Float64" Name="Phase1::Mobility" format="appended" RangeMin="33.188196388"         RangeMax="10000"                offset="1776"                />
13
 
        <DataArray type="Float64" Name="Phase1::TotalMobility" format="appended" RangeMin="5305.7840303"         RangeMax="10000"                offset="2120"                />
14
 
        <DataArray type="Float64" Name="Phase1::FractionalFlow" format="appended" RangeMin="0.0037230738064"      RangeMax="1"                    offset="2468"                />
15
 
        <DataArray type="Float64" Name="Phase1::SumSaturation" format="appended" RangeMin="1"                    RangeMax="1"                    offset="2792"                />
16
 
        <DataArray type="Float64" Name="Phase1::AveragePressure" format="appended" RangeMin="0"                    RangeMax="1299606.9094"         offset="2840"                />
17
 
        <DataArray type="Float64" Name="Phase1::RegressionSaturation" format="appended" RangeMin="0.052800192277"       RangeMax="1"                    offset="3216"                />
18
 
        <DataArray type="Float64" Name="Phase1::RegressionSaturationError" format="appended" RangeMin="0"                    RangeMax="0.090035911708"       offset="3540"                />
19
 
        <DataArray type="Float64" Name="Phase1::RegressionPressure" format="appended" RangeMin="0"                    RangeMax="1296600.7429"         offset="3852"                />
20
 
        <DataArray type="Float64" Name="Phase1::RegressionPressureError" format="appended" RangeMin="0"                    RangeMax="13249.516045"         offset="4224"                />
21
 
        <DataArray type="Float64" Name="Phase1::Time" format="appended" RangeMin="0.2"                  RangeMax="0.2"                  offset="4584"                />
22
 
        <DataArray type="Float64" Name="Phase2::Pressure" format="appended" RangeMin="0"                    RangeMax="1299606.9094"         offset="4632"                />
23
 
        <DataArray type="Float64" Name="Phase2::Saturation" format="appended" RangeMin="0"                    RangeMax="0.94357644925"        offset="5008"                />
24
 
        <DataArray type="Float64" Name="Phase2::RelativePermeability" format="appended" RangeMin="0"                    RangeMax="0.89033651557"        offset="5400"                />
25
 
        <DataArray type="Float64" Name="Phase2::Density" format="appended" RangeMin="1"                    RangeMax="1"                    offset="5784"                />
26
 
        <DataArray type="Float64" Name="Phase2::DarcyVelocityOverPorosityCFL" format="appended" RangeMin="0"                    RangeMax="0.31727346258"        offset="5832"                />
27
 
        <DataArray type="Float64" Name="Phase2::Mobility" format="appended" RangeMin="0"                    RangeMax="8881.004246"          offset="6204"                />
28
 
        <DataArray type="Float64" Name="Phase2::FractionalFlow" format="appended" RangeMin="0"                    RangeMax="0.99627692619"        offset="6596"                />
29
 
        <DataArray type="Float64" Name="Phase2::LinearInterpolatedAnalyticalSaturationSolution" format="appended" RangeMin="0"                    RangeMax="1"                    offset="6988"                />
30
 
        <DataArray type="Float64" Name="Phase2::AnalyticSaturationError" format="appended" RangeMin="0"                    RangeMax="0.59948015647"        offset="7092"                />
31
 
        <DataArray type="Float64" Name="Phase2::LinearInterpolatedAnalyticalPressureSolution" format="appended" RangeMin="0"                    RangeMax="1192176.597"          offset="7484"                />
32
 
        <DataArray type="Float64" Name="Phase2::AnalyticPressureError" format="appended" RangeMin="0"                    RangeMax="109590.86369"         offset="7616"                />
33
 
        <DataArray type="Float64" Name="Phase1::DarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0.0037138237804"      RangeMax="1.0109750632"         offset="7992"                />
34
 
        <DataArray type="Float64" Name="Phase1::TotalDarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0.97313071762"        RangeMax="1.0978039715"         offset="8796"                />
35
 
        <DataArray type="Float64" Name="Phase2::DarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="0.99380166303"        offset="9600"                />
 
12
        <DataArray type="Float64" Name="Phase1::Mobility" format="appended" RangeMin="33.188196388"         RangeMax="10000"                offset="1784"                />
 
13
        <DataArray type="Float64" Name="Phase1::TotalMobility" format="appended" RangeMin="5305.7840303"         RangeMax="10000"                offset="2128"                />
 
14
        <DataArray type="Float64" Name="Phase1::FractionalFlow" format="appended" RangeMin="0.0037230738064"      RangeMax="1"                    offset="2480"                />
 
15
        <DataArray type="Float64" Name="Phase1::SumSaturation" format="appended" RangeMin="1"                    RangeMax="1"                    offset="2808"                />
 
16
        <DataArray type="Float64" Name="Phase1::AveragePressure" format="appended" RangeMin="0"                    RangeMax="1299606.9094"         offset="2856"                />
 
17
        <DataArray type="Float64" Name="Phase1::RegressionSaturation" format="appended" RangeMin="0.056423550754"       RangeMax="1"                    offset="3232"                />
 
18
        <DataArray type="Float64" Name="Phase1::RegressionSaturationError" format="appended" RangeMin="0"                    RangeMax="7.3274719625e-15"     offset="3548"                />
 
19
        <DataArray type="Float64" Name="Phase1::RegressionPressure" format="appended" RangeMin="0"                    RangeMax="1299606.9094"         offset="3668"                />
 
20
        <DataArray type="Float64" Name="Phase1::RegressionPressureError" format="appended" RangeMin="0"                    RangeMax="9.5460563898e-09"     offset="4044"                />
 
21
        <DataArray type="Float64" Name="Phase1::Time" format="appended" RangeMin="0.2"                  RangeMax="0.2"                  offset="4188"                />
 
22
        <DataArray type="Float64" Name="Phase2::Pressure" format="appended" RangeMin="0"                    RangeMax="1299606.9094"         offset="4236"                />
 
23
        <DataArray type="Float64" Name="Phase2::Saturation" format="appended" RangeMin="0"                    RangeMax="0.94357644925"        offset="4612"                />
 
24
        <DataArray type="Float64" Name="Phase2::RelativePermeability" format="appended" RangeMin="0"                    RangeMax="0.89033651557"        offset="5004"                />
 
25
        <DataArray type="Float64" Name="Phase2::Density" format="appended" RangeMin="1"                    RangeMax="1"                    offset="5380"                />
 
26
        <DataArray type="Float64" Name="Phase2::DarcyVelocityOverPorosityCFL" format="appended" RangeMin="0"                    RangeMax="0.31727346258"        offset="5428"                />
 
27
        <DataArray type="Float64" Name="Phase2::Mobility" format="appended" RangeMin="0"                    RangeMax="8881.004246"          offset="5800"                />
 
28
        <DataArray type="Float64" Name="Phase2::FractionalFlow" format="appended" RangeMin="0"                    RangeMax="0.99627692619"        offset="6192"                />
 
29
        <DataArray type="Float64" Name="Phase2::LinearInterpolatedAnalyticalSaturationSolution" format="appended" RangeMin="0"                    RangeMax="1"                    offset="6584"                />
 
30
        <DataArray type="Float64" Name="Phase2::AnalyticSaturationError" format="appended" RangeMin="0"                    RangeMax="0.59948015647"        offset="6688"                />
 
31
        <DataArray type="Float64" Name="Phase2::LinearInterpolatedAnalyticalPressureSolution" format="appended" RangeMin="0"                    RangeMax="1192176.597"          offset="7080"                />
 
32
        <DataArray type="Float64" Name="Phase2::AnalyticPressureError" format="appended" RangeMin="0"                    RangeMax="109590.86369"         offset="7212"                />
 
33
        <DataArray type="Float64" Name="Phase1::DarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0.0037138237804"      RangeMax="1.0109750632"         offset="7588"                />
 
34
        <DataArray type="Float64" Name="Phase1::TotalDarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0.97313071762"        RangeMax="1.0978039715"         offset="8392"                />
 
35
        <DataArray type="Float64" Name="Phase1::BulkDarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0.51728996111"        RangeMax="1.0109748906"         offset="9192"                />
 
36
        <DataArray type="Float64" Name="Phase2::DarcyVelocity" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="0.99380166303"        offset="10000"               />
36
37
      </PointData>
37
38
      <CellData Scalars="Phase2::Viscosity">
38
 
        <DataArray type="Float64" Name="Phase1::Viscosity" format="appended" RangeMin="0.0001"               RangeMax="0.0001"               offset="10404"               />
39
 
        <DataArray type="Float64" Name="Phase1::Porosity" format="appended" RangeMin="0.5"                  RangeMax="0.5"                  offset="10456"               />
40
 
        <DataArray type="Float64" Name="Phase1::AbsolutePermeability" format="appended" RangeMin="1e-10"                RangeMax="1e-10"                offset="10504"               />
41
 
        <DataArray type="Float64" Name="Phase2::Viscosity" format="appended" RangeMin="0.0001"               RangeMax="0.0001"               offset="10556"               />
 
39
        <DataArray type="Float64" Name="Phase1::Viscosity" format="appended" RangeMin="0.0001"               RangeMax="0.0001"               offset="10804"               />
 
40
        <DataArray type="Float64" Name="Phase1::Porosity" format="appended" RangeMin="0.5"                  RangeMax="0.5"                  offset="10856"               />
 
41
        <DataArray type="Float64" Name="Phase1::AbsolutePermeability" format="appended" RangeMin="1e-10"                RangeMax="1e-10"                offset="10904"               />
 
42
        <DataArray type="Float64" Name="Phase2::Viscosity" format="appended" RangeMin="0.0001"               RangeMax="0.0001"               offset="10956"               />
42
43
      </CellData>
43
44
      <Points>
44
 
        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="10608"               />
 
45
        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="11008"               />
45
46
      </Points>
46
47
      <Cells>
47
 
        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="10816"               />
48
 
        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="11044"               />
49
 
        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="11192"               />
 
48
        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="11216"               />
 
49
        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="11444"               />
 
50
        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="11592"               />
50
51
      </Cells>
51
52
    </Piece>
52
53
  </UnstructuredGrid>
53
54
  <AppendedData encoding="base64">
54
 
   _AQAAAACAAAAIAQAABgEAAA==eJzLm9vWm7Xd2JEBCioyLryYdgXBl+puVtO+a+i45fBRsVft+o6fD4noBr/RctzAyRkpelDVcYbftemSL+Ucjf8tfsDcIeEYJXZeozJTyLH5xtblgf0cjkcY48w3T/jh8GauolWTvJHjgn/791xVMnCcfUZ+JtcLbcf4CzM3HjBUd5w4sSt2OoOSY/7BQLEj68QdOf0XXw5aL+goGilU0rSM3XHBpZVeNzZ8d5B/d76wfAfCfV/vb1D//N7QkftPxrWA4/qOp/7kii9113a0MtrJXBSv5ijy5Wna/0sKjn8euMu9uybhaH5aSmD2JCHHNbqyn0R3cji+PvzfzmrzDwcAVb5s6Q==AQAAAACAAAAIAQAA2wAAAA==eJzzSJMRLe/Ybs8ABh/sFcRzjXoer4HzW1K3NjxecNheYuue3SKBp+ybPgkmvZ94wV6gySPxT+wN+zeMs9uuyLy1Zy9fdDfm/3v7I5H/geA9XD+MFtjwZ/rJnH32exyS7lmzHbM/Yqx24tPB0/YsCe/3NmtdtC/6FPpbdtFN+5V8j8x2rnhr32ZzqmMF0Jzke6jmrVhbrm9QuwnOn7FbcWrPlgP2Gyo8svZoHLfX2XrzrNynM/anuqb+1s66ZO/z3rCnWf+xPbt2CIfuiff2EmcilmBzHwCzWXP1AQAAAACAAAAIAQAA3AAAAA==eJyry3RXKfRrtGcAgw/2u3o3pLeLZMH5elfjVfJEZthPbRBT6l+y1L7x3acrgpob7V9tXTGHr++QffWfdN3omlf2PlNa9u38996+Y/O/////v4frh9GvN2+8lP22x/7u/xBB9cB59oF7GgwfCq+w3z+57vWJrE323hrJGi0dR+wVmLMCPGtf2//pahYKBJpzbM9/FPP4hFze1MeWwvmZqkZdfIUT7e+7qG85dGiBPZuiy627eqvsc+/61S423GpvcIq9LZXxuv2zPq/i3RPe25vP3OiBzX0Ad4Nuhg==AQAAAACAAAAIAQAAEQAAAA==eJxjYACBD/YMI5wGAM+kJxA=AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/q0tla9zY2I/RhyiEvGB1D9VzxOG9mxTP6zwfcdafss/CZ6+wWzzdD8t+lyIcUWEP3AFcHreAJM/roAtWG3prT/fRnxuSw/DP5LINfh8Q8I/cJdfS2lBwj/+FlNhCUHCP93gUZbDQMI/ZktFHrnnZz/DvvD5Bqh6Pwr3qtOZ7oY/V8HVcTmOlD9M1jFsSd6wP3gGUVUwi8I/ax9gUyyQwT9ocjdXbKvBPxliQVsFvcE/naB2lb1mWj/NZqTJyvnPPybV7phgQm8/8zl+TAaEfz9NcKQS/ByLPy39DrGg0po/zmqc0HCOwD/XJmxNMYTCPxaMPKYETMI/54Xdk2FAwj/nDlPgaDrCP/m8fVg=AQAAAACAAAAIAQAA4AAAAA==eJxjYAAC98+2DEhgQckPW2GxjkSmG4/tQfyAOuZ9jEA65c4XsLoe/esg2sFX9/NeoHhTC6ewHSND7OfPhor7Gp7W2wUZCtiB1IXrfrcF6SvX/rMXxH8d/RlMv/j/CEynv/sBNu9s/0eQuoYWNpA9DvqHVn2x3aO/y+mIKvO+Nvm70R4nXtkyAdU9XcYHNldE9ZPtwd2Lnr/zWLMfyH3Q1cWyDyQuGncFZI7Due88IHUKd6w59gH1JZw59B8ofuOg8F72fb6RhQW8Bt/3Rjgb2Bcu5LYDua/txVtbADkRUxI=AQAAAACAAAAIAQAA8AAAAA==eJwzbd8SOdUv1OE/CLQfdlh+fcdFsRkODgwg0HHYISUlPfvuoUyHip7zSqE3KhyCvgQZus1scfjpkLH+Xet0h4vGsT8X/trlUBu6yMHJ/rDDr+VCcr+B5lyAmgcz1zBkI3uUTrSDwN0TW4WaChyMmCOqszdWOzB5KzUpiLc5HNhxUvVFxCyHI9Oa3GY/3O3we2LOsyT/ww4bY/bu+4tkzsPtvBvYz3s5/IPy99yc/m+acpLDquatf9rdix26TgiVBs6pd8gN4LIWW9rl8PC+dbmD3RqH/93PtJq1DzqYMsis/tp82EGlSBLFfQBCAow6AQAAAACAAAAIAQAA8wAAAA==eJxbrlHWUOx0wOE/CLQfdlhz+/NCicyDDgwg0HHY4eAU8+iCpXsdtHRWnpNZtsuBI+f9nh2WOxw+pi3p4fu2xeGPatZ0nlV7HS4sKPhh5HDYIfy6sNxvoDkXoObBzBW5t3wnN98Bh1Wr15qUN+9xWHbWSd9OcZfDmVcSFlyPtzsY/TTZdmLnFoeW3Zbe8Wb7HK7+0Kv/4H/YwXvu3n1/kcwJrn+zZuOjAw7/oPwzf5/Vt/Hud3ji9Xnb00u7HaLjjpzXqNvpkPhx/5JDDNsdfB66sGrP2+ZwYOt9prWJBx2miryJ+9Z82EGlSBLFfQCDoJTDAQAAAACAAAAIAQAA4AAAAA==eJwrjFThavzeYs8ABh/slVYuY79Vnwfnf2p94ld2brb9TpX3dXVX19rP23ymY7rpbnuu2Z2MS75ctBd9/OZdidRb+0fbJ9zJ//fe/gffv////7+H64fREwpW6k4R6bbffirlxPm2Rfb/q5Z/5fy6wf7bJr13D6bvtT+vpqP/J+yKPccjrUNv2t/aHwmakSAKNEdj2n8U8yxN+mRl0yvh/L2WIg5WGyfZ57wtexKks8ze5+wcru/XNtvr7b13oLbjkP0Z3o3hO389ti/bYdvmMee9/RkxjZv/sLgPACgKdyE=AQAAAACAAAAIAQAAEQAAAA==eJxjYACBD/YMI5wGAM+kJxA=AQAAAACAAAAIAQAABgEAAA==eJzLm9vWm7Xd2JEBCioyLryYdgXBl+puVtO+a+i45fBRsVft+o6fD4noBr/RctzAyRkpelDVcYbftemSL+Ucjf8tfsDcIeEYJXZeozJTyLH5xtblgf0cjkcY48w3T/jh8GauolWTvJHjgn/791xVMnCcfUZ+JtcLbcf4CzM3HjBUd5w4sSt2OoOSY/7BQLEj68QdOf0XXw5aL+goGilU0rSM3XHBpZVeNzZ8d5B/d76wfAfCfV/vb1D//N7QkftPxrWA4/qOp/7kii9113a0MtrJXBSv5ijy5Wna/0sKjn8euMu9uybhaH5aSmD2JCHHNbqyn0R3cji+PvzfzmrzDwcAVb5s6Q==AQAAAACAAAAIAQAA3wAAAA==eJzzUJNK+R+zzZ4BDD7Yf/CXKlvAsRrOd1nSPT/q2CH7A6Fnyq+nnrS/8MNczvbdefsvV6ubeT2v2af1def3pb62/5AUYKL077196cHf////fw/XD6MnNIgvn3hxj/2ehAbJnwJH7Zc8Co9NOXbKnkXb+Vld4wX7bT8P6rySv26/fIfFZUfNN/YnZErMWIHmXGL6j2LehE3XLCuzNsL5L5zKXqlt2W8fcdbrpLDZMfsG29Pzk3+ctrcLqz4QsvCivcgdzb4Pvg/se5SbF1pOfG9/sKV16T8s7gMA1MR0BA==AQAAAACAAAAIAQAA1gAAAA==eJxjYFARqN5eZM8AA8U6R3evy0XwNSxUZQKr7Rlec9etSK+1dzi8L04qqtn+QULqu7LlC+wTLL4XJJastmcocFScukoZrK/i0B5bBnTwI4ivfnMlUN5mTsSKOnsGlVORZwvr7RkWtH2dpdBq/+Do/JJ5gkvsFRKmWpV9X27PcOD7Ya8sEbB5CkWz4eY1bOSPytcuhLuvweXyH2V/oPuyt8qLhNbZN3w9dzDVvN6+YaJb8nOpTvsDMx0+TOfZbt+wtzljcyVQHcODVxm35tuhOw8AqVZRUA==AQAAAACAAAAIAQAABAEAAA==eJyrnvy5YNY6Y0cGKLhfqbPvxgkEn8slUTr/uqHjpcWmDXzN+o5fv8od2/xKy/FI6u5Arv2qjpI88ZEPH8o5erbHXwwplHAsE3+6SDJGyLGEJ/PVnAYOx5RrczX3NP9wWMXBefO3qJGjGZucvKmsgWN19YQzzx9pO87lf5b1VlfdcS7nLhejn4qOV+weGVg1iTvqtjdu2dQi6GiQ3BxoWszuqOUcoLWx4bvD3FNi2RG7EO5LeGGcuPyroaPIlXinVQ/1Ha3XHLy3vkTbMTqsrfvNYjXHEBbZ/CYrRUe19k6mHyySjvaH5NTmrBNy/DGf65HYLQ7HB1pcYScv/nAAAONbZP0=AQAAAACAAAAIAQAA/AAAAA==eJxjeCasZrlgoQMDDBjddQ+tWY7gO8h06nyY5sDQUe8QldfswPDgfEz+t0wHhgT5XKuIAgeGP5x3fz6ocWB4VebCvx6o7sTfy9z7uxwY5GbMrMqsdWA4ENYg/Q6ovkXL8BDPYgeGFQ/by5qA6hpcbI26WxwYBNjmmc6qc2C4IO+UPKXXocH34PZPLsscGN49nsAfsNSh4VOV+8RlMx0Ydh/XmSDcAbT/eO4U/jqE+9Yo8jy7ADSP88a37MhVDg07rq9Mbt3m4MDh+Xmh3EGg+Ud4Dz046cAgelnhEwvQXxzfmcUY5zgwqO38Vt4wwYHhaxqXeH+dAwBkjWJ2AQAAAACAAAAIAQAAEQAAAA==eJxbNBMETtovGuE0AJk7mYg=AQAAAACAAAAIAQAABgEAAA==eJzLm9vWm7Xd2JEBCioyLryYdgXBl+puVtO+a+i45fBRsVft+o6fD4noBr/RctzAyRkpelDVcYbftemSL+Ucjf8tfsDcIeEYJXZeozJTyLH5xtblgf0cjkcY48w3T/jh8GauolWTvJHjgn/791xVMnCcfUZ+JtcLbcf4CzM3HjBUd5w4sSt2OoOSY/7BQLEj68QdOf0XXw5aL+goGilU0rSM3XHBpZVeNzZ8d5B/d76wfAfCfV/vb1D//N7QkftPxrWA4/qOp/7kii9113a0MtrJXBSv5ijy5Wna/0sKjn8euMu9uybhaH5aSmD2JCHHNbqyn0R3cji+PvzfzmrzDwcAVb5s6Q==AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/jdzXB3xDu0/AAAAAAAAAACOLtk8xzHuP1FSb+eQoCofn6bSH8cX6z+60hDRumvpP78G904It+c/+L5bzwHR4z+g8CfLWxm3P9QoH7FLZBQ/NKZ354PHlD3qg1haeb1rOELWulN29sUi/mkAzWYy7D/Rb2cIcX7qPw9z9k2DD+k//k8IIb5q5z/HBlUC8S7jP7aK70g2urI/2mbecM3dBT+PAqOMXM5wPSVoFu+XjIM3SwoR+lmw7T+tULZmx9gKCRqRt9rc0us/1OFt5dA16j+1kolMeMPoP5s6NQLqyuY/ZyGc5/ig2T90fOrVe+l7P770+vnM6YY+aMCeQD2k0TvIGBCUgYkPMSxEgqI=AQAAAACAAAAIAQAADQEAAA==eJwTYfH6L538yp4BClbxSjEurn1jr8J0lAnEf+59YoPhh2f2tbzOU/YZPrH3y+65YTT5of3VsMPdfN037E9r1PF0rGmw3/s/8p3hX0s7FW+OFY1/rKzzTuxZLcfzwiDF4ojoadV5rGdmeXs233hhv7A90jP7/VP7xx9Z1rsseGxfUtQQf1Xpob3j+lZZNobr9gLPPwUVvi61N137mNfukazd8QNNq+SWf7Qy+ufvsOeRuP6uOb8+53W/hrs3fvMBTSfDF/bfzOW+nS1/al+Vvr3snuZj+ycBIiUcNg/s/0s2rFrXdsT+pJBZ34oIDvtPiZd41TIUbFl7XwjOKN5s7r4u/91OEX0lAAZic9A=AQAAAACAAAAIAQAAEQAAAA==eJxjYACBD/YMI5wGAM+kJxA=AQAAAACAAAAIAQAAAwEAAA==eJzLf7WaaeL1U/YMULDfNivK1O8KnC+6LjDDOeqgvdqGS+8ebj5gPy/ML+RE6H5795cxTxQjt9r3/L/1ezuvgb1F+M1Ja9dutF3InSxhtbrd4vbtekY+hmtKMHN+lGWdvjXtoP2D5bYKF+wO2lsVX9rPt/iAfW7JxB2aT/bZ714zXU/r1SZ7rrVHVs1lVbI/05pu7sDTa9vdtp9R684C8+pnPTImbw3ll5+4GFTXdB7uPmGFt2vWTT5oz82lb/Zd96B9hM3lg7sbDtgv25H06dfxvfYH2iK3CT3psL9lXiv8vmCu3XT1IFuFhufW/3ssjbL3Sxk2p7lumbyJixMAuDpv5Q==AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/pm/wlXWF8BAAAAAAAAAAADUIiKLgFjBQOtWfOJqLxIgnlncb1nXvEBjM3x8lxi5QH4NxfaPprVAZa2UwSUprkBakWfwVHuFQLdoasn5Fv4/d7vjycr1gj4UT7abgGpZOe+TGUK7GrQj2PYwC2avv0DZHVEJVnu7QHNNjXcnZrhAbKHUJyYgtUAFFoVZHUesQAMmuSUgpoJAh03bTBiz8T9XNKQzgiZgPtFlsDqVpXI4m3E8peLHwEC7mqtTVC7wCf4+8sdpfL5AqfInN22eukDSMXO4Y4S3QFOnXt0frLNAuOJNkxX+oECe4JT86xRbQLE49IsCenY/3kU7NxNOwjzSYraf6msAMjp3ea4=AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/ptu1/shrO8/AAAAAAAAAABWWfglgOHvPzPBxKFvyz0f0NuITYwh7z+1DREYqCLuP41p5g5Nmew/ejI7/60F5z9Y55iIWCy3P13qgfQ2Aik/KsbmvnsQrz3TwnNyINKEOPfmNNIzeOAiPlpJr62T7z9Vs3kDk9fuP6EQS8Fe4e0/oSk64gNN7D+a7GnogdTlP87vqO6ZwLM/TZun9exTHT/sOqym/XWKPaWF4csJjZ43l+PExTHN7z8EaL3gyoIaCTNe/y12au8/KJG40Tqd7j9WZqweIYXtP7RQyI/gXes/aOSdUI0K2D+V4pFw3u2IPw9CN5N+baI+dTvHV6r97TvqfwNhsucqMYdphLc=AQAAAACAAAAIAQAAPAAAAA==eJxjYACBD/YMKADBN+r72rL6+xv79T6votKfvbKvuNdvfcXmpX3E+fgnsseeo+nDDahlDq3cBwAATzsgAQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/khmHBV3iLc/AAAAAAAAAAAgF20yjOOsP1FSb+eQoCofMHkuUkb+rT9Qn5edyKqnP5B7jcm+XKg/AIMgqOCovz+g8CfLWxm3P9QoH7FLZBQ/NKZ354PHlD3qg1haeb1rOELWulN29sUigIak/paomD+AN7eglP2JP4C0tcx2qHY/gNbf1fDXhj/HBlUC8S7jP7aK70g2urI/2mbecM3dBT+PAqOMXM5wPSVoFu+XjIM3qK13LzB9sj+tULZmx9gKCYDR36PqTKI/YFuNr84Slj/A8JLB+1aOP6CXUkV8Zp8/ZyGc5/ig2T90fOrVe+l7P770+vnM6YY+aMCeQD2k0TvIGBCUgYkPMYgpgow=AQAAAACAAAAIAQAATwAAAA==eJz78vfKjA8GRo4MUPAFje9dGzXnTI2B4/M8kytKrLqOrGp7207P13Cco5jyfdJNZbC6hhY5iPosCQjfXwjK54DSPxyoZQ6t3AcAJ7FIWw==AQAAAACAAAAIAQAACAEAAA==eJxb8F3Se376DwcGKHBYse1/itUvOP/DnQmLXrB+c8hY1iViJfbFoeD5HiPblE8ODo1+UmoV7xwYlk3+mqq204HB+N/iB8zzHBiixM5r/Jzp0OD45tblFZMcGB4xrJ8d2ewgwFV0J1L7l8OJohgtV5ZfDhtMRCf/cPzlIFG29mLuwV8OHT0hb3Z9+uqgcOm4aW36cYcHctO6l7IecUi4u6V2Y/4OB4Y1R1d9U1nusIFpfX9uBcK9C7TStmYYfXdQyFwiWCD4zUFjW1vrGbuvDh9aViT90Pvi4OAYn/eN+6UDg/2Pie3Wux0OnN3EduTYBgeGbO/2mqOLHQ78+rC/32+SAwB4fHbvAQAAAACAAAAYAwAASQIAAA==eJzTXxu/Tn1vi73FBl/JtnPM+xmgoE1c+HSM5Ad7BjQwV8z8ybnsPHvJtAzfN7yccPU7NBa6yF5+b/8o+/GBGJ0cuL4J4ex/dn2ebX9Y/PGe3EJDuPr3FiL3d6xbaz+ZWbJ/YYInXPyKaHeCoctu+7kvlT31XiXAxdcpLs2NLrtkn9HQrWU3swUu/mf5dp+uje/s01QMNtfoL4SLq6ntsdml9sH+ffZHLqOJDXBxwatKvPOlP9inrizasaeiEC6eH/Nq/l+ZD/aTl7zYvvcjwl6FrfXTN8p+sH/x/6rNXDs3uPh2rfIZF9q77eNVVx08cVwfLh62xDQ/Y+Iie/2L+md8PyD8dWv2Hc2cNxvsu70nn3I9FA0Xl9hQWf2oZa+9xITPipPO58DFU6bN6LCqu2I/Madqb/C8Drg4W/Zmlom739l/3NKbHrS1Bh7Ozle2NK+Se2/P5mhi/CZrAVyckcfgsY/Fe/vl91e6/uiaBhc/bnHBZWfge/svNzLufHHvgIsfLYye6OVbaf+jZMq3vq88cHvP/0kJeaD+wb77qcWEOff0EPF43qhsm/xk+253jkdvvyLi92xsk/u/9GX2K953n0z+hQiHpMV7Lk9j22LvpXOv+px3DFx8S0mK9O+bh+xTIuOcznci4sVpq909Ro7H9st/bbEMPjIBLs4i18/uqPrBfur63GXuV8zh7l/0+I3JaZ0P9ukxZz/NNmmCiyd7rD56R+GDvdfzB3P52kvh4p5r+RbKSn2wnznrdVTi1xi4OABYEAi7AQAAAACAAAAYAwAARwIAAA==eJzTNE/62LLovb3f6ig+4Yjc/QxQ0CYufDpG8oM9Axpo6G/K+PHqvb1Z+LW4Ws8uuPodGgtdZC+/t3+U/fhAjE4OXN+U49NNzee+t1+Y++74+r1NcPVPy1y7Dsx8b895+eTcjnSEOT/S1EXVDry3bzuYl535YBJcnE3S/NmBm+/tRU68jd/5ez5cPCAhO2JV70f7Z5p7F2kHLUGIL9QL61D/YB8xeTqj4aQGuLjlOWXe+dIf7CN/FO/YU1EIF8+PeTX/r8wH+8lLXmzf+zEBLq6wtX76RtkP9i/+X7WZa+cGF7d2cZLobnpv7952f9dq8Ylw8cei7S5Mh97bl8U+uRpxaipcPEPLWXAp0F+fuG86JhycCRe/8t/E6/6y9/auN5l+2rbMg4vH/l7LOsPvvb3bYrv89MRFcPF7Hu4ssnZA/x5o5jFenwEP52Sd30198u/tn/TkXa/IXgAXn1do8NjH4r29gdMq1x9d0+Dixy0uuOwMfG//5UbGnS/uHXDxb59ttjy9+t7+rfKmmwrHEfF1/k9KyANgeHY/tZgw554eXDwxWPBax6339pUvTROvPkbE491PdzsPM32wD85+9X3lfUT4SBsm9n3j+GAva6v26Kf0NLg4j37IowcCH+wv/kgu5Xs2By6+9tvtcx6339sv0LtxSHMhInyijze2vUv5YH+HR+vM5pWicPdn7NE+c1rng/2pW3OWzjFpQoSPx+qjdxQ+2Hs9fzCXr70ULu65lm+hrNQH+5mzXkclfo2BiwMAdK0U+Q==AQAAAACAAAAYAwAARwIAAA==eJwr2bEiaIL/e/vekq9fb3Dk7mcgADIj5s3sOfPe/sK7jqJ8gy64ehtL5+bKlzbyTxlelSZ1nJeDid/o2R47c+87+8ch7FEqRk1w9WHGIX+u7Htrf8HpRkvdxQ64uB2L2Mv54W/sg1y8FikfngAX1+9I0Z8875n9A/YVXlLrp8LFJyfVTNBcftD+v3ZGTbloJVz8ll+qv95hTfss6bdHrtzbvw8m/nubgAOb6nrbdfFF81zem+yFiU8oCq4Xft1i8ab02P6HdX+2w8SVvuYoFU14oFRhcX224SKnRXD16Q+uZou8t/9jGKMpcwPhznm7J+99PeOd/dt7TenXfk2Bi7+dxt4jvvitveV1/afWrDPg4hdPHv3l9O21vVP6kS3er2bBxbU3JB6s4n9q7yJw4P+SYoT6HaY35TxY99nXXfF+uCO2AC7+7B/z6Rf1Mva7t1V7tk+5aAcT/9j5KiowvdP2WWtIhcUiC1uYuMdJ6amV3XPN3y8PZGXVczSHiQvZSsl7L35vX9G/20t7FSK+NG9c1fW/JMU5a/v7op5/uR0w8Ud74pW6Hd/b2y5X/aQagkgP9/Y767PPf2d/jXvezWv8E+Hip856OUlNfGsfEZXxyyYcET6vHc7/dcp+bV+RdGbJrGaEf681WT3sXXbT/nO6m9Dand1w8cDVIZX7r/Xbx+kkb76VaQQXP/bghMTt54vscv9mnXkfaQYPh29R6Y8lW99ZN/ZEa++bX24NEzc6f/Gwn6u2YcmPFsUbLosNYOIAhaUKKg==AQAAAACAAABAAQAAFQAAAA==eJzTdZZ5/chMyl53lCaLBgDoYHRBAQAAAACAAABAAQAAEgAAAA==eJxjYACBB/YMozRZNACmTSzZAQAAAACAAABAAQAAFQAAAA==eJzbvff6zfs1t213j9Jk0QDqnOA5AQAAAACAAABAAQAAFQAAAA==eJzTdZZ5/chMyl53lCaLBgDoYHRBAQAAAACAAAAYAwAAiAAAAA==eJx9kcEJgDAUQ93JAdzGnXSTjlDQk6feCgWh4MERpJVSkvCbW98Pofl/mkZ6l+FY5u3d+b4VOcn5uRc+V12GPxj/iUZOEr5W3UZ+Fn4eRY/hd0ZfLxx7sT8Ix16ck4RjL87PwrGXdWfea/NF477MvXDcA/uDcHxzThKOe+D8LBz30PkH5r6CWA==AQAAAACAAADAAwAAlwAAAA==eJxdjTUWw1AQxH6Ymcz2/U+ZIqMi2kZPKmZL+d0yrMr/4W+RflLH59pdhbU6Xon0szrOH3bXYaOO1yL9oo7zh91N2KrjjUi/quP8YXcbdup4K9Jv6jh/2N2FvTreifS7Os4fdvfhoI73Iv2hjvOH3UM4quODSH+q4/xh9xhO6vgo0l/qOH/YnYUfdXwS6Qt1nD9fDwgJHA==AQAAAACAAABAAQAAXQAAAA==eJwtxdEGwgAAAMCIGDFijIgYEWOMMUZExBgjYsQYY4yI2Of30N3LrVd/GwfeOvTOkWPvffDRiU8+O3Xm3IVLV7746pvvrt249cNPd3659+DRk2e//fHXi381vAmdAQAAAACAAAAoAAAADAAAAA==eJxjZSUOAAAQLADJ
 
55
   _AQAAAACAAAAIAQAABgEAAA==eJyrmdvWm7Xd2JEBCkoyLryYdgXB1+xuVtO+a+h46fBRsVft+o6Ch0V0g99oOR7n5IwUPajquNfv2nTJl3KOvv8WP2DukHDMEzuvUZkp5Djnxtblgf0cjncY48w3T/jh8GiuolWTvJFj37/9e64qGTimnpGfyfVC25H3wsyNBwzVHSMndsVOZ1By1DsYKHZknbjjQ7/Fl4PWCzo+jRAqaVrG7lh4aaXXjQ3fHbTenS8s34FwH9ODDeqf3xs66v/JuBZwXN/x259c8aXu2o6BRjuZi+LVHGW/PE37f0nBkfehu9y7axKOjqelBGZPEnLcriv7SXQnh+PXw//trDb/cAAABCVr7g==AQAAAACAAAAIAQAA2wAAAA==eJyzSJMRLe/Ybs8ABh/sBaRyjXoer4HzJ6RubXi84LA9x9Y9u0UCT9lXfBJMej/xgv2dRo/EP7E37C8zzm67IvPWnr180d2Y/+/tj0T+B4L3cP0w2mDDn+knc/bZL3BIumfNdsx+irHaiU8HT9srJbzf26x10f7bp9Dfsotu2p/ie2S2c8Vb+zabUx0rgOYk30M1r2Btub5B7SY4P2K34tSeLQfsN1R4ZO3ROG4vsPXmWblPZ+y3dE39rZ11yd7mvWFPs/5jezbtEA7dE+/tJc5ELMHmPgAzznRYAQAAAACAAAAIAQAA3AAAAA==eJxLz3RXKfRrtGcAgw/2HJM3pLeLZMH5PlfjVfJEZthXN4gp9S9Zap/z7tMVQc2N9rO3rpjD13fI3u9Pum50zSt7nykt+3b+e2/fsfnf////38P1w2iNLRsvZb/tsZ/wP0RQPXCe/b/dDYYPhVfY/59c9/pE1iZ7Oc1kjZaOI/axzFkBnrWv7f90NQsFAs05tuc/inm9gi5v6mNL4fwnKkZdfIUT7e+7qG85dGiB/WUFl1t39VbZW9/1q11suNWe6xR7WyrjdfsnfV7Fuye8tzefudEDm/sAlYhujg==AQAAAACAAAAIAQAAEQAAAA==eJxjYACBD/YMI5wGAM+kJxA=AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/uMula9zY2I/ph2iEvGB1D+X1BOG9mxTP3byfcdafss/d5++wWzzdD9U+1yIcUWEP7sGcHreAJM/p4ItWG3prT8qSHxuSw/DP87JNfh8Q8I/qphfS2lBwj84GFNhCUHCPxXiUZbDQMI/EU1FHrnnZz+XwPD5Bqh6P5j4qtOZ7oY/wsLVcTmOlD8d2DFsSd6wP4YHUVUwi8I/YiBgUyyQwT9rczdXbKvBPx9jQVsFvcE/zqF2lb1mWj/raKTJyvnPP2DW7phgQm8/rDt+TAaEfz/ycaQS/ByLP/b+DrGg0po/+Wuc0HCOwD8PKGxNMYTCP0eNPKYETMI/E4fdk2FAwj8TEFPgaDrCP6DBfSs=AQAAAACAAAAIAQAA5QAAAA==eJxjYGBgsHX7bMuABMxKf9i2S3UkMt14bA/kHgiuY97HyMDQUHrnC1jdD8PrYPqy7ue9QHGNPk5hO8YrMyf8N1TcF/S0vi7JUMAOKN0QovvdFqSvTvvPXpD6j1GfwfQvxsdget67H2Bzevo/gtQxVLCB7GkQmbvqi+1zI6nFZ1SZ9wUq341uOfHKlgkof3sZH8hchhsqn2wbNyx6/s5jzX6QPfO7WPaBxAsTr4DMeXDuOw9I3YH/1hz7gPocnh/6DxT/YaS8l32fyHOlD5IG3/dGOBvI5y3ktgO5b/qLt7YA8IpX+w==AQAAAACAAAAIAQAA8AAAAA==eJxza98SOdUv1IEBBDoOOzy5teOi2AwHB0YoPzclPfvuoUyHtJ7zSqE3KhzsvwQZus1scXjnkLH+Xet0h+vGsT8X/trlUBu6yMHJ/rDDr+VCcr/bDztc/A8EQPo/lC4I2cgepRPt8O/Oia1CTQUOf5kiqrM3Vjvweys1KYi3OXzeeVL1RcQsB9HpTW6zH+52YJyU8yzJ/7DDppi9+/4imbNlO+8G9vNecH7Vzen/piknOfQ3b/3T7l7sIH1CqDRwTr1DbgCXtdjSLgf9B9blDnZrHDh7nmk1ax90MGWQWf21+bCDWpEkWD/M3wB2q4BtAQAAAACAAAAIAQAA9AAAAA==eJxboVHWUOx0wIEBBDoOO1Td/rxQIvOgAyOUf2iKeXTB0r0Oejorz8ks2+XAl/N+zw7LHQ5f0pb08H3b4sCgljWdZ9VehwsLCn4YORx2CL8uLPe7/bDDxf9AAKT/Q2nue8t3cvMdcFi3eq1JefMeh41nnfTtFHc5BL+SsOB6vN1B/afJthM7tzhs3W3pHW+2z+H2D736D/6HHXzm7t33F8mc0Po3azY+OgDn3/z7rL6Nd7/Da6/P255e2u2gEnfkvEbdToekj/uXHGLY7hDz0IVVe942h5Nb7zOtTTzoMFXkTdy35sMOakWSYP0wfwMAZNeJaQ==AQAAAACAAAAIAQAA5AAAAA==eJxriFThavzeYs8ABh/sc9cvY79Vnwfn/2594ld2brb9fJX3dXVX19rXbj7TMd10t/3/WZ2MS75ctJd4/OZdidRb+0fbJ9zJ//fe/gffv////7+H64fR9wpW6k4R6baffyrlxPm2RfbXqpZ/5fy6wX7eZr13D6bvtT+mrqP/J+yKve4jrUNv2t/aHw6akSAKNEdj2n8U87hN+mRl0yvh/EJLEQerjZPsPd6WPQnSWWZvcHYO1/drm+119947UNtxyP4r78bwnb8e21fssG3zmPPe/oyYxs1/QPP+/0c1FwAxF32BAQAAAACAAAAIAQAAEQAAAA==eJxjYACBD/YMI5wGAM+kJxA=AQAAAACAAAAIAQAABgEAAA==eJyrmdvWm7Xd2JEBCkoyLryYdgXB1+xuVtO+a+h46fBRsVft+o6Ch0V0g99oOR7n5IwUPajquNfv2nTJl3KOvv8WP2DukHDMEzuvUZkp5Djnxtblgf0cjncY48w3T/jh8GiuolWTvJFj37/9e64qGTimnpGfyfVC25H3wsyNBwzVHSMndsVOZ1By1DsYKHZknbjjQ7/Fl4PWCzo+jRAqaVrG7lh4aaXXjQ3fHbTenS8s34FwH9ODDeqf3xs66v/JuBZwXN/x259c8aXu2o6BRjuZi+LVHGW/PE37f0nBkfehu9y7axKOjqelBGZPEnLcriv7SXQnh+PXw//trDb/cAAABCVr7g==AQAAAACAAAAIAQAA2wAAAA==eJzzSJMRLe/Ybs8ABh/sFcRzjXoer4HzW1K3NjxecNheYuue3SKBp+ybPgkmvZ94wV6gySPxT+wN+zeMs9uuyLy1Zy9fdDfm/3v7I5H/geA9XD+MFtjwZ/rJnH32exyS7lmzHbM/Yqx24tPB0/YsCe/3NmtdtC/6FPpbdtFN+5V8j8x2rnhr32ZzqmMF0Jzke6jmrVhbrm9QuwnOn7FbcWrPlgP2Gyo8svZoHLfX2XrzrNynM/anuqb+1s66ZO/z3rCnWf+xPbt2CIfuiff2EmcilmBzHwCzWXP1AQAAAACAAAAIAQAARgAAAA==eJxjYACBDTYMyKDhOyqfYQeUfwBKH4HSr6D0SzT1uABM/xkofQNK34PQDQy2EPoDAfPOoMk/sMEuf80GVX4BTnMBcvYUrA==AQAAAACAAAAIAQAABgEAAA==eJzLm9vWm7Xd2JEBCioyLryYdgXBl+puVtO+a+i45fBRsVft+o6fD4noBr/RctzAyRkpelDVcYbftemSL+Ucjf8tfsDcIeEYJXZeozJTyLH5xtblgf0cjkcY48w3T/jh8GauolWTvJHjgn/791xVMnCcfUZ+JtcLbcf4CzM3HjBUd5w4sSt2OoOSY/7BQLEj68QdOf0XXw5aL+goGilU0rSM3XHBpZVeNzZ8d5B/d76wfAfCfV/vb1D//N7QkftPxrWA4/qOp/7kii9113a0MtrJXBSv5ijy5Wna/0sKjn8euMu9uybhaH5aSmD2JCHHNbqyn0R3cji+PvzfzmrzDwcAVb5s6Q==AQAAAACAAAAIAQAAWgAAAA==eJxjYAABHTsGFCCAxtezw06rQ+gGJShfCkqLQGlOKP3DFkKrQPlGUNoaqt8FyreB0A4GdqjqJaDq2KF8NTT3aaGZa4bqPgYhKC0IpVmg9Deou57YAgAf7g+kAQAAAACAAAAIAQAAEQAAAA==eJxbNBMETtovGuE0AJk7mYg=AQAAAACAAAAIAQAABgEAAA==eJyrmdvWm7Xd2JEBCkoyLryYdgXB1+xuVtO+a+h46fBRsVft+o6Ch0V0g99oOR7n5IwUPajquNfv2nTJl3KOvv8WP2DukHDMEzuvUZkp5Djnxtblgf0cjncY48w3T/jh8GiuolWTvJFj37/9e64qGTimnpGfyfVC25H3wsyNBwzVHSMndsVOZ1By1DsYKHZknbjjQ7/Fl4PWCzo+jRAqaVrG7lh4aaXXjQ3fHbTenS8s34FwH9ODDeqf3xs66v/JuBZwXN/x259c8aXu2o6BRjuZi+LVHGW/PE37f0nBkfehu9y7axKOjqelBGZPEnLcriv7SXQnh+PXw//trDb/cAAABCVr7g==AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/jlzXB3xDu0/AAAAAAAAAABfLtk8xzHuP1y7bueQoCofnKbSH8cX6z++0hDRumvpP8QG904It+c/Er9bzwHR4z9n8SfLWxm3P0QqH7FLZBQ/A6h354PHlD0ajFhaeb1rOAHfulN29sUi+mkAzWYy7D/Yb2cIcX7qPxtz9k2DD+k/708IIb5q5z+FBlUC8S7jP66J70g2urI/vGLecM3dBT9y/qKMXM5wPdtEFu+XjIM3UgoR+lmw7T8ubbZmx9gKCSqRt9rc0us/1OFt5dA16j+8kolMeMPoP6Y6NQLqyuY/iCGc5/ig2T/nfOrVe+l7P6/1+vnM6YY+AMGeQD2k0TtGHBCUgYkPMaFjgIA=AQAAAACAAAAIAQAACAEAAA==eJyTYPH6L538yp4BCgJ5pRgX176B8x95n9hg+OGZfQuv85R9hk/sw7J7bhhNfmgvGn64m6/7hv0rzTqejjUN9unMUe8M/1raXffnWNH4x8p6wbU9q+V4XhikBxwRPa06j/XoLG/P5hsv7Ne0R3pmv39q/+0jy3qXBY/t44oa4q8qPbRfvK5Vlo3huv2yJ5+CCl+X2n9Z+JjX7pGsnfKOplVyyz9afVvu77Dnkbj+8Tm/Pud1v4a7r3rzAU0nwxf238zlvp0tf2rfmr697J7mY/vPASIlHDYP7E2kGlataztiP1HYrG9FBIe9d/IlXrUMBduQvheCM4o3m2tuzX+3U0RfCQCrN3H5AQAAAACAAAAIAQAAEQAAAA==eJxjYACBD/YMI5wGAM+kJxA=AQAAAACAAAAIAQAAAwEAAA==eJzjfrOaaeL1U/YMUMBjnxVl6ncFzpdbH5jhHHXQ3mjjpXcPNx+wfxrhF3IidL/94VcxTxQjt9ofY7j9ezuvgT1r5M1Ja9dutG0RSZawWt1uMftJPSMfwzUlmDl2FVmnb007aC+60lbhgt1B+8aSS/v5Fh+wX106cYfmk332u9dM19N6tcn+58ojq+ayKtnvrkg3d+DpteWM3c+odWeB+eOGHhmTt4byC09dDKprOg93n6ri2zXrJh+0l+XWN/uue9A+1Pbywd0NB+wP7kr69Ov4XvtF7ZHbhJ502C+yqBV+XzDXLlgjyFah4bn10X5Lo+z9UoZXmly3TN7ExQkA5l1tdQ==AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/pm/wlXWF8BAAAAAAAAAAACfIiKLgFjBQM/ue+JqLxIgn1ncb1nXvEBoM3x8lxi5QIYNxfaPprVAca2UwSUprkBMkWfwVHuFQAZpasn5Fv4/2bvjycr1gj7fVLabgGpZOZSaGUK7GrQjxfYwC2avv0DeHVEJVnu7QIFNjXcnZrhA8aDUJyYgtUBVFYVZHUesQP0kuSUgpoJAsUnbTBiz8T+QL6QzgiZgPiNDsDqVpXI4nXE8peLHwEDbqqtTVC7wCQ0/8sdpfL5AsvInN22eukCiMXO4Y4S3QFSnXt0frLNAieJNkxX+oEDi35T86xRbQC449IsCenY/SUU7NxNOwjyTY7af6msAMgVZeV0=AQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/plu1/shrO8/AAAAAAAAAABQWfglgOHvP54WxKFvyz0f0duITYwh7z+3DREYqCLuP5Bp5g5Nmew/gTI7/60F5z9G55iIWCy3P57qgfQ2Aik/y8bmvnsQrz2Rx3NyINKEOGjsNNIzeOAiPVpJr62T7z9Vs3kDk9fuP6QQS8Fe4e0/jSk64gNN7D8e7GnogdTlP5fuqO6ZwLM/6JSn9exTHT8WM6ym/XWKPdNM4csJjZ43l+PExTHN7z9ugr3gyoIaCTVe/y12au8/LJG40Tqd7j9ZZqweIYXtP7VQyI/gXes/FeSdUI0K2D/b4ZFw3u2IP6RBN5N+baI+fTrHV6r97TsmgQNhsucqMRJgg8U=AQAAAACAAAAIAQAAPAAAAA==eJxjYACBD/YMKADBN+r72rL6+xv79T6votKfvbKvuNdvfcXmpX3E+fgnsseeo+nDDahlDq3cBwAATzsgAQAAAACAAAAIAQAAEwEAAA==eJwBCAH3/jhmHBV3iLc/AAAAAAAAAAAQGm0yjOOsP1y7bueQoCofYHkuUkb+rT8Qn5edyKqnP0B7jcm+XKg/MIIgqOCovz9n8SfLWxm3P0QqH7FLZBQ/A6h354PHlD0ajFhaeb1rOAHfulN29sUiAIek/paomD/ANbeglP2JP4Cutcx2qHY/QNrf1fDXhj+FBlUC8S7jP66J70g2urI/vGLecM3dBT9y/qKMXM5wPdtEFu+XjIM3cK13LzB9sj8ubbZmx9gKCYDQ36PqTKI/YFuNr84Slj8A75LB+1aOP0CWUkV8Zp8/iCGc5/ig2T/nfOrVe+l7P6/1+vnM6YY+AMGeQD2k0TtGHBCUgYkPMe49fyE=AQAAAACAAAAIAQAATwAAAA==eJz78vfKjA8GRo4MUPAFje9dGzXnTI2B4/M8kytKrLqOrGp7207P13Cco5jyfdJNZbC6hhY5iPosCQjfXwjK54DSPxyoZQ6t3AcAJ7FIWw==AQAAAACAAAAIAQAABgEAAA==eJxr+CHpPT/9hwMDDKzY9j/F6hec/+DuhEUvWL85RCzvErES++KQ8GKPkW3KJ4cNTX5SahXvHBzWT/6aqrbTgcH33+IHzPMcGPLEzmv8nOnA4Pfm1uUVkxwY3jGsnx3Z7FDAWXQnUvuXw4rCGC1XFqD5xqKTfzj+cugoXnsx9+AvhxNdIW92ffrq4HDruGlt+nGHB8rTupeyHnFIeLyldmP+DocHG4+u+qay3CGBeX1/bgXCvQXaaVszjL47OGQtESwQ/ObQsb2t9YzdV4cVrSuSfuh9cbjgGJ/3jfulg4Pzj4nt1rsdHC5sYjtybIPDgVzv9pqjix0c/n7Y3+83yQEA7954jQ==AQAAAACAAAAYAwAASAIAAA==eJzTXRu/Tn1vi/3xlb6SbeeY9zNAwVRx4dMxkh/sGdDAdxnzJ+ey8+znx2b4vuHlhKtfobHQRfbye3uWoscHYnRy4Po6wtn/7Po8255N4vGe3EJDuPozFiL3d6xba7+NTbJ/YYInXPy5aHeCoctu+5WvlD31XiXAxfcoLs2NLrtkz9jQrWU3swUurrRiu0/Xxnf2l5QNNtfoL4SLW6vtsdml9sG+O/Mjl9HEBoT6q0q886U/2J9YUrRjT0UhXLw+5tX8vzIf7A/Ne7F970eEvXpb66dvlP1gz/P9qs1cOze4uL12+YwL7d32VmKrDp44rg8XL1limp8xcZH9lUP6Z3w/IPz1c/YdzZw3G+znWk0+5XooGi5uvqGy+lHLXvsnDZ8VJ53PgYvLTJ/RYVV3xV4+vWpv8LwOuLh49maWibvf2fte700P2loDD2fmK1uaV8m9txdzNTF+k7UALn6O2+Cxj8V7+54nK11/dE2Di8+0uOCyM/C9veG9jDtf3Dvg4ksKoyd6+Vbar02Z8q3vKw/c3pt/UkIeqH+wd7xpMWHOPT24+JTzRmXb5Cfbs9pzPHr7FRG/M2Ob3P+lL7Pf9KT7ZPIvRDi4L95zeRrbFntupXvV57xj4OLXSlKkf988ZL8qNM7pfCciXiq22t1j5HhsL/pji2XwkQlwcSG5fnZH1Q/20bfylrlfMYe7f9XjNyandT7YmySf/TTbpAkunuWx+ugdhQ/2894/mMvXXgoX913Lt1BW6oP9zRWvoxK/xsDFAUL2BSE=AQAAAACAAAAYAwAARgIAAA==eJzjMU/62LLovb3Xwig+4Yjc/QxQMFVc+HSM5Ad7BjQwtb8p48er9/bPPa7F1Xp2wdWv0FjoInv5vT1L0eMDMTo5cH3Vx6ebms99b++R++74+r1NcPX3yly7Dsx8b+9+7eTcjnSEOQbp6qJqB97bsx7Oy858MAkubipp/uzAzff2m46/jd/5ez5cPCUhO2JV70d7K829i7SDlsDFUxfqhXWof7D/MmE6o+GkBri41zll3vnSH+z3fC7esaeiEC5eH/Nq/l+ZD/aH5r3YvvdjAlxcb2v99I2yH+x5vl+1mWvnhlDv4iTR3fTe/n7t/V2rxSfCxa3F2l2YDr233xL45GrEqalw8UdazoJLgf7q/nPDMeHgTLh4zH8Tr/vL3ttrnWD6adsyDy7+5dda1hl+7+1nzbXLT09cBBc/5uHOImv30d6Ov4XHeH0GPJyVdX439cm/B/o373pF9gK4eHahwWMfi/f2ou6rXH90TYOLz7S44LIz8L294b2MO1/cO+Di7F9stjy9+t6+UnLTTYXjiPi6+Scl5AEwPB1vWkyYc08PLh4ULHit49Z7e+G7polXHyPi8eKnu52HmT7Y/0p+9X3lfUT4CBkm9n3j+GCvaar26Kf0NEQ464c8eiDwwd7lU3Ip37M5cPHz326f87j93p5D68YhzYWI8Ek/3tj2LuWDvcxpvTObV4rC3V+wR/vMaZ0P9tMfzlk6x6QJLp7lsfroHYUP9vPeP5jL114KF/ddy7dQVuqD/c0Vr6MSv8bAxQG1GhC+AQAAAACAAAAYAwAASwIAAA==eJxbKXlmR2TZG/sJdQzbb6dk7WeAgqniwqdjJD/YM6ABvwlOeZ+Y3tkLSNXN6t/eAVe/QmOhi+zl9/YsRY8PxOjkwPV9erW3XFrzlb2phqNgya86uPqujasv32N6Yf/o8uPqzVwtcHH3pOyGf4+e2q+YzqPf/aETLv7CmvH54taH9lr+bWIlkRPh4sdYdyyoDXpj33Bi4c3HV+fDxZ2uWL5KV/tg/2bm4ll3JzTAxUsblXjnS3+wZ/Uq2rGnohAuXh/zav5fmQ/2h+a92L73YwJcXG9r/fSNsh/seb5ftZlr5wYXj3nUMvFF/mv714/sW2P29cLFDY81iLp0vLRvFd3XmNGEcOem6h8PLFY/t39pcK+Yx2AyXPztTg/OxcVP7XkS263KV06Biysv/bdv464H9raGP+4tmoZQ//lG2XqfQ2/s72yMNzfKrISHc5jUtbJQuff2fQE3Q/ZnLYCLM7/Wf+xj8d7+3pmVrj+6psHFZ1pccNkZ+N7e8F7GnS/uHXDxB79sj/yLfWu/LSBBSDu4EW7vzT8pIQ/UP9g73rSYMOeeHlxcPnLyQSvH1/YOdcZ9dt/bEOG/09FoUuhL+//WfU0hdj1wcddy28jIouf2P5p/N/OrTICL39Faee7Aoyf261+m/X60YxJcXPeg/+Mz3Q/si5b8ZDI07YeLu+2MWtLK+cF+8dwNMdO6zOHur/4/g+20zgf7JStX7Jtt0gQXz/JYffSOwgf7ee8fzOVrL4WL+67lWygr9cH+5orXUYlfY+DiAD9OEpk=AQAAAACAAAAYAwAARwIAAA==eJwL37EiaIL/e/vFWV+/3uDI3c9AAFREzJvZc+a9feeDjqJ8gy64+sbJTs2VL23k7YJfliZ1nJeDie/v2R47c+87+44Q9igVoya4+kDjkD9X9r21v+N6o6XuYgdcvJhF7OX88Df2G1y9FikfngAXD+9I0Z8875n9RvYVXlLrp8LFpyTVTNBcftDeUTujply0Ei7u5Z/qr3dY095U4u2RK/f274OJ794u4MCmut6WNbponst7k70w8Zzy4Hrh1y0WXMXH9j+s+7MdJr72V45S0YQHSqbm12cbLnJaBBO/lP7garbIe/vNmjGaMjcQ7vyye/Le1zPe2c+93JR+7dcUuHjydPYe8cVv7R+e0H9qzToDLh568ugvp2+v7XkjjmzxfjULLp62PvFgFf9T+6XsB/4vKUaoX2F8U86DdZ/9uvPeD3fEFsDFXX4wn35RL2N/cXu1Z/uUi3Ywca2mV1GB6Z22XC0hFRaLLGxh4mKTpKdWds81X9EdyMqq52gOE1exlZL3XvzeXqR1t5f2KkR8ffx0Vdf/khRn38n3RT3/cjtg4lf3xCt1O763l56t+kk1BJEeTu531mef/87+Csu8m9f4J8LFd5/1cpKa+NaeIzjjl004InzUHc//dcp+bX8r5sySWc0I/65rsnrYu+ym/dNkN6G1O7vh4jWrQir3X+u3n6KevPlWphFcPOfBCYnbzxfZHWPNPvM+0gweDsJR6Y8lW99Z50yN1t43v9waJl5/4eJhP1dtw0lMrYo3XBYbwMQBJ1kBDg==AQAAAACAAABAAQAAFQAAAA==eJzTdZZ5/chMyl53lCaLBgDoYHRBAQAAAACAAABAAQAAEgAAAA==eJxjYACBB/YMozRZNACmTSzZAQAAAACAAABAAQAAFQAAAA==eJzbvff6zfs1t213j9Jk0QDqnOA5AQAAAACAAABAAQAAFQAAAA==eJzTdZZ5/chMyl53lCaLBgDoYHRBAQAAAACAAAAYAwAAiAAAAA==eJx9kcEJgDAUQ93JAdzGnXSTjlDQk6feCgWh4MERpJVSkvCbW98Pofl/mkZ6l+FY5u3d+b4VOcn5uRc+V12GPxj/iUZOEr5W3UZ+Fn4eRY/hd0ZfLxx7sT8Ix16ck4RjL87PwrGXdWfea/NF477MvXDcA/uDcHxzThKOe+D8LBz30PkH5r6CWA==AQAAAACAAADAAwAAlwAAAA==eJxdjTUWw1AQxH6Ymcz2/U+ZIqMi2kZPKmZL+d0yrMr/4W+RflLH59pdhbU6Xon0szrOH3bXYaOO1yL9oo7zh91N2KrjjUi/quP8YXcbdup4K9Jv6jh/2N2FvTreifS7Os4fdvfhoI73Iv2hjvOH3UM4quODSH+q4/xh9xhO6vgo0l/qOH/YnYUfdXwS6Qt1nD9fDwgJHA==AQAAAACAAABAAQAAXQAAAA==eJwtxdEGwgAAAMCIGDFijIgYEWOMMUZExBgjYsQYY4yI2Of30N3LrVd/GwfeOvTOkWPvffDRiU8+O3Xm3IVLV7746pvvrt249cNPd3659+DRk2e//fHXi381vAmdAQAAAACAAAAoAAAADAAAAA==eJxjZSUOAAAQLADJ
55
56
  </AppendedData>
56
57
</VTKFile>