Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def _NumberDensity(fields,data) :
- try:
- return data["nele"]+data["nion"]
- except:
- pass
- return data['pres']/(data['temp']*Na*kboltz*mh)
- add_field("NumberDensity", function=_NumberDensity,
- units=r'\rm{cm}^{-3}')
Advertisement
Add Comment
Please, Sign In to add comment