There are a few programs with which orbits can be initialized:
% ccdmath out=ccd1 fie=%x+%y size=2,4 Generating a map from scratch % tsf ccd1 set Image set Parameters int Nx 2 int Ny 4 int Nz 1 double Xmin 0.00000 double Ymin 0.00000 double Zmin 0.00000 double Dx 1.00000 double Dy 1.00000 double Dz 1.00000 double MapMin -4.00000 double MapMax 0.00000 int BeamType 0 double Beamx 0.00000 double Beamy 0.00000 double Beamz 0.00000 double Time 0.00000 char Storage[5] "CDef" tes set Map double MapValues[2][4] -4.00000 -3.00000 -2.00000 -1.00000 -3.00000 -2.00000 -1.00000 0.00000 tes tes % ccdprint ccd1 x= y= label=x,y Y\X 0 1 3 -1 0 2 -2 -1 1 -3 -2 0 -4 -3
nx,ny -> data[nx][ny] e.g. ccdmath out=ccd1 nx=10 ny=5 gives double MapValues[10][5] ccdmath "" - %x 3,2 | tsf - margin=100 | grep MapVal MapValues[3][2] -2.00000 -2.00000 -1.00000 -1.00000 0.00000 0.00000 ccdmath "" - %y 3,2 | tsf - margin=100 | grep MapVal MapValues[3][2] -2.00000 -1.00000 -2.00000 -1.00000 -2.00000 -1.00000