Hello zhaozc,
I failed to simulate this for several months, so I'm finding (however not found) is there some code public.
Please have a look at the parfiles and code that goes along with:
https://arxiv.org/abs/1509.08804
as well as the followup papers by the same group.
The code and parfiles are available here:
http://www.fis.unipr.it/gravity/Research/BNS2015.html
Yours, Roland