#!/usr/bin/octave -qf
output_precision=15;

UTFIL= ['flange_t28.0_EN13715.wheel'];
flange_thickness= 28.0;

de= 32.5 - flange_thickness;


D1y= 35+de;     D1z= 6.867;
Dmy=26.211+de;  Dmz=16.446;     RD= 13;
E1y= 38.427+de; E1z= 12.;

%       Flange Height 32
F1y= 41.497+de; F1z= 20.434;
Rfa= 23;
G1y= 46.153+de;
%       Flange Height 30
F1y= 40.530+de; F1z= 17.779;
Rfa= 21.5;
G1y= 47.758+de;
%
%       Flange Height 28
F1y= 39.765+de; F1z= 15.675;
Fmy= 58.558+de; Fmz=  8.835;    Rfa= 20;
G1y= 49.663+de; G1z= 26.748;
REmy= 55.000+de; REmz= 16.000;    RE= 12;
a = 15.-de*.5;
Sy= 70-a;       Sz= 28.;
RImy= 55.000;   RImz= 16.000;    RI= 12;
H1y= 62.765;    H1z= 25.149;
Hmy= 49.500;    Hmz=  9.519;    RH= 20.5;
H2y= 70.000;    H2z=  9.519;
% ===================================================== End Of Input Data

disp("#!/bin/sh")
disp("")
%disp("cat <<+")
disp("#")
disp("# Radius 13 from D1 to E1")
disp("# ----------------------------------------")
 disp("")
 disp("$gensys/exe/misc/func <<+")
 disp(" FUNC = circle")
 disp(" DX   = .02")
 disp("#           x0,      y0,      R,       x1,      x2,     side")
printf(" FIN  =   %14.10f, %14.10f, %14.10f, %14.10f, %14.10f, -1.,\n",
                  Dmy, Dmz,  RD,  D1y,E1y);
 disp("")
 disp(" DATA = SNGL")
 disp(' UTFIL= func_out.sngl  FORMUT = SNGL  TYPE_UTFIL = CREATE')
 disp(" eof")
 disp("+")
 disp("# ----------------------------------------------------------------------------------------------------")
 disp('if [ $? -ne 0 ]; then exit $?; fi')

 disp("")
 disp("#                               Create a strait line from E1 to F1")
 disp("$gensys/exe/misc/func <<+")
 disp(" FUNC = pline_1")
 disp(" DX   = .02")
 disp("#          x1,              y1,                x2,              y2")
printf(" FIN  =   %14.10f,  %14.10f,    %14.10f,  %14.10f,\n",
                  E1y, E1z,   F1y,  F1z);
 disp(" DATA = SNGL")
 disp(' UTFIL= func_out.sngl  FORMUT = SNGL  TYPE_UTFIL = APPEND')
 disp(" eof")
 disp("+")
 disp("# ----------------------------------------------------------------------------------------------------")
 disp('if [ $? -ne 0 ]; then exit $?; fi')

disp("#")
disp("# Radius Rfa from F1 to G1")
disp("# ----------------------------------------")
 disp("")
 disp("$gensys/exe/misc/func <<+")
 disp(" FUNC = circle")
 disp(" DX   = .02")
 disp("#           x0,      y0,      R,       x1,      x2,     side")
printf(" FIN  =   %14.10f, %14.10f, %14.10f, %14.10f, %14.10f,  1.,\n",
                  Fmy, Fmz,  Rfa,  F1y,G1y);
 disp("")
 disp(" DATA = SNGL")
 disp(' UTFIL= func_out.sngl  FORMUT = SNGL  TYPE_UTFIL = APPEND')
 disp(" eof")
 disp("+")
 disp("# ----------------------------------------------------------------------------------------------------")
 disp('if [ $? -ne 0 ]; then exit $?; fi')

disp("#")
disp("# Radius RE  from G1 to S")
disp("# ----------------------------------------")
 disp("")
 disp("$gensys/exe/misc/func <<+")
 disp(" FUNC = circle")
 disp(" DX   = .02")
 disp("#           x0,      y0,      R,       x1,      x2,     side")
printf(" FIN  =   %14.10f, %14.10f, %14.10f, %14.10f, %14.10f,  1.,\n",
                 REmy,REmz,  RE,   G1y, Sy);
 disp("")
 disp(" DATA = SNGL")
 disp(' UTFIL= func_out.sngl  FORMUT = SNGL  TYPE_UTFIL = APPEND')
 disp(" eof")
 disp("+")
 disp("# ----------------------------------------------------------------------------------------------------")
 disp('if [ $? -ne 0 ]; then exit $?; fi')

disp("#")
disp("# Radius RI  from S  to H1")
disp("# ----------------------------------------")
 disp("")
 disp("$gensys/exe/misc/func <<+")
 disp(" FUNC = circle")
 disp(" DX   = .02")
 disp("#           x0,      y0,      R,       x1,      x2,     side")
printf(" FIN  =   %14.10f, %14.10f, %14.10f, %14.10f, %14.10f,  1.,\n",
                 RImy,RImz,  RI,    Sy,H1y);
 disp("")
 disp(" DATA = SNGL")
 disp(' UTFIL= func_out.sngl  FORMUT = SNGL  TYPE_UTFIL = APPEND')
 disp(" eof")
 disp("+")
 disp("# ----------------------------------------------------------------------------------------------------")
 disp('if [ $? -ne 0 ]; then exit $?; fi')

disp("#")
disp("# Radius RH  from H1 to H2")
disp("# ----------------------------------------")
 disp("")
 disp("$gensys/exe/misc/func <<+")
 disp(" FUNC = circle")
 disp(" DX   = .02")
 disp("#           x0,      y0,      R,       x1,      x2,     side")
printf(" FIN  =   %14.10f, %14.10f, %14.10f, %14.10f, %14.10f,  1.,\n",
                  Hmy, Hmz,  RH,   H1y,H2y);
 disp("")
 disp(" DATA = SNGL")
 disp(' UTFIL= func_out.sngl  FORMUT = SNGL  TYPE_UTFIL = APPEND')
 disp(" eof")
 disp("+")
 disp("# ----------------------------------------------------------------------------------------------------")
 disp('if [ $? -ne 0 ]; then exit $?; fi')

 disp(['mv func_out.sngl ',UTFIL])
 disp(['echo File ',UTFIL,' has been created.'])