Boundaries

region 1 { define the domain boundary }

start(-far, -far) { start at the bottom left }

value(psi) = psi_far line to (far, -far) { walk the boundary Counter-Clockwise }

line to (far, far) line to (-far, far) line to close { return to close }

plots { write hardcopy files at termination }

contour(psi) zoom (-4, -4, 8, 8) as " line current numerically " { show the flow vectors: } contour(psi1) zoom (-4, -4, 8, 8) as " line current analytic "

vector(-dy(psi), dx(psi)) zoom (-4, -4, 8, 8)as " speed numerically " norm

vector(-dy(psi1), dx(psi1)) zoom (-4, -4, 8, 8)as " analytic speed " norm

contour(Pt) zoom (-4, -4, 8, 8) as " field of the pressures, analytic” painted

contour(P) zoom (-4, -4, 8, 8) as " field of the pressures, numerical " painted

end