| 1 | #Test Name: SquareSheetShelfTranSSA3dAggressiveRegionalOutput
|
|---|
| 2 | import numpy as np
|
|---|
| 3 | from model import *
|
|---|
| 4 | from socket import gethostname
|
|---|
| 5 | from triangle import triangle
|
|---|
| 6 | from setmask import setmask
|
|---|
| 7 | from parameterize import parameterize
|
|---|
| 8 | from setflowequation import setflowequation
|
|---|
| 9 | from solve import solve
|
|---|
| 10 | from frictioncoulomb import frictioncoulomb
|
|---|
| 11 | from ContourToMesh import ContourToMesh
|
|---|
| 12 | from regionaloutput import regionaloutput
|
|---|
| 13 |
|
|---|
| 14 | md = triangle(model(), '../Exp/Square.exp', 150000.)
|
|---|
| 15 | md = setmask(md, '../Exp/SquareShelf.exp', '')
|
|---|
| 16 | md = parameterize(md, '../Par/SquareSheetShelf.py')
|
|---|
| 17 | md.geometry.bed = copy.deepcopy(md.geometry.base)
|
|---|
| 18 | pos = np.nonzero(md.mask.ocean_levelset < 0.)
|
|---|
| 19 | md.geometry.bed[pos] = md.geometry.bed[pos] - 10
|
|---|
| 20 | md.friction = frictioncoulomb()
|
|---|
| 21 | md.friction.coefficient = 20 * np.ones(md.mesh.numberofvertices)
|
|---|
| 22 | md.friction.p = 1 * np.ones(md.mesh.numberofelements)
|
|---|
| 23 | md.friction.q = 1 * np.ones(md.mesh.numberofelements)
|
|---|
| 24 | md.friction.coefficientcoulomb = 0.02 * np.ones(md.mesh.numberofvertices)
|
|---|
| 25 | md.transient.isthermal = False
|
|---|
| 26 | md.transient.isgroundingline = True
|
|---|
| 27 | md.groundingline.migration = 'AggressiveMigration'
|
|---|
| 28 | md.mesh.scale_factor = 1.1 * np.ones((md.mesh.numberofvertices))
|
|---|
| 29 |
|
|---|
| 30 | md.settings.output_frequency = 2
|
|---|
| 31 | md = setflowequation(md, 'SSA', 'all')
|
|---|
| 32 | md.cluster = generic('name', gethostname(), 'np', 3)
|
|---|
| 33 |
|
|---|
| 34 | regionalmask = np.zeros((md.mesh.numberofvertices))
|
|---|
| 35 | inflag = ContourToMesh(md.mesh.elements, md.mesh.x, md.mesh.y, '../Exp/SquareHalfRight.exp', 'node', 1)
|
|---|
| 36 | regionalmask[np.nonzero(inflag)[0]] = 1.
|
|---|
| 37 | md.transient.requested_outputs = ['default', 'GroundedArea1', 'FloatingArea1', 'TotalFloatingBmb1', 'TotalGroundedBmb1', 'TotalSmb1', 'IceMass1', 'IceVolume1', 'IceVolumeAboveFloatation1', 'IceVolumeAboveFloatation2', 'IceVolumeAboveFloatation', 'IceMassScaled1', 'IceVolumeScaled1', 'IceVolumeAboveFloatationScaled1', 'IceVolumeAboveFloatationScaled2']
|
|---|
| 38 |
|
|---|
| 39 | md.outputdefinition.definitions = [regionaloutput('name', 'GroundedArea1', 'outputnamestring', 'GroundedArea', 'mask', regionalmask, 'definitionstring', 'Outputdefinition1'),
|
|---|
| 40 | regionaloutput('name', 'FloatingArea1', 'outputnamestring', 'FloatingArea', 'mask', regionalmask, 'definitionstring', 'Outputdefinition2'),
|
|---|
| 41 | regionaloutput('name', 'TotalFloatingBmb1', 'outputnamestring', 'TotalFloatingBmb', 'mask', regionalmask, 'definitionstring', 'Outputdefinition3'),
|
|---|
| 42 | regionaloutput('name', 'TotalGroundedBmb1', 'outputnamestring', 'TotalGroundedBmb', 'mask', regionalmask, 'definitionstring', 'Outputdefinition4'),
|
|---|
| 43 | regionaloutput('name', 'IceMass1', 'outputnamestring', 'IceMass', 'mask', regionalmask, 'definitionstring', 'Outputdefinition5'),
|
|---|
| 44 | regionaloutput('name', 'IceVolume1', 'outputnamestring', 'IceVolume', 'mask', regionalmask, 'definitionstring', 'Outputdefinition6'),
|
|---|
| 45 | regionaloutput('name', 'IceVolumeAboveFloatation1', 'outputnamestring', 'IceVolumeAboveFloatation', 'mask', regionalmask, 'definitionstring', 'Outputdefinition7'),
|
|---|
| 46 | regionaloutput('name', 'TotalSmb1', 'outputnamestring', 'TotalSmb', 'mask', regionalmask, 'definitionstring', 'Outputdefinition8'),
|
|---|
| 47 | regionaloutput('name', 'IceVolumeAboveFloatation2', 'outputnamestring', 'IceVolumeAboveFloatation', 'maskexpstring', '../Exp/SquareHalfRight.exp', 'definitionstring', 'Outputdefinition9', 'model', md),
|
|---|
| 48 | regionaloutput('name', 'IceMassScaled1', 'outputnamestring', 'IceMassScaled', 'mask', regionalmask, 'definitionstring', 'Outputdefinition10'),
|
|---|
| 49 | regionaloutput('name', 'IceVolumeScaled1', 'outputnamestring', 'IceVolumeScaled', 'mask', regionalmask, 'definitionstring', 'Outputdefinition11'),
|
|---|
| 50 | regionaloutput('name', 'IceVolumeAboveFloatationScaled1', 'outputnamestring', 'IceVolumeAboveFloatationScaled', 'mask', regionalmask, 'definitionstring', 'Outputdefinition12'),
|
|---|
| 51 | regionaloutput('name', 'IceVolumeAboveFloatationScaled2', 'outputnamestring', 'IceVolumeAboveFloatationScaled', 'maskexpstring', '../Exp/SquareHalfRight.exp', 'definitionstring', 'Outputdefinition13', 'model', md)]
|
|---|
| 52 |
|
|---|
| 53 | md.extrude(3, 1.)
|
|---|
| 54 | md2 = copy.deepcopy(md)
|
|---|
| 55 | md2.collapse()
|
|---|
| 56 | md = solve(md, 'Transient')
|
|---|
| 57 |
|
|---|
| 58 | #Fields and tolerances to track changes
|
|---|
| 59 | field_names = ['IceMass1', 'IceVolume1', 'IceVolumeAboveFloatation1', 'IceVolumeAboveFloatation21', 'Thickness1', 'GroundedArea1', 'FloatingArea1', 'TotalFloatingBmb1', 'TotalGroundedBmb1', 'TotalSmb1', 'IceMassScaled1', 'IceVolumeScaled1', 'IceVolumeAboveFloatationScaled1', 'IceVolumeAboveFloatationScaled21', 'IceMass3', 'IceVolume3', 'IceVolumeAboveFloatation3', 'IceVolumeAboveFloatation23', 'Thickness3', 'GroundedArea3', 'FloatingArea3', 'TotalFloatingBmb3', 'TotalGroundedBmb3', 'TotalSmb3', 'IceMassScaled3', 'IceVolumeScaled3', 'IceVolumeAboveFloatationScaled3', 'IceVolumeAboveFloatationScaled23', 'ExtrudedMask', 'CollapsedMask']
|
|---|
| 60 | field_tolerances = [1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13]
|
|---|
| 61 | field_values = [md.results.TransientSolution[0].IceMass1,
|
|---|
| 62 | md.results.TransientSolution[0].IceVolume1,
|
|---|
| 63 | md.results.TransientSolution[0].IceVolumeAboveFloatation1,
|
|---|
| 64 | md.results.TransientSolution[0].IceVolumeAboveFloatation2,
|
|---|
| 65 | md.results.TransientSolution[0].Thickness,
|
|---|
| 66 | md.results.TransientSolution[0].GroundedArea1,
|
|---|
| 67 | md.results.TransientSolution[0].FloatingArea1,
|
|---|
| 68 | md.results.TransientSolution[0].TotalFloatingBmb1,
|
|---|
| 69 | md.results.TransientSolution[0].TotalGroundedBmb1,
|
|---|
| 70 | md.results.TransientSolution[0].TotalSmb1,
|
|---|
| 71 | md.results.TransientSolution[0].IceMassScaled1,
|
|---|
| 72 | md.results.TransientSolution[0].IceVolumeScaled1,
|
|---|
| 73 | md.results.TransientSolution[0].IceVolumeAboveFloatationScaled1,
|
|---|
| 74 | md.results.TransientSolution[0].IceVolumeAboveFloatationScaled2,
|
|---|
| 75 | md.results.TransientSolution[2].IceMass1,
|
|---|
| 76 | md.results.TransientSolution[2].IceVolume1,
|
|---|
| 77 | md.results.TransientSolution[2].IceVolumeAboveFloatation1,
|
|---|
| 78 | md.results.TransientSolution[2].IceVolumeAboveFloatation2,
|
|---|
| 79 | md.results.TransientSolution[2].Thickness,
|
|---|
| 80 | md.results.TransientSolution[2].GroundedArea1,
|
|---|
| 81 | md.results.TransientSolution[2].FloatingArea1,
|
|---|
| 82 | md.results.TransientSolution[2].TotalFloatingBmb1,
|
|---|
| 83 | md.results.TransientSolution[2].TotalGroundedBmb1,
|
|---|
| 84 | md.results.TransientSolution[2].TotalSmb1,
|
|---|
| 85 | md.results.TransientSolution[2].IceMassScaled1,
|
|---|
| 86 | md.results.TransientSolution[2].IceVolumeScaled1,
|
|---|
| 87 | md.results.TransientSolution[2].IceVolumeAboveFloatationScaled1,
|
|---|
| 88 | md.results.TransientSolution[2].IceVolumeAboveFloatationScaled2,
|
|---|
| 89 | len(md.outputdefinition.definitions[0].mask),
|
|---|
| 90 | len(md2.outputdefinition.definitions[0].mask)]
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
