/[escript]/trunk/modellib/test/python/run_convection.py
ViewVC logotype

Diff of /trunk/modellib/test/python/run_convection.py

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

revision 727 by gross, Wed Apr 12 23:58:02 2006 UTC revision 728 by gross, Mon May 15 02:01:18 2006 UTC
# Line 20  from esys.modellib.temperature import Te Line 20  from esys.modellib.temperature import Te
20  from esys.modellib.materials import SimpleEarthModel,GravityForce  from esys.modellib.materials import SimpleEarthModel,GravityForce
21  from esys.modellib.visualization import WriteVTK  from esys.modellib.visualization import WriteVTK
22    
23    try:
24       WORKDIR=os.environ['MODELLIB_WORKDIR']
25    except KeyError:
26       WORKDIR='.'
27    
28    
29  dom=RectangularDomain()  dom=RectangularDomain()
30  dom.order=2  dom.order=2
31    
# Line 88  vis=WriteVTK() Line 94  vis=WriteVTK()
94  vis.t=Link(sq)  vis.t=Link(sq)
95  vis.scalar=Link(temp,"temperature")  vis.scalar=Link(temp,"temperature")
96  vis.vector=Link(vel,"velocity")  vis.vector=Link(vel,"velocity")
97  vis.stride=5  vis.dt=0.0005
98  vis.filename=os.environ['MODELLIB_WORKDIR']+"/temp.xml"  vis.filename=WORKDIR+"/temp.xml"
99    
100  per=GaussianProfile()  per=GaussianProfile()
101  per.domain=Link(dom)  per.domain=Link(dom)
# Line 105  lc.f1=1. Line 111  lc.f1=1.
111  lc.v1=Link(temp_val,"out")  lc.v1=Link(temp_val,"out")
112  temp.temperature=Link(lc,"out")  temp.temperature=Link(lc,"out")
113    
114  # s=Simulation([dom,sq,temp_constraints,vel_constraints,temp_val,per,lc,mat,grav,Simulation([vel],debug=True),temp,vis],debug=True)  s=Simulation([sq,Simulation([vel],debug=True),temp,vis],debug=True)
 s=Simulation([dom,sq,Simulation([vel],debug=True),temp,vis],debug=True)  
115  s.writeXML()  s.writeXML()
116  s.run()  s.run()
Legend:
Removed from v.727  
changed lines
  Added in v.728
  ViewVC Help
Powered by ViewVC 1.1.26