[13280] | 1 | from model import *
|
---|
| 2 | from triangle import *
|
---|
| 3 | from setmask import *
|
---|
| 4 | from parameterize import *
|
---|
| 5 | from setflowequation import *
|
---|
| 6 | from EnumDefinitions import *
|
---|
| 7 | from solve import *
|
---|
| 8 | from MatlabFuncs import *
|
---|
| 9 |
|
---|
| 10 | md=triangle(model(),'../Exp/Square.exp',200000)
|
---|
| 11 | md=setmask(md,'all','')
|
---|
| 12 | md=parameterize(md,'../Par/SquareShelfConstrained.py')
|
---|
| 13 | md.extrude(3,1)
|
---|
| 14 | md=setflowequation(md,'pattyn','all')
|
---|
| 15 | md.cluster=generic('name',oshostname(),'np',3)
|
---|
| 16 | md.transient.requested_outputs=IceVolumeEnum()
|
---|
| 17 | md=solve(md,TransientSolutionEnum())
|
---|
| 18 |
|
---|
| 19 | #Fields and tolerances to track changes
|
---|
| 20 | field_names=[\
|
---|
| 21 | 'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface2','Thickness1','Temperature1','BasalforcingsMeltingRate1','Volume1', \
|
---|
| 22 | 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsMeltingRate2','Volume2', \
|
---|
| 23 | 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3','Volume3']
|
---|
| 24 | field_tolerances=[\
|
---|
| 25 | 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-13,\
|
---|
| 26 | 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,2e-12,\
|
---|
| 27 | 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-11]
|
---|
| 28 | field_values=[\
|
---|
| 29 | md.results['TransientSolution'][1]['Vx'],\
|
---|
| 30 | md.results['TransientSolution'][1]['Vy'],\
|
---|
| 31 | md.results['TransientSolution'][1]['Vz'],\
|
---|
| 32 | md.results['TransientSolution'][1]['Vel'],\
|
---|
| 33 | md.results['TransientSolution'][1]['Pressure'],\
|
---|
| 34 | md.results['TransientSolution'][1]['Bed'],\
|
---|
| 35 | md.results['TransientSolution'][1]['Surface'],\
|
---|
| 36 | md.results['TransientSolution'][1]['Thickness'],\
|
---|
| 37 | md.results['TransientSolution'][1]['Temperature'],\
|
---|
| 38 | md.results['TransientSolution'][1]['BasalforcingsMeltingRate'],\
|
---|
| 39 | md.results['TransientSolution'][1]['IceVolume'],\
|
---|
| 40 | md.results['TransientSolution'][2]['Vx'],\
|
---|
| 41 | md.results['TransientSolution'][2]['Vy'],\
|
---|
| 42 | md.results['TransientSolution'][2]['Vz'],\
|
---|
| 43 | md.results['TransientSolution'][2]['Vel'],\
|
---|
| 44 | md.results['TransientSolution'][2]['Pressure'],\
|
---|
| 45 | md.results['TransientSolution'][2]['Bed'],\
|
---|
| 46 | md.results['TransientSolution'][2]['Surface'],\
|
---|
| 47 | md.results['TransientSolution'][2]['Thickness'],\
|
---|
| 48 | md.results['TransientSolution'][2]['Temperature'],\
|
---|
| 49 | md.results['TransientSolution'][2]['BasalforcingsMeltingRate'],\
|
---|
| 50 | md.results['TransientSolution'][2]['IceVolume'],\
|
---|
| 51 | md.results['TransientSolution'][3]['Vx'],\
|
---|
| 52 | md.results['TransientSolution'][3]['Vy'],\
|
---|
| 53 | md.results['TransientSolution'][3]['Vz'],\
|
---|
| 54 | md.results['TransientSolution'][3]['Vel'],\
|
---|
| 55 | md.results['TransientSolution'][3]['Pressure'],\
|
---|
| 56 | md.results['TransientSolution'][3]['Bed'],\
|
---|
| 57 | md.results['TransientSolution'][3]['Surface'],\
|
---|
| 58 | md.results['TransientSolution'][3]['Thickness'],\
|
---|
| 59 | md.results['TransientSolution'][3]['Temperature'],\
|
---|
| 60 | md.results['TransientSolution'][3]['BasalforcingsMeltingRate'],\
|
---|
| 61 | md.results['TransientSolution'][3]['IceVolume'],\
|
---|
| 62 | ]
|
---|