If I want to simulate the situation of the binary neutron star with the gas surrounding, and the gas can be in different density and different velocity distribution, how could I do ?(Should I insert the thorns "GRHydro_InitData"and "EOS_Polytrope" into the parameter file of binary neutron stars--"NsNsToHMNS.par" ?)

Best,

Shih-Jie