#include "../h/which.h"
#ifdef PUMA
#include "../h/pumadata.h"
#endif
#ifdef STAN
#include "../h/standata.h"
#endif

#define RD      57.29577951308232100
#define CHK     20
main()
{
	int i;
	short e[6];
	double r[6], rr[6];


	r[0] = JCAL1; r[1] = JCAL2; r[2] = JCAL3;
	r[3] = JCAL4; r[4] = JCAL5; r[5] = JCAL6;
printf(">>>>> cal      %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f\n",
r[0] * RD, r[1] * RD, r[2] * RD, r[3] * RD, r[4] * RD, r[5] * RD);
		rngtoenc(e, r);
printf("rng > enc     %-6d %-6d %-6d %-6d %-6d %-6d\n",
e[0] & 0xffff, e[1] & 0xffff, e[2] & 0xffff, e[3] & 0xffff, e[4] & 0xffff, e[5]
		enctorng(rr, e);
printf("enc > rng      %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f\n",
rr[0] * RD, rr[1] * RD, rr[2] * RD, rr[3] * RD, rr[4] * RD, rr[5] * RD);
printf("off            %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f\n",
(r[0] - rr[0]) * RD, (r[1] - rr[1]) * RD, (r[2] - rr[2]) * RD,
(r[3] - rr[3]) * RD, (r[4] - rr[4]) * RD, (r[5] - rr[5]) * RD);
	r[0] = 0.; r[1] = 0.; r[2] = 0.;
	r[3] = 0.; r[4] = 0.; r[5] = 0.;

	for (i = 0; i < CHK + 1; ++i) {
printf(">>>>> rng      %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f\n",
r[0] * RD, r[1] * RD, r[2] * RD, r[3] * RD, r[4] * RD, r[5] * RD);
		rngtoenc(e, r);
printf("rng > enc     %-6d %-6d %-6d %-6d %-6d %-6d\n",
e[0] & 0xffff, e[1] & 0xffff, e[2] & 0xffff, e[3] & 0xffff, e[4] & 0xffff, e[5]
		enctorng(rr, e);
printf("enc > rng      %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f\n",
rr[0] * RD, rr[1] * RD, rr[2] * RD, rr[3] * RD, rr[4] * RD, rr[5] * RD);
printf("off            %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f  %-4.1f\n",
(r[0] - rr[0]) * RD, (r[1] - rr[1]) * RD, (r[2] - rr[2]) * RD,
(r[3] - rr[3]) * RD, (r[4] - rr[4]) * RD, (r[5] - rr[5]) * RD);

	r[0] += JRNG1 / CHK;
	r[1] += JRNG2 / CHK;
	r[2] += JRNG3 / CHK;
	r[3] += JRNG4 / CHK;
	r[4] += JRNG5 / CHK;
	r[5] += JRNG6 / CHK;

	}
}
