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