It works well after I reinstalled fipy.
Thank you so much Zebo 2015-09-10 15:33 GMT-05:00 Keller, Trevor <[email protected]>: > >>> print C_V((1.0,1.0)) > > [0. 0.] > > >>> print C_V_xy((20,20)) > > 0.0 > > > Works for me: the problem comes after the initialization. > > Please provide a minimal working example, or dig into why C_V_xy gives a > scalar value instead of the array you're expecting. > > Trevor Keller, Ph.D. > Materials Science and Engineering Division > National Institute of Standards and Technology > 100 Bureau Dr. MS 8550; Gaithersburg, MD 20899 > Office: 223/A131 or (301) 975-2889 > > > > > ------------------------------ > *From:* [email protected] <[email protected]> on behalf of Zebo > LI <[email protected]> > *Sent:* Thursday, September 10, 2015 4:07 PM > *To:* FIPY > *Subject:* Extract data from fipy variables > > I defined two meshes in the following way, > (1) > nLnR = 100 > nCta = 50 > dLnR = 0.02 > dCta = 2*np.pi/nCta > PL = nLnR*dLnR > ##Create mesh in Polar Coordinate system (LnR,Cta) > Pmesh = PeriodicGrid2DTopBottom(dx=dLnR, dy=dCta, nx=nLnR, ny=nCta) > (2) > ##Create mesh in the Cartesian Coordinate system (X,Y) > nx = 100 > ny = nx > dx = 1.0 > dy = dx > CL = nx*dx > Cmesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny) > > I also defined two variables on these two meshes: > C_V = CellVariable(name = "Vacancy Concentration", mesh = Pmesh, value = > 0.0) > C_V_xy = CellVariable(name = "C_V in XY plane",mesh=Cmesh,value = 0.0) > > After solving the equation, I want to check the values at certain points, > If I use > C_V((1.0,1.0)) > The out put is array([ 0., 0.]). > But when I use > C_V_xy((20,20)) > I got an error: > IndexError: 0-d arrays can't be indexed > > I wonder why this happens and can I extract data at a certain point using > the coordinate values for the second one? > > best, > Zebo > > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > >
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
