1 0 GfsAxi GfsBox GfsGEdge {} {
    Time { end = 0 }
    PhysicalParams { L = 50 }
    AdvectionParams { scheme = none }
    ApproxProjectionParams { tolerance = 1e-10 }
    Refine 4
    Refine (LEVEL + 1./50.*(x*x + y*y)*(4. - LEVEL))
    Global {
	#define A0 0.5
	#define U0 1.
    }
    Solid (ellipse (0., 0., A0, A0))
    Init {} {
	U = U0
	Phi = {
	    double r = sqrt (cx*cx + cy*cy);
	    return U0*A0*A0*A0*cx/(2.*r*r*r);
	}
    }
    OutputErrorNorm { start = end } { awk '{ print LEVEL " " $7 " " $9}' } { v = U } {
 	s = (dx("Phi") + 1.)
    }
    OutputSimulation { start = end } sim-LEVEL.gfs
}
GfsBox {
    left = Boundary { BcDirichlet U U0 }
    right = Boundary { BcDirichlet U U0 }
}