/* ./funcalc -f wcs.in $S foo.fits */
local
double dval1, dval2;
void *wcs;
int init=0;
end

if( !init++ ) FunInfoGet(fun, FUN_WCS, &wcs, 0);
pix2wcs(wcs, (double)cur->x, (double)cur->y, &dval1, &dval2);
cur->RA:D = dval1;
cur->Dec:D = dval2;