Paste #178
Welcome To LodgeIt
Welcome to the LodgeIt pastebin. In order to use the notification feature a 31 day cookie with an unique ID was created for you. The lodgeit database does not store any information about you, it's just used for an advanced pastebin experience :-). Read more on the about lodgeit page. Have fun :-)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | import yt import numpy as np loc='/share/Part1/guido/harris_iso/test/md_hdf5_plt_cnt_0000' ds =yt.load(loc) ray =ds.ortho_ray(0,[0,0]) by =ray[('flash','magy')] P =ray[('flash','pres')] #gas pressure P_B =by**2/(8*np.pi) # magnetic pressure (cgs) P_t =P+P_B # total pressure print("P: {}".format(P)) print("P_B: {}".format(P_B)) print("P_t: {}".format(P_t)) |