Rumen,
As Steve said, fix gravity won’t work for SRD, there are some minor changes to the code that we had to make for forced flow.
Attached is the version of fix srd that I used for those papers. Unfortunately there are MANY other changes/additions, some of which aren’t used or working, so it’s a bit messy. If you search for “FPOISS” in the code you’ll find the relevant parts. Basically, the SRD advection step is modified slightly to simulate the effect of a force acting on SRD particles (I believe only in the x-direction, but it’d be easy to modify for others).
I’m leaving for vacation in a few hours, so I’ll be out of touch, but if you have questions feel free to send them my way and I can try to help when I come back (6/21).
-Dan Bolintineanu
fix_srd_mod.tar (250 KB)