[9023] | 1 | md=mesh(model,'../Exp/Square.exp',180000);
|
---|
| 2 | md=geography(md,'all','');
|
---|
| 3 | md=parameterize(md,'../Par/SquareShelf.par');
|
---|
| 4 | md=extrude(md,3,1);
|
---|
| 5 | md=setelementstype(md,'macayeal','all');
|
---|
| 6 | md.cluster=none;
|
---|
| 7 | md.spctemperature=[md.spctemperature, md.spctemperature+5, md.spctemperature+10, md.spctemperature+15; 1.5 2.5 3.5 4];
|
---|
| 8 | md.dt=1;
|
---|
| 9 | md.ndt=4;
|
---|
| 10 | md=solve(md,ThermalSolutionEnum);
|
---|
| 11 |
|
---|
| 12 | %Fields and tolerances to track changes
|
---|
[9611] | 13 | field_names ={'Temperature1','BasalforcingsMeltingRate1','Temperature2','BasalforcingsMeltingRate2','Temperature3','BasalforcingsMeltingRate3','Temperature4','BasalforcingsMeltingRate4'};
|
---|
[9023] | 14 | field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
|
---|
| 15 | field_values={...
|
---|
| 16 | PatchToVec(md.results.ThermalSolution(1).Temperature),...
|
---|
[9611] | 17 | PatchToVec(md.results.ThermalSolution(1).BasalforcingsMeltingRate),...
|
---|
[9023] | 18 | PatchToVec(md.results.ThermalSolution(2).Temperature),...
|
---|
[9611] | 19 | PatchToVec(md.results.ThermalSolution(2).BasalforcingsMeltingRate),...
|
---|
[9023] | 20 | PatchToVec(md.results.ThermalSolution(3).Temperature),...
|
---|
[9611] | 21 | PatchToVec(md.results.ThermalSolution(3).BasalforcingsMeltingRate),...
|
---|
[9023] | 22 | PatchToVec(md.results.ThermalSolution(4).Temperature),...
|
---|
[9611] | 23 | PatchToVec(md.results.ThermalSolution(4).BasalforcingsMeltingRate),...
|
---|
[9023] | 24 | };
|
---|