! 
! draw a map of amplitudes or phases using the eoutput file from SOLVE HOLO 
! 
! @ holo file amp|pha number 
! 
! 
if "&1*".eq."*" then 
  say "Usage : @ PR:HOLO file AMP|PHA [NUMBER]" 
  return 
endif 
if (exist(a)) then 
  del /var a 
endif 
if (exist(rms)) then 
  del /var rms rmsw 
endif 
def real rms rmsw /global 
gtvl\clear whole 
sic\define image a &1 read /global 
! 
! plot first map 
greg1\set expan 0.9 
greg1\set box_location 6 26 1 21 
if "&2*".eq."AMP*" then 
  levels -30 to 0 by 3 
  greg2\rgdata a[3] /variable 
  greg1\limits -7.5 7.5 -7.5 7.5 
  greg2\plot a[3] /scaling lin -30 0 
else 
  level -pi to pi by pi|18 
  greg2\rgdata a[2] /variable 
  greg1\limits -7.5 7.5 -7.5 7.5 
  greg2\plot a[2] /scaling lin -pi|2 pi|2 
  @ holo_panels.graphic &3 
endif 
rgmap 
greg1\box n n 
@ holo_rms.graphic a 
draw text 0 -0.1 "&1 "'0.001*nint(1000*rms)'" (unw) "'0.001*nint(1000*rmsw)'" (w) rad." 2 /box 2 
!