以下我们使用一个很简单的共振器范例模型进行功能示范 此范例将会示范如何使用GLAD。我们将使用resonator.inp这个档案示范共振器。即使使用者对共振器不感兴趣,这个范例也将会展示在GALD中解决问题的几个重要步骤: 初始化计算机数组及单位 选择波长 定义初始分布 使用宏进行重复运算 建立数据显示计算结果 此共振器将使用半对称的结构进行计算,由半径50cm的球面镜及平面镜所组成。共振器的长度为46cm。输出将由平面镜输出。下表1显示其结构参数: 图1-稳定的共振器结构。其光腰将会在平面镜上形成,及其相曲率将与曲面镜的理想模态曲率相同 为了简化讨论,我们将忽略增益及形成bare-cavity分析。我们开始分析从准备一个命令档案如下: variab/dec/int pass macro/def reson/o pass = pass + 1 # increment pass counter prop 45 # propagate 45 cm. mirror/sph 1 -50 # mirror of 50 cm. radius clap/c/n 1 .14 # .14 cm. radius aperture prop 45 # propagate 45 cm. along beam mirror/flat 1 # flat mirror variab/set Energy 1 energy # set variable to energy value Energy = Energy - 1 # calculate energy difference udata/set pass pass Energy # store energy differences energy/norm 1 1 # renormalize energy plot/l 1 xrad=.15 # make a plot at each pass macro/end array/set 1 64 # set array size wavelength/set 0 1.064 # set wavelengths units/set 1 .005 # set .005 cm sample spacing resonator/name reson # set name of resonator macro resonator/eigen/test 1 # find resonator properties resonator/eigen/set 1 # set surrogate beam to eigen mode clear 1 0 # clear the array noise 1 1 # start from noise energy/norm 1 1 # normalize energy pass = 0 # initialize pass counter reson/run 100 # run resonator 100 times title Energy loss per pass plot/watch plot1.plt # set plot name plot/udata min=-.05 max=.0 # plot summary of eigenvalues title diffraction mode shape set/density 32 # set plot grid to 32 x 32 set/window/abs -.05 .05 -.05 .05 # set plot window plot/watch plot2.plt # set plot name plot/iso 1 # make an isometric plot 以下就对每一项指令来做介绍: variab/dec/int pass 此行定义一整数变量叫pass。我们将使用pass来储存数据,变量如果不清楚的定义为整数,将会被定义为实数变量。 macro/def reson/o 此行开始定义宏,就像是子程序或函式一样。所有介于macro/def与macro/end之间的指令都将定义为宏。这些指令暂时不会被执行。这些指令列将被放在MACLIB中留待以后使用。这些宏指令列将不需缩排。但使用缩排将会使这些指令更容易阅读。 |