zemax中宏编程调用程序求助

cathyzhyq

大家好~~现在要求通过ZPL产生一个随机位置的光源（在一个固定圆柱体内的随机光源），我用macro编了一个程序，要求重复调用几万到几十万次产生随机位置，但在程序测试时发现仅仅第一条和第二条是不同的，后面的调用都重复了第二次的结果，求教各位大神如何解决~


下附程序问题集中点：
!This is the main program


For i,1,100,1
GOSUB source_position

PRINT
PRINT i
PRINT "Reading on detector is", x,  y,  z
NEXT

END   # This is the end of the main program


! This is the sub-routine
SUB source_position
r=RAND(9.9)
angle=RAND(360)
z0=RAND(9.9)
pi=3.1415926


y=r*SINE(angle*pi/180)
x=y/TANG(angle*pi/180)
z=z0+0.1
RETURN   # This is the end of the sub-routine
谢谢大家^_^

yinlovetian

FOR i=1,100,1看看这样行不，我在我电脑上100次都不一样

cathyzhyq

yinlovetian 发表于 2017-5-7 19:31
FOR i=1,100,1看看这样行不，我在我电脑上100次都不一样

多谢啦，我刚刚试了一下，这个是可以了，但是这个是一个测试程序，主程序里面其实还有一个调用过程，用x，y，z这三个随机数作为source的坐标，用到了
SETNSCPOSITION 1,4,1,x
SETNSCPOSITION 1,4,2,y
SETNSCPOSITION 1,4,3,z
ray trace之后再重新调用下一组x，y，z
但是执行到第二次就失败了。
求问如何将光源的坐标归零，我猜想是不是这个原因
但是写完
SETNSCPOSITION 1,4,1,0
SETNSCPOSITION 1,4,2,0
SETNSCPOSITION 1,4,3,0     
之后依然不行。  
求问如何多次设置光源坐标啊，每次都采用一组随机数，只改变x positon，y position 和z position，其他的parameter直接在NSC editor中设置，每次循环均不改变。
多谢多谢~~

cathyzhyq

下面是ZPL macro的全部程序段，麻烦大家帮我挑一挑错。多谢啦~~NSC editor中只有一个光源和一个探测器
分别是 object 4 ——source point；
和 object 5 ——detector rectangle；
ZPL旨在产生几万~几十万个随机位置作为光源坐标，带入NSC系统，进行ray trace， 最终在探测器 detector rectangle上得到能量值并且一条一条输出
坐标-能量 对应数据：x, y, z, energy。

!This is the main program

For i=0,10,1
GOSUB source_position
NEXT
END   # This is the end of the main program

! This is the sub-routine
SUB source_position

!create 3 randoms
r=RAND(9.9)
angle=RAND(360)
z0=RAND(9.9)
pi=3.1415926

!create the coordination parameter
y=r*SINE(angle*pi/180)
x=y/TANG(angle*pi/180)
z=z0+0.1

!set the source_position
SETNSCPOSITION 1,4,1,x
SETNSCPOSITION 1,4,2,y
SETNSCPOSITION 1,4,3,z

!ray tracing
surf=1
source=4
split=1
scatt=1
pol=1
ignore_err=1
random_seed=1
save=0
object=0
obj=5
pix=0
data=0

temp=NSDD(surf, object, pix, data)   #clear the detector
NSTR surf, source, split, scatt, pol, ignore_err, random_seed, save

!print
PRINT
PRINT i
PRINT "Reading on detector is", x,  y,  z,  NSDD(surf, obj, pix, data)

!recover source_position
SETNSCPOSITION 1,4,1,0
SETNSCPOSITION 1,4,2,0
SETNSCPOSITION 1,4,3,0

RETURN  # This is the end of the sub-routine