local
double dval1, dval2;
FITSHead wcs;
int offscl;
int init=0;
end

if( !init++ ) FunInfoGet(fun, FUN_WCS, &wcs, 0);
wcs2pix(wcs, (double)cur->ra, (double)cur->dec, &dval1, &dval2, &offscl);
cur->nx:J:1:1024 = (int)(dval1+0.5);
cur->ny:J:1:1024 = (int)(dval2+0.5);

