-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On 28/10/15 14:21, Buesing, Henrik wrote:
[Buesing, Henrik] Not sure that we are on the same page here. I attached the code of calc_temp(p,h). DATA_temp is an array with 146,888 (602*244) entries.
Of course, when doing
gradtemp = (temp('+') - temp('-'))/Delta_h_facet
Ah, I think I see. What space does this resulting temp field live in, is it just DG0. In which case, you could do: temp_coeff = Function(DG0) temp_coeff.dat.data[:] = calc_temp(p, h) then use temp_coeff later in gradtemp. Does that make sense? Lawrence -----BEGIN PGP SIGNATURE----- Version: GnuPG v2.0.22 (GNU/Linux) iQEcBAEBAgAGBQJWMNtFAAoJECOc1kQ8PEYvFJAH/RYtWVeEj3/awCoVDRwWJoh0 QE1VLdqo9vlYcEkvvDGa018BmF+EUGs8sIRqPBtGDKEWeVP5u9L+h1Zw/dfvspdv dd6EkJ2KI/xm2gIioO4422qBQTgTVonKqke3s1l7fRaOR67FeWQJPxlbllqyQJFP Zejqb7Y8HbOyWzCbhKJLJrSJVFyV308kyQJ783jB5Q7nxJ56c6y0Jxr9196jdbao mm1YkXoqvCq4azkdKr00zOGvA2nj7RAcbixltmfx/Kk/kfiDE2HByPIrzVv2czlU azKfTq4qsHp4mlqxnTMapO0v89+L8BxGBaDnE9KWP7y69xoGyxqqDgBFz130TVY= =RcxF -----END PGP SIGNATURE-----