procedure tapcheck (imageName, coordFile, apDiam, innerSky, outerSky, frm) string imageName {prompt="Full image name of strip mosaic"} string coordFile {prompt="Name of coordinate file"} int apDiam {prompt="Aperture diameter"} int innerSky {prompt="Inner radius of sky annulus"} int outerSky {prompt="Outer radius of sky annulus"} int frm {1, prompt="Ximtool frame to display in"} begin # Variable declarations: int apRad # aperture radius # Make sure the IMAGES and IMGTOOLS packages are loaded: if (! defpac ("images")) { print("IMAGES package not loaded") bye } print (" ") print ("Displaying apertures in frame "//frm//" ...") print (" ") display (image=imageName,fra=frm,bpdispl="none",erase+,bord-, sele+,repe-,fill-,zscal+,ztrans="linear") apRad = int(apDiam/2) tvmark(frame=frm,coords=coordFile,logfile="",autolog-,outimag="",mark="circle", radii=apRad,lengths=0,font="raster",color=0,label-,number-,pointsi=3,txsize=2, interac-,flpar-,mode="ql") tvmark(frame=frm,coords=coordFile,logfile="",autolog-,outimag="",mark="circle", radii=innerSky,lengths=0,font="raster",color=0,label-,number-,pointsi=3,txsize=2, interac-,flpar-,mode="ql") tvmark(frame=frm,coords=coordFile,logfile="",autolog-,outimag="",mark="circle", radii=outerSky,lengths=0,font="raster",color=0,label-,number-,pointsi=3,txsize=2, interac-,flpar-,mode="ql") print (" ") print ("Done.") print (" ") end