I think the more direct approach is to edit the numpy array qb.dat.data. E.g. something like qb.dat.data[np.isnan(qb.dat.data)] = 0 On 17 June 2016 at 03:52, Francis Poulin <fpoulin@uwaterloo.ca> wrote:

Hello,

I want to define an initial condition for a field that I am pretty sure will have nan's. It does in python so I'm pretty sure that it will in firedrake as well.

In python I know how to remove that nan's by doing the following.

qb[np.isnan(qb)] = 0.

Is there a way of doing this in firedrake?

Cheers, Francis

------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo

email: fpoulin@uwaterloo.ca Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637