# wings1 probleme avec uniquement la circulation
# fev 2011
# gerris2D -m wing1.gfs | gfsview2D w.gfv  -s
Define L0 100.


1 0 GfsPoisson GfsBox GfsGEdge {} {
    Time { iend = 2 }
    Refine 8
    RefineSolid 15
    PhysicalParams { L = L0 }
    Solid aile.gts
    SurfaceBc P Dirichlet 1
     Init { istep = 1 } { U =  dy("P")
                          V = -dx("P")
                          press = 0.5-0.5*(U*U+V*V) }
    AdaptGradient { istep = 1 } { cmax = 1e-2 maxlevel = 13 } P
    ApproxProjectionParams { nitermin = 1 }
    OutputProjectionStats { istep = 1 } stderr
    OutputSimulation { istep = 1 } stdout
    OutputLocation { istep = 1 } ailepress1.dat ailepos.dat
    OutputLocation { istep = 1 } vbdfp1.dat bdfp.dat
    OutputLocation { istep = 1 } vbdfm1.dat bdfm.dat  
}
GfsBox {
    top = Boundary {
	BcDirichlet P 0
    }
    left = Boundary {
	BcDirichlet P 0
    }
    right = Boundary {
	BcDirichlet P 0 
    }
    bottom = Boundary {
	BcDirichlet P 0
    }
}