On Sun, Jun 5, 2011 at 12:33 PM, Burak Atakan <herz...@sonnenkinder.org> wrote: > > Dear all, > > we need the cellVolumes of a 2D cylindrical mesh and tried the following: > > from fipy import * > mesh=CylindricalGrid2D(dr=1.,dz=1.,nr=10,nz=1) > a=mesh.getCellVolumes() > print a > > The output is: > > [ 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5] > > Could it be that a factor of pi is missing here? Assuming that the half > cylinder > is modelled...
Hi Burak, The assumption is a 1 Rad wedge. For the inner working of FiPy the factor is not important. Hope this helps. -- Daniel Wheeler
