#
#  This makes some nice graphics for a display
#  ON BEAM - OFF BEAM
#
#
#

include 'image.g'

makeexample := function(rad=150, dx=30, dy=90) {

dc.copy('PLAY.0', 'IM.dummy1');
im := image('IM.dummy1');

mysh :=  im.shape();
#  should be [2048, 2048, 1, 1, 1];
cent1 := [1500, 1000, 1, 1, 1];
blc := cent1; blc[1] -:= rad; blc[2] -:= rad
trc := cent1; trc[1] +:= rad; trc[2] +:= rad
myarr := im.getchunk(blc=blc , trc=trc)
im.set(pixels=0.0);
#
mypos2 := 0 * myarr;
iy := 1:(2*rad+1)
for (ix in [1:(2*rad+1)]) {
   mypos2[ix,,1,1,1] := (ix - rad)^2 + (iy - rad)^2;
}
myarr[(mypos2 > rad^2)] := 0;

im.putchunk(myarr, blc=blc);
im.done();

mypos := F
myarr := F




dc.copy('PLAY.0', 'IM.dummy2');
im := image('IM.dummy2');
mysh :=  im.shape();
#  should be [2048, 2048, 1, 1, 1];
cent2 := [(1500-dx), (1000-dy), 1, 1, 1];
blc := cent2; blc[1] -:= rad; blc[2] -:= rad
trc := cent2; trc[1] +:= rad; trc[2] +:= rad
myarr := im.getchunk(blc=blc , trc=trc)
im.set(pixels=0.0);
#
mypos2 := 0 * myarr;
iy := 1:(2*rad+1)
for (ix in [1:(2*rad+1)]) {
   mypos2[ix,,1,1,1] := (ix - rad)^2 + (iy - rad)^2;
}
myarr[(mypos2 > rad^2)] := 0;

im.putchunk(myarr, blc=blc);
im.done();


mypos := F
myarr := F

}