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