; Source: region1 ;Units: SI: meter-kilogram-second config tpflow grid 3655,84 gen 0.0,0.0 0.0,402.5 7254.6147,472.4137 7254.6147,0.0 i=1,265 j=1,22 gen 7254.6147,0.0 7254.6147,472.4137 9055.931,598.16144 9055.931,0.0 i=265,331 j=1,22 gen 9055.931,0.0 9055.931,598.16144 9092.72,598.16943 9092.72,0.0 i=332,335 j=1,22 gen 9092.72,0.0 9092.72,598.16943 9373.119,598.22925 9373.119,0.0 i=336,341 j=1,22 gen 9373.119,0.0 9373.119,598.22925 9748.182,598.30927 9748.182,0.0 i=341,347 j=1,22 gen 9748.182,0.0 9748.182,598.30927 10163.381,598.39795 10163.381,0.0 i=347,353 j=1,22 gen 10163.381,0.0 10163.381,598.39795 10484.994,598.4667 10484.994,0.0 i=353,359 j=1,22 gen 10484.994,0.0 10484.994,598.4667 10884.735,598.55194 10884.735,0.0 i=359,374 j=1,22 gen 10884.735,0.0 10884.735,598.55194 11108.172,598.5996 11108.172,0.0 i=374,381 j=1,22 gen 11108.172,0.0 11108.172,598.5996 11254.615,598.6309 11254.615,0.0 i=381,386 j=1,22 gen 11254.615,0.0 11254.615,598.6309 16668.37,599.78674 16668.367,0.0 i=386,583 j=1,22 gen 16668.367,0.0 16668.37,599.78674 22000.0,600.9251 22000.0,0.0 i=583,777 j=1,22 gen 0.0,402.5 0.0,893.0 7254.6147,893.0 7254.6147,472.4137 i=778,1042 j=1,5 gen 7254.6147,472.4137 7254.6147,893.0 9055.931,893.0 9055.931,598.16144 i=1042,1108 j=1,5 gen 9055.931,598.16144 9055.931,893.0 9276.447,893.5458 9311.266,855.53204 i=1108,1117 j=1,5 gen 9055.931,598.16144 9311.266,855.53204 9331.556,833.37946 9092.72,598.16943 i=1122,1125 j=1,10 gen 9092.72,598.16943 9331.556,833.37946 9563.43,835.3389 9373.119,598.22925 i=336,341 j=22,31 gen 9373.119,598.22925 9563.43,835.3389 9874.733,837.9697 9748.182,598.30927 i=341,347 j=22,31 gen 9748.182,598.30927 9874.733,837.9697 10206.734,840.7753 10163.381,598.39795 i=347,353 j=22,31 gen 10163.381,598.39795 10206.734,840.7753 10484.994,843.12683 10484.994,598.4667 i=353,359 j=22,31 gen 10484.994,598.4667 10484.994,843.12683 10884.735,846.50494 10884.735,598.55194 i=359,374 j=22,31 gen 10884.735,598.55194 10884.735,846.50494 11108.172,848.3932 11108.172,598.5996 i=374,381 j=22,31 gen 11108.172,598.5996 11108.172,848.3932 11254.615,849.63055 11254.615,598.6309 i=381,386 j=22,31 gen 11254.615,598.6309 11254.615,849.63055 16668.371,864.0277 16668.37,599.78674 i=386,583 j=22,31 gen 16668.37,599.78674 16668.371,864.0277 22000.0,878.2064 22000.0,600.9251 i=583,777 j=22,31 gen 9331.556,833.37946 9311.266,855.53204 9681.219,1231.71 9702.608,1210.6748 i=1126,1141 j=1,4 gen 9331.556,833.37946 9702.608,1210.6748 9864.686,1210.6746 9563.43,835.3389 i=336,341 j=31,46 gen 9563.43,835.3389 9864.686,1210.6746 10071.54,1210.6746 9874.733,837.9697 i=341,347 j=31,46 gen 9874.733,837.9697 10071.54,1210.6746 10272.896,1210.6746 10206.734,840.7753 i=347,353 j=31,46 gen 10206.734,840.7753 10272.896,1210.6746 10484.994,1210.6746 10484.994,843.12683 i=353,359 j=31,46 gen 10484.994,843.12683 10484.994,1210.6746 10884.735,1210.6746 10884.735,846.50494 i=359,374 j=31,46 gen 9311.266,855.53204 9311.266,1084.8022 9681.807,1328.5612 9681.219,1231.71 i=1126,1141 j=4,10 gen 9276.447,893.5458 9276.447,1061.8972 9311.266,1084.8022 9311.266,855.53204 i=1117,1121 j=5,11 gen 10884.735,846.50494 10884.735,1210.6746 11108.172,1210.6746 11108.172,848.3932 i=374,381 j=31,46 gen 11108.172,848.3932 11108.172,1210.6746 11254.615,1210.6746 11254.615,849.63055 i=381,386 j=31,46 gen 11254.615,849.63055 11254.615,1210.6746 16668.371,1226.2603 16668.371,864.0277 i=386,583 j=31,46 gen 16668.371,864.0277 16668.371,1226.2603 22000.0,1241.6094 22000.0,878.2064 i=583,777 j=31,46 gen 0.0,893.0 0.0,1031.5 7254.6147,1043.8618 7254.6147,893.0 i=1142,1406 j=1,7 gen 7254.6147,893.0 7254.6147,1043.8618 9055.931,1059.9302 9055.931,893.0 i=1406,1472 j=1,7 gen 9055.931,893.0 9055.931,1059.9302 9276.447,1061.8972 9276.447,893.5458 i=1108,1117 j=5,11 gen 0.0,1031.5 0.0,1234.0 7254.6147,1235.0245 7254.6147,1043.8618 i=1473,1737 j=1,8 gen 7254.6147,1043.8618 7254.6147,1235.0245 9055.931,1230.3433 9055.931,1059.9302 i=1406,1472 j=7,14 gen 9055.931,1059.9302 9055.931,1230.3433 9276.447,1230.0217 9276.447,1061.8972 i=1108,1117 j=11,18 gen 9276.447,1061.8972 9276.447,1230.0217 9311.266,1251.7979 9311.266,1084.8022 i=1117,1121 j=11,18 gen 9311.266,1084.8022 9311.266,1251.7979 9682.752,1484.1327 9681.807,1328.5612 i=1126,1141 j=10,17 gen 9702.608,1210.6748 9865.459,1376.2646 9993.941,1371.7136 9864.686,1210.6746 i=336,341 j=46,53 gen 9681.219,1231.71 9824.798,1377.7048 9865.459,1376.2646 9702.608,1210.6748 i=1738,1741 j=1,8 gen 9864.686,1210.6746 9993.941,1371.7136 10153.59,1366.0588 10071.54,1210.6746 i=341,347 j=46,53 gen 10071.54,1210.6746 10153.59,1366.0588 10299.764,1360.8813 10272.896,1210.6746 i=347,353 j=46,53 gen 10272.896,1210.6746 10299.764,1360.8813 10484.994,1360.0 10484.994,1210.6746 i=353,359 j=46,53 gen 10484.994,1210.6746 10484.994,1360.0 10884.735,1340.0 10884.735,1210.6746 i=359,374 j=46,53 gen 10884.735,1210.6746 10884.735,1340.0 11108.172,1340.0 11108.172,1210.6746 i=374,381 j=46,53 gen 11108.172,1210.6746 11108.172,1340.0 11254.615,1340.6123 11254.615,1210.6746 i=381,386 j=46,53 gen 11254.615,1210.6746 11254.615,1340.6123 16668.371,1457.5435 16668.371,1226.2603 i=386,583 j=46,53 gen 0.0,1234.0 0.0,1430.0 7255.112,1430.0 7254.6147,1235.0245 i=1742,2006 j=1,9 gen 7254.6147,1235.0245 7255.112,1430.0 9055.931,1430.0 9055.931,1230.3433 i=2006,2072 j=1,9 gen 9055.931,1230.3433 9055.931,1430.0 9276.447,1463.552 9276.447,1230.0217 i=2073,2082 j=1,9 gen 9276.447,1230.0217 9276.447,1463.552 9311.266,1475.6217 9311.266,1251.7979 i=2082,2086 j=1,9 gen 9681.219,1231.71 9681.807,1328.5612 9760.0,1380.0 9824.798,1377.7048 i=2087,2094 j=1,7 gen 16668.371,1226.2603 16668.371,1457.5435 22000.0,1572.7007 22000.0,1241.6094 i=583,777 j=46,53 gen 9311.266,1251.7979 9311.266,1475.6217 9683.484,1604.6503 9682.752,1484.1327 i=1126,1141 j=17,25 gen 9681.807,1328.5612 9682.752,1484.1327 9837.205,1533.5361 9760.0,1380.0 i=2095,2102 j=1,8 gen 10884.735,1340.0 10884.735,1532.8855 11108.172,1565.2068 11108.172,1340.0 i=374,381 j=53,60 gen 10484.994,1360.0 10484.994,1499.9999 10884.735,1532.8855 10884.735,1340.0 i=359,374 j=53,60 gen 11108.172,1340.0 11108.172,1565.2068 11254.615,1583.7418 11254.615,1340.6123 i=381,386 j=53,60 gen 11254.615,1340.6123 11254.615,1583.7418 16668.37,1731.465 16668.371,1457.5435 i=386,583 j=53,60 gen 9760.0,1380.0 9837.205,1533.5361 9973.167,1528.5701 9824.798,1377.7048 i=2103,2109 j=1,8 gen 9824.798,1377.7048 9973.167,1528.5701 10013.785,1527.0867 9865.459,1376.2646 i=2110,2113 j=1,8 gen 9993.941,1371.7136 10111.795,1518.5466 10228.731,1508.3574 10153.59,1366.0588 i=341,347 j=53,60 gen 9865.459,1376.2646 10013.785,1527.0867 10111.795,1518.5466 9993.941,1371.7136 i=2113,2118 j=1,8 gen 10153.59,1366.0588 10228.731,1508.3574 10324.647,1500.0 10299.764,1360.8813 i=347,353 j=53,60 gen 10299.764,1360.8813 10324.647,1500.0 10484.994,1499.9999 10484.994,1360.0 i=353,359 j=53,60 gen 0.0,1430.0 0.0,1580.0 7255.1123,1580.0 7255.112,1430.0 i=1742,2006 j=9,13 gen 7255.112,1430.0 7255.1123,1580.0 9055.931,1580.0 9055.931,1430.0 i=2006,2072 j=9,13 gen 9055.931,1430.0 9055.931,1580.0 9276.447,1580.0 9276.447,1463.552 i=2073,2082 j=9,13 gen 9276.447,1463.552 9276.447,1580.0 9311.266,1589.6936 9311.266,1475.6217 i=2082,2086 j=9,13 gen 16668.371,1457.5435 16668.37,1731.465 22000.0,1700.7007 22000.0,1572.7007 i=583,777 j=53,60 gen 9311.266,1475.6217 9311.266,1589.6936 9684.024,1693.471 9683.484,1604.6503 i=1126,1141 j=25,29 gen 9682.752,1484.1327 9683.484,1604.6503 9912.966,1684.1991 9837.205,1533.5361 i=2095,2102 j=8,16 gen 10324.647,1500.0 10404.821,1706.3105 10484.994,1722.6211 10484.994,1499.9999 i=353,359 j=60,68 gen 10228.731,1508.3574 10324.647,1690.0 10404.821,1706.3105 10324.647,1500.0 i=2125,2131 j=1,9 gen 10111.795,1518.5466 10249.408,1690.0 10324.647,1690.0 10228.731,1508.3574 i=2119,2125 j=1,9 gen 10484.994,1499.9999 10484.994,1722.6211 10884.735,1803.9453 10884.735,1532.8855 i=359,374 j=60,68 gen 9973.167,1528.5701 10131.926,1690.0 10174.002,1690.0 10013.785,1527.0867 i=2132,2135 j=1,9 gen 9837.205,1533.5361 9912.966,1684.1991 10131.926,1690.0 9973.167,1528.5701 i=2103,2109 j=8,16 gen 10013.785,1527.0867 10174.002,1690.0 10249.408,1690.0 10111.795,1518.5466 i=2135,2140 j=1,9 gen 10884.735,1532.8855 10884.735,1803.9453 11108.172,1849.4019 11108.172,1565.2068 i=374,381 j=60,68 gen 11108.172,1565.2068 11108.172,1849.4019 11254.615,1879.1945 11254.615,1583.7418 i=381,386 j=60,68 gen 0.0,1580.0 0.0,1690.0 7255.112,1690.0 7255.1123,1580.0 i=2141,2405 j=1,4 gen 7255.1123,1580.0 7255.112,1690.0 9055.931,1802.9473 9055.931,1580.0 i=2405,2471 j=1,4 gen 9055.931,1580.0 9055.931,1802.9473 9276.447,1815.4193 9276.447,1580.0 i=2472,2481 j=1,4 gen 9276.447,1580.0 9276.447,1815.4193 9311.266,1817.1143 9311.266,1589.6936 i=2481,2485 j=1,4 gen 9311.266,1589.6936 9311.266,1817.1143 9684.886,1835.3027 9684.024,1693.471 i=1126,1141 j=29,32 gen 11254.615,1583.7418 11254.615,1879.1945 16668.37,1879.1945 16668.37,1731.465 i=386,583 j=60,68 gen 9683.484,1604.6503 9684.024,1693.471 9912.966,1757.2096 9912.966,1684.1991 i=2486,2493 j=1,5 gen 9912.966,1684.1991 9912.966,1757.2096 10131.926,1818.1688 10131.926,1690.0 i=2103,2109 j=16,20 gen 0.0,1690.0 0.0,1884.0 7255.1123,1930.7445 7255.112,1690.0 i=2141,2405 j=4,8 gen 7255.112,1690.0 7255.1123,1930.7445 9055.931,1942.347 9055.931,1802.9473 i=2405,2471 j=4,8 gen 9684.024,1693.471 9684.886,1835.3027 9912.966,1846.406 9912.966,1757.2096 i=2486,2493 j=5,8 gen 10131.926,1690.0 10131.926,1818.1688 10174.168,1820.6124 10174.002,1690.0 i=2132,2135 j=9,13 gen 10174.002,1690.0 10174.168,1820.6124 10249.408,1824.9651 10249.408,1690.0 i=2135,2140 j=9,13 gen 10249.408,1690.0 10249.408,1824.9651 10324.647,1829.3176 10324.647,1690.0 i=2119,2125 j=9,13 gen 10324.647,1690.0 10324.647,1829.3176 10404.821,1833.9554 10404.821,1706.3105 i=2125,2131 j=9,13 gen 10404.821,1706.3105 10404.821,1833.9554 10484.994,1838.5933 10484.994,1722.6211 i=2494,2500 j=1,5 gen 16668.37,1731.465 16668.37,1879.1945 22000.0,1879.1945 22000.0,1700.7007 i=583,777 j=60,68 gen 10484.994,1722.6211 10484.994,1838.5933 10884.735,1861.7175 10884.735,1803.9453 i=2501,2516 j=1,5 gen 9912.966,1757.2096 9912.966,1846.406 10131.926,1857.0654 10131.926,1818.1688 i=2103,2109 j=20,23 gen 9055.931,1802.9473 9055.931,1942.347 9276.447,1943.7678 9276.447,1815.4193 i=2472,2481 j=4,8 gen 9276.447,1815.4193 9276.447,1943.7678 9311.266,1943.9478 9311.266,1817.1143 i=2481,2485 j=4,8 gen 9311.266,1817.1143 9311.266,1943.9478 9685.558,1945.8816 9684.886,1835.3027 i=1126,1141 j=32,36 gen 10131.926,1818.1688 10131.926,1857.0654 10174.168,1859.1217 10174.168,1820.6124 i=2517,2520 j=1,4 gen 10174.168,1820.6124 10174.168,1859.1217 10249.408,1862.7845 10249.408,1824.9651 i=2520,2525 j=1,4 gen 10884.735,1803.9453 10884.735,1861.7175 11036.958,1881.6385 11108.172,1849.4019 i=2526,2533 j=1,5 gen 9684.886,1835.3027 9685.558,1945.8816 9912.966,1947.8689 9912.966,1846.406 i=2534,2541 j=1,5 gen 10324.647,1829.3176 10324.647,1866.4474 10404.821,1870.3503 10404.821,1833.9554 i=2548,2554 j=1,4 gen 10404.821,1833.9554 10404.821,1870.3503 10484.994,1874.2534 10484.994,1838.5933 i=2494,2500 j=5,8 gen 10249.408,1824.9651 10249.408,1862.7845 10324.647,1866.4474 10324.647,1829.3176 i=2542,2548 j=1,4 gen 10484.994,1838.5933 10484.994,1874.2534 10884.735,1893.7133 10884.735,1861.7175 i=2501,2516 j=5,8 gen 9912.966,1846.406 9912.966,1947.8689 10131.926,1949.2795 10131.926,1857.0654 i=2103,2109 j=23,27 gen 10131.926,1857.0654 10131.926,1949.2795 10174.168,1949.5518 10174.168,1859.1217 i=2517,2520 j=4,8 gen 10174.168,1859.1217 10174.168,1949.5518 10249.408,1950.0365 10249.408,1862.7845 i=2520,2525 j=4,8 gen 10249.408,1862.7845 10249.408,1950.0365 10324.647,1950.5212 10324.647,1866.4474 i=2542,2548 j=4,8 gen 10324.647,1866.4474 10324.647,1950.5212 10404.821,1951.0378 10404.821,1870.3503 i=2548,2554 j=4,8 gen 10884.735,1861.7175 10884.735,1893.7133 11039.38,1901.2416 11036.958,1881.6385 i=2526,2533 j=5,8 gen 11108.172,1849.4019 11108.172,1904.5906 11254.615,1911.7195 11254.615,1879.1945 i=381,386 j=68,71 gen 11036.958,1881.6385 11039.38,1901.2416 11108.172,1904.5906 11108.172,1849.4019 i=2555,2559 j=1,4 gen 0.0,1884.0 0.0,2001.0 7255.112,2001.0 7255.1123,1930.7445 i=2560,2824 j=1,4 gen 10404.821,1870.3503 10404.821,1951.0378 10484.994,1951.5544 10484.994,1874.2534 i=2891,2897 j=1,5 gen 10484.994,1874.2534 10484.994,1951.5544 10884.734,1954.1299 10884.735,1893.7133 i=2898,2913 j=1,5 gen 11254.615,1879.1945 11254.615,1911.7195 16668.37,1911.7195 16668.37,1879.1945 i=386,583 j=68,71 gen 16668.37,1879.1945 16668.37,1911.7195 22000.0,1911.7196 22000.0,1879.1945 i=583,777 j=68,71 gen 10884.735,1893.7133 10884.734,1954.1299 11039.38,1955.1263 11039.38,1901.2416 i=2526,2533 j=8,12 gen 11039.38,1901.2416 11039.38,1955.1263 11108.172,1955.5695 11108.172,1904.5906 i=2555,2559 j=4,8 gen 11108.172,1904.5906 11108.172,1955.5695 11254.615,1956.5131 11254.615,1911.7195 i=381,386 j=71,75 gen 7255.1123,1930.7445 7255.112,2001.0 9055.931,2001.0 9055.931,1942.347 i=2824,2890 j=1,4 gen 11254.615,1911.7195 11254.615,1956.5131 16668.37,1956.5131 16668.37,1911.7195 i=386,583 j=71,75 gen 16668.37,1911.7195 16668.37,1956.5131 22000.0,1956.5131 22000.0,1911.7196 i=2914,3108 j=1,5 gen 9055.931,1942.347 9055.931,2001.0 9276.447,2001.0 9276.447,1943.7678 i=3109,3118 j=1,4 gen 9276.447,1943.7678 9276.447,2001.0 9311.266,2001.0 9311.266,1943.9478 i=3118,3122 j=1,4 gen 9311.266,1943.9478 9311.266,2001.0 9685.893,2001.0 9685.558,1945.8816 i=1126,1141 j=36,39 gen 9685.558,1945.8816 9685.893,2001.0 9912.966,2001.0 9912.966,1947.8689 i=2534,2541 j=5,8 gen 9912.966,1947.8689 9912.966,2001.0 10131.926,2001.0 10131.926,1949.2795 i=2103,2109 j=27,30 gen 10174.168,1949.5518 10174.168,2001.0 10249.408,2001.0 10249.408,1950.0365 i=3126,3131 j=1,4 gen 10131.926,1949.2795 10131.926,2001.0 10174.168,2001.0 10174.168,1949.5518 i=3123,3126 j=1,4 gen 10249.408,1950.0365 10249.408,2001.0 10324.646,2001.0 10324.647,1950.5212 i=3132,3138 j=1,4 gen 10324.647,1950.5212 10324.646,2001.0 10404.821,2001.0 10404.821,1951.0378 i=3138,3144 j=1,4 gen 10404.821,1951.0378 10404.821,2001.0 10484.994,2001.0 10484.994,1951.5544 i=2891,2897 j=5,8 gen 10484.994,1951.5544 10484.994,2001.0 10884.735,2001.0 10884.734,1954.1299 i=2898,2913 j=5,8 gen 10884.734,1954.1299 10884.735,2001.0 11039.38,2001.0 11039.38,1955.1263 i=2526,2533 j=12,15 gen 11039.38,1955.1263 11039.38,2001.0 11108.172,2001.0 11108.172,1955.5695 i=2555,2559 j=8,11 gen 11108.172,1955.5695 11108.172,2001.0 11254.615,2001.0 11254.615,1956.5131 i=381,386 j=75,78 gen 11254.615,1956.5131 11254.615,2001.0 16668.367,2001.0 16668.37,1956.5131 i=386,583 j=75,78 gen 16668.37,1956.5131 16668.367,2001.0 22000.0,2001.0 22000.0,1956.5131 i=2914,3108 j=5,8 gen 0.0,2001.0 0.0,2200.0 7255.1113,2200.0 7255.112,2001.0 i=3145,3409 j=1,8 gen 7255.112,2001.0 7255.1113,2200.0 9055.931,2200.0 9055.931,2001.0 i=2824,2890 j=4,11 gen 9055.931,2001.0 9055.931,2200.0 9276.447,2200.0 9276.447,2001.0 i=3410,3419 j=1,8 gen 9276.447,2001.0 9276.447,2200.0 9311.266,2200.0 9311.266,2001.0 i=3118,3122 j=4,11 gen 9311.266,2001.0 9311.266,2200.0 9687.102,2200.0 9685.893,2001.0 i=1126,1141 j=39,46 gen 9685.893,2001.0 9687.102,2200.0 9912.967,2200.0 9912.966,2001.0 i=3420,3427 j=1,8 gen 9912.966,2001.0 9912.967,2200.0 10131.926,2200.0 10131.926,2001.0 i=2103,2109 j=30,37 gen 10131.926,2001.0 10131.926,2200.0 10174.168,2200.0 10174.168,2001.0 i=3428,3431 j=1,8 gen 10174.168,2001.0 10174.168,2200.0 10249.408,2200.0 10249.408,2001.0 i=3126,3131 j=4,11 gen 10404.821,2001.0 10404.821,2200.0 10484.994,2200.0 10484.994,2001.0 i=3432,3438 j=1,8 gen 10324.646,2001.0 10324.647,2200.0 10404.821,2200.0 10404.821,2001.0 i=3138,3144 j=4,11 gen 10249.408,2001.0 10249.408,2200.0 10324.647,2200.0 10324.646,2001.0 i=3439,3445 j=1,8 gen 10484.994,2001.0 10484.994,2200.0 10884.735,2200.0 10884.735,2001.0 i=3446,3461 j=1,8 gen 10884.735,2001.0 10884.735,2200.0 11039.379,2200.0 11039.38,2001.0 i=2526,2533 j=15,22 gen 11039.38,2001.0 11039.379,2200.0 11108.171,2200.0 11108.172,2001.0 i=2555,2559 j=11,18 gen 11108.172,2001.0 11108.171,2200.0 11254.615,2200.0 11254.615,2001.0 i=381,386 j=78,85 gen 11254.615,2001.0 11254.615,2200.0 16668.37,2200.0 16668.367,2001.0 i=386,583 j=78,85 gen 16668.367,2001.0 16668.37,2200.0 22000.0,2200.0 22000.0,2001.0 i=3462,3656 j=1,8 ;Define material property ;;;;;;NORTLAND ;;;;;; def prop_NORTLAND e_NORTLAND = 400e6 ;E mod p_NORTLAND = 0.45 ;porsion bul_NORTLAND = e_NORTLAND/(3*(1-2*p_NORTLAND)) ;bulk mod she_NORTLAND = e_NORTLAND/(2*(1+p_NORTLAND)) ;shear mod coh_NORTLAND = 21e3 ;cohesion fric_NORTLAND = 24.51 ;firction dens_NORTLAND = 2138.3 ;density k_NORTLAND=1.02e-21 ;perm coef 1.02e-11 eo_NORTLAND=0.01 ;porosity end prop_NORTLAND ; ;;;;;;SELE ;;;;;; def prop_SELE e_SELE = 1230e6 ;E mod p_SELE = 0.39 ;porsion bul_SELE = e_SELE/(3*(1-2*p_SELE)) ;bulk mod she_SELE = e_SELE/(2*(1+p_SELE)) ;shear mod coh_SELE = 21e3 ;cohesion fric_SELE = 26.96 ;firction dens_SELE = 2185 ;density k_SELE=1.02e-21 ;perm coef eo_SELE=0.01 ;porosity end prop_SELE ;;;;;;LISTA ;;;;;; def prop_LISTA e_LISTA = 2570e6 ;E mod p_LISTA = 0.35 ;porsion bul_LISTA = e_LISTA/(3*(1-2*p_LISTA)) ;bulk mod she_LISTA = e_LISTA/(2*(1+p_LISTA)) ;shear mod coh_LISTA = 21e3 ;cohesion fric_LISTA = 26.19 ;firction dens_LISTA = 2305.3 ;density k_LISTA=1.02e-21 ;perm coef eo_LISTA=0.01 ;porosity end prop_LISTA ;;;;;;SHETLAND ;;;;;; def prop_SHETLAND e_SHETLAND =1440e6 ;E mod p_SHETLAND = 0.39 ;porsion bul_SHETLAND = e_SHETLAND/(3*(1-2*p_SHETLAND)) ;bulk mod she_SHETLAND = e_SHETLAND/(2*(1+p_SHETLAND)) ;shear mod coh_SHETLAND = 21e3 ;cohesion fric_SHETLAND = 26.19 ;firction dens_SHETLAND = 2239.2 ;density k_SHETLAND=1.02e-21 ;perm coef eo_SHETLAND=0.01 ;porosity end prop_SHETLAND ;;;;;;SVARTE ;;;;;; def prop_SVARTE e_SVARTE = 1810e6 ;E mod p_SVARTE = 0.37 ;porsion bul_SVARTE = e_SVARTE/(3*(1-2*p_SVARTE)) ;bulk mod she_SVARTE = e_SVARTE/(2*(1+p_SVARTE)) ;shear mod coh_SVARTE = 21e3 ;cohesion fric_SVARTE = 26.19 ;firction dens_SVARTE = 2270 ;density k_SVARTE=1.02e-21 ;perm coef eo_SVARTE=0.08 ;porosity 0.0835 end prop_SVARTE ;;;;;;ASGARD ;;;;;; def prop_ASGARD e_ASGARD = 1440e6 ;E mod p_ASGARD = 0.39 ;porsion bul_ASGARD = e_ASGARD/(3*(1-2*p_ASGARD)) ;bulk mod she_ASGARD = e_ASGARD/(2*(1+p_ASGARD)) ;shear mod coh_ASGARD = 21e3 ;cohesion fric_ASGARD = 25.2 ;firction dens_ASGARD = 2254.7 ;density k_ASGARD=1.02e-21 ;perm coef eo_ASGARD=0.01 ;porosity end prop_ASGARD ;;;;;;CAPROCK ;;;;;; def prop_CAPROCK e_CAPROCK = 2150e6 ;E mod p_CAPROCK = 0.37 ;porsion bul_CAPROCK = e_CAPROCK/(3*(1-2*p_CAPROCK)) ;bulk mod she_CAPROCK = e_CAPROCK/(2*(1+p_CAPROCK)) ;shear mod coh_CAPROCK = 21e3 ;cohesion fric_CAPROCK =25.43 ;firction dens_CAPROCK = 2314.9 ;density k_CAPROCK=1.02e-21 ;perm coef eo_CAPROCK=0.01 ;porosity end prop_CAPROCK ;;;;;;SOGNEFJORD ;;;;;; def prop_SOGNEFJORD e_SOGNEFJORD = 2790e6 ;E mod p_SOGNEFJORD = 0.335 ;porsion bul_SOGNEFJORD = e_SOGNEFJORD/(3*(1-2*p_SOGNEFJORD)) ;bulk mod she_SOGNEFJORD = e_SOGNEFJORD/(2*(1+p_SOGNEFJORD)) ;shear mod coh_SOGNEFJORD = 21e3 ;cohesion fric_SOGNEFJORD = 27.97 ;firction dens_SOGNEFJORD = 2240 ;density k_SOGNEFJORD=4.7e-13 ;perm coef ; k_SOGNEFJORD=(7e-12)*(1.02e-6)*e-12 eo_SOGNEFJORD=0.16 ;porosity end prop_SOGNEFJORD ;;;;;;FENSSCROSS ;;;;;; def prop_FENSSCROSS e_FENSSCROSS = 3700e6 ;E mod p_FENSSCROSS = 0.315 ;porsion bul_FENSSCROSS = e_FENSSCROSS/(3*(1-2*p_FENSSCROSS)) ;bulk mod she_FENSSCROSS = e_FENSSCROSS/(2*(1+p_FENSSCROSS)) ;shear mod coh_FENSSCROSS = 21e3 ;cohesion fric_FENSSCROSS = 27.83 ;firction dens_FENSSCROSS = 2310 ;density k_FENSSCROSS=2.3e-13 ;perm coef ; k_FENSSCROSS=(7e-12)*(1.02e-6) ;perm coef eo_FENSSCROSS=0.178 ;porosity end prop_FENSSCROSS ;;;;;;UNDERBURDEN ;;;;;; def prop_UNDERBURDEN e_UNDERBURDEN = 4480e6 ;E mod p_UNDERBURDEN = 0.3 ;porsion bul_UNDERBURDEN = e_UNDERBURDEN/(3*(1-2*p_UNDERBURDEN)) ;bulk mod she_UNDERBURDEN = e_UNDERBURDEN/(2*(1+p_UNDERBURDEN)) ;shear mod coh_UNDERBURDEN = 21e3 ;cohesion fric_UNDERBURDEN = 26.78 ;firction dens_UNDERBURDEN = 2437 ;density k_UNDERBURDEN=1.0e-13 ;perm coef ; k_UNDERBURDEN=(7e-12)*(1.02e-6) perm coef eo_UNDERBURDEN=0.11 ;porosity end prop_UNDERBURDEN ;;;;;; ;;;;;; def prop_FAULT e_FAULT = 1933e6 ;E mod p_FAULT = 0.3621 ;porsion bul_FAULT = e_FAULT/(3*(1-2*p_FAULT)) ;bulk mod she_FAULT = e_FAULT/(2*(1+p_FAULT)) ;shear mod coh_FAULT = 21e3 ;cohesion fric_FAULT = 21.47 ;firction dens_FAULT = 2246.5 ;density k_FAULT=8e-11 ;perm coef ; k_FAULT=(7e-12)*(1.02e-6) ;perm coef eo_FAULT=0.25 ;porosity end prop_FAULT ; Define material models group 'UNDERBURDEN' i=1,264 j=1,21 group 'UNDERBURDEN' i=265,330 j=1,21 group 'UNDERBURDEN' i=332,334 j=1,21 group 'UNDERBURDEN' i=336,340 j=1,21 group 'UNDERBURDEN' i=341,346 j=1,21 group 'UNDERBURDEN' i=347,352 j=1,21 group 'UNDERBURDEN' i=353,358 j=1,21 group 'UNDERBURDEN' i=359,373 j=1,21 group 'UNDERBURDEN' i=374,380 j=1,21 group 'UNDERBURDEN' i=381,385 j=1,21 group 'UNDERBURDEN' i=386,582 j=1,21 group 'UNDERBURDEN' i=583,776 j=1,21 group 'UNDERBURDEN' i=336,340 j=22,30 group 'UNDERBURDEN' i=341,346 j=22,30 group 'UNDERBURDEN' i=347,352 j=22,30 group 'UNDERBURDEN' i=353,358 j=22,30 group 'UNDERBURDEN' i=359,373 j=22,30 group 'UNDERBURDEN' i=374,380 j=22,30 group 'UNDERBURDEN' i=381,385 j=22,30 group 'UNDERBURDEN' i=386,582 j=22,30 group 'UNDERBURDEN' i=583,776 j=22,30 group 'FENSSCROSS' i=778,1041 j=1,4 group 'FENSSCROSS' i=1042,1107 j=1,4 group 'FENSSCROSS' i=1108,1116 j=1,4 group 'FENSSCROSS' i=336,340 j=31,45 group 'FENSSCROSS' i=341,346 j=31,45 group 'FENSSCROSS' i=347,352 j=31,45 group 'FENSSCROSS' i=353,358 j=31,45 group 'FENSSCROSS' i=359,373 j=31,45 group 'FENSSCROSS' i=374,380 j=31,45 group 'FENSSCROSS' i=381,385 j=31,45 group 'FENSSCROSS' i=386,582 j=31,45 group 'FENSSCROSS' i=583,776 j=31,45 group 'CAPROCK' i=1126,1140 j=4,9 group 'CAPROCK' i=1117,1120 j=5,10 group 'CAPROCK' i=1142,1405 j=1,6 group 'CAPROCK' i=1406,1471 j=1,6 group 'CAPROCK' i=1108,1116 j=5,10 group 'CAPROCK' i=336,340 j=46,52 group 'CAPROCK' i=341,346 j=46,52 group 'CAPROCK' i=347,352 j=46,52 group 'CAPROCK' i=353,358 j=46,52 group 'CAPROCK' i=359,373 j=46,52 group 'CAPROCK' i=374,380 j=46,52 group 'CAPROCK' i=381,385 j=46,52 group 'CAPROCK' i=386,582 j=46,52 group 'CAPROCK' i=2087,2093 j=1,6 group 'CAPROCK' i=583,776 j=46,52 group 'ASGARD' i=1473,1736 j=1,7 group 'ASGARD' i=1406,1471 j=7,13 group 'ASGARD' i=1108,1116 j=11,17 group 'ASGARD' i=1117,1120 j=11,17 group 'ASGARD' i=1126,1140 j=10,16 group 'ASGARD' i=2095,2101 j=1,7 group 'ASGARD' i=374,380 j=53,59 group 'ASGARD' i=359,373 j=53,59 group 'ASGARD' i=381,385 j=53,59 group 'ASGARD' i=386,582 j=53,59 group 'ASGARD' i=2103,2108 j=1,7 group 'ASGARD' i=341,346 j=53,59 group 'ASGARD' i=2113,2117 j=1,7 group 'ASGARD' i=347,352 j=53,59 group 'ASGARD' i=353,358 j=53,59 group 'ASGARD' i=583,776 j=53,59 group 'SVARTE' i=1742,2005 j=1,8 group 'SVARTE' i=2006,2071 j=1,8 group 'SVARTE' i=2073,2081 j=1,8 group 'SVARTE' i=2082,2085 j=1,8 group 'SVARTE' i=1126,1140 j=17,24 group 'SVARTE' i=2095,2101 j=8,15 group 'SVARTE' i=353,358 j=60,67 group 'SVARTE' i=2125,2130 j=1,8 group 'SVARTE' i=2119,2124 j=1,8 group 'SVARTE' i=359,373 j=60,67 group 'SVARTE' i=2103,2108 j=8,15 group 'SVARTE' i=2135,2139 j=1,8 group 'SVARTE' i=374,380 j=60,67 group 'SVARTE' i=381,385 j=60,67 group 'SVARTE' i=386,582 j=60,67 group 'SVARTE' i=583,776 j=60,67 group 'SHETLAND' i=1742,2005 j=9,12 group 'SHETLAND' i=2006,2071 j=9,12 group 'SHETLAND' i=2073,2081 j=9,12 group 'SHETLAND' i=2082,2085 j=9,12 group 'SHETLAND' i=1126,1140 j=25,28 group 'SHETLAND' i=2141,2404 j=1,3 group 'SHETLAND' i=2405,2470 j=1,3 group 'SHETLAND' i=2472,2480 j=1,3 group 'SHETLAND' i=2481,2484 j=1,3 group 'SHETLAND' i=1126,1140 j=29,31 group 'SHETLAND' i=2486,2492 j=1,4 group 'SHETLAND' i=2103,2108 j=16,19 group 'SHETLAND' i=2486,2492 j=5,7 group 'SHETLAND' i=2132,2134 j=9,12 group 'SHETLAND' i=2135,2139 j=9,12 group 'SHETLAND' i=2119,2124 j=9,12 group 'SHETLAND' i=2125,2130 j=9,12 group 'SHETLAND' i=2494,2499 j=1,4 group 'SHETLAND' i=2501,2515 j=1,4 group 'SHETLAND' i=2103,2108 j=20,22 group 'SHETLAND' i=2517,2519 j=1,3 group 'SHETLAND' i=2520,2524 j=1,3 group 'SHETLAND' i=2526,2532 j=1,4 group 'SHETLAND' i=2548,2553 j=1,3 group 'SHETLAND' i=2494,2499 j=5,7 group 'SHETLAND' i=2542,2547 j=1,3 group 'SHETLAND' i=2501,2515 j=5,7 group 'SHETLAND' i=2526,2532 j=5,7 group 'SHETLAND' i=381,385 j=68,70 group 'SHETLAND' i=2555,2558 j=1,3 group 'LISTA' i=2141,2404 j=4,7 group 'LISTA' i=2405,2470 j=4,7 group 'LISTA' i=2472,2480 j=4,7 group 'LISTA' i=2481,2484 j=4,7 group 'LISTA' i=1126,1140 j=32,35 group 'LISTA' i=2534,2540 j=1,4 group 'LISTA' i=2103,2108 j=23,26 group 'LISTA' i=2517,2519 j=4,7 group 'LISTA' i=2520,2524 j=4,7 group 'LISTA' i=2542,2547 j=4,7 group 'LISTA' i=2548,2553 j=4,7 group 'LISTA' i=2891,2896 j=1,4 group 'LISTA' i=2898,2912 j=1,4 group 'LISTA' i=2526,2532 j=8,11 group 'LISTA' i=2555,2558 j=4,7 group 'LISTA' i=381,385 j=71,74 group 'SELE' i=2560,2823 j=1,3 group 'SELE' i=2824,2889 j=1,3 group 'SELE' i=3109,3117 j=1,3 group 'SELE' i=3118,3121 j=1,3 group 'SELE' i=1126,1140 j=36,38 group 'SELE' i=2534,2540 j=5,7 group 'SELE' i=2103,2108 j=27,29 group 'SELE' i=3126,3130 j=1,3 group 'SELE' i=3123,3125 j=1,3 group 'SELE' i=3132,3137 j=1,3 group 'SELE' i=3138,3143 j=1,3 group 'SELE' i=2891,2896 j=5,7 group 'SELE' i=2898,2912 j=5,7 group 'SELE' i=2526,2532 j=12,14 group 'SELE' i=2555,2558 j=8,10 group 'SELE' i=381,385 j=75,77 group 'NORTLAND' i=386,582 j=68,70 group 'NORTLAND' i=583,776 j=68,70 group 'NORTLAND' i=386,582 j=71,74 group 'NORTLAND' i=2914,3107 j=1,4 group 'NORTLAND' i=386,582 j=75,77 group 'NORTLAND' i=2914,3107 j=5,7 group 'NORTLAND' i=3145,3408 j=1,7 group 'NORTLAND' i=2824,2889 j=4,10 group 'NORTLAND' i=3410,3418 j=1,7 group 'NORTLAND' i=3118,3121 j=4,10 group 'NORTLAND' i=1126,1140 j=39,45 group 'NORTLAND' i=3420,3426 j=1,7 group 'NORTLAND' i=2103,2108 j=30,36 group 'NORTLAND' i=3428,3430 j=1,7 group 'NORTLAND' i=3126,3130 j=4,10 group 'NORTLAND' i=3432,3437 j=1,7 group 'NORTLAND' i=3138,3143 j=4,10 group 'NORTLAND' i=3439,3444 j=1,7 group 'NORTLAND' i=3446,3460 j=1,7 group 'NORTLAND' i=2526,2532 j=15,21 group 'NORTLAND' i=2555,2558 j=11,17 group 'NORTLAND' i=381,385 j=78,84 group 'NORTLAND' i=386,582 j=78,84 group 'NORTLAND' i=3462,3655 j=1,7 group 'FAULT' i=1122,1124 j=1,9 group 'FAULT' i=1126,1140 j=1,3 group 'FAULT' i=1738,1740 j=1,7 group 'FAULT' i=2110,2112 j=1,7 group 'FAULT' i=2132,2134 j=1,8 ;apply material prop model mohr group 'UNDERBURDEN' prop density=dens_UNDERBURDEN bulk=bul_UNDERBURDEN shear=she_UNDERBURDEN cohesion=coh_UNDERBURDEN friction=fric_UNDERBURDEN dilation=0.0 tension=0.0 group 'UNDERBURDEN' initial wperm=k_UNDERBURDEN poro=eo_UNDERBURDEN rsat=0.3 group 'UNDERBURDEN' model mohr group 'FENSSCROSS' prop density=dens_FENSSCROSS bulk=bul_FENSSCROSS shear=she_FENSSCROSS cohesion=coh_FENSSCROSS friction=fric_FENSSCROSS dilation=0.0 tension=0.0 group 'FENSSCROSS' initial wperm=k_FENSSCROSS poro=eo_FENSSCROSS rsat=0.3 group 'FENSSCROSS' model mohr group 'CAPROCK' prop density=dens_CAPROCK bulk=bul_CAPROCK shear=she_CAPROCK cohesion=coh_CAPROCK friction=fric_CAPROCK dilation=0.0 tension=0.0 group 'CAPROCK' initial wperm=k_CAPROCK poro=eo_CAPROCK rsat=0.3 group 'CAPROCK' model mohr group 'ASGARD' prop density=dens_ASGARD bulk=bul_ASGARD shear=she_ASGARD cohesion=coh_ASGARD friction=fric_ASGARD dilation=0.0 tension=0.0 group 'ASGARD' initial wperm=k_ASGARD poro=eo_ASGARD rsat=0.3 group 'ASGARD' model mohr group 'SVARTE' prop density=dens_SVARTE bulk=bul_SVARTE shear=she_SVARTE cohesion=coh_SVARTE friction=fric_SVARTE dilation=0.0 tension=0.0 group 'SVARTE' initial wperm=k_SVARTE poro=eo_SVARTE rsat=0.3 group 'SVARTE' model mohr group 'SHETLAND' prop density=dens_SHETLAND bulk=bul_SHETLAND shear=she_SHETLAND cohesion=coh_SHETLAND friction=fric_SHETLAND dilation=0.0 tension=0.0 group 'SHETLAND' initial wperm=k_SHETLAND poro=eo_SHETLAND rsat=0.3 group 'SHETLAND' model mohr group 'LISTA' prop density=dens_LISTA bulk=bul_LISTA shear=she_LISTA cohesion=coh_LISTA friction=fric_LISTA dilation=0.0 tension=0.0 group 'LISTA' initial wperm=k_LISTA poro=eo_LISTA rsat=0.3 group 'LISTA' model mohr group 'SELE' prop density=dens_SELE bulk=bul_SELE shear=she_SELE cohesion=coh_SELE friction=fric_SELE dilation=0.0 tension=0.0 group 'SELE' initial wperm=k_SELE visrat=1.0 poro=eo_SELE rsat=0.3 group 'SELE' model mohr group 'NORTLAND' prop density=dens_NORTLAND bulk=bul_NORTLAND shear=she_NORTLAND cohesion=coh_NORTLAND friction=fric_NORTLAND dilation=0.0 tension=0.0 group 'NORTLAND' initial wperm=k_NORTLAND poro=eo_NORTLAND rsat=0.3 group 'NORTLAND' model mohr group 'FAULT' prop density=dens_FAULT bulk=bul_FAULT shear=she_FAULT cohesion=coh_FAULT friction=fric_FAULT dilation=0.0 tension=0.0 group 'FAULT' initial wperm=k_FAULT poro=eo_FAULT rsat=0.3 group 'FAULT' ; Attach grids attach aside from 778,1 to 1042,1 bside from 1,22 to 265,22 attach aside from 1042,1 to 1108,1 bside from 265,22 to 331,22 attach aside from 332,1 to 332,22 bside from 331,1 to 331,22 attach aside from 1122,1 to 1125,1 bside from 332,22 to 335,22 attach aside from 336,1 to 336,22 bside from 335,1 to 335,22 attach aside from 1142,1 to 1406,1 bside from 778,5 to 1042,5 attach aside from 1406,1 to 1472,1 bside from 1042,5 to 1108,5 attach aside from 1117,5 to 1121,5 bside from 1117,5 to 1117,1 attach aside from 1122,1 to 1122,10 bside from 1108,1 to 1117,1 attach aside from 1126,4 to 1126,1 bside from 1122,10 to 1125,10 attach aside from 336,22 to 336,31 bside from 1125,1 to 1125,10 attach aside from 1738,1 to 1741,1 bside from 1141,4 to 1141,1 attach aside from 336,31 to 336,46 bside from 1126,1 to 1141,1 attach aside from 1121,11 to 1121,5 bside from 1126,10 to 1126,4 attach aside from 2087,7 to 2087,1 bside from 1141,10 to 1141,4 attach aside from 1473,1 to 1737,1 bside from 1142,7 to 1406,7 attach aside from 1108,11 to 1108,5 bside from 1472,7 to 1472,1 attach aside from 1742,1 to 2006,1 bside from 1473,8 to 1737,8 attach aside from 1406,7 to 1406,14 bside from 1737,1 to 1737,8 attach aside from 2006,1 to 2072,1 bside from 1406,14 to 1472,14 attach aside from 1108,11 to 1108,18 bside from 1472,7 to 1472,14 attach aside from 2073,1 to 2082,1 bside from 1108,18 to 1117,18 attach aside from 2082,1 to 2086,1 bside from 1117,18 to 1121,18 attach aside from 1126,10 to 1126,17 bside from 1121,11 to 1121,18 attach aside from 2095,1 to 2095,8 bside from 1141,10 to 1141,17 attach aside from 1741,1 to 1741,8 bside from 336,46 to 336,53 attach aside from 2113,1 to 2118,1 bside from 336,53 to 341,53 attach aside from 2087,1 to 2094,1 bside from 1738,1 to 1738,8 attach aside from 2110,1 to 2113,1 bside from 1738,8 to 1741,8 attach aside from 2073,1 to 2073,9 bside from 2072,1 to 2072,9 attach aside from 1126,17 to 1126,25 bside from 2086,1 to 2086,9 attach aside from 2095,1 to 2102,1 bside from 2087,7 to 2094,7 attach aside from 2103,1 to 2109,1 bside from 2094,7 to 2094,1 attach aside from 2095,8 to 2095,16 bside from 1141,17 to 1141,25 attach aside from 2103,1 to 2103,8 bside from 2102,1 to 2102,8 attach aside from 2110,1 to 2110,8 bside from 2109,1 to 2109,8 attach aside from 2132,1 to 2135,1 bside from 2110,8 to 2113,8 attach aside from 2118,1 to 2118,8 bside from 341,53 to 341,60 attach aside from 2119,1 to 2125,1 bside from 341,60 to 347,60 attach aside from 2135,1 to 2140,1 bside from 2113,8 to 2118,8 attach aside from 2125,1 to 2131,1 bside from 347,60 to 353,60 attach aside from 2141,1 to 2405,1 bside from 1742,13 to 2006,13 attach aside from 2405,1 to 2471,1 bside from 2006,13 to 2072,13 attach aside from 2073,13 to 2073,9 bside from 2072,13 to 2072,9 attach aside from 2472,1 to 2481,1 bside from 2073,13 to 2082,13 attach aside from 2481,1 to 2485,1 bside from 2082,13 to 2086,13 attach aside from 1126,29 to 1126,25 bside from 2086,13 to 2086,9 attach aside from 2486,5 to 2486,1 bside from 1141,29 to 1141,25 attach aside from 2486,1 to 2493,1 bside from 2095,16 to 2102,16 attach aside from 2103,8 to 2103,16 bside from 2102,8 to 2102,16 attach aside from 2131,1 to 2131,9 bside from 353,60 to 353,68 attach aside from 2494,1 to 2500,1 bside from 353,68 to 359,68 attach aside from 2140,1 to 2140,9 bside from 2119,1 to 2119,9 attach aside from 2501,1 to 2516,1 bside from 359,68 to 374,68 attach aside from 2109,8 to 2109,16 bside from 2132,1 to 2132,9 attach aside from 2526,1 to 2533,1 bside from 374,68 to 381,68 attach aside from 2472,1 to 2472,4 bside from 2471,1 to 2471,4 attach aside from 1126,29 to 1126,32 bside from 2485,1 to 2485,4 attach aside from 2486,5 to 2486,8 bside from 1141,29 to 1141,32 attach aside from 2103,20 to 2103,16 bside from 2493,5 to 2493,1 attach aside from 2132,13 to 2132,9 bside from 2109,20 to 2109,16 attach aside from 2560,1 to 2824,1 bside from 2141,8 to 2405,8 attach aside from 2824,1 to 2890,1 bside from 2405,8 to 2471,8 attach aside from 2472,4 to 2472,8 bside from 2471,4 to 2471,8 attach aside from 2534,1 to 2541,1 bside from 2486,8 to 2493,8 attach aside from 2103,20 to 2103,23 bside from 2493,5 to 2493,8 attach aside from 2517,1 to 2520,1 bside from 2132,13 to 2135,13 attach aside from 2520,1 to 2525,1 bside from 2135,13 to 2140,13 attach aside from 2119,13 to 2119,9 bside from 2140,13 to 2140,9 attach aside from 2542,1 to 2548,1 bside from 2119,13 to 2125,13 attach aside from 2548,1 to 2554,1 bside from 2125,13 to 2131,13 attach aside from 2494,5 to 2494,1 bside from 2131,13 to 2131,9 attach aside from 2501,5 to 2501,1 bside from 2500,5 to 2500,1 attach aside from 2526,5 to 2526,1 bside from 2516,5 to 2516,1 attach aside from 2517,1 to 2517,4 bside from 2109,20 to 2109,23 attach aside from 3109,1 to 3118,1 bside from 2472,8 to 2481,8 attach aside from 3118,1 to 3122,1 bside from 2481,8 to 2485,8 attach aside from 1126,32 to 1126,36 bside from 2485,4 to 2485,8 attach aside from 2534,1 to 2534,5 bside from 1141,32 to 1141,36 attach aside from 2542,1 to 2542,4 bside from 2525,1 to 2525,4 attach aside from 2555,1 to 2559,1 bside from 2533,5 to 2533,1 attach aside from 2103,23 to 2103,27 bside from 2541,1 to 2541,5 attach aside from 2494,5 to 2494,8 bside from 2554,1 to 2554,4 attach aside from 2891,1 to 2897,1 bside from 2494,8 to 2500,8 attach aside from 2501,5 to 2501,8 bside from 2500,5 to 2500,8 attach aside from 2898,1 to 2913,1 bside from 2501,8 to 2516,8 attach aside from 2526,5 to 2526,8 bside from 2516,5 to 2516,8 attach aside from 2517,4 to 2517,8 bside from 2109,23 to 2109,27 attach aside from 3123,1 to 3126,1 bside from 2517,8 to 2520,8 attach aside from 3126,1 to 3131,1 bside from 2520,8 to 2525,8 attach aside from 2542,4 to 2542,8 bside from 2525,4 to 2525,8 attach aside from 3132,1 to 3138,1 bside from 2542,8 to 2548,8 attach aside from 3138,1 to 3144,1 bside from 2548,8 to 2554,8 attach aside from 2891,1 to 2891,5 bside from 2554,4 to 2554,8 attach aside from 2555,1 to 2555,4 bside from 2533,5 to 2533,8 attach aside from 2559,1 to 2559,4 bside from 381,68 to 381,71 attach aside from 3145,1 to 3409,1 bside from 2560,4 to 2824,4 attach aside from 2898,1 to 2898,5 bside from 2897,1 to 2897,5 attach aside from 2526,8 to 2526,12 bside from 2913,1 to 2913,5 attach aside from 2914,1 to 3108,1 bside from 583,71 to 777,71 attach aside from 2555,4 to 2555,8 bside from 2533,8 to 2533,12 attach aside from 381,71 to 381,75 bside from 2559,4 to 2559,8 attach aside from 3109,1 to 3109,4 bside from 2890,1 to 2890,4 attach aside from 2914,1 to 2914,5 bside from 583,71 to 583,75 attach aside from 3410,1 to 3419,1 bside from 3109,4 to 3118,4 attach aside from 1126,36 to 1126,39 bside from 3122,1 to 3122,4 attach aside from 2534,5 to 2534,8 bside from 1141,36 to 1141,39 attach aside from 3420,1 to 3427,1 bside from 2534,8 to 2541,8 attach aside from 2103,27 to 2103,30 bside from 2541,5 to 2541,8 attach aside from 3123,1 to 3123,4 bside from 2109,27 to 2109,30 attach aside from 3132,1 to 3132,4 bside from 3131,1 to 3131,4 attach aside from 3428,1 to 3431,1 bside from 3123,4 to 3126,4 attach aside from 3439,1 to 3445,1 bside from 3132,4 to 3138,4 attach aside from 2891,5 to 2891,8 bside from 3144,1 to 3144,4 attach aside from 3432,1 to 3438,1 bside from 2891,8 to 2897,8 attach aside from 2898,5 to 2898,8 bside from 2897,5 to 2897,8 attach aside from 3446,1 to 3461,1 bside from 2898,8 to 2913,8 attach aside from 2526,12 to 2526,15 bside from 2913,5 to 2913,8 attach aside from 2555,8 to 2555,11 bside from 2533,12 to 2533,15 attach aside from 381,75 to 381,78 bside from 2559,8 to 2559,11 attach aside from 2914,5 to 2914,8 bside from 583,75 to 583,78 attach aside from 3462,1 to 3656,1 bside from 2914,8 to 3108,8 attach aside from 2824,4 to 2824,11 bside from 3409,1 to 3409,8 attach aside from 3410,1 to 3410,8 bside from 2890,4 to 2890,11 attach aside from 3118,4 to 3118,11 bside from 3419,1 to 3419,8 attach aside from 1126,39 to 1126,46 bside from 3122,4 to 3122,11 attach aside from 3420,1 to 3420,8 bside from 1141,39 to 1141,46 attach aside from 2103,30 to 2103,37 bside from 3427,1 to 3427,8 attach aside from 3428,1 to 3428,8 bside from 2109,30 to 2109,37 attach aside from 3126,4 to 3126,11 bside from 3431,1 to 3431,8 attach aside from 3439,1 to 3439,8 bside from 3131,4 to 3131,11 attach aside from 3144,4 to 3144,11 bside from 3432,1 to 3432,8 attach aside from 3446,1 to 3446,8 bside from 3438,1 to 3438,8 attach aside from 3445,1 to 3445,8 bside from 3138,4 to 3138,11 attach aside from 2526,15 to 2526,22 bside from 3461,1 to 3461,8 attach aside from 2555,11 to 2555,18 bside from 2533,15 to 2533,22 attach aside from 381,78 to 381,85 bside from 2559,11 to 2559,18 attach aside from 3462,1 to 3462,8 bside from 583,78 to 583,85 ; Fixed boundary conditions fix x i=1 j=1,22 fix x y i=1,265 j=1 fix x y i=265,331 j=1 fix x y i=332,335 j=1 fix x y i=336,341 j=1 fix x y i=341,347 j=1 fix x y i=347,353 j=1 fix x y i=353,359 j=1 fix x y i=359,374 j=1 fix x y i=374,381 j=1 fix x y i=381,386 j=1 fix x y i=386,583 j=1 fix x i=777 j=1,22 fix x y i=583,777 j=1 fix x i=778 j=1,5 fix x i=777 j=22,31 fix x i=777 j=31,46 fix x i=1142 j=1,7 fix x i=1473 j=1,8 fix x i=1742 j=1,9 fix x i=777 j=46,53 fix x i=1742 j=9,13 fix x i=777 j=53,60 fix x i=2141 j=1,4 fix x i=2141 j=4,8 fix x i=777 j=60,68 fix x i=2560 j=1,4 fix x i=777 j=68,71 fix x i=3108 j=1,5 fix x i=3108 j=5,8 fix x i=3145 j=1,8 fix x i=3656 j=1,8 save new_S0_model.sav