# Boundaries T = Transmissive_boundary(domain) domain.set_boundary( {'left': T, 'right': T, 'bottom': T, 'top': T, \ 'ghost': None} ) # Set Evolve parameters domain.set_default_order(2) domain.set_timestepping_method('rk2') print domain.get_timestepping_method() #domain.use_edge_limiter = True #domain.tight_slope_limiters = True #domain.use_centroid_velocities = False domain.CFL = 1.0 domain.set_beta(0.8) # Ensure that the domain definitions make sense domain.check_integrity() # Set the inititial conditions