refactored to change implicit real*8 to implicit none
modified advection_3d_inflow to smaller diameter inflow circle and cfl = 0.7 for stability