- Timestamp:
- 06/07/17 10:50:54 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/branches/trunk-larour-NatGeoScience2016/test/NightlyRun/test273.py
r21056 r21759 1 1 #Test Name: SquareShelfStressSSA2dDamageUpdate 2 import numpy 3 from model import model 2 import numpy as np 3 from model import * 4 from socket import gethostname 4 5 from triangle import triangle 5 6 from setmask import setmask … … 9 10 from matdamageice import matdamageice 10 11 from generic import generic 11 from socket import gethostname as oshostname12 from socket import gethostname 12 13 13 14 md=triangle(model(),'../Exp/Square.exp',150000.) … … 16 17 md=parameterize(md,'../Par/SquareShelf.py') 17 18 md.damage.isdamage=1 18 md.damage.D=0.*n umpy.ones(md.mesh.numberofvertices)19 md.damage.spcdamage=n umpy.nan*numpy.ones(md.mesh.numberofvertices)19 md.damage.D=0.*np.ones(md.mesh.numberofvertices) 20 md.damage.spcdamage=np.nan*np.ones(md.mesh.numberofvertices) 20 21 md=setflowequation(md,'SSA','all') 21 md.cluster=generic('name', oshostname(),'np',3)22 md.cluster=generic('name',gethostname(),'np',3) 22 23 23 24 md.stressbalance.requested_outputs=['default','NewDamage'] … … 30 31 field_names =['Vx','Vy','Vel','Pressure','NewDamage'] 31 32 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13] 32 field_values=[\ 33 md.results.StressbalanceSolution.Vx,\ 34 md.results.StressbalanceSolution.Vy,\ 35 md.results.StressbalanceSolution.Vel,\ 36 md.results.StressbalanceSolution.Pressure,\ 37 md.results.StressbalanceSolution.NewDamage,\ 38 ] 33 field_values=[md.results.StressbalanceSolution.Vx, 34 md.results.StressbalanceSolution.Vy, 35 md.results.StressbalanceSolution.Vel, 36 md.results.StressbalanceSolution.Pressure, 37 md.results.StressbalanceSolution.NewDamage]
Note:
See TracChangeset
for help on using the changeset viewer.