Hi Yuan,
I don’t have any advice to give on shock-finding in yt or Enzo, but I thought I’d point out this paper:
http://arxiv.org/abs/1407.4117
You might find Section 5 interesting. Basically, a common assumption of shock finders is to filter using the dot product of the temperature and entropy gradients. When you do this, the computed dissipated energy can be sensitive to your assumptions for the other jump requirements, and depending on what you assume, can identify spurious shocks. This paper instead of using the entropy gradient, uses the density gradient, and gives somewhat different results. More knowledgable people than me can probably give you better advice, but I thought you might find this interesting, especially if you’re going to write something yourself.
Christine
On Jun 23, 2015, at 5:51 PM, Yuan Li

Dear all,

I am trying to figure out a way to mask out cells that contain shocks. One way I can think of is to create a derived field by identifying shocks. It seems that there are different ways of finding shocks: Enzo has a refinement criterion (refine by shocks) and there is also shock finding written by Sam.

Has anybody tried to find shocks using yt? I would appreciate it if you would like to share your code/experience/advice. Thank you!

Yuan _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org