[求助] zemax鳞甲反光杯建模

突然要做一个鳞甲反光杯，从0开始学zemax非序列模式，跟着南京光研的教程，但是我用他的程序生成出来是这样的：


程序是这样的：
N=40
F=50
TH=1
R2=50


PI=ACOS(-1)
DFI=(2*PI/N)*(2*PI/N)
Z1=F/2
DZ1=0


FOR j,1,3,1
Y1=SQRT(4*F*Z1)
A1=1+(4*F*F/(Y1*Y1))*(1-DFI)
B1=-4*F*DFI
C1=-(Y1*Y1*DFI)
DZ1=(-B1+SQRT(B1*B1-4*A1*C1))/(2*A1)
R1=2*PI*Y1/N
ARFA=(ATAN(2*F/Y1))*57.29578
gosub drawlens1
Z1=Z1+DZ1
next

end

sub drawlens1
k=0
m=0
for q,1,N,1
INSERTOBJECT 1,2
SETNSCPROPERTY 1,2,0,0,"NSC_SLEN"
SETNSCPROPERTY 1,2,4,0,"PMMA"
SETNSCPROPERTY 1,2,2,0,-1
SETNSCPOSITION 1,2,3,Z1
SETNSCPOSITION 1,2,2,Y1*cosi(m)
SETNSCPOSITION 1,2,1,Y1*sine(m)

SETNSCPOSITION 1,2,4,90
SETNSCPOSITION 1,2,5,k
SETNSCPARAMETER 1,2,1,0
SETNSCPARAMETER 1,2,2,0
SETNSCPARAMETER 1,2,3,R1
SETNSCPARAMETER 1,2,4,R1
SETNSCPARAMETER 1,2,5,TH
SETNSCPARAMETER 1,2,6,R2
SETNSCPARAMETER 1,2,7,0
SETNSCPARAMETER 1,2,8,R1
SETNSCPARAMETER 1,2,9,R1
SETNSCPROPERTY 1,2,5,2,"HEAR1"

SETNSCPOSITION 1,2,4,-ARFA
k=k-(360/N)
m=m+((2*360*3.14159)/(360*N))
next
Return



如果不要SETNSCPOSITION 1,2,4,-ARFA这句话生成出来的就是这样的：


怎么样才能搞成他教程这样的：


各位大佬有什么建议吗？

各个公司有不同的做法，你可以先做TIR形式的光面，然后用3D软件往上贴鳞片；或者自建鳞片数据，在非序列里调一列鳞片然后圆形阵列，再输出到3D软件里面裁剪组合。

那可以帮忙看看我这个程序哪里有问题嘛

首先我们公司不用这种方法做；

其次个人不太看好这种方法，这种批量阵列出来的，感觉对光线不够细分，缺少控制；
做出来肯定有问题

好吧 才搞这个两周感觉还有很多不懂的