/[escript]/trunk/modellib/py_src/crustal/newcastle.py
ViewVC logotype

Diff of /trunk/modellib/py_src/crustal/newcastle.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 996 by gross, Tue Feb 20 06:19:30 2007 UTC revision 997 by gross, Mon Feb 26 06:31:45 2007 UTC
# Line 19  __version__="$Revision$" Line 19  __version__="$Revision$"
19  __date__="$Date$"  __date__="$Date$"
20    
21  d=True  d=True
22  from setups import MiningHistory  from setups import MiningHistory, DensityChange, LinearElasticStressChange
23  from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox  from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
24  from esys.modellib.input import Sequencer  from esys.modellib.input import Sequencer
25  from esys.escript.modelframe import Link,Simulation, DataSource  from esys.escript.modelframe import Link,Simulation, DataSource
26    
27  dom=FinleyReader(debug=d)  dom=FinleyReader(debug=d)
28  dom.source=DataSource("./newcastle_mines.msh","gmsh")  dom.source=DataSource("./newcastle_mines.msh","gmsh")
29  dom.region_tag_map_source=DataSource("./vtag.xml", "ESysXML")  dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
 dom.surface_tag_map_source=DataSource("./vtag.xml", "ESysXML")  
30    
31  sq=Sequencer(debug=d)  sq=Sequencer(debug=d)
32  sq.t=1840.  sq.t=1840.
# Line 38  hist=MiningHistory(debug=d) Line 37  hist=MiningHistory(debug=d)
37  hist.history=DataSource("./newcastle_mining.xml")  hist.history=DataSource("./newcastle_mining.xml")
38  hist.t=Link(sq,"t")  hist.t=Link(sq,"t")
39    
40  # hist.mine_locations=Link(dom,"region_tag_map")  dens_dot=DensityChange(debug=d)
41  # hist.domain=Link(dom,"domain")  dens_dot.domain=Link(dom,"domain")
42    dens_dot.tag_map=Link(dom,"tag_map")
43    dens_dot.mass_rate=Link(hist,"mass_changes")
44    
45  fix=VectorConstrainerOverBox(debug=d)  fix=VectorConstrainerOverBox(debug=d)
46  fix.domain=Link(dom,"domain")  fix.domain=Link(dom,"domain")
# Line 50  fix.front=False Line 51  fix.front=False
51  fix.back=False  fix.back=False
52  fix.left=False  fix.left=False
53  fix.right=False  fix.right=False
54    # dens_dot.density_rate=Link(hist,"mass_changes")
55    
56    el=LinearElasticStressChange(debug=d)
57    el.domain=Link(dom,"domain")
58    el.tag_map=Link(dom,"tag_map")
59    el.density=1.
60    el.lame_lambda=2.
61    el.lame_mu=1.
62    el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
63    el.density_rate=Link(dens_dot,"density_rate")
64    
65    # hist.domain=Link(dom,"domain")
66    
67    
68  s=Simulation([sq, hist, fix], debug=d)  s=Simulation([sq, hist, dens_dot, fix, el], debug=d)
69  s.run()  s.run()

Legend:
Removed from v.996  
changed lines
  Added in v.997

  ViewVC Help
Powered by ViewVC 1.1.26