Index: /issm/workshop/07UQ_PIG/Exp_Par/ErrorContour.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/ErrorContour.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/ErrorContour.exp (revision 18013) @@ -0,0 +1,24 @@ +## Name:ErrorContour +## Icon:0 +# Points Count Value +18 1.000000 +# X pos Y pos +-1642090.6657102096 -328038.7002928721 +-1665417.0018035804 -319119.8070807009 +-1681882.6508106657 -304712.3641995014 +-1666103.0705122089 -245710.4552574461 +-1632485.7037894099 -232675.1497935036 +-1616020.0547823247 -224442.3252899610 +-1598182.2683579824 -200429.9204879618 +-1596810.1309407253 -181219.9966463624 +-1615333.9860736963 -160637.9353875059 +-1559762.4206747836 -147602.6299235635 +-1527517.1913692418 -188080.6837326479 +-1566623.1077610692 -218953.7756209326 +-1533691.8097468987 -229244.8062503609 +-1523400.7791174706 -309514.8451599012 +-1533691.8097468987 -350678.9676776142 +-1594751.9248148398 -371261.0289364707 +-1636602.1160411814 -349306.8302603571 +-1642090.6657102096 -328038.7002928721 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux1.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux1.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux1.exp (revision 18013) @@ -0,0 +1,22 @@ +## Name:MassFlux1 +## Icon:0 +# Points Count Value +16 1.000000 +# X pos Y pos +-1638344.3538415846 -328022.1178486745 +-1636175.0592654669 -328022.1178486745 +-1634647.1933393457 -327710.7373731936 +-1631836.4701132313 -326087.2311287920 +-1629847.9500851235 -324428.7567974643 +-1626413.2336729371 -321756.7703747695 +-1622074.6445207016 -320098.2960434417 +-1618278.3790124957 -319637.6087291840 +-1614210.9516822749 -319729.7461920356 +-1609601.2007080249 -320466.8458948479 +-1604539.5133637502 -322954.5573918396 +-1602370.2187876324 -324613.0317231673 +-1599206.6641974607 -326824.3308316044 +-1595681.5605112694 -327008.6057573075 +-1593150.7168391321 -327653.5679972683 +-1590710.2604409996 -330786.2417342208 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux10.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux10.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux10.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux11 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1539536.1089629673 -171232.8272031187 +-1540384.1546699132 -200614.0179257037 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux11.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux11.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux11.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux12 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1559431.5986180613 -236969.7272119786 +-1549112.1094185982 -261719.2194756830 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux12.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux12.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux12.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux13 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1572656.1032135738 -292391.9615161452 +-1573651.1418990945 -315407.7936937447 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux13.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux13.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux13.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux14 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1579674.9916810417 -321516.2892665596 +-1587436.3185859676 -337140.9687259690 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux14.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux14.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux14.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux14 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1579674.9916810417 -321516.2892665596 +-1587436.3185859676 -337140.9687259690 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux15.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux15.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux15.exp (revision 18013) @@ -0,0 +1,127 @@ +## Name: +## Icon:0 +# Points Count Value +121 25000.000000 +# X pos Y pos +-1712012.0818103175 -348661.1275967784 +-1711738.9446965172 -348380.8584401013 +-1709339.1962288069 -348358.9741029009 +-1707566.1624406802 -348048.0535193763 +-1706142.6423542583 -346463.0619396989 +-1703644.8906504354 -344421.7875470470 +-1701428.0245438961 -343211.6412474416 +-1699098.6802015656 -343202.0651453367 +-1698767.2671226400 -341251.8958583454 +-1694089.8466347910 -339325.0408948531 +-1691249.4280928215 -338644.2293668062 +-1688271.0930389061 -338423.0784381799 +-1686131.0738047434 -339115.2406933784 +-1682853.5065837551 -340300.2702040364 +-1681648.3224467863 -342890.2203795200 +-1681160.5378740367 -343956.3892846339 +-1681343.2453770319 -345059.2736684235 +-1680476.2008784728 -345129.3389066600 +-1679673.6094558411 -346405.5378093903 +-1679898.5797234357 -347818.5803631411 +-1677215.4777375804 -347397.3820658240 +-1674583.5118423651 -349115.8291018623 +-1675803.3920552148 -352985.6713725598 +-1672935.8509424068 -352801.1873251720 +-1668529.5198238282 -353048.3463822120 +-1666370.5689502691 -354473.4893832008 +-1665423.6161039397 -356733.4691418489 +-1662720.0686237188 -357334.8975107344 +-1660276.4961961196 -357869.6209268258 +-1657124.6115999895 -361080.6428832713 +-1655993.6074009922 -363193.2597271361 +-1653688.2093358042 -364017.9494011777 +-1654435.7188006472 -361506.6539926908 +-1652050.9793420888 -355295.4585798742 +-1649364.4081354018 -349945.5944371450 +-1645761.7354861272 -350839.5079576963 +-1645996.5509334437 -351783.7951824765 +-1645345.5744557772 -351250.5634369315 +-1644300.2406060672 -352221.4198542684 +-1644633.1139884400 -348521.7302174515 +-1647057.7191411837 -348406.3266729190 +-1645415.0137059502 -341247.9400814587 +-1644395.7076754270 -340417.8557739360 +-1644169.8943326229 -336477.7151767056 +-1643216.0658811626 -334392.3056475453 +-1641434.0471865973 -334447.2336099123 +-1641190.7738541283 -334332.2252321915 +-1641050.9552025823 -333796.1902735520 +-1639787.1689887117 -331073.8439645328 +-1637807.1154673321 -329711.1825533187 +-1635805.1587214153 -327827.0308139678 +-1632065.9795259857 -326063.4189332099 +-1625929.4636559838 -324999.2470617331 +-1619584.9279456623 -324892.7557348280 +-1615193.9710026928 -325097.1546172730 +-1611296.7956007582 -324995.1262699560 +-1608304.4173934418 -325211.0183405159 +-1603316.9037486394 -326424.5673739603 +-1600245.1111844976 -328836.3685103582 +-1595818.7639296530 -331148.7367850624 +-1591649.1976614639 -334787.7632693424 +-1588371.5600852505 -337840.7799059978 +-1587648.7608976758 -343365.9382903088 +-1587132.8425864622 -345049.3782468919 +-1586012.7810109658 -347625.9740446581 +-1588260.6470941922 -349356.1732009141 +-1584626.3352827206 -352498.3947002917 +-1584121.7278491815 -354575.2387267592 +-1584906.9380743054 -356537.0120279597 +-1585446.5634405916 -359133.6460951389 +-1586894.3150289655 -361098.6217195836 +-1585768.2162689434 -363527.8170254802 +-1584807.6500740163 -365694.8871387623 +-1585744.0045738907 -368358.0113804646 +-1585152.9531241022 -369694.9777909298 +-1584499.3876485692 -371017.3807505022 +-1585456.9434211247 -372944.8898924028 +-1429208.6248950495 -363616.6734221856 +-1415050.5145378630 -326851.1264040734 +-1413047.7406149551 -325411.1419014899 +-1405442.9810691341 -295790.2533492565 +-1388007.7873070226 -266002.2718567938 +-1375120.7873420289 -286146.6923060659 +-1373971.2210923850 -285345.4403176561 +-1346034.9931278315 -322527.8743193840 +-1351274.5360716779 -327267.6902373450 +-1382747.3419384742 -297522.4999399224 +-1382639.2845126297 -296126.1700367680 +-1382479.2357768724 -295804.4343703665 +-1397758.7640077863 -316499.0520702488 +-1408929.4248104726 -362684.2359936420 +-1202573.2098436793 -350320.7957165864 +-1192371.9994205239 -349816.1301302491 +-1233144.4335198263 -292276.3631183850 +-1232225.6096719459 -288608.8406540032 +-1360692.9373680721 -108664.8182541504 +-1424223.8165527256 -81531.7934136631 +-1403445.8481024427 -45960.5859752788 +-1478312.1378749039 59339.3203815984 +-1513118.6444103392 3225.4273728387 +-1524169.2804273698 11860.8708815790 +-1504759.9852699745 36302.7657860445 +-1528911.0294425790 62751.8076344339 +-1622627.2364271644 61937.4929792271 +-1647011.3496706884 32864.6807476235 +-1637527.3859760596 19912.7371761968 +-1650069.2631712460 30042.0704865309 +-1724734.6043381831 -61378.0233897057 +-1718614.9048978437 -69045.3133313945 +-1758091.3164752426 -113926.3236823874 +-1754430.7541207569 -116894.8054778867 +-1760233.4203404333 -126776.6398906443 +-1757170.3120773330 -130521.7753713153 +-1762186.6328915418 -136423.2047200473 +-1765126.7323215308 -132420.4516091067 +-1769785.4415677825 -135199.8713571306 +-1774818.6171028516 -132429.3647437230 +-1793280.1507095622 -153411.6536596067 +-1711636.5537849038 -224740.3773715630 +-1729008.9305659365 -335286.5679692774 +-1712012.0818103175 -348661.1275967784 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux2.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux2.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux2.exp (revision 18013) @@ -0,0 +1,171 @@ +## Name:MassFlux2 +## Icon:0 +# Points Count Value +165 1.000000 +# X pos Y pos +-1609993.8450579999 -266189.3103270000 +-1609790.4465310001 -265898.8271340000 +-1609587.0480040000 -265608.3439410000 +-1609383.6494770001 -265317.8607480000 +-1609136.7086169999 -265274.3184150000 +-1608933.3100900000 -264983.8352220000 +-1608686.3692300001 -264940.2928890000 +-1608236.0298430000 -264606.2673630000 +-1607989.0889830000 -264562.7250310000 +-1607785.6904559999 -264272.2418380000 +-1607538.7495960000 -264228.6995050000 +-1607291.8087370000 -264185.1571720000 +-1607044.8678770000 -264141.6148390000 +-1606797.9270170000 -264098.0725060000 +-1606304.0452970001 -264010.9878410000 +-1606057.1044369999 -263967.4455080000 +-1605810.1635769999 -263923.9031750000 +-1605563.2227179999 -263880.3608420000 +-1605316.2818580000 -263836.8185100000 +-1605069.3409980000 -263793.2761770000 +-1604778.8578049999 -263996.6747040000 +-1604531.9169449999 -263953.1323710000 +-1604038.0352250000 -263866.0477060000 +-1603791.0943660000 -263822.5053730000 +-1603500.6111730000 -264025.9039000000 +-1603253.6703130000 -263982.3615680000 +-1602963.1871199999 -264185.7600950000 +-1602716.2462599999 -264142.2177620000 +-1602425.7630680001 -264345.6162900000 +-1601888.3390150000 -264505.4724840000 +-1601641.3981550001 -264461.9301510000 +-1601350.9149630000 -264665.3286790000 +-1601060.4317699999 -264868.7272060000 +-1600479.4653840000 -265275.5242610000 +-1600188.9821919999 -265478.9227880000 +-1599898.4989990001 -265682.3213160000 +-1599608.0158060000 -265885.7198430000 +-1599027.0494210001 -266292.5168980000 +-1598736.5662280000 -266495.9154250000 +-1598446.0830349999 -266699.3139520000 +-1598402.5407030000 -266946.2548130000 +-1598605.9392299999 -267236.7380050000 +-1598518.8545639999 -267730.6197260000 +-1598722.2530910000 -268021.1029190000 +-1598678.7107579999 -268268.0437790000 +-1598635.1684250000 -268514.9846390000 +-1598591.6260919999 -268761.9254990000 +-1598795.0246190000 -269052.4086920000 +-1598751.4822869999 -269299.3495520000 +-1598707.9399540001 -269546.2904130000 +-1598620.8552880001 -270040.1721330000 +-1598577.3129550000 -270287.1129930000 +-1598533.7706220001 -270534.0538530000 +-1598737.1691490000 -270824.5370460000 +-1598693.6268170001 -271071.4779060000 +-1598650.0844840000 -271318.4187660000 +-1598606.5421509999 -271565.3596270000 +-1598316.0589580000 -271768.7581540000 +-1598228.9742930001 -272262.6398740000 +-1597938.4911000000 -272466.0384020000 +-1597648.0079069999 -272669.4369290000 +-1597067.0415220000 -273076.2339840000 +-1597023.4991890001 -273323.1748440000 +-1596979.9568560000 -273570.1157040000 +-1596936.4145229999 -273817.0565640000 +-1596892.8721900000 -274063.9974240000 +-1596805.7875250001 -274557.8791450000 +-1596762.2451919999 -274804.8200050000 +-1596718.7028590001 -275051.7608650000 +-1596675.1605260000 -275298.7017250000 +-1596631.6181930001 -275545.6425850000 +-1596341.1350010000 -275749.0411130000 +-1596297.5926679999 -275995.9819730000 +-1595963.5671420000 -276446.3213600000 +-1595920.0248090001 -276693.2622210000 +-1595629.5416170000 -276896.6607480000 +-1595339.0584239999 -277100.0592750000 +-1595048.5752310001 -277303.4578030000 +-1594467.6088459999 -277710.2548570000 +-1594177.1256530001 -277913.6533850000 +-1593886.6424600000 -278117.0519120000 +-1593392.7607410001 -278029.9672470000 +-1593102.2775480000 -278233.3657740000 +-1592898.8790209999 -277942.8825810000 +-1592695.4804940000 -277652.3993880000 +-1592201.5987740001 -277565.3147220000 +-1591998.2002470000 -277274.8315290000 +-1592041.7425800001 -277027.8906690000 +-1591838.3440530000 -276737.4074760000 +-1591881.8863860001 -276490.4666160000 +-1591722.0301910001 -275953.0425630000 +-1591765.5725240000 -275706.1017030000 +-1591562.1739970001 -275415.6185100000 +-1591605.7163300000 -275168.6776500000 +-1591402.3178030001 -274878.1944570000 +-1591198.9192760000 -274587.7112640000 +-1591242.4616090001 -274340.7704040000 +-1591039.0630820000 -274050.2872110000 +-1591082.6054150001 -273803.3463510000 +-1590922.7492209999 -273265.9222970000 +-1590966.2915530000 -273018.9814370000 +-1590762.8930259999 -272728.4982440000 +-1590806.4353590000 -272481.5573840000 +-1590603.0368319999 -272191.0741910000 +-1590399.6383050000 -271900.5909980000 +-1590443.1806379999 -271653.6501380000 +-1590239.7821110000 -271363.1669450000 +-1590036.3835839999 -271072.6837520000 +-1590123.4682499999 -270578.8020320000 +-1589920.0697230001 -270288.3188390000 +-1589716.6711960000 -269997.8356460000 +-1589760.2135280001 -269750.8947860000 +-1589556.8150010000 -269460.4115930000 +-1589600.3573340001 -269213.4707330000 +-1589440.5011400001 -268676.0466800000 +-1589484.0434729999 -268429.1058190000 +-1589280.6449460001 -268138.6226260000 +-1589077.2464190000 -267848.1394340000 +-1588873.8478920001 -267557.6562410000 +-1588423.5085050000 -267223.6307150000 +-1588220.1099779999 -266933.1475220000 +-1587973.1691180000 -266889.6051890000 +-1587726.2282580000 -266846.0628560000 +-1587479.2873980000 -266802.5205230000 +-1586985.4056780001 -266715.4358580000 +-1586738.4648190001 -266671.8935250000 +-1586491.5239589999 -266628.3511920000 +-1586244.5830989999 -266584.8088590000 +-1585997.6422389999 -266541.2665270000 +-1585750.7013790000 -266497.7241940000 +-1585256.8196590000 -266410.6395280000 +-1585009.8787990001 -266367.0971960000 +-1584762.9379400001 -266323.5548630000 +-1584472.4547470000 -266526.9533900000 +-1584428.9124139999 -266773.8942500000 +-1584138.4292210001 -266977.2927780000 +-1583804.4036960001 -267427.6321650000 +-1583760.8613630000 -267674.5730250000 +-1583470.3781699999 -267877.9715530000 +-1583179.8949770001 -268081.3700800000 +-1583136.3526450000 -268328.3109400000 +-1582845.8694519999 -268531.7094680000 +-1582802.3271190000 -268778.6503280000 +-1582468.3015940001 -269228.9897150000 +-1582424.7592610000 -269475.9305750000 +-1582134.2760679999 -269679.3291030000 +-1582090.7337350000 -269926.2699630000 +-1581800.2505419999 -270129.6684900000 +-1581509.7673500001 -270333.0670180000 +-1581219.2841570000 -270536.4655450000 +-1581132.1994910000 -271030.3472650000 +-1580841.7162990000 -271233.7457930000 +-1580551.2331060001 -271437.1443200000 +-1580507.6907730000 -271684.0851800000 +-1580217.2075799999 -271887.4837080000 +-1580173.6652470001 -272134.4245680000 +-1579839.6397220001 -272584.7639550000 +-1579796.0973890000 -272831.7048160000 +-1579505.6141959999 -273035.1033430000 +-1579462.0718640001 -273282.0442030000 +-1579418.5295309999 -273528.9850630000 +-1579128.0463380001 -273732.3835910000 +-1579084.5040050000 -273979.3244510000 +-1578750.4784800000 -274429.6638380000 +-1578706.9361469999 -274676.6046980000 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux3.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux3.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux3.exp (revision 18013) @@ -0,0 +1,10 @@ +## Name:MassFlux4 +## Icon:0 +# Points Count Value +4 1.000000 +# X pos Y pos +-1605390.0162510793 -209047.7381093690 +-1592454.5133393104 -208141.0538114636 +-1574086.0992045982 -208918.2117810968 +-1560891.8862345938 -210472.5277203633 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux4.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux4.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux4.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux5 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1660144.8928561222 -319887.0521198768 +-1630434.8468677206 -298387.4931837295 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux5.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux5.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux5.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux6 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1635925.9455644884 -157610.4872330402 +-1625288.0602165582 -120861.2411445094 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux6.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux6.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux6.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux7 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1642052.5574172530 -95522.7999921589 +-1610982.9267521810 -78949.7704229522 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux7.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux7.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux7.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux8 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1605014.5365509845 -81707.1563440522 +-1566756.1417461073 -82763.2940976564 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux8.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux8.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux8.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux9 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1555378.1885604258 -45179.9925785796 +-1550583.9034971579 -68415.0231578706 + Index: /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux9.exp =================================================================== --- /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux9.exp (revision 18013) +++ /issm/workshop/07UQ_PIG/Exp_Par/MassFluxes/MassFlux9.exp (revision 18013) @@ -0,0 +1,8 @@ +## Name:MassFlux10 +## Icon:0 +# Points Count Value +2 1.000000 +# X pos Y pos +-1565010.6613481115 -111795.8673475993 +-1537474.3225972611 -138442.9656094175 + Index: /issm/workshop/07UQ_PIG/runme.m =================================================================== --- /issm/workshop/07UQ_PIG/runme.m (revision 18013) +++ /issm/workshop/07UQ_PIG/runme.m (revision 18013) @@ -0,0 +1,290 @@ +%PIG Uncertainty Quantification Application +steps=[1]; + +if any(steps==1) + disp(' Step 1: plot flux gates'); + + md = loadmodel('../06PIG_Application/Models/PIG.Control_drag_default'); + + texts=cell(1,13); + textpositions=cell(1,13); + + for i=1:13, + contour=expread(['./Exp_Par/MassFluxes/MassFlux' num2str(i) '.exp']); + textpositions{i}=[contour.x(end) contour.y(end)]; + end + + plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'log',10,'expdisp',... + {'Exp_Par/MassFluxes/MassFlux1.exp','Exp_Par/MassFluxes/MassFlux2.exp',... + 'Exp_Par/MassFluxes/MassFlux3.exp','Exp_Par/MassFluxes/MassFlux4.exp',... + 'Exp_Par/MassFluxes/MassFlux5.exp','Exp_Par/MassFluxes/MassFlux6.exp',... + 'Exp_Par/MassFluxes/MassFlux7.exp','Exp_Par/MassFluxes/MassFlux8.exp',... + 'Exp_Par/MassFluxes/MassFlux9.exp','Exp_Par/MassFluxes/MassFlux10.exp',... + 'Exp_Par/MassFluxes/MassFlux11.exp','Exp_Par/MassFluxes/MassFlux12.exp',... + 'Exp_Par/MassFluxes/MassFlux13.exp'},... + 'expstyle',{'k-','k-','k-','k-','k-','k-','k-',... + 'k-','k-','k-','k-','k-','k-'},'linewidth',2,... + 'text',{'1','2','3','4','5','6','7',... + '8','9','10','11','12','13'},... + 'textposition',textpositions); + +end +if any(steps==2) + disp(' Step 2: compute cross overs from CRESIS'); + + md = loadmodel('../06PIG_Application/Models/PIG.Control_drag_default'); + + %load cross overs: CRESIS McCord Antarctica, 2009 (courtesy of John Paden) + load('./Data/CrossOvers2009.mat'); + + %interpolate cross over errors over our mesh vertices + DeltaHH=InterpFromMeshToMesh2d(index,x,y,dhh,md.mesh.x,md.mesh.y); + + %avoid NaN values + pos=find(isnan(DeltaHH)); DeltaHH(pos)=0; + + %filter out unrealistic error ranges + flags=ContourToNodes(md.mesh.x,md.mesh.y,'Exp_Par/ErrorContour.exp',1); + pos=find(~flags); DeltaHH(pos)=0; + + %avoid large unrealistic values + pos=find(DeltaHH>1); DeltaHH(pos)=1; + pos=find(DeltaHH<-1); DeltaHH(pos)=-1; + + %transform into absolute errors and setup a minimum error everywhere. + DeltaHH=abs(DeltaHH); + pos=find(DeltaHH==0); pos2=find(DeltaHH~=0); + DeltaHH(pos)=min(DeltaHH(pos2)); + + save Models/PIG.CrossOvers DeltaHH + +end +if any(steps==3) + disp(' Step 3: sampling analysis'); + + %load model and cross over errors + md = loadmodel('../06PIG_Application/Models/PIG.Control_drag_default'); + load -mat Models/PIG.CrossOvers + + %partition the mesh + md.qmu.numberofpartitions=50; + md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,... + 'weighting','on'); + md.qmu.partition=md.qmu.partition-1; %switch partition to c-indexing + + %make DeltaHH into our 3 sigma deviation + DeltaHH=DeltaHH/6; %2 (to transform DeltaHH into a radius) x 3 (for 3 sigma) + DeltaHH_on_partition=AreaAverageOntoPartition(md,DeltaHH); + DeltaHH_on_grids=DeltaHH_on_partition(md.qmu.partition+1); %just to check in case + + md.qmu.variables.thickness=normal_uncertain('scaled_Thickness',1,1); + md.qmu.variables.thickness.stddev=DeltaHH_on_partition; + + %responses + md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); %grounding line + md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux6=response_function('indexed_MassFlux_6',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux7=response_function('indexed_MassFlux_7',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux8=response_function('indexed_MassFlux_8',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux9=response_function('indexed_MassFlux_9',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux10=response_function('indexed_MassFlux_10',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux11=response_function('indexed_MassFlux_11',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux12=response_function('indexed_MassFlux_12',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux13=response_function('indexed_MassFlux_13',[],... + [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + + %mass flux profiles + md.qmu.mass_flux_profiles={'MassFlux1.exp',... + 'MassFlux2.exp',... + 'MassFlux3.exp',... + 'MassFlux4.exp',... + 'MassFlux5.exp',... + 'MassFlux6.exp',... + 'MassFlux7.exp',... + 'MassFlux8.exp',... + 'MassFlux9.exp',... + 'MassFlux10.exp',... + 'MassFlux11.exp',... + 'MassFlux12.exp',... + 'MassFlux13.exp'}; + md.qmu.mass_flux_profile_directory='../Exp_Par/MassFluxes/'; + + %% sampling analysis + md.qmu.method =dakota_method('nond_samp'); + md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),... + 'seed',1234,... + 'samples',30,... + 'sample_type','lhs'); %random or lhs + + %% a variety of parameters + md.qmu.params.evaluation_concurrency=1; + md.qmu.params.analysis_driver=''; + md.qmu.params.analysis_components=''; + md.qmu.params.tabular_graphics_data=true; + md.verbose=verbose('qmu',true); + + md.stressbalance.restol=10^-5; %tighten tolerances for UQ analyses + + %solve + md.qmu.isdakota=1; md.inversion.iscontrol=0; + md.cluster=generic('name',oshostname,'np',2); + md=solve(md,StressbalanceSolutionEnum,'overwrite','y'); + + save ./Models/PIG.Sampling md; +end +if any(steps==4) + disp(' Step 4: sensitivity analysis'); + + %load model + md = loadmodel('../06PIG_Application/Models/PIG.Control_drag_default'); + + %partition the mesh + md.qmu.numberofpartitions=10; + md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); + md.qmu.partition=md.qmu.partition-1; %switch partition to c-indexing + + %all types of variables and responses: scaled_Thickness, indexed_MassFlux_i,MaxVel,nodal_DragCoefficient_i. scaled variables are expanded. + + %variables + md.qmu.variables.DragCoefficient=normal_uncertain(... + 'scaled_FrictionCoefficient',1,0.05); + md.qmu.variables.rheology_B=normal_uncertain(... + 'scaled_MaterialsRheologyB',1,0.05); + md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05); + + %responses + md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux7=response_function('indexed_MassFlux_7',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux8=response_function('indexed_MassFlux_8',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux9=response_function('indexed_MassFlux_9',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux10=response_function('indexed_MassFlux_10',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux11=response_function('indexed_MassFlux_11',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux12=response_function('indexed_MassFlux_12',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + md.qmu.responses.MassFlux13=response_function('indexed_MassFlux_13',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); + + %mass flux profiles + md.qmu.mass_flux_profiles={'MassFlux1.exp',... + 'MassFlux2.exp',... + 'MassFlux3.exp',... + 'MassFlux4.exp',... + 'MassFlux5.exp',... + 'MassFlux6.exp',... + 'MassFlux7.exp',... + 'MassFlux8.exp',... + 'MassFlux9.exp',... + 'MassFlux10.exp',... + 'MassFlux11.exp',... + 'MassFlux12.exp',... + 'MassFlux13.exp'}; + md.qmu.mass_flux_profile_directory='../Exp_Par/MassFluxes/'; + + %method: local reliability + md.qmu.method =dakota_method('nond_l'); + md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),... + 'output','quiet'); + + %parameters + md.qmu.params.evaluation_concurrency=1; + md.qmu.params.analysis_driver=''; + md.qmu.params.analysis_components=''; + md.qmu.params.tabular_graphics_data=false; + + md.stressbalance.restol=10^-5; %tighten for qmu analyses + + %solve + md.qmu.isdakota=1; + md.inversion.iscontrol=0; + md.cluster=generic('name',oshostname,'np',2); + md.verbose=verbose('qmu',true); + md=solve(md,StressbalanceSolutionEnum,'overwrite','y'); + + save ./Models/PIG.Sensitivity md; +end +if any(steps==5) + disp(' Step 5: plot partition'); + + %load model + md = loadmodel('./Models/PIG.Sampling'); + + plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,... + 'linewidth',2, 'axis#all','image','unit','km','colorbar','off',... + 'title','','meshcolor','b','grid','on'); + +end +if any(steps==6) + disp(' Step 6: plot histogram'); + + %load model + md = loadmodel('./Models/PIG.Sampling'); + + %which profile are we looking at? + index=1; + + %retrieve results for the specific profile, mass flux in m^3 water equiv/s + result=md.results.dakota.dresp_dat(md.qmu.numberofpartitions+index); + result.sample=result.sample/1e12*60*60*24*365; + + %plot histogram + plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off',... + 'xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor',... + 'none','EdgeColor','red'); + +end +if any(steps==7) + disp(' Step 7: plot sensitivity'); + + %load model + md = loadmodel('./Models/PIG.Sensitivity'); + %copy dakota results into model qmu + md.qmu.results=md.results.dakota; + + %which profile are we looking at? + index=1; + + %To plot sensitivities + sa=md.results.dakota.dresp_out(index).sens(1:10); sa=sa(md.qmu.partition+1)/1e12*60*60*24*365; + sb=md.results.dakota.dresp_out(index).sens(11:20); sb=sb(md.qmu.partition+1)/1e12*60*60*24*365; + sh=md.results.dakota.dresp_out(index).sens(21:30); sh=sh(md.qmu.partition+1)/1e12*60*60*24*365; + + plotmodel(md,'data',sh,'data',sa,'data',sb,'expdisp#all',... + ['Exp_Par/MassFluxes/MassFlux' num2str(index) '.exp'],... + 'expstyle#all','b-','linewidth#all',2,... + 'nlines',3,'ncols',1, 'axis#all','image',... + 'colorbar#all','on','colorbarfontsize#all',10,... + 'colorbartitle#1','S_{H}', 'colorbartitle#2','S_{\alpha}',... + 'colorbartitle#3','S_{B}','unit#all','km','figure',1); + + %To plot importance factors + ifa=importancefactors(md,'scaled_FrictionCoefficient',['indexed_MassFlux_' num2str(index)]); + ifb=importancefactors(md,'scaled_MaterialsRheologyB',['indexed_MassFlux_' num2str(index)]); + ifh=importancefactors(md,'scaled_Thickness',['indexed_MassFlux_' num2str(index)]); + + plotmodel(md,'data',ifh,'data',ifa,'data',ifb,'expdisp#all',... + ['Exp_Par/MassFluxes/MassFlux' num2str(index) '.exp'],... + 'expstyle#all','b-','linewidth#all',2,'log#all',10,... + 'nlines',3,'ncols',1, 'axis#all','image','caxis#all',[1e-10 1],... + 'colorbar#all','on','colorbarfontsize#all',10,... + 'colorbartitle#1','If_{H}', 'colorbartitle#2','If_{\alpha}',... + 'colorbartitle#3','If_{B}','unit#all','km','figure',2); + +end Index: /issm/workshop/09SeaRISE_GrIS/Exp_Par/DomainOutline.exp =================================================================== --- /issm/workshop/09SeaRISE_GrIS/Exp_Par/DomainOutline.exp (revision 18013) +++ /issm/workshop/09SeaRISE_GrIS/Exp_Par/DomainOutline.exp (revision 18013) @@ -0,0 +1,733 @@ +## Name:DomainOutline +## Icon:0 +# Points Count Value +727 1.000000 +# X pos Y pos +47662.8073760000 -2621047.2585999998 +39869.9055420000 -2638077.6292949999 +26037.5047870000 -2644088.3483640002 +24868.5695120000 -2659115.1460360000 +25647.8596950000 -2674141.9437079998 +16880.8451320000 -2687666.0616130000 +515.7512810000 -2685662.4885900002 +-7666.7956450000 -2672639.2639410002 +-22083.6640370000 -2677147.3032430001 +-30071.3884170000 -2690170.5278920000 +-44683.0793550000 -2694678.5671930001 +-41760.7411680000 -2709705.3648660001 +-37669.4677050000 -2724231.2692820001 +-44488.2568100000 -2738757.1736989999 +-55787.9644690000 -2753283.0781149999 +-67672.1397650000 -2742263.4264890002 +-87154.3943500000 -2741762.5332330000 +-105078.0685680000 -2741011.1933490001 +-97090.3441880000 -2760546.0303230002 +-99817.8598300000 -2780581.7605530000 +-82478.6532500000 -2789096.9459000002 +-77218.4445120000 -2804123.7435730002 +-77023.6219660000 -2820152.3277560002 +-95336.9412760000 -2829669.2996149999 +-112870.9704020000 -2832674.6591500002 +-123196.5653320000 -2818649.6479890002 +-122417.2751480000 -2834177.3389170002 +-120274.2271440000 -2849204.1365890000 +-105662.5362060000 -2854713.9624020001 +-89492.2649000000 -2850706.8163560000 +-77997.7346950000 -2860724.6814710000 +-92024.9579960000 -2867737.1870519998 +-84816.5238000000 -2890277.3835600000 +-94752.4736380000 -2901797.9284419999 +-114234.7282230000 -2894284.5296060001 +-127677.4838860000 -2883765.7712349999 +-146964.9159250000 -2888774.7037930000 +-139951.3042750000 -2905805.0744880000 +-154173.3501210000 -2916824.7261140002 +-165667.8803260000 -2904302.3947210000 +-162550.7195930000 -2920330.9789040000 +-152809.5923000000 -2934355.9900650000 +-147744.2061080000 -2949883.6809930000 +-163719.6548680000 -2950384.5742489998 +-175603.8301650000 -2963908.6921540000 +-189825.8760110000 -2956896.1865730002 +-195670.5523870000 -2971422.0909899999 +-210282.2433250000 -2966413.1584319999 +-222750.8862590000 -2978301.3022220000 +-214763.1618800000 -2995247.5385030000 +-222750.8862590000 -3011386.8111520000 +-206970.2600460000 -3016228.5929470002 +-217295.8549760000 -3027122.6019839998 +-232686.8360980000 -3031964.3837790000 +-212035.6462380000 -3035192.2383090002 +-212814.9364210000 -3050928.0291410000 +-197229.1327540000 -3048910.6200600001 +-185344.9574570000 -3058190.7018329999 +-174824.5399810000 -3071102.1199520002 +-182812.2643610000 -3084820.5017030002 +-192163.7465610000 -3098942.3652710002 +-176577.9428940000 -3092486.6562120002 +-184955.3123650000 -3105801.5561469998 +-182617.4418150000 -3121537.3469790001 +-199372.1807580000 -3120326.9015310002 +-218464.7902510000 -3115888.6015519998 +-205022.0345870000 -3126782.6105900002 +-189825.8760110000 -3128800.0196710001 +-198592.8905740000 -3146149.7377680000 +-213009.7589670000 -3140904.4741580002 +-225283.5793550000 -3131220.9105679998 +-222750.8862590000 -3146149.7377680000 +-208334.0178670000 -3151395.0013790000 +-195670.5523870000 -3161078.5649680002 +-213009.7589670000 -3171569.0921900002 +-227621.4499060000 -3161078.5649680002 +-242622.7859360000 -3157850.7104389998 +-240090.0928400000 -3174393.4649040001 +-220023.3706180000 -3185690.9557579998 +-233076.4811890000 -3193357.1102660000 +-258013.7670580000 -3188111.8466550000 +-275158.1510920000 -3196181.4829790001 +-258208.5896040000 -3198198.8920600000 +-252753.5583200000 -3213127.7192600002 +-270482.4099920000 -3215952.0919740000 +-258793.0572410000 -3226442.6191960000 +-275547.7961840000 -3232969.1157640000 +-276521.9089130000 -3248275.5695210001 +-293276.6478560000 -3241868.2167850002 +-282561.4078340000 -3256818.7065010001 +-288211.2616640000 -3273549.0164200002 +-302433.3075110000 -3267141.6636850000 +-305745.2907900000 -3251479.2458879999 +-313733.0151700000 -3237596.6482950002 +-310810.6769820000 -3218018.6260489998 +-325812.0130120000 -3213747.0575589999 +-338280.6559470000 -3203780.0644160002 +-336527.2530340000 -3187761.6825779998 +-321720.7395500000 -3185269.9342919998 +-325812.0130120000 -3165691.9120459999 +-313538.1926240000 -3155368.9548610002 +-329708.4639290000 -3158928.5952699999 +-338085.8334010000 -3173167.1569030001 +-351723.4116100000 -3164268.0558819999 +-342566.7519550000 -3151809.3144530002 +-325812.0130120000 -3153233.1706160000 +-310615.8544360000 -3148249.6740450002 +-329123.9962920000 -3150029.4942490002 +-343151.2195930000 -3144334.0695949998 +-361854.1839940000 -3145401.9617180000 +-377634.8102080000 -3142910.2134320000 +-388350.0502290000 -3154657.0267800000 +-397701.5324300000 -3178150.6534750001 +-417573.4321060000 -3179930.4736790000 +-432720.8850460000 -3186308.0579260001 +-451813.4945390000 -3187438.8412219998 +-467788.9432980000 -3177261.7915590000 +-482985.1018740000 -3173492.5139070000 +-490583.1811620000 -3158792.3310610000 +-492726.2291670000 -3141453.6538590002 +-490193.5360710000 -3119214.9157079998 +-505389.6946470000 -3113560.9992289999 +-489803.8909790000 -3110922.5048719998 +-478309.3607740000 -3100368.5274439999 +-495648.5673540000 -3101122.3829750000 +-512792.9513890000 -3108660.9382799999 +-501103.5986380000 -3097353.1053220001 +-493895.1644420000 -3083783.7057719999 +-486491.9077000000 -3068329.6673960001 +-503636.2917340000 -3066068.1008040002 +-522728.9012270000 -3064183.4619780001 +-521559.9659520000 -3046467.8570099999 +-531301.0932440000 -3029506.1075729998 +-543769.7361790000 -3012167.4303710000 +-545912.7841830000 -2992944.1143419999 +-550004.0576459999 -2974474.6538439998 +-567538.0867720000 -2983897.8479749998 +-585851.4060810000 -2995959.5364640001 +-589942.6795440000 -2981259.3536189999 +-572213.8278720001 -2975982.3649050002 +-565200.2162220000 -2959397.5432330002 +-548055.8321870000 -2954120.5545190000 +-545717.9616370000 -2933389.5274290000 +-536561.3019820000 -2920197.0556450002 +-545523.1390910000 -2907758.4393910002 +-533444.1412490000 -2898335.2452590000 +-530568.0950890000 -2883304.2548420001 +-545624.1658510000 -2879383.1998239998 +-559750.8495300000 -2868110.1666489998 +-577223.3267109999 -2867129.9028940001 +-590606.5007220000 -2859287.7928590002 +-573319.9009570000 -2851445.6828239998 +-551572.2431889999 -2846054.2321750000 +-530939.8499220000 -2834781.1989989998 +-528337.5660870000 -2817626.5832969998 +-543765.3916830000 -2804883.1544900001 +-528151.6886699999 -2802432.4951040000 +-528151.6886699999 -2785768.0112790000 +-513467.3727410000 -2778416.0331210000 +-513839.1275750000 -2760771.2855420001 +-525735.2822510001 -2749498.2523670001 +-509563.9469880000 -2738225.2191909999 +-517184.9210780000 -2724501.5266300002 +-532240.9918400000 -2725971.9222610001 +-524248.2629170000 -2707837.0428050002 +-512537.9856570000 -2694113.3502429998 +-521088.3468310000 -2678429.1301730000 +-538189.0691780000 -2669606.7563829999 +-556033.3011930001 -2677448.8664190001 +-539676.0885130001 -2666175.8332429999 +-526664.6693360000 -2658333.7232080000 +-517556.6759110000 -2639708.7118739998 +-526107.0370850000 -2626965.2830670001 +-538932.5788460000 -2618142.9092779998 +-550456.9786890000 -2607360.0079790000 +-565698.9268680000 -2599517.8979440001 +-581126.7524640000 -2600988.2935759998 +-583110.5023000001 -2578346.1085480000 +-572583.8192940000 -2566324.6704099998 +-556315.3091930000 -2558811.2715739999 +-538132.8567270000 -2571333.6029679999 +-524735.2601730000 -2578346.1085480000 +-519759.0100240000 -2594875.5859869998 +-503107.7114500000 -2585859.5073839999 +-486456.4128760000 -2576342.5355250002 +-477843.6722340000 -2559312.1648300001 +-474015.7875040000 -2536771.9683220000 +-473441.6047950000 -2521244.2773940000 +-474972.7586870000 -2504714.7999539999 +-485690.8359300000 -2491691.5753049999 +-486073.6244030000 -2476664.7776330002 +-492581.0284430000 -2459634.4069380001 +-474207.1817410000 -2447112.0755440001 +-465594.4410990000 -2434589.7441509999 +-462340.7390790000 -2418060.2667109999 +-462914.9217880000 -2399026.3229930000 +-468465.3546460000 -2375484.3399729999 +-473058.8163220000 -2360457.5423010001 +-478417.8549440000 -2343427.1716060000 +-504256.0768690000 -2331405.7334679998 +-489901.5091330000 -2322389.6548649999 +-475164.1529230000 -2307863.7504480002 +-462577.4401970000 -2287577.5735909999 +-457317.2314590000 -2269044.5231289999 +-439783.2023330000 -2268042.7366169998 +-436666.0415990000 -2252014.1524330000 +-435886.7514160000 -2233481.1019710000 +-448550.2168960000 -2220457.8773220000 +-440172.8474240000 -2206933.7594170002 +-429262.7848570000 -2193409.6415120000 +-433743.7034120000 -2176880.1640720000 +-426730.0917610000 -2162354.2596560000 +-419716.4801110000 -2147327.4619840002 +-428678.3172200000 -2134805.1305900002 +-427314.5593990000 -2110762.2543150000 +-426535.2692150000 -2094733.6701310000 +-427509.3819440000 -2067184.5410650000 +-428093.8495820000 -2050655.0636260000 +-425366.3339400000 -2035628.2659539999 +-432379.9455910000 -2021102.3615369999 +-446017.5238000000 -2003571.0975860001 +-430821.3652240000 -2002569.3110750001 +-428678.3172200000 -1986039.8336360001 +-424197.3986650000 -1971013.0359630000 +-418060.4884710000 -1957238.4714309999 +-433451.4695930000 -1950225.9658500000 +-426243.0353960000 -1933696.4884110000 +-441634.0165180000 -1943714.3535259999 +-433841.1146840000 -1926683.9828300001 +-448647.6281690000 -1933696.4884110000 +-464428.2543820000 -1948222.3928270000 +-465402.3671120000 -1933195.5951550000 +-463843.7867450000 -1914662.5446929999 +-456245.7074570000 -1901138.4267879999 +-468714.3503910000 -1892122.3481840000 +-480208.8805960000 -1905646.4660890000 +-479429.5904130000 -1888616.0953939999 +-478455.4776830000 -1871585.7246989999 +-478650.3002290000 -1856558.9270270001 +-484494.9766050000 -1841031.2360990001 +-494430.9264430000 -1822498.1856370000 +-505925.4566480000 -1807972.2812200000 +-508847.7948350000 -1791442.8037810000 +-504172.0537350000 -1776416.0061079999 +-510990.8428400000 -1759886.5286689999 +-502223.8282770000 -1745861.5175079999 +-512939.0682980000 -1728330.2535570001 +-511964.9555690000 -1712802.5626300001 +-505340.9890100000 -1697775.7649570000 +-506315.1017390000 -1681747.1807740000 +-507289.2144690000 -1662713.2370559999 +-508446.9486300000 -1645056.7497910000 +-502180.9448550000 -1623518.3397939999 +-505662.0580630000 -1602480.8230530000 +-505836.1137240000 -1586953.1321250000 +-509317.2269320000 -1571425.4411970000 +-516105.3976870000 -1557400.4300360000 +-530552.0175000000 -1537865.5930620001 +-533336.9080660000 -1513321.8235309999 +-522719.5127820000 -1498295.0258589999 +-506532.3363650000 -1492284.3067900001 +-522893.5684430000 -1484770.9079539999 +-515757.2863660000 -1469744.1102819999 +-532118.5184440000 -1465236.0709800001 +-535947.7429720000 -1449708.3800520001 +-544476.4703320001 -1428670.8633109999 +-551960.8637290000 -1415146.7454059999 +-565363.1495800000 -1405629.7735470000 +-567277.7618440000 -1386094.9365729999 +-582420.6042990000 -1384091.3635499999 +-593212.0552430000 -1370066.3523899999 +-606266.2297730000 -1362052.0602980000 +-613793.7528780000 -1344395.5730330001 +-621618.1005840000 -1359923.2639609999 +-634079.0987820000 -1369440.2358200001 +-633354.6221430000 -1354413.4381480000 +-638860.6446030000 -1339386.6404760000 +-649438.0035390001 -1351908.9718690000 +-667405.0241970000 -1352409.8651250000 +-679431.3364110000 -1363429.5167510000 +-679721.1270670000 -1346900.0393119999 +-676533.4298540000 -1330370.5618720001 +-690008.6953470000 -1339386.6404760000 +-701165.6355950000 -1327365.2023380001 +-688994.4280520000 -1315844.6574560001 +-709279.7739560000 -1318349.1237339999 +-717249.0169900000 -1331372.3483839999 +-735216.0376480001 -1329869.6686160001 +-741420.6668040000 -1315969.8807699999 +-750947.3474350000 -1299440.4033299999 +-734597.5036490000 -1299440.4033299999 +-722496.0444700000 -1284413.6056580001 +-707948.5456690000 -1288420.7517039999 +-689023.9233340001 -1290424.3247270000 +-678209.8534290000 -1276399.3135660000 +-672545.3406220000 -1262374.3024060000 +-693401.0468670001 -1264377.8754280000 +-709493.4127980000 -1260871.6226379999 +-725971.9955110000 -1250853.7575240000 +-725714.5176560000 -1234825.1733400000 +-720951.1773400001 -1215290.3363659999 +-737172.2821980000 -1218295.6958999999 +-749402.4803050000 -1205272.4712510000 +-736399.8486340001 -1196757.2859040000 +-720822.4384130000 -1191247.4600899999 +-710652.0631450000 -1204771.5779949999 +-697906.9093280000 -1218796.5891559999 +-684131.8440910000 -1227812.6677590001 +-668940.6506530000 -1229816.2407819999 +-666237.1331770000 -1253358.2238020001 +-653620.7182870000 -1262374.3024060000 +-648471.1611900000 -1247347.5047330000 +-634696.0959530000 -1258868.0496149999 +-632625.0227419999 -1243541.0164699999 +-618416.8083510000 -1233878.1370890001 +-601000.2874830000 -1228080.4094600000 +-615361.2783740000 -1217451.2421410000 +-598861.4165000001 -1216968.0981719999 +-610625.2069099999 -1201990.6351310001 +-617652.9258560000 -1185080.5962139999 +-631708.3637490000 -1178799.7246160000 +-648666.5551200000 -1183631.1643070001 +-656458.1565610000 -1169136.8452349999 +-653708.1795820000 -1154159.3821950001 +-665166.4169950000 -1166237.9814210001 +-676166.3249110000 -1153193.0942569999 +-676166.3249110000 -1133384.1915249999 +-689305.1038110000 -1118406.7284850001 +-709013.2721609999 -1111159.5689490000 +-711610.4726409999 -1091833.8101870001 +-716346.5441050000 -1074440.6273010001 +-701527.2237180000 -1070092.3315790000 +-680291.2903790000 -1068642.8996720000 +-664096.9815030000 -1078788.9230219999 +-650958.2026030000 -1100530.4016300000 +-628041.7277770001 -1095698.9619390001 +-611236.3129050001 -1091350.6662170000 +-588930.9440750000 -1094732.6740009999 +-577319.9301640000 -1078305.7790530000 +-565250.5867559999 -1066710.3237960001 +-550431.2663680000 -1071541.7634860000 +-540042.4644469999 -1060429.4521979999 +-531203.6819710000 -1047613.0696300000 +-510357.6695660000 -1037937.0529810000 +-498473.4942690000 -1027821.2173940000 +-486199.6738810000 -1015066.4681760001 +-487953.0767930000 -994834.7970020000 +-464964.0163830000 -990876.4265550000 +-449378.2127160000 -998353.3485110000 +-430285.6032230000 -1001432.0810810000 +-417622.1377430000 -1012427.5545450000 +-399308.8184330000 -1004950.6325890000 +-384697.1274945000 -1000332.5337340001 +-370085.4365560000 -995714.4348789999 +-375540.4678400000 -980760.5909680000 +-380021.3863940000 -966246.5659950000 +-383528.1922200000 -949973.2652690000 +-383528.1922200000 -930181.4130330000 +-384112.6598570000 -911709.0176140000 +-380216.2089400000 -883560.6055450001 +-366919.5701860000 -892846.9764360000 +-353087.1694310000 -880338.0618260000 +-336527.2530340000 -873692.7009390000 +-322694.8522790000 -888547.0370390000 +-307888.3387940000 -897146.9158330000 +-294835.2282230000 -913564.8662580000 +-275937.4412760000 -919037.5164000000 +-252169.0906820000 -915128.4805850000 +-245350.3015780000 -901446.8552300000 +-244181.3663030000 -918646.6128190000 +-231615.3120950000 -930563.0925760000 +-215445.0407900000 -924017.9382920000 +-206872.8487730000 -937108.2468600000 +-190507.7549220000 -933426.5975750000 +-186611.3040050000 -912972.9904390000 +-185961.8955186667 -896882.8194913333 +-185312.4870323333 -880792.6485436666 +-184663.0785460000 -864702.4775960000 +-175311.5963460000 -883519.7961620000 +-174532.3061620000 -902337.1147280000 +-172584.0807040000 -918700.0004370000 +-168882.4523330000 -933835.6697180000 +-152322.5359360000 -933835.6697180000 +-142581.4086430000 -921563.5054359999 +-127969.7177050000 -912154.8461530000 +-131086.8784390000 -930563.0925760000 +-115501.0747710000 -924427.0104349999 +-112189.0914910000 -905609.6918690000 +-114916.6071330000 -887610.5175890001 +-101473.8514700000 -899882.6818710000 +-84329.4674350000 -896610.1047290000 +-100110.0936490000 -892110.3111590000 +-87836.2732610000 -879838.1468770000 +-74685.7514160000 -869305.2030530000 +-57541.3673810000 -863251.9975030000 +-49553.6430020000 -850388.9357080000 +-34357.4844260000 -846227.3568920000 +-47995.0626350000 -825797.7881600000 +-30071.3884170000 -835255.9218320000 +-18771.6807580000 -850767.2610550000 +-6303.0378240000 -859847.0693810000 +7139.7178400000 -852658.8877899999 +22141.0538700000 -862873.6721560000 +31297.7135250000 -877628.3606850000 +47662.8073760000 -874601.7579100000 +49611.0328340000 -893518.0252550000 +66365.7717770000 -886708.1690110000 +81123.5796250000 -882220.1151260000 +97683.4960220000 -884202.5846670000 +80149.4668960000 -888960.5115630000 +65537.7759570000 -896493.8958170000 +68460.1141450000 -914732.6155880000 +84435.5629050000 -917508.0729450000 +100605.8342100000 -907595.7252430000 +120282.9113400000 -911957.1582320000 +129049.9259030000 -927420.4206470000 +153402.7441340000 -925437.9511060000 +169767.8379850000 -920680.0242090001 +180677.9005530000 -910371.1825990001 +192367.2533040000 -899269.3531730000 +201134.2678670000 -912353.6521400000 +182626.1260110000 -920680.0242090001 +169183.3703480000 -936539.7805320000 +163338.6939720000 -950417.0673150000 +182820.9485570000 -955174.9942120000 +187301.8671120000 -970241.7627180000 +200744.6227750000 -980550.6043280000 +198796.3973170000 -995617.3728350000 +211118.9233410000 -983385.4589040000 +211508.5684330000 -999969.0227530000 +225146.1466420000 -1008484.9068910000 +240731.9503100000 -1009829.5201760000 +238588.9023060000 -1032239.7415930000 +253785.0608820000 -1029550.5150230000 +273072.4929210000 -1036273.5814480000 +281255.0398460000 -1056890.9851510001 +286320.4260380000 -1071681.7312860000 +305413.0355310000 -1086024.2729930000 +316907.5657360000 -1095884.7704169999 +326453.8704830000 -1081990.4331380001 +341260.3839670000 -1076611.9799980000 +355677.2523600000 -1064958.6648609999 +372319.0245690000 -1056006.2810410000 +387981.0391590000 -1068027.7191780000 +383779.0352450000 -1084056.3033620000 +366207.0188750000 -1085558.9831290001 +366971.0195860000 -1102088.4605690001 +352646.0062410000 -1113108.1121950001 +362960.0158500000 -1126131.3368440000 +380914.0325760000 -1133143.8424249999 +394093.0448530000 -1123626.8705660000 +409182.0589100000 -1113108.1121950001 +409182.0589100000 -1131641.1626569999 +404407.0544620000 -1158689.3984670001 +386262.0375580000 -1153179.5726540000 +374229.0263480000 -1166202.7973030000 +364488.0172730000 -1182732.2747430000 +377285.0291950000 -1191748.3533460000 +388172.0393370000 -1202768.0049719999 +392947.0437860000 -1219798.3756680000 +380341.0320420000 -1228814.4542710001 +383588.0350670000 -1251354.6507790000 +389509.0405830000 -1265880.5551960000 +389509.0405830000 -1282910.9258910001 +399632.0500130000 -1295934.1505400001 +408609.0583760000 -1313966.3077469999 +410137.0598000000 -1333000.2514650000 +413526.9681570000 -1350406.2921020000 +416526.2265730000 -1367436.6627970000 +415776.4119690000 -1385468.8200040001 +419900.3922910000 -1406005.4434890000 +430960.1577010000 -1421032.2411610000 +431897.4259560000 -1436059.0388330000 +448018.4399430000 -1436559.9320890000 +448580.8008960000 -1456094.7690630001 +438458.3037410000 -1469117.9937120001 +432459.7869090000 -1486148.3644069999 +449705.5228020000 -1494663.5497550000 +448393.3472450000 -1513196.6002170001 +463202.1856750000 -1521210.8923090000 +447830.9862920000 -1525218.0383550001 +458140.9370980000 -1539243.0495160001 +462639.8247220000 -1554269.8471880001 +473887.0437820000 -1564287.7123030000 +458703.2980510000 -1557275.2067219999 +463014.7320240000 -1574305.5774170000 +478948.2923600000 -1567293.0718370001 +484009.5409370000 -1584824.3357879999 +475011.7656880000 -1597346.6671819999 +484009.5409370000 -1614877.9311319999 +488508.4285610000 -1630906.5153160000 +470512.8780640000 -1618384.1839230000 +461702.5564670000 -1631407.4085720000 +474074.4974330000 -1641926.1669419999 +488883.3358630000 -1648938.6725230000 +495661.8324810000 -1662963.6836830000 +478923.6439970000 -1667471.7229849999 +463095.1396700000 -1663965.4701950001 +449449.8773190000 -1652945.8185690001 +455817.6664170000 -1670477.0825199999 +441626.5935720000 -1662963.6836830000 +433439.4361610000 -1648437.7792670000 +409690.0603430000 -1652315.3206390000 +394504.9542530000 -1655951.1781029999 +411243.1427370000 -1662462.7904280000 +426707.7734010000 -1676487.8015880000 +437260.1096190000 -1691514.5992610001 +455999.6032480000 -1702534.2508870000 +441626.5935720000 -1710047.6497229999 +428891.0153770000 -1692516.3857720001 +435076.8676430000 -1708544.9699560001 +419976.1106410000 -1704036.9306540000 +408696.0270980000 -1692516.3857720001 +395050.7647470000 -1701532.4643750000 +388501.0388190000 -1718562.8350710000 +378494.5130940000 -1731085.1664640000 +376129.3342870000 -1713553.9025129999 +361938.2614420000 -1720566.4080940001 +373582.2186480000 -1732086.9529760000 +358481.4616460000 -1738598.5652999999 +369397.6715270000 -1754627.1494839999 +352659.4830430000 -1751621.7899490001 +335805.3526830000 -1754501.9261700001 +338338.0457790000 -1777042.1226780000 +354897.9621760000 -1775539.4429110000 +349053.2858010000 -1792569.8136060000 +330739.9664910000 -1791568.0270950000 +344182.7221550000 -1806093.9315110000 +341455.2065130000 -1823625.1954620001 +351780.8014430000 -1809600.1843010001 +344377.5447010000 -1826129.6617409999 +356066.8974510000 -1836147.5268550001 +345546.4799760000 -1853678.7908059999 +363080.5091020000 -1853678.7908059999 +371847.5236650000 -1841156.4594129999 +370288.9432980000 -1859689.5098750000 +386848.8596950000 -1868204.6952229999 +404772.5339130000 -1856183.2570849999 +415098.1288430000 -1867703.8019669999 +428735.7070520000 -1860691.2963870000 +421137.6277640000 -1874215.4142920000 +413929.1935680000 -1888741.3187080000 +418799.7572140000 -1903768.1163800000 +428735.7070520000 -1916791.3410300000 +415098.1288430000 -1935825.2847480001 +398148.5673540000 -1926809.2061439999 +380030.0705910000 -1929814.5656790000 +364444.2669230000 -1929814.5656790000 +348468.8181630000 -1924304.7398659999 +333077.8370420000 -1919796.7005640001 +348079.1730720000 -1935825.2847480001 +332883.0144960000 -1939832.4307939999 +318855.7911950000 -1951853.8689309999 +331908.9017660000 -1964376.2003250001 +349442.9308930000 -1959868.1610230000 +358209.9454560000 -1943338.6835840000 +358404.7680020000 -1958365.4812560000 +345936.1250670000 -1971388.7059050000 +345351.6574300000 -1987918.1833450000 +330739.9664910000 -1991925.3293910001 +337753.5781420000 -2005449.4472960001 +353534.2043550000 -2007453.0203180001 +365223.5571060000 -1997435.1552040000 +385679.9244200000 -1986415.5035770000 +399317.5026290000 -1977399.4249740001 +399122.6800840000 -1993928.9024139999 +381393.8284120000 -1998436.9417150000 +399512.3251750000 -1999939.6214820000 +395615.8742580000 -2018973.5652000001 +387822.9724240000 -2032998.5763610001 +383975.2271440000 -2048275.8206610000 +398976.5631740000 -2062801.7250780000 +397223.1602620000 -2078830.3092610000 +379299.4860440000 -2080833.8822840001 +362739.5696470000 -2085842.8148419999 +374428.9223980000 -2097363.3597240001 +390599.1937030000 -2089849.9608880000 +379299.4860440000 -2103374.0787930000 +364492.9725590000 -2111388.3708839999 +355531.1354500000 -2097363.3597240001 +343257.3150620000 -2111889.2641400001 +351829.5070790000 -2127917.8483239999 +365661.9078340000 -2119402.6629760000 +357284.5383630000 -2133928.5673929998 +374039.2773060000 -2136933.9269269998 +385533.8075110000 -2147452.6852980000 +375792.6802180000 -2160976.8032030002 +390988.8387940000 -2162479.4829699998 +378520.1958600000 -2174000.0278520002 +394105.9995280000 -2164983.9492489998 +409107.3355580000 -2161978.5897140000 +421186.3334010000 -2149957.1515759998 +437941.0723440000 -2150958.9380879998 +450409.7152780000 -2171495.5615730002 +454695.8112860000 -2156969.6571570002 +471645.3727750000 -2159975.0166910002 +459566.3749330000 -2147452.6852980000 +465605.8738540000 -2132425.8876260002 +477490.0491500000 -2120404.4494880000 +485867.4186220000 -2133427.6741370000 +493264.3211700000 -2119402.6629760000 +495978.8328800000 -2102873.1855370002 +509098.9728100000 -2110386.5843730001 +523877.9810080000 -2104876.7585600000 +541371.5009160000 -2109885.6911169998 +529608.6168400000 -2096862.4664680001 +544990.8498620000 -2099366.9327469999 +561579.5325330000 -2107381.2248390000 +576509.3469370001 -2104876.7585600000 +596868.1847600000 -2102873.1855370002 +611948.8053700000 -2108883.9046060001 +597622.2157910001 -2119402.6629760000 +580430.3082950000 -2125413.3820449999 +599130.2778520000 -2126916.0618119999 +583144.8200050000 -2134930.3539040000 +568818.2304260000 -2148454.4718089998 +582843.2075930000 -2161477.6964580002 +566706.9435400000 -2162980.3762260000 +552983.5787850000 -2170994.6683180002 +543331.9815950000 -2188025.0390130002 +547554.5553660000 -2207058.9827310001 +538053.7643820000 -2193033.9715700001 +525084.4306570000 -2209062.5557539999 +516036.0582910000 -2221083.9938920001 +529910.2292520000 -2214572.3815669999 +541220.6947100000 -2225091.1399369999 +526140.0741000000 -2223087.5669140001 +526743.2989240000 -2239617.0443540001 +521615.8879170000 -2256647.4150490002 +506233.6548950000 -2255144.7352820002 +508194.1355740000 -2275681.3587670000 +489343.3598110000 -2268167.9599310001 +488740.1349870000 -2283695.6508590002 +475770.8012620000 -2293212.6227179999 +460539.3744460000 -2281692.0778359999 +449228.9089890000 -2296217.9822519999 +448022.4593400000 -2314250.1394590000 +440180.5366230000 -2298722.4485309999 +441085.3738600000 -2277184.0385340001 +426000.8914760000 -2286450.5637659999 +431125.9775380000 -2301477.3614380001 +426190.7094780000 -2317505.9456210001 +424672.1654600000 -2332532.7432940002 +421445.2594210000 -2313498.7995750001 +409296.9072740000 -2298472.0019029998 +411764.5413040000 -2318006.8388769999 +403222.7312010000 -2304983.6142279999 +411005.2692950000 -2289956.8165560002 +393542.0130840000 -2293463.0693460000 +378356.5729000000 -2295967.5356250000 +388416.9270220000 -2308489.8670180002 +405880.1832330000 -2317005.0523660001 +417838.7173770000 -2328525.5972480001 +402463.4591910000 -2342049.7151529999 +385000.2029800000 -2341047.9286409998 +375699.1208680000 -2327022.9174799998 +384240.9309710000 -2314500.5860870001 +368865.6727850000 -2321012.1984120002 +373421.3048400000 -2337541.6758510000 +363360.9507190000 -2353069.3667790000 +351402.4165740000 -2336539.8893400002 +340392.9724410000 -2324518.4512020000 +329573.3463100000 -2313999.6928309998 +335267.8863790000 -2329527.3837589999 +336216.9763910000 -2347058.6477100002 +337166.0664020000 -2362085.4453819999 +318374.0841750000 -2366593.4846839998 +303758.0979980000 -2362586.3386380002 +290314.2882280000 -2354697.2698599999 +292067.6911410000 -2339670.4721880001 +275702.5972900000 -2347684.7642799998 +277880.2580262116 -2321542.8592617917 +269420.4169677348 -2331176.2138948478 +260960.5759092579 -2323294.3782859840 +252500.7348507810 -2343436.8470641924 +254834.4841082919 -2357157.0794203631 +258918.5453089359 -2372336.9109633607 +266740.7601810000 -2359706.2024170002 +271754.1662252456 -2369125.7927523423 +266545.9376350000 -2379241.0393909998 +291483.2235030000 -2386253.5449720002 +275897.4198350000 -2396271.4100859999 +288755.7078610000 -2404786.5954339998 +276871.5325650000 -2415806.2470600000 +263818.4219930000 -2423820.5391520001 +259532.3259840000 -2445358.9491490000 +250570.4888750000 -2432836.6177559998 +248817.0859630000 -2452371.4547290001 +247258.5055960000 -2470403.6119360002 +229919.2990160000 -2474911.6512380000 +226022.8480990000 -2490439.3421660000 +239660.4263080000 -2499956.3140250002 +227776.2510110000 -2510975.9656509999 +220178.1717230000 -2523999.1902999999 +205761.3033310000 -2533015.2689029998 +210631.8669770000 -2560063.5047129998 +203228.6102350000 -2545036.7070410000 +194266.7731260000 -2558560.8249459998 +199624.3931360000 -2581101.0214539999 +185791.9923810000 -2568077.7968049999 +189493.6207520000 -2583104.5944770002 +186571.2825650000 -2601136.7516839998 +185207.5247440000 -2584106.3809890002 +169816.5436220000 -2595126.0326149999 +156568.6105040000 -2586610.8472679998 +134748.4853690000 -2581101.0214539999 +149749.8214000000 -2579598.3416869999 +136112.2431900000 -2592120.6730809999 +146437.8381200000 -2612657.2965660002 +130462.3893610000 -2620671.5886579999 +137281.1784650000 -2635197.4930739999 +125007.3580770000 -2621172.4819140001 +106888.8613130000 -2617666.2291230001 +91887.5252830000 -2609151.0437759999 +109031.9093180000 -2607648.3640089999 +92277.1703750000 -2603641.2179629998 +97537.3791130000 -2587111.7405229998 +88380.7194580000 -2573587.6226180000 +74743.1412490000 -2567076.0102940002 +61105.5630390000 -2577093.8754090001 +66950.2394150000 -2594124.2461040001 +57598.7572140000 -2581101.0214539999 +53117.8386600000 -2598131.3921500002 +41818.1310010000 -2586610.8472679998 +30810.6571600000 -2599133.1786610000 +44253.4128240000 -2608650.1505200001 +47662.8073760000 -2621047.2585999998 + Index: /issm/workshop/09SeaRISE_GrIS/Exp_Par/Greenland.par =================================================================== --- /issm/workshop/09SeaRISE_GrIS/Exp_Par/Greenland.par (revision 18013) +++ /issm/workshop/09SeaRISE_GrIS/Exp_Par/Greenland.par (revision 18013) @@ -0,0 +1,81 @@ +%Name and hemisphere +md.miscellaneous.name='SeaRISEgreenland'; +md.mesh.hemisphere='n'; + +disp(' Loading SeaRISE data from NetCDF'); +ncdata='Data/Greenland_5km_dev1.2.nc'; +x1 = ncread(ncdata,'x1'); +y1 = ncread(ncdata,'y1'); +usrf = ncread(ncdata,'usrf')'; +topg = ncread(ncdata,'topg')'; +velx = ncread(ncdata,'surfvelx')'; +vely = ncread(ncdata,'surfvely')'; +temp = ncread(ncdata,'airtemp2m')'; +smb = ncread(ncdata,'smb')'; +gflux = ncread(ncdata,'bheatflx')'; + +disp(' Interpolating surface and bedrock'); +md.geometry.base = InterpFromGridToMesh(x1,y1,topg,md.mesh.x,md.mesh.y,0); +md.geometry.surface = InterpFromGridToMesh(x1,y1,usrf,md.mesh.x,md.mesh.y,0); + +disp(' Constructing thickness'); +md.geometry.thickness=md.geometry.surface-md.geometry.base; + +%Set min thickness to 1 meter +pos0=find(md.geometry.thickness<=0); +md.geometry.thickness(pos0)=1; +md.geometry.surface=md.geometry.thickness+md.geometry.base; + +disp(' Interpolating velocities '); +md.inversion.vx_obs = InterpFromGridToMesh(x1,y1,velx,md.mesh.x,md.mesh.y,0); +md.inversion.vy_obs = InterpFromGridToMesh(x1,y1,vely,md.mesh.x,md.mesh.y,0); +md.inversion.vel_obs = sqrt(md.inversion.vx_obs.^2+md.inversion.vy_obs.^2); +md.initialization.vx = md.inversion.vx_obs; +md.initialization.vy = md.inversion.vy_obs; +md.initialization.vz = zeros(md.mesh.numberofvertices,1); +md.initialization.vel= md.inversion.vel_obs; + +disp(' Interpolating temperatures'); +md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15; + +disp(' Interpolating surface mass balance'); +md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0); +md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice; + +disp(' Construct basal friction parameters'); +md.friction.coefficient=30*ones(md.mesh.numberofvertices,1); +pos=find(md.mask.groundedice_levelset<0); +md.friction.coefficient(pos)=0; %no friction applied on floating ice +md.friction.p=ones(md.mesh.numberofelements,1); +md.friction.q=ones(md.mesh.numberofelements,1); + +disp(' Construct ice rheological properties'); +md.materials.rheology_n=3*ones(md.mesh.numberofelements,1); +md.materials.rheology_B=paterson(md.initialization.temperature); +md.friction.q=ones(md.mesh.numberofelements,1); +md.friction.p=ones(md.mesh.numberofelements,1); + +disp(' Set other boundary conditions'); +md.mask.ice_levelset(md.mesh.vertexonboundary==1)=0; +md.basalforcings.melting_rate = zeros(md.mesh.numberofvertices,1); +md.thermal.spctemperature = [md.initialization.temperature;1]; %impose observed temperature on surface +md.masstransport.spcthickness = NaN*ones(md.mesh.numberofvertices,1); + +disp(' Set geothermal heat flux'); +md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0); + +disp(' Set Pressure'); +md.initialization.pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness; + +disp(' Single point constraint for continental model'); +%Set at least one vertex to velocity 0 so as to not get a singular problem (point on the wet peninsula) +md.stressbalance.referential=NaN*ones(md.mesh.numberofvertices,6); +md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices,1); +md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices,1); +md.stressbalance.spcvz = NaN*ones(md.mesh.numberofvertices,1); +location = 1.0e+06 *[.32011 -2.2039]; +[dist pos]=min(sqrt((md.mesh.x - location(1)).^2 + (md.mesh.y - location(2)).^2)); +md.stressbalance.spcvx(pos) = 0; +md.stressbalance.spcvy(pos) = 0; +md.stressbalance.spcvz(pos) = 0; + Index: /issm/workshop/09SeaRISE_GrIS/Exp_Par/Jak_grounded.exp =================================================================== --- /issm/workshop/09SeaRISE_GrIS/Exp_Par/Jak_grounded.exp (revision 18013) +++ /issm/workshop/09SeaRISE_GrIS/Exp_Par/Jak_grounded.exp (revision 18013) @@ -0,0 +1,19 @@ +## Name:Jak_grouded +## Icon:0 +# Points Count Value +13 1.000000 +# X pos Y pos +-414397.6663804040 -2243466.6708395295 +-421486.9037621063 -2215731.3059818177 +-342860.8164377716 -2178965.8223332223 +-275835.2993744044 -2209281.2211311869 +-226210.6377024883 -2262816.9253914217 +-228144.0660793162 -2295067.3496445753 +-333193.6745536321 -2328607.7908678548 +-379595.9555975018 -2325382.7484425395 +-422775.8560133249 -2285392.2223686292 +-427287.1888925899 -2265396.9593316740 +-411819.7618779668 -2254431.8150856020 +-413108.7141291854 -2243466.6708395295 +-414397.6663804040 -2243466.6708395295 + Index: /issm/workshop/09SeaRISE_GrIS/Exp_Par/Jak_outline.exp =================================================================== --- /issm/workshop/09SeaRISE_GrIS/Exp_Par/Jak_outline.exp (revision 18013) +++ /issm/workshop/09SeaRISE_GrIS/Exp_Par/Jak_outline.exp (revision 18013) @@ -0,0 +1,17 @@ +## Name:Jak +## Icon:0 +# Points Count Value +11 1.000000 +# X pos Y pos +-421561.7247949215 -2216716.3488182197 +-386019.6416099078 -2203459.4896082953 +-351446.8879663036 -2185999.2360147359 +-289086.6874689616 -2214129.6445821370 +-255160.1535196303 -2282353.9688088228 +-340138.0433165266 -2317597.8140254519 +-383306.1449952667 -2302506.6581698526 +-416391.9672407377 -2279120.5885137189 +-448056.7322601135 -2246463.4475331730 +-429126.8716234095 -2231167.6696134652 +-421561.7247949215 -2216716.3488182197 + Index: /issm/workshop/09SeaRISE_GrIS/runme.m =================================================================== --- /issm/workshop/09SeaRISE_GrIS/runme.m (revision 18013) +++ /issm/workshop/09SeaRISE_GrIS/runme.m (revision 18013) @@ -0,0 +1,272 @@ +clear all; +steps=[ 1 2 3 4 5 ]; + +%Location of SeaRISE dataset +ncdata='./Data/Greenland_5km_dev1.2.nc'; + +if any(steps==1) + disp(' Step 1: Mesh creation'); + + %Generate initial uniform mesh (resolution = 20000 m) + md=triangle(model,'./Exp_Par/DomainOutline.exp',20000); + + % Get velocities (Note: You can use ncdisp('file') to see an ncdump) + x1 = ncread(ncdata,'x1'); + y1 = ncread(ncdata,'y1'); + velx = ncread(ncdata,'surfvelx'); + vely = ncread(ncdata,'surfvely'); + vx = InterpFromGridToMesh(x1,y1,velx',md.mesh.x,md.mesh.y,0); + vy = InterpFromGridToMesh(x1,y1,vely',md.mesh.x,md.mesh.y,0); + vel = sqrt(vx.^2+vy.^2); + + %Mesh Greenland + md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8); + + %convert x,y coordinates (Polar stereo) to lat/lon + [md.mesh.lat,md.mesh.long]=xy2ll(md.mesh.x,md.mesh.y,+1,39,71); + + save ./Models/Greenland.Mesh_generation md; + + plotmodel (md,'data','mesh'); +end + +if any(steps==2) + disp(' Step 2: Parameterization'); + md = loadmodel('./Models/Greenland.Mesh_generation'); + + md = setmask(md,'',''); + md = parameterize(md,'./Exp_Par/Greenland.par'); + md = setflowequation(md,'SSA','all'); + + save ./Models/Greenland.Parameterization md; +end + +if any(steps==3) + disp(' Step 3: Control method friction'); + md = loadmodel('./Models/Greenland.Parameterization'); + + %Control general + md.inversion.iscontrol=1; + md.inversion.nsteps=30; + md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1); + md.inversion.maxiter_per_step=5*ones(md.inversion.nsteps,1); + + %Cost functions + md.inversion.cost_functions=[101 103 501]; + md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,3); + md.inversion.cost_functions_coefficients(:,1)=350; + md.inversion.cost_functions_coefficients(:,2)=0.6; + md.inversion.cost_functions_coefficients(:,3)=2e-6; + + %Controls + md.inversion.control_parameters={'FrictionCoefficient'}; + md.inversion.gradient_scaling(1:md.inversion.nsteps)=50; + md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1); + md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1); + + %Additional parameters + md.stressbalance.restol=0.01; md.stressbalance.reltol=0.1; + md.stressbalance.abstol=NaN; + + %Go solve + md.cluster=generic('name',oshostname,'np',2); + md.toolkits=toolkits; + md.verbose=verbose('solution',true,'control',true); + md=solve(md,StressbalanceSolutionEnum); + + %Update model friction fields accordingly + md.friction.coefficient=md.results.StressbalanceSolution.FrictionCoefficient; + + save ./Models/Greenland.Control_drag md; +end + +if any(steps==4) + disp(' Step 4: Transient run'); + md = loadmodel('./Models/Greenland.Control_drag'); + + %Set surface mass balance + x1 = ncread(ncdata,'x1'); + y1 = ncread(ncdata,'y1'); + smb = ncread(ncdata,'smb'); + smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0); + smb = smb*md.materials.rho_freshwater/md.materials.rho_ice; + smb = [smb smb smb-1.0]; + md.surfaceforcings.mass_balance = [smb;1 10 20]; + + %Set transient options, run for 20 years, saving every timestep + md.timestepping.time_step=0.2; + md.timestepping.final_time=20; + md.settings.output_frequency=1; + + %Additional options + md.inversion.iscontrol=0; + md.transient.requested_outputs={'IceVolume','TotalSmb', ... + 'SurfaceforcingsMassBalance'}; + md.verbose=verbose('solution',true,'module',true,'convergence',true); + + %Go solve + md.cluster=generic('name',oshostname,'np',2); + md=solve(md,TransientSolutionEnum); + + save ./Models/Greenland.Transient md; +end + +if any(steps==5) + disp(' Step 5: Plotting'); + md = loadmodel('./Models/Greenland.Transient'); + %md = loadmodel('./Models/Greenland.HistoricTransient'); + + %Planview plots + plotmodel(md,'data',md.results.TransientSolution(end).Vel,'caxis',[1e-1 6000],... + 'log', 10, 'title', 'Velocity (m/y)','gridded',1, ... + 'data', md.results.TransientSolution(1).SurfaceforcingsMassBalance, ... + 'title', 'Surface mass balance (m/y)', 'gridded',1,... + 'data',md.results.TransientSolution(end).Thickness,... + 'title','Thickness (m)', 'gridded',1, ... + 'data',md.results.TransientSolution(end).Surface, ... + 'title', 'Surface (m)', 'gridded',1); + + %Line Plots + figure + + %Plot surface mass balance + surfmb=[]; for i=1:100; surfmb=[surfmb ... + md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end + subplot(3,1,1); plot([0.2:0.2:20],mean(surfmb)); title('Mean Surface mass balance'); + + %Plot velocity + vel=[]; for i=1:100; vel=[vel md.results.TransientSolution(i).Vel]; end + subplot(3,1,2); plot([0.2:0.2:20],mean(vel)); title('Mean Velocity'); + + %Plot Volume + volume=[]; for i=1:100; volume=[volume md.results.TransientSolution(i).IceVolume]; end + subplot(3,1,3); plot([0.2:0.2:20],volume); title('Ice Volume'); + xlabel('years') +end + +if any(steps==6) + disp(' Step 6: Historical Relaxation run'); + md = loadmodel('./Models/Greenland.Control_drag'); + + ncbox='./Data/Box_Greenland_SMB_monthly_1840-2012_5km_ver20140421.nc'; + + %convert mesh x,y into the Box projection + [md.mesh.lat,md.mesh.long] = xy2ll(md.mesh.x,md.mesh.y,+1,39,71); + [xi,yi]= ll2xy(md.mesh.lat,md.mesh.long,+1,45,70); + + %Set surface mass balance + lat = ncread(ncbox,'lat'); + lon = ncread(ncbox,'lon'); + smbbox = ncread(ncbox,'SMB'); + [x1 y1]=ll2xy(lat,lon,+1,45,70); + + %surface mass balance + smbmean = mean(mean(smbbox,3),4); + save -v7.3 ./Data/smbbox smbbox x1 y1; + clear smbbox + + %Interpolate and set surface mass balance + index = BamgTriangulate(x1(:),y1(:)); + smb_mo = InterpFromMeshToMesh2d(index,x1(:),y1(:),smbmean(:),xi,yi); + smb = smb_mo*12/1000*md.materials.rho_freshwater/md.materials.rho_ice; + md.surfaceforcings.mass_balance = [smb;1 ]; + + %Set transient options, run for 20 years, saving every timestep + md.timestepping.time_step=0.2; + md.timestepping.final_time=20; + md.settings.output_frequency=1; + + %Additional options + md.inversion.iscontrol=0; + md.transient.requested_outputs={'IceVolume','TotalSmb', ... + 'SurfaceforcingsMassBalance'}; + md.verbose=verbose('solution',true,'module',true); + + %Go solve + md.cluster=generic('name',oshostname,'np',2); + md=solve(md,TransientSolutionEnum); + + save ./Models/Greenland.HistoricTransient md; +end + +if any(steps==7) + disp(' Step 7: Box Transient run'); + md = loadmodel('./Models/Greenland.HistoricTransient_1000yr'); + + %load past transient results + md.geometry.base=md.results.TransientSolution(end).Base; + md.geometry.thickness=md.results.TransientSolution(end).Thickness; + md.geometry.surface=md.geometry.base+md.geometry.thickness; + md.initialization.vx=(md.results.TransientSolution(end).Vx); + md.initialization.vy=(md.results.TransientSolution(end).Vy); + md.results=[]; + + %convert mesh x,y into the Box projection + [md.mesh.lat,md.mesh.long] = xy2ll(md.mesh.x,md.mesh.y,+1,39,71); + [xi,yi]= ll2xy(md.mesh.lat,md.mesh.long,+1,45,70); + + %Set surface mass balance + load ./Data/smbbox + index = BamgTriangulate(x1(:),y1(:)); + + %Set years to run + years_of_simulation = 2003:2012; + + %initialize surface mass balance matrix + smb = nan*ones(md.mesh.numberofvertices,length(years_of_simulation)*12); + + %Interpolate and set surface mass balance + for year=years_of_simulation + for month=1:12 + smb_mo = griddata(double(x1),double(y1),... + double(squeeze(smbbox(:,:,month,year-1839))),xi,yi,'nearest'); + smb(:,(year-years_of_simulation(1))*12+month) = smb_mo; + end + end + md.surfaceforcings.mass_balance = ... + [smb*12/1000*md.materials.rho_freshwater/md.materials.rho_ice; ... + [1/24:1/12:length(years_of_simulation)]]; + + %Set transient options, monthly timestep, saving every month + md.timestepping.time_step=1/12; + md.timestepping.final_time=length(years_of_simulation); + md.settings.output_frequency=1; + + %Additional options + md.inversion.iscontrol=0; + md.transient.requested_outputs={'IceVolume','TotalSmb', ... + 'SurfaceforcingsMassBalance'}; + md.verbose=verbose('solution',true,'module',true); + + %Go solve + md.cluster=generic('name',oshostname,'np',2); + md=solve(md,TransientSolutionEnum); + + save ./Models/Greenland.BoxTransient md; +end + +if any(steps==8) + disp(' Step 8: Plot Box Transient'); + md = loadmodel('./Models/Greenland.BoxTransient'); + + %Set years run + years_of_simulation = 2003:2012; + t = [years_of_simulation(1):1/12:years_of_simulation(end)+11/12]; + + %Line Plots + figure + + %Plot surface mass balance + surfmb=[]; for i=1:length(t); surfmb=[surfmb ... + md.results.TransientSolution(i).TotalSmb]; end + subplot(3,1,1); plot(t,surfmb); title('Total Surface mass balance'); + + %Plot velocity + vel=[]; for i=1:length(t); vel=[vel md.results.TransientSolution(i).Vel]; end + subplot(3,1,2); plot(t,max(vel)); title('Max Velocity'); + + %Plot Volume + volume=[]; for i=1:length(t); volume=[volume md.results.TransientSolution(i).IceVolume]; end + subplot(3,1,3); plot(t,volume); title('Ice Volume'); + xlabel('years') +end Index: /issm/workshop/10IceBridge/Exp_Par/DomainOutline.exp =================================================================== --- /issm/workshop/10IceBridge/Exp_Par/DomainOutline.exp (revision 18013) +++ /issm/workshop/10IceBridge/Exp_Par/DomainOutline.exp (revision 18013) @@ -0,0 +1,733 @@ +## Name:DomainOutline +## Icon:0 +# Points Count Value +727 1.000000 +# X pos Y pos +47662.8073760000 -2621047.2585999998 +39869.9055420000 -2638077.6292949999 +26037.5047870000 -2644088.3483640002 +24868.5695120000 -2659115.1460360000 +25647.8596950000 -2674141.9437079998 +16880.8451320000 -2687666.0616130000 +515.7512810000 -2685662.4885900002 +-7666.7956450000 -2672639.2639410002 +-22083.6640370000 -2677147.3032430001 +-30071.3884170000 -2690170.5278920000 +-44683.0793550000 -2694678.5671930001 +-41760.7411680000 -2709705.3648660001 +-37669.4677050000 -2724231.2692820001 +-44488.2568100000 -2738757.1736989999 +-55787.9644690000 -2753283.0781149999 +-67672.1397650000 -2742263.4264890002 +-87154.3943500000 -2741762.5332330000 +-105078.0685680000 -2741011.1933490001 +-97090.3441880000 -2760546.0303230002 +-99817.8598300000 -2780581.7605530000 +-82478.6532500000 -2789096.9459000002 +-77218.4445120000 -2804123.7435730002 +-77023.6219660000 -2820152.3277560002 +-95336.9412760000 -2829669.2996149999 +-112870.9704020000 -2832674.6591500002 +-123196.5653320000 -2818649.6479890002 +-122417.2751480000 -2834177.3389170002 +-120274.2271440000 -2849204.1365890000 +-105662.5362060000 -2854713.9624020001 +-89492.2649000000 -2850706.8163560000 +-77997.7346950000 -2860724.6814710000 +-92024.9579960000 -2867737.1870519998 +-84816.5238000000 -2890277.3835600000 +-94752.4736380000 -2901797.9284419999 +-114234.7282230000 -2894284.5296060001 +-127677.4838860000 -2883765.7712349999 +-146964.9159250000 -2888774.7037930000 +-139951.3042750000 -2905805.0744880000 +-154173.3501210000 -2916824.7261140002 +-165667.8803260000 -2904302.3947210000 +-162550.7195930000 -2920330.9789040000 +-152809.5923000000 -2934355.9900650000 +-147744.2061080000 -2949883.6809930000 +-163719.6548680000 -2950384.5742489998 +-175603.8301650000 -2963908.6921540000 +-189825.8760110000 -2956896.1865730002 +-195670.5523870000 -2971422.0909899999 +-210282.2433250000 -2966413.1584319999 +-222750.8862590000 -2978301.3022220000 +-214763.1618800000 -2995247.5385030000 +-222750.8862590000 -3011386.8111520000 +-206970.2600460000 -3016228.5929470002 +-217295.8549760000 -3027122.6019839998 +-232686.8360980000 -3031964.3837790000 +-212035.6462380000 -3035192.2383090002 +-212814.9364210000 -3050928.0291410000 +-197229.1327540000 -3048910.6200600001 +-185344.9574570000 -3058190.7018329999 +-174824.5399810000 -3071102.1199520002 +-182812.2643610000 -3084820.5017030002 +-192163.7465610000 -3098942.3652710002 +-176577.9428940000 -3092486.6562120002 +-184955.3123650000 -3105801.5561469998 +-182617.4418150000 -3121537.3469790001 +-199372.1807580000 -3120326.9015310002 +-218464.7902510000 -3115888.6015519998 +-205022.0345870000 -3126782.6105900002 +-189825.8760110000 -3128800.0196710001 +-198592.8905740000 -3146149.7377680000 +-213009.7589670000 -3140904.4741580002 +-225283.5793550000 -3131220.9105679998 +-222750.8862590000 -3146149.7377680000 +-208334.0178670000 -3151395.0013790000 +-195670.5523870000 -3161078.5649680002 +-213009.7589670000 -3171569.0921900002 +-227621.4499060000 -3161078.5649680002 +-242622.7859360000 -3157850.7104389998 +-240090.0928400000 -3174393.4649040001 +-220023.3706180000 -3185690.9557579998 +-233076.4811890000 -3193357.1102660000 +-258013.7670580000 -3188111.8466550000 +-275158.1510920000 -3196181.4829790001 +-258208.5896040000 -3198198.8920600000 +-252753.5583200000 -3213127.7192600002 +-270482.4099920000 -3215952.0919740000 +-258793.0572410000 -3226442.6191960000 +-275547.7961840000 -3232969.1157640000 +-276521.9089130000 -3248275.5695210001 +-293276.6478560000 -3241868.2167850002 +-282561.4078340000 -3256818.7065010001 +-288211.2616640000 -3273549.0164200002 +-302433.3075110000 -3267141.6636850000 +-305745.2907900000 -3251479.2458879999 +-313733.0151700000 -3237596.6482950002 +-310810.6769820000 -3218018.6260489998 +-325812.0130120000 -3213747.0575589999 +-338280.6559470000 -3203780.0644160002 +-336527.2530340000 -3187761.6825779998 +-321720.7395500000 -3185269.9342919998 +-325812.0130120000 -3165691.9120459999 +-313538.1926240000 -3155368.9548610002 +-329708.4639290000 -3158928.5952699999 +-338085.8334010000 -3173167.1569030001 +-351723.4116100000 -3164268.0558819999 +-342566.7519550000 -3151809.3144530002 +-325812.0130120000 -3153233.1706160000 +-310615.8544360000 -3148249.6740450002 +-329123.9962920000 -3150029.4942490002 +-343151.2195930000 -3144334.0695949998 +-361854.1839940000 -3145401.9617180000 +-377634.8102080000 -3142910.2134320000 +-388350.0502290000 -3154657.0267800000 +-397701.5324300000 -3178150.6534750001 +-417573.4321060000 -3179930.4736790000 +-432720.8850460000 -3186308.0579260001 +-451813.4945390000 -3187438.8412219998 +-467788.9432980000 -3177261.7915590000 +-482985.1018740000 -3173492.5139070000 +-490583.1811620000 -3158792.3310610000 +-492726.2291670000 -3141453.6538590002 +-490193.5360710000 -3119214.9157079998 +-505389.6946470000 -3113560.9992289999 +-489803.8909790000 -3110922.5048719998 +-478309.3607740000 -3100368.5274439999 +-495648.5673540000 -3101122.3829750000 +-512792.9513890000 -3108660.9382799999 +-501103.5986380000 -3097353.1053220001 +-493895.1644420000 -3083783.7057719999 +-486491.9077000000 -3068329.6673960001 +-503636.2917340000 -3066068.1008040002 +-522728.9012270000 -3064183.4619780001 +-521559.9659520000 -3046467.8570099999 +-531301.0932440000 -3029506.1075729998 +-543769.7361790000 -3012167.4303710000 +-545912.7841830000 -2992944.1143419999 +-550004.0576459999 -2974474.6538439998 +-567538.0867720000 -2983897.8479749998 +-585851.4060810000 -2995959.5364640001 +-589942.6795440000 -2981259.3536189999 +-572213.8278720001 -2975982.3649050002 +-565200.2162220000 -2959397.5432330002 +-548055.8321870000 -2954120.5545190000 +-545717.9616370000 -2933389.5274290000 +-536561.3019820000 -2920197.0556450002 +-545523.1390910000 -2907758.4393910002 +-533444.1412490000 -2898335.2452590000 +-530568.0950890000 -2883304.2548420001 +-545624.1658510000 -2879383.1998239998 +-559750.8495300000 -2868110.1666489998 +-577223.3267109999 -2867129.9028940001 +-590606.5007220000 -2859287.7928590002 +-573319.9009570000 -2851445.6828239998 +-551572.2431889999 -2846054.2321750000 +-530939.8499220000 -2834781.1989989998 +-528337.5660870000 -2817626.5832969998 +-543765.3916830000 -2804883.1544900001 +-528151.6886699999 -2802432.4951040000 +-528151.6886699999 -2785768.0112790000 +-513467.3727410000 -2778416.0331210000 +-513839.1275750000 -2760771.2855420001 +-525735.2822510001 -2749498.2523670001 +-509563.9469880000 -2738225.2191909999 +-517184.9210780000 -2724501.5266300002 +-532240.9918400000 -2725971.9222610001 +-524248.2629170000 -2707837.0428050002 +-512537.9856570000 -2694113.3502429998 +-521088.3468310000 -2678429.1301730000 +-538189.0691780000 -2669606.7563829999 +-556033.3011930001 -2677448.8664190001 +-539676.0885130001 -2666175.8332429999 +-526664.6693360000 -2658333.7232080000 +-517556.6759110000 -2639708.7118739998 +-526107.0370850000 -2626965.2830670001 +-538932.5788460000 -2618142.9092779998 +-550456.9786890000 -2607360.0079790000 +-565698.9268680000 -2599517.8979440001 +-581126.7524640000 -2600988.2935759998 +-583110.5023000001 -2578346.1085480000 +-572583.8192940000 -2566324.6704099998 +-556315.3091930000 -2558811.2715739999 +-538132.8567270000 -2571333.6029679999 +-524735.2601730000 -2578346.1085480000 +-519759.0100240000 -2594875.5859869998 +-503107.7114500000 -2585859.5073839999 +-486456.4128760000 -2576342.5355250002 +-477843.6722340000 -2559312.1648300001 +-474015.7875040000 -2536771.9683220000 +-473441.6047950000 -2521244.2773940000 +-474972.7586870000 -2504714.7999539999 +-485690.8359300000 -2491691.5753049999 +-486073.6244030000 -2476664.7776330002 +-492581.0284430000 -2459634.4069380001 +-474207.1817410000 -2447112.0755440001 +-465594.4410990000 -2434589.7441509999 +-462340.7390790000 -2418060.2667109999 +-462914.9217880000 -2399026.3229930000 +-468465.3546460000 -2375484.3399729999 +-473058.8163220000 -2360457.5423010001 +-478417.8549440000 -2343427.1716060000 +-504256.0768690000 -2331405.7334679998 +-489901.5091330000 -2322389.6548649999 +-475164.1529230000 -2307863.7504480002 +-462577.4401970000 -2287577.5735909999 +-457317.2314590000 -2269044.5231289999 +-439783.2023330000 -2268042.7366169998 +-436666.0415990000 -2252014.1524330000 +-435886.7514160000 -2233481.1019710000 +-448550.2168960000 -2220457.8773220000 +-440172.8474240000 -2206933.7594170002 +-429262.7848570000 -2193409.6415120000 +-433743.7034120000 -2176880.1640720000 +-426730.0917610000 -2162354.2596560000 +-419716.4801110000 -2147327.4619840002 +-428678.3172200000 -2134805.1305900002 +-427314.5593990000 -2110762.2543150000 +-426535.2692150000 -2094733.6701310000 +-427509.3819440000 -2067184.5410650000 +-428093.8495820000 -2050655.0636260000 +-425366.3339400000 -2035628.2659539999 +-432379.9455910000 -2021102.3615369999 +-446017.5238000000 -2003571.0975860001 +-430821.3652240000 -2002569.3110750001 +-428678.3172200000 -1986039.8336360001 +-424197.3986650000 -1971013.0359630000 +-418060.4884710000 -1957238.4714309999 +-433451.4695930000 -1950225.9658500000 +-426243.0353960000 -1933696.4884110000 +-441634.0165180000 -1943714.3535259999 +-433841.1146840000 -1926683.9828300001 +-448647.6281690000 -1933696.4884110000 +-464428.2543820000 -1948222.3928270000 +-465402.3671120000 -1933195.5951550000 +-463843.7867450000 -1914662.5446929999 +-456245.7074570000 -1901138.4267879999 +-468714.3503910000 -1892122.3481840000 +-480208.8805960000 -1905646.4660890000 +-479429.5904130000 -1888616.0953939999 +-478455.4776830000 -1871585.7246989999 +-478650.3002290000 -1856558.9270270001 +-484494.9766050000 -1841031.2360990001 +-494430.9264430000 -1822498.1856370000 +-505925.4566480000 -1807972.2812200000 +-508847.7948350000 -1791442.8037810000 +-504172.0537350000 -1776416.0061079999 +-510990.8428400000 -1759886.5286689999 +-502223.8282770000 -1745861.5175079999 +-512939.0682980000 -1728330.2535570001 +-511964.9555690000 -1712802.5626300001 +-505340.9890100000 -1697775.7649570000 +-506315.1017390000 -1681747.1807740000 +-507289.2144690000 -1662713.2370559999 +-508446.9486300000 -1645056.7497910000 +-502180.9448550000 -1623518.3397939999 +-505662.0580630000 -1602480.8230530000 +-505836.1137240000 -1586953.1321250000 +-509317.2269320000 -1571425.4411970000 +-516105.3976870000 -1557400.4300360000 +-530552.0175000000 -1537865.5930620001 +-533336.9080660000 -1513321.8235309999 +-522719.5127820000 -1498295.0258589999 +-506532.3363650000 -1492284.3067900001 +-522893.5684430000 -1484770.9079539999 +-515757.2863660000 -1469744.1102819999 +-532118.5184440000 -1465236.0709800001 +-535947.7429720000 -1449708.3800520001 +-544476.4703320001 -1428670.8633109999 +-551960.8637290000 -1415146.7454059999 +-565363.1495800000 -1405629.7735470000 +-567277.7618440000 -1386094.9365729999 +-582420.6042990000 -1384091.3635499999 +-593212.0552430000 -1370066.3523899999 +-606266.2297730000 -1362052.0602980000 +-613793.7528780000 -1344395.5730330001 +-621618.1005840000 -1359923.2639609999 +-634079.0987820000 -1369440.2358200001 +-633354.6221430000 -1354413.4381480000 +-638860.6446030000 -1339386.6404760000 +-649438.0035390001 -1351908.9718690000 +-667405.0241970000 -1352409.8651250000 +-679431.3364110000 -1363429.5167510000 +-679721.1270670000 -1346900.0393119999 +-676533.4298540000 -1330370.5618720001 +-690008.6953470000 -1339386.6404760000 +-701165.6355950000 -1327365.2023380001 +-688994.4280520000 -1315844.6574560001 +-709279.7739560000 -1318349.1237339999 +-717249.0169900000 -1331372.3483839999 +-735216.0376480001 -1329869.6686160001 +-741420.6668040000 -1315969.8807699999 +-750947.3474350000 -1299440.4033299999 +-734597.5036490000 -1299440.4033299999 +-722496.0444700000 -1284413.6056580001 +-707948.5456690000 -1288420.7517039999 +-689023.9233340001 -1290424.3247270000 +-678209.8534290000 -1276399.3135660000 +-672545.3406220000 -1262374.3024060000 +-693401.0468670001 -1264377.8754280000 +-709493.4127980000 -1260871.6226379999 +-725971.9955110000 -1250853.7575240000 +-725714.5176560000 -1234825.1733400000 +-720951.1773400001 -1215290.3363659999 +-737172.2821980000 -1218295.6958999999 +-749402.4803050000 -1205272.4712510000 +-736399.8486340001 -1196757.2859040000 +-720822.4384130000 -1191247.4600899999 +-710652.0631450000 -1204771.5779949999 +-697906.9093280000 -1218796.5891559999 +-684131.8440910000 -1227812.6677590001 +-668940.6506530000 -1229816.2407819999 +-666237.1331770000 -1253358.2238020001 +-653620.7182870000 -1262374.3024060000 +-648471.1611900000 -1247347.5047330000 +-634696.0959530000 -1258868.0496149999 +-632625.0227419999 -1243541.0164699999 +-618416.8083510000 -1233878.1370890001 +-601000.2874830000 -1228080.4094600000 +-615361.2783740000 -1217451.2421410000 +-598861.4165000001 -1216968.0981719999 +-610625.2069099999 -1201990.6351310001 +-617652.9258560000 -1185080.5962139999 +-631708.3637490000 -1178799.7246160000 +-648666.5551200000 -1183631.1643070001 +-656458.1565610000 -1169136.8452349999 +-653708.1795820000 -1154159.3821950001 +-665166.4169950000 -1166237.9814210001 +-676166.3249110000 -1153193.0942569999 +-676166.3249110000 -1133384.1915249999 +-689305.1038110000 -1118406.7284850001 +-709013.2721609999 -1111159.5689490000 +-711610.4726409999 -1091833.8101870001 +-716346.5441050000 -1074440.6273010001 +-701527.2237180000 -1070092.3315790000 +-680291.2903790000 -1068642.8996720000 +-664096.9815030000 -1078788.9230219999 +-650958.2026030000 -1100530.4016300000 +-628041.7277770001 -1095698.9619390001 +-611236.3129050001 -1091350.6662170000 +-588930.9440750000 -1094732.6740009999 +-577319.9301640000 -1078305.7790530000 +-565250.5867559999 -1066710.3237960001 +-550431.2663680000 -1071541.7634860000 +-540042.4644469999 -1060429.4521979999 +-531203.6819710000 -1047613.0696300000 +-510357.6695660000 -1037937.0529810000 +-498473.4942690000 -1027821.2173940000 +-486199.6738810000 -1015066.4681760001 +-487953.0767930000 -994834.7970020000 +-464964.0163830000 -990876.4265550000 +-449378.2127160000 -998353.3485110000 +-430285.6032230000 -1001432.0810810000 +-417622.1377430000 -1012427.5545450000 +-399308.8184330000 -1004950.6325890000 +-384697.1274945000 -1000332.5337340001 +-370085.4365560000 -995714.4348789999 +-375540.4678400000 -980760.5909680000 +-380021.3863940000 -966246.5659950000 +-383528.1922200000 -949973.2652690000 +-383528.1922200000 -930181.4130330000 +-384112.6598570000 -911709.0176140000 +-380216.2089400000 -883560.6055450001 +-366919.5701860000 -892846.9764360000 +-353087.1694310000 -880338.0618260000 +-336527.2530340000 -873692.7009390000 +-322694.8522790000 -888547.0370390000 +-307888.3387940000 -897146.9158330000 +-294835.2282230000 -913564.8662580000 +-275937.4412760000 -919037.5164000000 +-252169.0906820000 -915128.4805850000 +-245350.3015780000 -901446.8552300000 +-244181.3663030000 -918646.6128190000 +-231615.3120950000 -930563.0925760000 +-215445.0407900000 -924017.9382920000 +-206872.8487730000 -937108.2468600000 +-190507.7549220000 -933426.5975750000 +-186611.3040050000 -912972.9904390000 +-185961.8955186667 -896882.8194913333 +-185312.4870323333 -880792.6485436666 +-184663.0785460000 -864702.4775960000 +-175311.5963460000 -883519.7961620000 +-174532.3061620000 -902337.1147280000 +-172584.0807040000 -918700.0004370000 +-168882.4523330000 -933835.6697180000 +-152322.5359360000 -933835.6697180000 +-142581.4086430000 -921563.5054359999 +-127969.7177050000 -912154.8461530000 +-131086.8784390000 -930563.0925760000 +-115501.0747710000 -924427.0104349999 +-112189.0914910000 -905609.6918690000 +-114916.6071330000 -887610.5175890001 +-101473.8514700000 -899882.6818710000 +-84329.4674350000 -896610.1047290000 +-100110.0936490000 -892110.3111590000 +-87836.2732610000 -879838.1468770000 +-74685.7514160000 -869305.2030530000 +-57541.3673810000 -863251.9975030000 +-49553.6430020000 -850388.9357080000 +-34357.4844260000 -846227.3568920000 +-47995.0626350000 -825797.7881600000 +-30071.3884170000 -835255.9218320000 +-18771.6807580000 -850767.2610550000 +-6303.0378240000 -859847.0693810000 +7139.7178400000 -852658.8877899999 +22141.0538700000 -862873.6721560000 +31297.7135250000 -877628.3606850000 +47662.8073760000 -874601.7579100000 +49611.0328340000 -893518.0252550000 +66365.7717770000 -886708.1690110000 +81123.5796250000 -882220.1151260000 +97683.4960220000 -884202.5846670000 +80149.4668960000 -888960.5115630000 +65537.7759570000 -896493.8958170000 +68460.1141450000 -914732.6155880000 +84435.5629050000 -917508.0729450000 +100605.8342100000 -907595.7252430000 +120282.9113400000 -911957.1582320000 +129049.9259030000 -927420.4206470000 +153402.7441340000 -925437.9511060000 +169767.8379850000 -920680.0242090001 +180677.9005530000 -910371.1825990001 +192367.2533040000 -899269.3531730000 +201134.2678670000 -912353.6521400000 +182626.1260110000 -920680.0242090001 +169183.3703480000 -936539.7805320000 +163338.6939720000 -950417.0673150000 +182820.9485570000 -955174.9942120000 +187301.8671120000 -970241.7627180000 +200744.6227750000 -980550.6043280000 +198796.3973170000 -995617.3728350000 +211118.9233410000 -983385.4589040000 +211508.5684330000 -999969.0227530000 +225146.1466420000 -1008484.9068910000 +240731.9503100000 -1009829.5201760000 +238588.9023060000 -1032239.7415930000 +253785.0608820000 -1029550.5150230000 +273072.4929210000 -1036273.5814480000 +281255.0398460000 -1056890.9851510001 +286320.4260380000 -1071681.7312860000 +305413.0355310000 -1086024.2729930000 +316907.5657360000 -1095884.7704169999 +326453.8704830000 -1081990.4331380001 +341260.3839670000 -1076611.9799980000 +355677.2523600000 -1064958.6648609999 +372319.0245690000 -1056006.2810410000 +387981.0391590000 -1068027.7191780000 +383779.0352450000 -1084056.3033620000 +366207.0188750000 -1085558.9831290001 +366971.0195860000 -1102088.4605690001 +352646.0062410000 -1113108.1121950001 +362960.0158500000 -1126131.3368440000 +380914.0325760000 -1133143.8424249999 +394093.0448530000 -1123626.8705660000 +409182.0589100000 -1113108.1121950001 +409182.0589100000 -1131641.1626569999 +404407.0544620000 -1158689.3984670001 +386262.0375580000 -1153179.5726540000 +374229.0263480000 -1166202.7973030000 +364488.0172730000 -1182732.2747430000 +377285.0291950000 -1191748.3533460000 +388172.0393370000 -1202768.0049719999 +392947.0437860000 -1219798.3756680000 +380341.0320420000 -1228814.4542710001 +383588.0350670000 -1251354.6507790000 +389509.0405830000 -1265880.5551960000 +389509.0405830000 -1282910.9258910001 +399632.0500130000 -1295934.1505400001 +408609.0583760000 -1313966.3077469999 +410137.0598000000 -1333000.2514650000 +413526.9681570000 -1350406.2921020000 +416526.2265730000 -1367436.6627970000 +415776.4119690000 -1385468.8200040001 +419900.3922910000 -1406005.4434890000 +430960.1577010000 -1421032.2411610000 +431897.4259560000 -1436059.0388330000 +448018.4399430000 -1436559.9320890000 +448580.8008960000 -1456094.7690630001 +438458.3037410000 -1469117.9937120001 +432459.7869090000 -1486148.3644069999 +449705.5228020000 -1494663.5497550000 +448393.3472450000 -1513196.6002170001 +463202.1856750000 -1521210.8923090000 +447830.9862920000 -1525218.0383550001 +458140.9370980000 -1539243.0495160001 +462639.8247220000 -1554269.8471880001 +473887.0437820000 -1564287.7123030000 +458703.2980510000 -1557275.2067219999 +463014.7320240000 -1574305.5774170000 +478948.2923600000 -1567293.0718370001 +484009.5409370000 -1584824.3357879999 +475011.7656880000 -1597346.6671819999 +484009.5409370000 -1614877.9311319999 +488508.4285610000 -1630906.5153160000 +470512.8780640000 -1618384.1839230000 +461702.5564670000 -1631407.4085720000 +474074.4974330000 -1641926.1669419999 +488883.3358630000 -1648938.6725230000 +495661.8324810000 -1662963.6836830000 +478923.6439970000 -1667471.7229849999 +463095.1396700000 -1663965.4701950001 +449449.8773190000 -1652945.8185690001 +455817.6664170000 -1670477.0825199999 +441626.5935720000 -1662963.6836830000 +433439.4361610000 -1648437.7792670000 +409690.0603430000 -1652315.3206390000 +394504.9542530000 -1655951.1781029999 +411243.1427370000 -1662462.7904280000 +426707.7734010000 -1676487.8015880000 +437260.1096190000 -1691514.5992610001 +455999.6032480000 -1702534.2508870000 +441626.5935720000 -1710047.6497229999 +428891.0153770000 -1692516.3857720001 +435076.8676430000 -1708544.9699560001 +419976.1106410000 -1704036.9306540000 +408696.0270980000 -1692516.3857720001 +395050.7647470000 -1701532.4643750000 +388501.0388190000 -1718562.8350710000 +378494.5130940000 -1731085.1664640000 +376129.3342870000 -1713553.9025129999 +361938.2614420000 -1720566.4080940001 +373582.2186480000 -1732086.9529760000 +358481.4616460000 -1738598.5652999999 +369397.6715270000 -1754627.1494839999 +352659.4830430000 -1751621.7899490001 +335805.3526830000 -1754501.9261700001 +338338.0457790000 -1777042.1226780000 +354897.9621760000 -1775539.4429110000 +349053.2858010000 -1792569.8136060000 +330739.9664910000 -1791568.0270950000 +344182.7221550000 -1806093.9315110000 +341455.2065130000 -1823625.1954620001 +351780.8014430000 -1809600.1843010001 +344377.5447010000 -1826129.6617409999 +356066.8974510000 -1836147.5268550001 +345546.4799760000 -1853678.7908059999 +363080.5091020000 -1853678.7908059999 +371847.5236650000 -1841156.4594129999 +370288.9432980000 -1859689.5098750000 +386848.8596950000 -1868204.6952229999 +404772.5339130000 -1856183.2570849999 +415098.1288430000 -1867703.8019669999 +428735.7070520000 -1860691.2963870000 +421137.6277640000 -1874215.4142920000 +413929.1935680000 -1888741.3187080000 +418799.7572140000 -1903768.1163800000 +428735.7070520000 -1916791.3410300000 +415098.1288430000 -1935825.2847480001 +398148.5673540000 -1926809.2061439999 +380030.0705910000 -1929814.5656790000 +364444.2669230000 -1929814.5656790000 +348468.8181630000 -1924304.7398659999 +333077.8370420000 -1919796.7005640001 +348079.1730720000 -1935825.2847480001 +332883.0144960000 -1939832.4307939999 +318855.7911950000 -1951853.8689309999 +331908.9017660000 -1964376.2003250001 +349442.9308930000 -1959868.1610230000 +358209.9454560000 -1943338.6835840000 +358404.7680020000 -1958365.4812560000 +345936.1250670000 -1971388.7059050000 +345351.6574300000 -1987918.1833450000 +330739.9664910000 -1991925.3293910001 +337753.5781420000 -2005449.4472960001 +353534.2043550000 -2007453.0203180001 +365223.5571060000 -1997435.1552040000 +385679.9244200000 -1986415.5035770000 +399317.5026290000 -1977399.4249740001 +399122.6800840000 -1993928.9024139999 +381393.8284120000 -1998436.9417150000 +399512.3251750000 -1999939.6214820000 +395615.8742580000 -2018973.5652000001 +387822.9724240000 -2032998.5763610001 +383975.2271440000 -2048275.8206610000 +398976.5631740000 -2062801.7250780000 +397223.1602620000 -2078830.3092610000 +379299.4860440000 -2080833.8822840001 +362739.5696470000 -2085842.8148419999 +374428.9223980000 -2097363.3597240001 +390599.1937030000 -2089849.9608880000 +379299.4860440000 -2103374.0787930000 +364492.9725590000 -2111388.3708839999 +355531.1354500000 -2097363.3597240001 +343257.3150620000 -2111889.2641400001 +351829.5070790000 -2127917.8483239999 +365661.9078340000 -2119402.6629760000 +357284.5383630000 -2133928.5673929998 +374039.2773060000 -2136933.9269269998 +385533.8075110000 -2147452.6852980000 +375792.6802180000 -2160976.8032030002 +390988.8387940000 -2162479.4829699998 +378520.1958600000 -2174000.0278520002 +394105.9995280000 -2164983.9492489998 +409107.3355580000 -2161978.5897140000 +421186.3334010000 -2149957.1515759998 +437941.0723440000 -2150958.9380879998 +450409.7152780000 -2171495.5615730002 +454695.8112860000 -2156969.6571570002 +471645.3727750000 -2159975.0166910002 +459566.3749330000 -2147452.6852980000 +465605.8738540000 -2132425.8876260002 +477490.0491500000 -2120404.4494880000 +485867.4186220000 -2133427.6741370000 +493264.3211700000 -2119402.6629760000 +495978.8328800000 -2102873.1855370002 +509098.9728100000 -2110386.5843730001 +523877.9810080000 -2104876.7585600000 +541371.5009160000 -2109885.6911169998 +529608.6168400000 -2096862.4664680001 +544990.8498620000 -2099366.9327469999 +561579.5325330000 -2107381.2248390000 +576509.3469370001 -2104876.7585600000 +596868.1847600000 -2102873.1855370002 +611948.8053700000 -2108883.9046060001 +597622.2157910001 -2119402.6629760000 +580430.3082950000 -2125413.3820449999 +599130.2778520000 -2126916.0618119999 +583144.8200050000 -2134930.3539040000 +568818.2304260000 -2148454.4718089998 +582843.2075930000 -2161477.6964580002 +566706.9435400000 -2162980.3762260000 +552983.5787850000 -2170994.6683180002 +543331.9815950000 -2188025.0390130002 +547554.5553660000 -2207058.9827310001 +538053.7643820000 -2193033.9715700001 +525084.4306570000 -2209062.5557539999 +516036.0582910000 -2221083.9938920001 +529910.2292520000 -2214572.3815669999 +541220.6947100000 -2225091.1399369999 +526140.0741000000 -2223087.5669140001 +526743.2989240000 -2239617.0443540001 +521615.8879170000 -2256647.4150490002 +506233.6548950000 -2255144.7352820002 +508194.1355740000 -2275681.3587670000 +489343.3598110000 -2268167.9599310001 +488740.1349870000 -2283695.6508590002 +475770.8012620000 -2293212.6227179999 +460539.3744460000 -2281692.0778359999 +449228.9089890000 -2296217.9822519999 +448022.4593400000 -2314250.1394590000 +440180.5366230000 -2298722.4485309999 +441085.3738600000 -2277184.0385340001 +426000.8914760000 -2286450.5637659999 +431125.9775380000 -2301477.3614380001 +426190.7094780000 -2317505.9456210001 +424672.1654600000 -2332532.7432940002 +421445.2594210000 -2313498.7995750001 +409296.9072740000 -2298472.0019029998 +411764.5413040000 -2318006.8388769999 +403222.7312010000 -2304983.6142279999 +411005.2692950000 -2289956.8165560002 +393542.0130840000 -2293463.0693460000 +378356.5729000000 -2295967.5356250000 +388416.9270220000 -2308489.8670180002 +405880.1832330000 -2317005.0523660001 +417838.7173770000 -2328525.5972480001 +402463.4591910000 -2342049.7151529999 +385000.2029800000 -2341047.9286409998 +375699.1208680000 -2327022.9174799998 +384240.9309710000 -2314500.5860870001 +368865.6727850000 -2321012.1984120002 +373421.3048400000 -2337541.6758510000 +363360.9507190000 -2353069.3667790000 +351402.4165740000 -2336539.8893400002 +340392.9724410000 -2324518.4512020000 +329573.3463100000 -2313999.6928309998 +335267.8863790000 -2329527.3837589999 +336216.9763910000 -2347058.6477100002 +337166.0664020000 -2362085.4453819999 +318374.0841750000 -2366593.4846839998 +303758.0979980000 -2362586.3386380002 +290314.2882280000 -2354697.2698599999 +292067.6911410000 -2339670.4721880001 +275702.5972900000 -2347684.7642799998 +277880.2580262116 -2321542.8592617917 +269420.4169677348 -2331176.2138948478 +260960.5759092579 -2323294.3782859840 +252500.7348507810 -2343436.8470641924 +254834.4841082919 -2357157.0794203631 +258918.5453089359 -2372336.9109633607 +266740.7601810000 -2359706.2024170002 +271754.1662252456 -2369125.7927523423 +266545.9376350000 -2379241.0393909998 +291483.2235030000 -2386253.5449720002 +275897.4198350000 -2396271.4100859999 +288755.7078610000 -2404786.5954339998 +276871.5325650000 -2415806.2470600000 +263818.4219930000 -2423820.5391520001 +259532.3259840000 -2445358.9491490000 +250570.4888750000 -2432836.6177559998 +248817.0859630000 -2452371.4547290001 +247258.5055960000 -2470403.6119360002 +229919.2990160000 -2474911.6512380000 +226022.8480990000 -2490439.3421660000 +239660.4263080000 -2499956.3140250002 +227776.2510110000 -2510975.9656509999 +220178.1717230000 -2523999.1902999999 +205761.3033310000 -2533015.2689029998 +210631.8669770000 -2560063.5047129998 +203228.6102350000 -2545036.7070410000 +194266.7731260000 -2558560.8249459998 +199624.3931360000 -2581101.0214539999 +185791.9923810000 -2568077.7968049999 +189493.6207520000 -2583104.5944770002 +186571.2825650000 -2601136.7516839998 +185207.5247440000 -2584106.3809890002 +169816.5436220000 -2595126.0326149999 +156568.6105040000 -2586610.8472679998 +134748.4853690000 -2581101.0214539999 +149749.8214000000 -2579598.3416869999 +136112.2431900000 -2592120.6730809999 +146437.8381200000 -2612657.2965660002 +130462.3893610000 -2620671.5886579999 +137281.1784650000 -2635197.4930739999 +125007.3580770000 -2621172.4819140001 +106888.8613130000 -2617666.2291230001 +91887.5252830000 -2609151.0437759999 +109031.9093180000 -2607648.3640089999 +92277.1703750000 -2603641.2179629998 +97537.3791130000 -2587111.7405229998 +88380.7194580000 -2573587.6226180000 +74743.1412490000 -2567076.0102940002 +61105.5630390000 -2577093.8754090001 +66950.2394150000 -2594124.2461040001 +57598.7572140000 -2581101.0214539999 +53117.8386600000 -2598131.3921500002 +41818.1310010000 -2586610.8472679998 +30810.6571600000 -2599133.1786610000 +44253.4128240000 -2608650.1505200001 +47662.8073760000 -2621047.2585999998 + Index: /issm/workshop/10IceBridge/Exp_Par/Greenland.par =================================================================== --- /issm/workshop/10IceBridge/Exp_Par/Greenland.par (revision 18013) +++ /issm/workshop/10IceBridge/Exp_Par/Greenland.par (revision 18013) @@ -0,0 +1,103 @@ +%Name and hemisphere +md.miscellaneous.name='SeaRISEgreenland'; +md.mesh.hemisphere='n'; + +disp(' Loading SeaRISE data from NetCDF'); +ncdata='../Data/Greenland_5km_dev1.2.nc'; +x1 = ncread(ncdata,'x1'); +y1 = ncread(ncdata,'y1'); +usrf = ncread(ncdata,'usrf')'; +topg = ncread(ncdata,'topg')'; +velx = ncread(ncdata,'surfvelx')'; +vely = ncread(ncdata,'surfvely')'; +temp = ncread(ncdata,'airtemp2m')'; +smb = ncread(ncdata,'smb')'; +gflux = ncread(ncdata,'bheatflx')'; + +disp(' Interpolating surface and bedrock'); +md.geometry.base = InterpFromGridToMesh(x1,y1,topg,md.mesh.x,md.mesh.y,0); +md.geometry.surface = InterpFromGridToMesh(x1,y1,usrf,md.mesh.x,md.mesh.y,0); + +disp(' Constructing thickness'); +md.geometry.thickness=md.geometry.surface-md.geometry.base; + +%Set min thickness to 1 meter +pos0=find(md.geometry.thickness<=0); +md.geometry.thickness(pos0)=1; +md.geometry.surface=md.geometry.thickness+md.geometry.base; + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%Reading IceBridge data for Jakobshavn +disp(' reading IceBridge Jakobshavn bedrock'); +fid = fopen('../Data/Jakobshavn_2008_2011_Composite_XYZGrid.txt'); +data = fscanf(fid,'%g,%g,%g,%g,%g',[5 266400])'; +fclose(fid); + +[xi,yi]= ll2xy(md.mesh.lat,md.mesh.long,+1,45,70); +bed = flipud(reshape(data(:,5),[360 740])); bed(find(bed==-9999))=NaN; +bedy = flipud(reshape(data(:,1),[360 740])); +bedx = flipud(reshape(data(:,2),[360 740])); + +%Insert Icebridge bed and recalculate thickness +bed_jks=InterpFromGridToMesh(bedx(1,:)',bedy(:,1),bed,xi,yi,NaN); +in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,... + './Exp_Par/Jak_grounded.exp','node',1); +bed_jks(~in)=NaN; +pos=find(~isnan(bed_jks)); +md.geometry.base(pos)=bed_jks(pos); +md.geometry.thickness=md.geometry.surface-md.geometry.base; +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +disp(' Interpolating velocities '); +md.inversion.vx_obs = InterpFromGridToMesh(x1,y1,velx,md.mesh.x,md.mesh.y,0); +md.inversion.vy_obs = InterpFromGridToMesh(x1,y1,vely,md.mesh.x,md.mesh.y,0); +md.inversion.vel_obs = sqrt(md.inversion.vx_obs.^2+md.inversion.vy_obs.^2); +md.initialization.vx = md.inversion.vx_obs; +md.initialization.vy = md.inversion.vy_obs; +md.initialization.vz = zeros(md.mesh.numberofvertices,1); +md.initialization.vel= md.inversion.vel_obs; + +disp(' Interpolating temperatures'); +md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15; + +disp(' Interpolating surface mass balance'); +md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0); +md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice; + +disp(' Construct basal friction parameters'); +md.friction.coefficient=30*ones(md.mesh.numberofvertices,1); +pos=find(md.mask.groundedice_levelset<0); +md.friction.coefficient(pos)=0; %no friction applied on floating ice +md.friction.p=ones(md.mesh.numberofelements,1); +md.friction.q=ones(md.mesh.numberofelements,1); + +disp(' Construct ice rheological properties'); +md.materials.rheology_n=3*ones(md.mesh.numberofelements,1); +md.materials.rheology_B=paterson(md.initialization.temperature); +md.friction.q=ones(md.mesh.numberofelements,1); +md.friction.p=ones(md.mesh.numberofelements,1); + +disp(' Set other boundary conditions'); +md.mask.ice_levelset(md.mesh.vertexonboundary==1)=0; +md.basalforcings.melting_rate = zeros(md.mesh.numberofvertices,1); +md.thermal.spctemperature = [md.initialization.temperature;1]; %impose observed temperature on surface +md.masstransport.spcthickness = NaN*ones(md.mesh.numberofvertices,1); + +disp(' Set geothermal heat flux'); +md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0); + +disp(' Set Pressure'); +md.initialization.pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness; + +disp(' Single point constraint for continental model'); +%Set at least one vertex to velocity 0 so as to not get a singular problem (point on the wet peninsula) +md.stressbalance.referential=NaN*ones(md.mesh.numberofvertices,6); +md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices,1); +md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices,1); +md.stressbalance.spcvz = NaN*ones(md.mesh.numberofvertices,1); +location = 1.0e+06 *[.32011 -2.2039]; +[dist pos]=min(sqrt((md.mesh.x - location(1)).^2 + (md.mesh.y - location(2)).^2)); +md.stressbalance.spcvx(pos) = 0; +md.stressbalance.spcvy(pos) = 0; +md.stressbalance.spcvz(pos) = 0; + Index: /issm/workshop/10IceBridge/Exp_Par/Greenland_solution.par =================================================================== --- /issm/workshop/10IceBridge/Exp_Par/Greenland_solution.par (revision 18013) +++ /issm/workshop/10IceBridge/Exp_Par/Greenland_solution.par (revision 18013) @@ -0,0 +1,107 @@ +%Name and hemisphere +md.miscellaneous.name='SeaRISEgreenland'; +md.mesh.hemisphere='n'; + +disp(' Loading SeaRISE data from NetCDF'); +ncdata='../Data/Greenland_5km_dev1.2.nc'; +x1 = ncread(ncdata,'x1'); +y1 = ncread(ncdata,'y1'); +usrf = ncread(ncdata,'usrf')'; +topg = ncread(ncdata,'topg')'; +velx = ncread(ncdata,'surfvelx')'; +vely = ncread(ncdata,'surfvely')'; +temp = ncread(ncdata,'airtemp2m')'; +smb = ncread(ncdata,'smb')'; +gflux = ncread(ncdata,'bheatflx')'; + +disp(' Interpolating surface and bedrock'); +md.geometry.base = InterpFromGridToMesh(x1,y1,topg,md.mesh.x,md.mesh.y,0); +md.geometry.surface = InterpFromGridToMesh(x1,y1,usrf,md.mesh.x,md.mesh.y,0); + +disp(' Constructing thickness'); +md.geometry.thickness=md.geometry.surface-md.geometry.base; + +%Set min thickness to 1 meter +pos0=find(md.geometry.thickness<=0); +md.geometry.thickness(pos0)=1; +md.geometry.surface=md.geometry.thickness+md.geometry.base; + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%Reading IceBridge data for Jakobshavn +disp(' reading IceBridge Jakobshavn bedrock'); +fid = fopen('../Data/Jakobshavn_2008_2011_Composite_XYZGrid.txt'); +data = fscanf(fid,'%g,%g,%g,%g,%g',[5 266400])'; +fclose(fid); + +[xi,yi]= ll2xy(md.mesh.lat,md.mesh.long,+1,45,70); +bed = flipud(reshape(data(:,5),[360 740])); bed(find(bed==-9999))=NaN; +surf = flipud(reshape(data(:,4),[360 740])); surf(find(surf==-9999))=NaN; +bedy = flipud(reshape(data(:,1),[360 740])); +bedx = flipud(reshape(data(:,2),[360 740])); + +%Insert Icebridge bed and recalculate thickness +bed_jks=InterpFromGridToMesh(bedx(1,:)',bedy(:,1),bed,xi,yi,NaN); +surf_jks=InterpFromGridToMesh(bedx(1,:)',bedy(:,1),surf,xi,yi,NaN); +in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,... + './Exp_Par/Jak_grounded.exp','node',1); +bed_jks(~in)=NaN; +surf_jks(~in)=NaN; +pos=find(~isnan(bed_jks)); +md.geometry.base(pos)=bed_jks(pos); +md.geometry.surface(pos)=surf_jks(pos); +md.geometry.thickness=md.geometry.surface-md.geometry.base; +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +disp(' Interpolating velocities '); +md.inversion.vx_obs = InterpFromGridToMesh(x1,y1,velx,md.mesh.x,md.mesh.y,0); +md.inversion.vy_obs = InterpFromGridToMesh(x1,y1,vely,md.mesh.x,md.mesh.y,0); +md.inversion.vel_obs = sqrt(md.inversion.vx_obs.^2+md.inversion.vy_obs.^2); +md.initialization.vx = md.inversion.vx_obs; +md.initialization.vy = md.inversion.vy_obs; +md.initialization.vz = zeros(md.mesh.numberofvertices,1); +md.initialization.vel= md.inversion.vel_obs; + +disp(' Interpolating temperatures'); +md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15; + +disp(' Interpolating surface mass balance'); +md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0); +md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice; + +disp(' Construct basal friction parameters'); +md.friction.coefficient=30*ones(md.mesh.numberofvertices,1); +pos=find(md.mask.groundedice_levelset<0); +md.friction.coefficient(pos)=0; %no friction applied on floating ice +md.friction.p=ones(md.mesh.numberofelements,1); +md.friction.q=ones(md.mesh.numberofelements,1); + +disp(' Construct ice rheological properties'); +md.materials.rheology_n=3*ones(md.mesh.numberofelements,1); +md.materials.rheology_B=paterson(md.initialization.temperature); +md.friction.q=ones(md.mesh.numberofelements,1); +md.friction.p=ones(md.mesh.numberofelements,1); + +disp(' Set other boundary conditions'); +md.mask.ice_levelset(md.mesh.vertexonboundary==1)=0; +md.basalforcings.melting_rate = zeros(md.mesh.numberofvertices,1); +md.thermal.spctemperature = [md.initialization.temperature;1]; %impose observed temperature on surface +md.masstransport.spcthickness = NaN*ones(md.mesh.numberofvertices,1); + +disp(' Set geothermal heat flux'); +md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0); + +disp(' Set Pressure'); +md.initialization.pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness; + +disp(' Single point constraint for continental model'); +%Set at least one vertex to velocity 0 so as to not get a singular problem (point on the wet peninsula) +md.stressbalance.referential=NaN*ones(md.mesh.numberofvertices,6); +md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices,1); +md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices,1); +md.stressbalance.spcvz = NaN*ones(md.mesh.numberofvertices,1); +location = 1.0e+06 *[.32011 -2.2039]; +[dist pos]=min(sqrt((md.mesh.x - location(1)).^2 + (md.mesh.y - location(2)).^2)); +md.stressbalance.spcvx(pos) = 0; +md.stressbalance.spcvy(pos) = 0; +md.stressbalance.spcvz(pos) = 0; + Index: /issm/workshop/10IceBridge/Exp_Par/Jak_grounded.exp =================================================================== --- /issm/workshop/10IceBridge/Exp_Par/Jak_grounded.exp (revision 18013) +++ /issm/workshop/10IceBridge/Exp_Par/Jak_grounded.exp (revision 18013) @@ -0,0 +1,19 @@ +## Name:Jak_grouded +## Icon:0 +# Points Count Value +13 1.000000 +# X pos Y pos +-414397.6663804040 -2243466.6708395295 +-421486.9037621063 -2215731.3059818177 +-342860.8164377716 -2178965.8223332223 +-275835.2993744044 -2209281.2211311869 +-226210.6377024883 -2262816.9253914217 +-228144.0660793162 -2295067.3496445753 +-333193.6745536321 -2328607.7908678548 +-379595.9555975018 -2325382.7484425395 +-422775.8560133249 -2285392.2223686292 +-427287.1888925899 -2265396.9593316740 +-411819.7618779668 -2254431.8150856020 +-413108.7141291854 -2243466.6708395295 +-414397.6663804040 -2243466.6708395295 + Index: /issm/workshop/10IceBridge/Exp_Par/Jak_outline.exp =================================================================== --- /issm/workshop/10IceBridge/Exp_Par/Jak_outline.exp (revision 18013) +++ /issm/workshop/10IceBridge/Exp_Par/Jak_outline.exp (revision 18013) @@ -0,0 +1,17 @@ +## Name:Jak +## Icon:0 +# Points Count Value +11 1.000000 +# X pos Y pos +-421561.7247949215 -2216716.3488182197 +-386019.6416099078 -2203459.4896082953 +-351446.8879663036 -2185999.2360147359 +-289086.6874689616 -2214129.6445821370 +-255160.1535196303 -2282353.9688088228 +-340138.0433165266 -2317597.8140254519 +-383306.1449952667 -2302506.6581698526 +-416391.9672407377 -2279120.5885137189 +-448056.7322601135 -2246463.4475331730 +-429126.8716234095 -2231167.6696134652 +-421561.7247949215 -2216716.3488182197 + Index: /issm/workshop/10IceBridge/Exp_Par/data_gaps.exp =================================================================== --- /issm/workshop/10IceBridge/Exp_Par/data_gaps.exp (revision 18013) +++ /issm/workshop/10IceBridge/Exp_Par/data_gaps.exp (revision 18013) @@ -0,0 +1,70 @@ +## Name:data_gaps +## Icon:0 +# Points Count Value +13 1.000000 +# X pos Y pos +-183076.0390091268 -1794119.5324780210 +-233886.6983670916 -1838953.5447257366 +-239864.4229974404 -2006333.8571172084 +-299641.6693009282 -2081057.2108634012 +-269753.0461491842 -2116924.4206615738 +-257797.5968884868 -2185669.9061080711 +-278719.6330947075 -2248437.5232548732 +-254808.7345733123 -2266371.1281539593 +-224920.1114215683 -2167736.3012089850 +-198020.3505849990 -2093012.9474627920 +-123298.7927056390 -1967477.7131691882 +-180087.1766939524 -1797108.4666278686 +-183076.0390091268 -1794119.5324780210 + +## Name:data_gaps +## Icon:0 +# Points Count Value +12 1.000000 +# X pos Y pos +-90421.3072387208 -1824008.8739764981 +-117321.0680752902 -1946555.1741202541 +-156176.2781725575 -2075079.3425637058 +-177098.3143787782 -2203603.5110071572 +-174109.4520636038 -2314194.0745515227 +-72488.1333476745 -2293271.5355025884 +-93410.1695538952 -2203603.5110071572 +-30644.0609352330 -2206592.4451570050 +-6733.1624138379 -2081057.2108634012 +44077.4969441269 -1967477.7131691882 +-93410.1695538952 -1824008.8739764981 +-90421.3072387208 -1824008.8739764981 + +## Name:data_gaps +## Icon:0 +# Points Count Value +15 1.000000 +# X pos Y pos +-87432.4449235464 -2365005.9550989335 +-287686.2200402305 -2639987.8968849229 +-284697.3577250564 -2879102.6288727401 +-470006.8212658686 -2983715.3241174100 +-538750.6545148797 -3157073.5048085772 +-356430.0532892416 -3267664.0683529424 +-189053.7636394757 -3270653.0025027902 +-272741.9084643587 -3243752.5951541606 +-218942.3867912197 -3076372.2827626891 +-236875.5606822660 -3019582.5339155826 +-117321.0680752902 -2834268.6166250245 +20166.5984227317 -2669877.2383834003 +-12710.8870441865 -2469618.6503436035 +-87432.4449235464 -2367994.8892487814 +-87432.4449235464 -2365005.9550989335 + +## Name:data_gaps +## Icon:0 +# Points Count Value +6 1.000000 +# X pos Y pos +-186064.9013243013 -1555004.8004902040 +-177098.3143787782 -1695484.7055330465 +-81454.7202931975 -1698473.6396828941 +-36621.7855655816 -1513159.7223923360 +-186064.9013243013 -1552015.8663403564 +-186064.9013243013 -1555004.8004902040 + Index: /issm/workshop/10IceBridge/readjak.m =================================================================== --- /issm/workshop/10IceBridge/readjak.m (revision 18013) +++ /issm/workshop/10IceBridge/readjak.m (revision 18013) @@ -0,0 +1,20 @@ +clear all +disp('Reading Cresis''ss Jakobshavn '); + +rows=266400; + +A=dlmread('../Data/Jakobshavn_2008_2011_Composite_XYZGrid.txt'); + +y=A(:,1); +x=A(:,2); +thick=A(:,3); +surf=A(:,4); +bed=A(:,5); + +bed(bed < -9000)=NaN; +thick(thick < -9000)=NaN; +surf(surf < -9000)=NaN; + +save CresisBedJakComposite x y bed +save CresisThicknessJakComposite x y thick +save CresisSurfaceJakComposite x y surf Index: /issm/workshop/10IceBridge/runme.m =================================================================== --- /issm/workshop/10IceBridge/runme.m (revision 18013) +++ /issm/workshop/10IceBridge/runme.m (revision 18013) @@ -0,0 +1,155 @@ +clear all; +steps=[1]; + +%Location of SeaRISE dataset +ncdata='../Data/Greenland_5km_dev1.2.nc'; + +if any(steps==1) + disp(' Step 1: Mesh creation'); + + %Generate initial uniform mesh (resolution = 20000 m) + md=triangle(model,'./Exp_Par/DomainOutline.exp',20000); + + % Get velocities (Note: You can use ncdisp('file') to see an ncdump) + x1 = ncread(ncdata,'x1'); + y1 = ncread(ncdata,'y1'); + velx = ncread(ncdata,'surfvelx'); + vely = ncread(ncdata,'surfvely'); + vx = InterpFromGridToMesh(x1,y1,velx',md.mesh.x,md.mesh.y,0); + vy = InterpFromGridToMesh(x1,y1,vely',md.mesh.x,md.mesh.y,0); + vel = sqrt(vx.^2+vy.^2); + + %Mesh greenland without refinement in Jak basin + md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8); + return; + + %Refine mesh in the region of Jakobshavn (resolution = 3000 m) + hmaxVertices=NaN*ones(md.mesh.numberofvertices,1); + in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,... + './Exp_Par/Jak_outline.exp','node',1); + hmaxVertices(find(in))=3000; + md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,... + 'err',8,'hmaxVertices',hmaxVertices); + + %convert x,y coordinates (Polar stereo) to lat/lon + [md.mesh.lat,md.mesh.long]=xy2ll(md.mesh.x,md.mesh.y,+1,39,71); + + save ./Models/Greenland.Mesh_generation md; +end + +if any(steps==2) + disp(' Step 2: Parameterization'); + md = loadmodel('./Models/Greenland.Mesh_generation'); + + md = setmask(md,'',''); + md = parameterize(md,'./Exp_Par/Greenland.par'); + md = setflowequation(md,'SSA','all'); + + save ./Models/Greenland.Parameterization2 md; +end + +if any(steps==3) + disp(' Step 3: Control method friction'); + md = loadmodel('./Models/Greenland.Parameterization2'); + + %Control general + md.inversion.iscontrol=1; + md.inversion.nsteps=30; + md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1); + md.inversion.maxiter_per_step=5*ones(md.inversion.nsteps,1); + md.verbose=verbose('solution',true,'control',true); + + %Cost functions + md.inversion.cost_functions=[101 103 501]; + md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,3); + md.inversion.cost_functions_coefficients(:,1)=350; + md.inversion.cost_functions_coefficients(:,2)=0.6; + md.inversion.cost_functions_coefficients(:,3)=2e-6; + + %Controls + md.inversion.control_parameters={'FrictionCoefficient'}; + md.inversion.gradient_scaling(1:md.inversion.nsteps)=50; + md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1); + md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1); + in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,'./Exp_Par/data_gaps.exp','node',1); + md.inversion.cost_functions_coefficients(find(in),1)=0.0; + md.inversion.cost_functions_coefficients(find(in),2)=0.0; + + %Additional parameters + md.stressbalance.restol=0.01; + md.stressbalance.reltol=0.1; + md.stressbalance.abstol=NaN; + + %Go solve + md.cluster=generic('name',oshostname,'np',2); + + + md.verbose=verbose('solution',true,'control',true); + md=solve(md,StressbalanceSolutionEnum); + + %Update model friction fields accordingly + md.friction.coefficient=md.results.StressbalanceSolution.FrictionCoefficient; + + save ./Models/Greenland.Control_drag md; +end + +if any(steps==4) + disp(' Step 4: Transient run'); + md = loadmodel('./Models/Greenland.Control_drag'); + + %Set surface mass balance + x1 = ncread(ncdata,'x1'); + y1 = ncread(ncdata,'y1'); + smb = ncread(ncdata,'smb'); + smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0)*1000/md.materials.rho_ice; + smb = [smb smb smb-1.0]; + md.surfaceforcings.mass_balance = [smb;1 10 20]; + + %Set transient options, run for 20 years, saving every year + md.timestepping.time_step=0.2; %This must be reduced for finer resolutions + md.timestepping.final_time=20; + md.settings.output_frequency=1; + + %Additional options + md.inversion.iscontrol=0; + md.transient.requested_outputs={'IceVolume','TotalSmb','SurfaceforcingsMassBalance'}; + md.verbose=verbose('solution',true,'module',true,'convergence',true); + + %Go solve + md.cluster=generic('name',oshostname,'np',2); + md=solve(md,TransientSolutionEnum); + + save ./Models/Greenland.Transient md; +end + +if any(steps==5) + disp(' Step 5: Plotting'); + md = loadmodel('./Models/Greenland.Transient'); + + %Planview plots + plotmodel(md,'data',md.results.TransientSolution(end).Vel,'caxis',[1e-1 6000],... + 'log', 10, 'title', 'Velocity (m/y)','gridded',1, ... + 'data', md.results.TransientSolution(end).SurfaceforcingsMassBalance, ... + 'title', 'Surface mass balance (m/y)','gridded',1, ... + 'data',md.results.TransientSolution(end).Thickness,... + 'title','Thickness (m)','gridded',1, ... + 'data',md.results.TransientSolution(end).Surface, ... + 'title', 'Surface (m)','gridded',1); + + %Line Plots + figure + + %Plot surface mass balance + surfmb=[]; for i=1:100; surfmb=[surfmb ... + md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end + subplot(3,1,1); plot([0.2:0.2:20],mean(surfmb)); title('Mean Surface mass balance'); + + %Plot velocity + vel=[]; for i=1:100; vel=[vel md.results.TransientSolution(i).Vel]; end + subplot(3,1,2); plot([0.2:0.2:20],mean(vel)); title('Mean Velocity'); + + %Plot Volume + volume=[]; for i=1:100; volume=[volume md.results.TransientSolution(i).IceVolume]; end + subplot(3,1,3); plot([0.2:0.2:20],volume); title('Ice Volume'); + xlabel('years') +end Index: /issm/workshop/Data/Download.sh =================================================================== --- /issm/workshop/Data/Download.sh (revision 18012) +++ /issm/workshop/Data/Download.sh (revision 18013) @@ -1,5 +1,5 @@ #!/bin/bash -echo "Downloading SeaRISE dataset" +echo "Downloading SeaRISE dataset - Antarctica" $ISSM_DIR/scripts/DownloadExternalPackage2.py 'http://websrv.cs.umt.edu/isis/index.php/Present_Day_Antarctica/Antarctica_5km_withshelves_v0.75.nc' 'Antarctica_5km_withshelves_v0.75.nc' @@ -7,3 +7,16 @@ $ISSM_DIR/scripts/DownloadExternalPackage2.py 'ftp://sidads.colorado.edu/pub/DATASETS/nsidc0484_MEASURES_antarc_vel_V01/900m/antarctica_ice_velocity.nc' 'Antarctica_ice_velocity.nc' +echo "Downloading PIG control model" +$ISSM_DIR/scripts/DownloadExternalPackage2.py 'http://issm.jpl.nasa.gov/files/workshop2014/PIG.Control_drag_default' 'PIG.Control_drag_default' + +echo "Downloading PIG errors" +$ISSM_DIR/scripts/DownloadExternalPackage2.py 'http://issm.jpl.nasa.gov/files/workshop2014/CrossOvers2009.mat' 'CrossOvers2009.mat' + +echo "Downloading SeaRISE dataset - Greenland" +$ISSM_DIR/scripts/DownloadExternalPackage2.py 'http://websrv.cs.umt.edu/isis/index.php/Present_Day_Greenland/Greenland_5km_dev1.2.nc' 'Greenland_5km_dev1.2.nc' + echo "Downloading Jason Box's SMB" +$ISSM_DIR/scripts/DownloadExternalPackage2.py 'http://issm.jpl.nasa.gov/files/workshop2014/Box_Greenland_SMB_monthly_1840-2012_5km_ver20140421.nc' 'Box_Greenland_SMB_monthly_1840-2012_5km_ver20140421.nc' + +echo "Downloading IceBridge Jakobshavn bedrock" +$ISSM_DIR/scripts/DownloadExternalPackage2.py 'https://data.cresis.ku.edu/data/grids/old_versions/Jakobshavn_2008_2011_Composite.zip' 'Jakobshavn_2008_2011_Composite.zip'