%begin
%include "edwin:specs.inc"
%include "edwin:consts.inc"
%include "utils:region"
%integer ref,size,i,j,dev,numlines,a,b,c,d
%bytearray epic(0:128*512)
%string(255) cubefile,linefile,plotfile
%constreal     xmult=1.264705882

printline("Hello there!")
prompt("Input Display File:");read(cubefile)
prompt("Input Line File:");read(linefile)
prompt("Input plot file:");read(plotfile)

accessfile(cubefile,0,ref,size)
readregion(ref,0,size,epic(0))
deaccessfile(ref)

openoutput(1,cubefile.".pdf")

dev=default device
initialise for(dev)
newframe
storeon(1)

setcolour(green)
moveabs(0,0)
lineabs(512*2,0)
lineabs(512*2,128*2)
lineabs(0,128*2)
lineabs(0,0)

moveabs(0,270)
lineabs(1056,270)
lineabs(1056,990)
lineabs(0,990)
lineabs(0,270)

setcolour(red)
%for i=0,1,128 %cycle
  %for j=0,1,512 %cycle
    %if epic(j+i*512)=1 %then markerabs(0,2*j,2*i)
  %repeat
%repeat

openinput(1,linefile)
selectinput(1)
setcolour(black)
read(numlines)
%for i=1,1,numlines %cycle
  read(a);read(b);read(c);read(d)
  a=int(a/xmult)
  c=int(c/xmult)
  moveabs(2*a,2*b)
  lineabs(2*c,2*d)
%repeat
closeinput
selectinput(0)

setcolour(red)
openinput(2,plotfile)
selectinput(2)
read(numlines)
%for i=1,1,numlines %cycle
  read(a);read(b)
  markerabs(0,2*a,270+2*b)
%repeat
read(numlines)
setcolour(black)
%for i=1,1,numlines %cycle
  read(a);read(b);read(c)
  moveabs(2*b,270+2*c)
  linerel(0,2*a)
%repeat
closeinput
selectinput(0)
terminate edwin
%endofprogram
