On Tue, Aug 18, 2015 at 3:18 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Hi,
I have a user defined field with the following definition.
def _epsilon(field, data): return (4.0/3.0)*(data['vShock']/data["VelocityMagnitude"]) add_field("epsilon", function=_epsilon, take_log=False, units=r" " )
Now I would like take ln(1 + data["epsilon"]). I am trying with the following definition. Is this the correct way of doing this ?
def _lnepsilon(field, data): return 1 + data["epsilon"] add_field("lnepsilon", function=_lnepsilon, *take_log=True*, units=r" " )
The take_log keyword argument only affects how the plot is displayed by default in e.g. . If you want the field to return the logarithm of 1+data['epsilon'], then you should do something like: def _lnepsilon(field, data): return np.log(1 + data['epsilon']) add_field("lnepsilon", function=_lnepsilon, units="")
Thank you -- Reju Sam John
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org