!    filename: spot_sz.zpl	date: 12/9/96	rev: 1 
!  Computes Gaussian Beam Diameter after each element.

n = NSUR()
INPUT "Enter surface# for start of first lens element: ",surfstart
pie=3.14159

FOR i=1,n,1
  GETSYSTEMDATA 1
  f[1]=196.72
  w[1]=VEC1(11)
  d[i]=SDIA(surfstart)
  wave=WAVL(PWAV())
    
  PRINT wave
  PRINT VEC1(11)

  z[i]=(pie*POWR(w[1],2))/wave
  w[i]=powr(powr(f[1]*w[1],2)/((powr(d[i]-f[1],2))+powr(z[i],2)),.5)
  spotsize=2*w[i]  

  PRINT z[i]
  PRINT w[i]
  PRINT spotsize
  PRINT "The Gaussian beam diameter after lens ",i, " is: ",spotsize, " microns"
  surfstart=surfstart+2
NEXT