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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 987 - (hide annotations)
Tue Feb 20 06:19:30 2007 UTC (12 years, 8 months ago) by gross
File MIME type: text/x-python
File size: 1733 byte(s)
file missing
1 gross 986 """
2     mining activities in modelframe
3    
4     @var __author__: name of author
5     @var __licence__: licence agreement
6     @var __url__: url entry point on documentation
7     @var __version__: version
8     @var __date__: date of the version
9     """
10    
11     __copyright__=""" Copyright (c) 2006 by ACcESS MNRF
12     http://www.access.edu.au
13     Primary Business: Queensland, Australia"""
14     __license__="""Licensed under the Open Software License version 3.0
15     http://www.opensource.org/licenses/osl-3.0.php"""
16     __author__="Lutz Gross, l.gross@uq.edu.au"
17     __url__="http://www.iservo.edu.au/esys/escript"
18     __version__="$Revision$"
19     __date__="$Date$"
20    
21     d=True
22     from setups import MiningHistory
23     from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
24     from esys.modellib.input import Sequencer
25     from esys.escript.modelframe import Link,Simulation, DataSource
26    
27     dom=FinleyReader(debug=d)
28     dom.source=DataSource("./newcastle_mines.msh","gmsh")
29     dom.region_tag_map_source=DataSource("./vtag.xml", "ESysXML")
30     dom.surface_tag_map_source=DataSource("./vtag.xml", "ESysXML")
31    
32     sq=Sequencer(debug=d)
33     sq.t=1840.
34     sq.t_end=2000.
35     sq.dt_max=100.
36    
37     hist=MiningHistory(debug=d)
38     hist.history=DataSource("./newcastle_mining.xml")
39     hist.t=Link(sq,"t")
40    
41     # hist.mine_locations=Link(dom,"region_tag_map")
42     # hist.domain=Link(dom,"domain")
43    
44     fix=VectorConstrainerOverBox(debug=d)
45     fix.domain=Link(dom,"domain")
46     fix.value=0.
47     fix.top=False
48     fix.bottom=True
49     fix.front=False
50     fix.back=False
51     fix.left=False
52     fix.right=False
53    
54    
55     s=Simulation([sq, hist, fix], debug=d)
56     s.run()

  ViewVC Help
Powered by ViewVC 1.1.26