Define THETA atan2(y,x) 
1 0 GfsSimulation GfsBox GfsGEdge { x = 0.5 y = -0.5 } {
    Time { dtmax = 1 }
    Refine LEVEL
    AdvectionParams { scheme = none }
    Solid (-difference(ellipse(0,0,0.7,0.7), ellipse(0,0,0.3,0.3)))
    SourceCoriolis 0 1
    EventStop { istep = 1 } U 1e-9
    Variable k
    Init {} { k = THETA>-M_PI/3 ? 1 : 0.5 }
    PhysicalParams { alpha = k }

    EventList { start = end } {
	OutputErrorNorm {} { awk '{print LEVEL, $5, $7, $9}' >> p } { v = P } {
	    s = THETA>-M_PI/3 ? THETA : 2*THETA+M_PI/3
	}
	OutputErrorNorm {} { awk '{print LEVEL, $5, $7, $9}' >> U } { v = U } {
	    s = y / (x*x + y*y)
	}
	OutputErrorNorm {} { awk '{print LEVEL, $5, $7, $9}' >> V } { v = V } {
	    s = -x / (x*x + y*y)
	}
	OutputSimulation {} solution.gfs	
    }
}
Box { 
    top = Boundary { 
	BcDirichlet P 0
	BcDirichlet U 0
	BcNeumann V 0
    } 
    left = Boundary {
	BcDirichlet P -2*M_PI/3
	BcNeumann U 0
	BcDirichlet V 0
    }
}