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() |