// 不依赖于Builtin内核
SetFactory("OpenCASCADE");

lc = 0.00025;

// 在x-y平面画点
Point(1) = {-0.01, 0.02, 0, lc};
Point(2) = {0, 0.02, 0, lc};
Point(3) = {0, 0.025, 0, lc};
Point(4) = {0.08, 0.025, 0, lc};
Point(5) = {0.09, 0.035, 0, lc};
Point(6) = {0.1, 0.035, 0, lc};
Point(7) = {0.125, 0.01, 0, lc};
Point(8) = {0.127, 0.01, 0, lc};
Point(9) = {0.129, 0.0105, 0, lc};
Point(10) = {0.129, 0, 0, lc};
Point(11) = {-0.01, 0, 0, lc};
Point(12) = {0, 0, 0, lc};
Point(13) = {0.08, 0, 0, lc};
Point(14) = {0.1, 0, 0, lc};

// 马上旋转以上所有点
Rotate {{1,0,0},{0,0,0},2.5*Pi/180.0}
{
	Point{1:14};
}

// 在旋转过的点上定义线
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 5};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 9};
Line(9) = {9, 10};
Line(10) = {10, 14};
Line(11) = {14, 13};
Line(12) = {13, 12};
Line(13) = {12, 11};
Line(14) = {11, 1};

// tune for structured mesh begain
Line(15) = {2, 12};
Line(16) = {4, 13};
Line(17) = {6, 14};
Transfinite Curve {14, 15} = 80 Using Progression 1;
Transfinite Curve {16, 17, 9} = 120 Using Progression 1;
Transfinite Curve {2} = 41 Using Progression 1;
Transfinite Curve {1, 13} = 20 Using Progression 1;
Transfinite Curve {3, 12} = 200 Using Progression 1;
Transfinite Curve {11} = 99 Using Progression 1;
Transfinite Curve {4, 5} = 50 Using Progression 1;
Transfinite Curve {10} = 250 Using Progression 1;
Transfinite Curve {6} = 192 Using Progression 1;
Transfinite Curve {7, 8} = 30 Using Progression 1;
Curve Loop(101) = {14, 1, 15, 13};
Plane Surface(201) = {101};
Curve Loop(102) = {15, -12, -16, -3, -2};
Plane Surface(202) = {102};
Curve Loop(103) = {16, -11, -17, -5, -4};
Plane Surface(203) = {103};
Curve Loop(104) = {17, -10, -9, -8, -7, -6};
Plane Surface(204) = {104};
Transfinite Surface {201} = {11, 1, 2, 12};
Transfinite Surface {202} = {12, 3, 4, 13};
Transfinite Surface {203} = {13, 4, 6, 14};
Transfinite Surface {204} = {14, 6, 9, 10};
Recombine Surface {201, 202, 203, 204};
// tune for structured mesh end

// 旋转扫略
fluid[] = Extrude {{1,0,0},{0,0,0},-5*Pi/180.0}
{
	Surface{201}; Surface{202}; Surface{203}; Surface{204};
	Layers{1};
	Recombine;
};

Physical Surface("wedge0") = {201, 202, 203, 204};
Physical Surface("wedge1") = {208, 212, 216, 221};
Physical Surface("inlet") = {205};
Physical Surface("inlet-wall") = {206, 211};
Physical Surface("propellant") = {210};
Physical Surface("after-burner-wall") = {215, 214, 220};
Physical Surface("nozzle-wall") = {219, 218};
Physical Surface("outlet") = {217};

Physical Volume("fluidVolume") = {1, 2, 3, 4};