The functional form of the profile is chosen by setting
where
and
and are input quantities.
The number of rational surfaces is set equal to and the corresponding -values and layer widths are defined through input arrays and . The degree of flattening is determined by , where has no effect on and will make at the rational surface.
The coefficients are obtained from a DATA file:
call link('unit11=(qpoly,open)//') read(11,2000)mpolyq 2000 format(i2) read(11,2100)(cq(i),i=1,mpolyq) 2100 format(1e12.5)
call link('unit21=(qdata,open)//') read(21,2000)nsurf 2000 format(i2) read(21,2100)(psidk(i),qdk(i),i=1,nsurf) 2100 format(2e13.5)