Index: /issm/trunk-jpl/examples/IceBridge/DomainOutline.exp =================================================================== --- /issm/trunk-jpl/examples/IceBridge/DomainOutline.exp (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/DomainOutline.exp (revision 18267) @@ -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/trunk-jpl/examples/IceBridge/Greenland.par =================================================================== --- /issm/trunk-jpl/examples/IceBridge/Greenland.par (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/Greenland.par (revision 18267) @@ -0,0 +1,104 @@ +%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'); +titles = fgets(fid); 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,... + './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.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1); +md.basalforcings.floatingice_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/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par =================================================================== --- /issm/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par (revision 18267) @@ -0,0 +1,108 @@ +%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'); +titles = fgets(fid); 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,... + './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.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1); +md.basalforcings.floatingice_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/trunk-jpl/examples/IceBridge/Greenland_noOIB.par =================================================================== --- /issm/trunk-jpl/examples/IceBridge/Greenland_noOIB.par (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/Greenland_noOIB.par (revision 18267) @@ -0,0 +1,82 @@ +%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.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1); +md.basalforcings.floatingice_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/trunk-jpl/examples/IceBridge/Jak_grounded.exp =================================================================== --- /issm/trunk-jpl/examples/IceBridge/Jak_grounded.exp (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/Jak_grounded.exp (revision 18267) @@ -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/trunk-jpl/examples/IceBridge/Jak_outline.exp =================================================================== --- /issm/trunk-jpl/examples/IceBridge/Jak_outline.exp (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/Jak_outline.exp (revision 18267) @@ -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/trunk-jpl/examples/IceBridge/data_gaps.exp =================================================================== --- /issm/trunk-jpl/examples/IceBridge/data_gaps.exp (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/data_gaps.exp (revision 18267) @@ -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/trunk-jpl/examples/IceBridge/runme.m =================================================================== --- /issm/trunk-jpl/examples/IceBridge/runme.m (revision 18267) +++ /issm/trunk-jpl/examples/IceBridge/runme.m (revision 18267) @@ -0,0 +1,159 @@ +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,'./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); + save ./Models/Greenland.Mesh_generation_old md; + 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,... + './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,'./Greenland_cheatsheet.par'); + %md = parameterize(md,'./Greenland_noOIB.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,... + './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/trunk-jpl/examples/Pig/runme.m =================================================================== --- /issm/trunk-jpl/examples/Pig/runme.m (revision 18266) +++ /issm/trunk-jpl/examples/Pig/runme.m (revision 18267) @@ -188,14 +188,18 @@ % Load Model - - % Disable inversion + md = loadmodel('./Models/PIG.Control_drag'); + md.inversion.iscontrol=0; + + disp(' Extruding mesh') + number_of_layers=3; + md=extrude(md,number_of_layers,0.9); % Extrude Mesh - % Set Flowequation - % Solve + md=solve(md,StressbalanceSolutionEnum); % Save Model + save ./Models/PIG.ModelHO md; end @@ -215,14 +219,14 @@ 'ylim#all',[min(mdHO.mesh.y) max(mdHO.mesh.y)]/10^3,... 'FontSize#all',12,... - 'data',mdHO.initialization.vel,'title','Observed velocity',... + 'data',mdSSA.initialization.vel,'title','Observed velocity',... 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdHO.initialization.vel(surf)),'title','(HO-observed) velocities',... 'data',mdSSA.results.StressbalanceSolution.Vel,'title','Modeled SSA Velocity',... 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdSSA.results.StressbalanceSolution.Vel),'title','(HO-SSA) velocities',... - 'data',mdHO.results.StressbalanceSolution.Vel,'title','Modeled HO surface Velocities',... + 'data',project2d(mdHO,mdHO.results.StressbalanceSolution.Vel,1),'title','Modeled HO surface Velocities',... 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdHO.results.StressbalanceSolution.Vel(basal)),'title','(HOsurf-HO base) velocities',... 'caxis#1',([1.5,4000]),'caxis#3',([1.5,4000]),'caxis#5',([1.5,4000]),... 'colorbar#all','on','view#all',2,... 'colorbartitle#all','[m/yr]',... - 'layer#5',1, 'log#1', 10,'log#3', 10,'log#5', 10); -end -% }}} + 'log#1', 10,'log#3', 10,'log#5', 10,'gridded#all',1); +end +% }}} Index: /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m =================================================================== --- /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m (revision 18266) +++ /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m (revision 18267) @@ -1,4 +1,4 @@ %PIG Uncertainty Quantification Application -steps=[1:3]; +steps=[7]; if any(steps==1) @@ -162,9 +162,9 @@ %variables + md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05); 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