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 domain.set_quantity('stage', Set_Stage(0.2, 0.4, 1.0)) # Let processor 0 output some timing information visualise = False if visualise:
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 domain.set_quantity('stage', Set_Stage(0.2,0.4,1.0)) # Let processor 0 output some timing information visualise = False