Changeset 24313 for issm/trunk/test/NightlyRun/test275.py
- Timestamp:
- 11/01/19 12:01:57 (5 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
-
issm/trunk/test
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test275.py
r21729 r24313 1 #Test Name: SquareShelfDamageEvolutionSSA2dPralong 1 #Test Name: SquareShelfDamageEvolutionSSA2dPralong 2 2 import numpy as np 3 3 from triangle import triangle … … 6 6 from setmask import setmask 7 7 from parameterize import parameterize 8 from verbose import verbose9 8 from setflowequation import setflowequation 10 9 from generic import generic 11 from socket import gethostname12 10 from solve import solve 13 11 from matdamageice import matdamageice 14 12 15 md =triangle(model(),'../Exp/Square.exp',150000.)16 md =setmask(md,'all','')17 md.materials =matdamageice()18 md =parameterize(md,'../Par/SquareShelf.py')19 md.damage.isdamage =120 md.damage.D =0.1*np.ones(md.mesh.numberofvertices)21 md.damage.spcdamage =np.nan*np.ones(md.mesh.numberofvertices)22 md.damage.law =113 md = triangle(model(), '../Exp/Square.exp', 150000.) 14 md = setmask(md, 'all', '') 15 md.materials = matdamageice() 16 md = parameterize(md, '../Par/SquareShelf.py') 17 md.damage.isdamage = 1 18 md.damage.D = 0.1 * np.ones(md.mesh.numberofvertices) 19 md.damage.spcdamage = np.nan * np.ones(md.mesh.numberofvertices) 20 md.damage.law = 1 23 21 24 md.damage.c1 =1.e-1125 md.damage.c2 =0.426 md.damage.c3 =1.e-327 md.damage.healing =0.428 md.damage.stress_threshold =1.e529 md.damage.stabilization =122 md.damage.c1 = 1.e-11 23 md.damage.c2 = 0.4 24 md.damage.c3 = 1.e-3 25 md.damage.healing = 0.4 26 md.damage.stress_threshold = 1.e5 27 md.damage.stabilization = 1 30 28 31 md.damage.requested_outputs =['default','DamageF']29 md.damage.requested_outputs = ['default', 'DamageF'] 32 30 33 md =setflowequation(md,'SSA','all')34 md.cluster =generic('name',gethostname(),'np',3)35 md =solve(md,'DamageEvolution')31 md = setflowequation(md, 'SSA', 'all') 32 md.cluster = generic('name', gethostname(), 'np', 3) 33 md = solve(md, 'DamageEvolution') 36 34 37 field_names =['D','F']38 field_tolerances =[1.e-13,1.e-13]39 field_values =[md.results.DamageEvolutionSolution.DamageDbar,40 35 field_names = ['D', 'F'] 36 field_tolerances = [1.e-13, 1.e-13] 37 field_values = [md.results.DamageEvolutionSolution.DamageDbar, 38 md.results.DamageEvolutionSolution.DamageF]
Note:
See TracChangeset
for help on using the changeset viewer.