Dear developers,
I am new to fipy and at the moment I am trying to implement it for a 2D
convection-diffusion problem. The package is great and rather easy to use.
Nonetheless, I have an issue that might sound trivial to you. The issue
regards the plot of the 2D mesh. When I implement a value for the length of
the x-axis that is much smaller than the value for the length of the
y-axis, the viewer does not scale the axes to obtain a more or less
"square" plot, but displays a very squashed plot instead. I was wondering
if you could give me a hand in solving this problem. I added a simplified
version of the code below.
Thanks in advance.
Kind regards,
Davide Cretti
N=100
Lx=0.1
Ly=1
dt=1
steps=1
alpha=1
mesh = Grid2D(dx=Lx/N, dy=Ly/N, nx=N, ny=N)
phi = CellVariable(mesh=mesh, value=1)
eq = TransientTerm() == DiffusionTerm(coeff=alpha)
viewer = Viewer(phi, mesh)
for steps in range(steps):
eq.solve(var = phi, dt=dt)
viewer.plot()
if __name__ == '__main__':
input()
_______________________________________________
fipy mailing list
fipy@nist.gov
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
