#
include 'image.g';
#
doubleimage := function(in1, out1) {
  im := image(in1);
  oldcsys := im.coordsys()
  newshape := im.shape()
  newshape[1] := newshape[1] * 2
  oldchunk := im.getchunk()
  im.done();
#
  blc1:= 0 * newshape + 1
  newim := imagefromshape(outfile=out1, shape=newshape, csys=oldcsys)
  newim.putchunk(pixels=oldchunk, blc=blc1, replicate=T);
  newim.view();
#  works great!
#
  newim.done()
  oldchunk := F;
  return T;
}