Guest User

Untitled

a guest
Jul 9th, 2013
169
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.24 KB | None | 0 0
  1. def _NumberDensity(fields,data) :
  2.     try:
  3.         return data["nele"]+data["nion"]
  4.     except:
  5.         pass
  6.     return data['pres']/(data['temp']*Na*kboltz*mh)
  7. add_field("NumberDensity", function=_NumberDensity,
  8.         units=r'\rm{cm}^{-3}')
Advertisement
Add Comment
Please, Sign In to add comment