--- trunk/esys2/doc/user/examples/diffusion.py	2004/12/15 07:08:39	102
+++ trunk/esys2/doc/user/examples/diffusion.py	2005/01/27 06:21:59	108
@@ -1,34 +1,34 @@
 # $Id$
-from mytools import *
-from esys.escript import *
-import esys.finley
+from mytools import Helmholtz
+from esys.escript import Lsup
+from esys.finley import Rectangle
 #... set some parameters ...
-x_c=[0.02,0.002]
+xc=[0.02,0.002]
 r=0.001
-q0=50.e6
+qc=50.e6
 Tref=0.
 rhocp=2.6e6
 eta=75.
 kappa=240.
-t_end=5.
-# ...time step size and counter ...
+tend=5.
+# ... time, time step size and counter ...
+t=0
 h=0.1
 i=0
-t=0
 #... generate domain ...
-mydomain = esys.finley.Rectangle(l0=0.05,l1=0.01,n0=250, n1=50)
+mydomain = Rectangle(l0=0.05,l1=0.01,n0=250, n1=50)
 #... open PDE ...
 mypde=Helmholtz(mydomain)
 # ... set heat source: ....
 x=mydomain.getX()
-q=q0*(length(x-x_c)-r).whereNegative()
+q=qc*(length(x-xc)-r).whereNegative()
 # ... set initial temperature ....
 T=Tref
 # ... start iteration:
-while t<t_end:
+while t<tend:
       i+=1
       t+=h
       print "time step :",t
       mypde.setValue(kappa=kappa,omega=rhocp/h,f=q+rhocp/h*T,eta=eta,g=eta*Tref)
       T=mypde.getSolution()
-      T.saveDX("T%d.dx"%i)
+      T.saveDX("T%d.dx"%i)
\ No newline at end of file