Hi Mo, I tried running your code but I got "name 'Pf' is not defined" -- a different error.
On Fri, May 6, 2016 at 12:41 AM, Mohammad Kazemi <kazem...@gmail.com> wrote: > > CODE: > ------------------------------------------------ > from __future__ import division > from numpy import * > from pylab import * > from matplotlib import rc > from fipy import * > > #################### INITIALLIZATION > ########################################### > L = 0.0508 #+0.0021+0.0014 > nx = 10.0 > dx = L / nx > Pdi=6894757.29 > Pui=7170547.58 > D=1e-6 > tmax = 1000 > sigma = 400 > k0 = 10e-16 > T = 300 > M = 0.0399 > R = 8.314 > muf=1.96e-5 > hf=30e-6 > phif=0.25 > A1=1.14 > A2=0.28 > wkf=4 > wkd=4 > wku=4 > Kn=sqrt((3.14*R*T/(2*M))*muf/(Pf*hf)) At this point "Pf" is not defined. -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]