| 7 |
|
|
| 8 |
from escript.modelframe import Link,Simulation |
from escript.modelframe import Link,Simulation |
| 9 |
from modellib.geometry import RectangularDomain,ScalarConstrainer,VectorConstrainer |
from modellib.geometry import RectangularDomain,ScalarConstrainer,VectorConstrainer |
| 10 |
from modellib.input import Sequencer,InterpolateOverBox,GausseanProfile,LinearCombination |
from modellib.input import Sequencer,InterpolateOverBox,GaussianProfile,LinearCombination |
| 11 |
from modellib.flow import SteadyIncompressibleFlow |
from modellib.flow import SteadyIncompressibleFlow |
| 12 |
from modellib.temperature import TemperatureAdvection |
from modellib.temperature import TemperatureAdvection |
| 13 |
from modellib.materials import SimpleEarthModel,GravityForce |
from modellib.materials import SimpleEarthModel,GravityForce |
| 89 |
vis.vector=Link(vel,"velocity") |
vis.vector=Link(vel,"velocity") |
| 90 |
vis.stride=5 |
vis.stride=5 |
| 91 |
|
|
| 92 |
per=GausseanProfile() |
per=GaussianProfile() |
| 93 |
per.domain=Link(dom) |
per.domain=Link(dom) |
| 94 |
per.x_c=[0.5,0.5,0.5] |
per.x_c=[0.5,0.5,0.5] |
| 95 |
per.A=0.0001 |
per.A=0.0001 |
| 103 |
lc.v1=Link(temp_val,"out") |
lc.v1=Link(temp_val,"out") |
| 104 |
temp.temperature=Link(lc,"out") |
temp.temperature=Link(lc,"out") |
| 105 |
|
|
| 106 |
s=Simulation([dom,sq,temp_constraints,vel_constraints,temp_val,per,lc,mat,grav,Simulation([vel],debug=True),temp,vis],debug=True) |
# s=Simulation([dom,sq,temp_constraints,vel_constraints,temp_val,per,lc,mat,grav,Simulation([vel],debug=True),temp,vis],debug=True) |
| 107 |
|
s=Simulation([dom,sq,Simulation([vel],debug=True),temp,vis],debug=True) |
| 108 |
s.writeXML() |
s.writeXML() |
| 109 |
s.run() |
s.run() |
Parent Directory
| 
Revision Log
| 
Patch