Nb 36 "Number of blades"
psi 2*pi[rad]/Nb "Sector angle"
rtip 100[mm] "Tip radius"
rroot 0.8*rtip "Root radius"
rrms sqrt((rroot^2+rtip^2)/2) "Root mean square radius"
S 2*rrms*sin(psi/2) "Blade spacing at rms radius"
B S "Blade length, S0=S/B=1"
alpha 30[deg] "Blade angle"
rh B*sin(alpha) "Row height"
wp11x rroot "Work plane 1, point 1, x-coordinate"
wp11y 0 "Work plane 1, point 1, y-coordinate"
wp11z 0 "Work plane 1, point 1, z-coordinate"
wp12x rtip "Work plane 1, point 2, x-coordinate"
wp12y 0 "Work plane 1, point 2, y-coordinate"
wp12z 0 "Work plane 1, point 2, z-coordinate"
wp13x rroot "Work plane 1, point 3, x-coordinate"
wp13y 1[cm] "Work plane 1, point 3, y-coordinate"
wp13z 1[cm]*tan(alpha) "Work plane 1, point 3, z-coordinate"
wp21x wp11x*cos(psi)-wp11y*sin(psi) "Work plane 2, point 1, x-coordinate"
wp21y wp11x*sin(psi)+wp11y*cos(psi) "Work plane 2, point 1, y-coordinate"
wp21z wp11z "Work plane 2, point 1, z-coordinate"
wp22x wp12x*cos(psi)-wp12y*sin(psi) "Work plane 2, point 2, x-coordinate"
wp22y wp12x*sin(psi)+wp12y*cos(psi) "Work plane 2, point 2, y-coordinate"
wp22z wp12z "Work plane 2, point 2, z-coordinate"
wp23x wp13x*cos(psi)-wp13y*sin(psi) "Work plane 2, point 3, x-coordinate"
wp23y wp13x*sin(psi)+wp13y*cos(psi) "Work plane 2, point 3, y-coordinate"
wp23z wp13z "Work plane 2, point 3, z-coordinate"