/[escript]/trunk/esys2/doc/user/examples/diffusion.py
ViewVC logotype

Diff of /trunk/esys2/doc/user/examples/diffusion.py

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

revision 112 by jgs, Thu Jan 27 06:21:59 2005 UTC revision 113 by jgs, Mon Feb 28 07:06:33 2005 UTC
# Line 1  Line 1 
1  # $Id$  # $Id$
2  from mytools import Helmholtz  from mytools import Helmholtz
3  from esys.escript import Lsup  from esys.escript import Lsup,length
4  from esys.finley import Rectangle  from esys.finley import Rectangle
5  #... set some parameters ...  #... set some parameters ...
6  xc=[0.02,0.002]  xc=[0.02,0.002]
# Line 31  while t<tend: Line 31  while t<tend:
31        print "time step :",t        print "time step :",t
32        mypde.setValue(kappa=kappa,omega=rhocp/h,f=q+rhocp/h*T,eta=eta,g=eta*Tref)        mypde.setValue(kappa=kappa,omega=rhocp/h,f=q+rhocp/h*T,eta=eta,g=eta*Tref)
33        T=mypde.getSolution()        T=mypde.getSolution()
       T.saveDX("T%d.dx"%i)  
34          T.saveDX("T%d.dx"%i)

Legend:
Removed from v.112  
changed lines
  Added in v.113

  ViewVC Help
Powered by ViewVC 1.1.26