def new_instance_of_athena_selfgravity_code(self):
result=Athena(mode = AthenaInterface.MODE_SELF_GRAVITY, redirection="none", number_of_workers=1)
result.initialize_code()
result.parameters.gamma = self.gamma
result.parameters.courant_number=0.3
result.set_grav_mean_rho(self.rho_mean.value_in(density))
return result
def new_instance_of_athena_selfgravity_code(self):
result = Athena(mode=AthenaInterface.MODE_SELF_GRAVITY,
redirection="none",
number_of_workers=1)
result.initialize_code()
result.parameters.gamma = self.gamma
result.parameters.courant_number = 0.3
result.set_grav_mean_rho(self.rho_mean)
return result
