# wings0 probleme sans portance P0=0, sauf final P0 -> bonne valeur
# fev 2011
# gerris2D -DP0=0 -Dalph=$alpha wing0.gfs | gfsview2D w.gfv -s
Define L0 100.
Define pis180 0.0174532925199433

1 0 GfsPoisson GfsBox GfsGEdge {} {
    Time { iend = 2 }
    Refine 8
    RefineSolid 15
    PhysicalParams { L = L0 }
    Solid aile.gts
    SurfaceBc P Dirichlet P0
     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 } ailepress0.dat ailepos.dat
    OutputLocation { istep = 1 } vbdfp0.dat bdfp.dat
    OutputLocation { istep = 1 } vbdfm0.dat bdfm.dat  
}
GfsBox {
    top = Boundary {
	BcDirichlet P L0/2.*cos(alph*pis180)-x*sin(alph*pis180) 
    }
    left = Boundary {
	BcDirichlet P y*cos(alph*pis180)+L0/2*sin(alph*pis180) 
    }
    right = Boundary {
	BcDirichlet P y*cos(alph*pis180)-L0/2*sin(alph*pis180) 
    }
    bottom = Boundary {
	BcDirichlet P -L0/2.*cos(alph*pis180)-x*sin(alph*pis180) 
    }
}