/*--------------------------------------------------------------------------*/
int ffiimgll(fitsfile *fptr, /* I - FITS file pointer */
int bitpix, /* I - bits per pixel */
int naxis, /* I - number of axes in the array */
LONGLONG *naxes, /* I - size of each axis */
int *status) /* IO - error status */
/*
insert an IMAGE extension following the current HDU
*/
{
int bytlen, nexthdu, maxhdu, ii, onaxis;
long nblocks;
LONGLONG npixels, newstart, datasize;
char errmsg[FLEN_ERRMSG], card[FLEN_CARD], naxiskey[FLEN_KEYWORD];
if (*status > 0)
return(*status);
if (fptr->HDUposition != (fptr->Fptr)->curhdu)
ffmahd(fptr, (fptr->HDUposition) + 1, NULL, status);
maxhdu = (fptr->Fptr)->maxhdu;
if (*status != PREPEND_PRIMARY)
{
/* if the current header is completely empty ... */
if (( (fptr->Fptr)->headend == (fptr->Fptr)->headstart[(fptr->Fptr)->curhdu])
/* or, if we are at the end of the file, ... */
|| ( (((fptr->Fptr)->curhdu) == maxhdu ) &&
((fptr->Fptr)->headstart[maxhdu + 1] >= (fptr->Fptr)->logfilesize ) ) )
{
/* then simply append new image extension */
ffcrimll(fptr, bitpix, naxis, naxes, status);
return(*status);
}
}
if (bitpix == 8)
bytlen = 1;
else if (bitpix == 16)
bytlen = 2;
else if (bitpix == 32 || bitpix == -32)
bytlen = 4;
else if (bitpix == 64 || bitpix == -64)
bytlen = 8;
else
{
sprintf(errmsg,
"Illegal value for BITPIX keyword: %d", bitpix);
ffpmsg(errmsg);
return(*status = BAD_BITPIX); /* illegal bitpix value */
}
if (naxis < 0 || naxis > 999)
{
sprintf(errmsg,
"Illegal value for NAXIS keyword: %d", naxis);
ffpmsg(errmsg);
return(*status = BAD_NAXIS);
}
for (ii = 0; ii < naxis; ii++)
{
if (naxes[ii] < 0)
{
sprintf(errmsg,
"Illegal value for NAXIS%d keyword: %ld", ii + 1, (long) naxes[ii]);
ffpmsg(errmsg);
return(*status = BAD_NAXES);
}
}
/* calculate number of pixels in the image */
if (naxis == 0)
npixels = 0;
else
npixels = naxes[0];
for (ii = 1; ii < naxis; ii++)
npixels = npixels * naxes[ii];
datasize = npixels * bytlen; /* size of image in bytes */
nblocks = (long) (((datasize + 2879) / 2880) + 1); /* +1 for the header */
if ((fptr->Fptr)->writemode == READWRITE) /* must have write access */
{ /* close the CHDU */
ffrdef(fptr, status); /* scan header to redefine structure */
ffpdfl(fptr, status); /* insure correct data file values */
}
else
return(*status = READONLY_FILE);
if (*status == PREPEND_PRIMARY)
{
/* inserting a new primary array; the current primary */
/* array must be transformed into an image extension. */
*status = 0;
ffmahd(fptr, 1, NULL, status); /* move to the primary array */
ffgidm(fptr, &onaxis, status);
if (onaxis > 0)
ffkeyn("NAXIS",onaxis, naxiskey, status);
else
strcpy(naxiskey, "NAXIS");
ffgcrd(fptr, naxiskey, card, status); /* read last NAXIS keyword */
ffikyj(fptr, "PCOUNT", 0, "required keyword", status); /* add PCOUNT and */
ffikyj(fptr, "GCOUNT", 1, "required keyword", status); /* GCOUNT keywords */
if (*status > 0)
return(*status);
if (ffdkey(fptr, "EXTEND", status) ) /* delete the EXTEND keyword */
*status = 0;
/* redefine internal structure for this HDU */
ffrdef(fptr, status);
/* insert space for the primary array */
if (ffiblk(fptr, nblocks, -1, status) > 0) /* insert the blocks */
return(*status);
nexthdu = 0; /* number of the new hdu */
newstart = 0; /* starting addr of HDU */
}
else
{
nexthdu = ((fptr->Fptr)->curhdu) + 1; /* number of the next (new) hdu */
newstart = (fptr->Fptr)->headstart[nexthdu]; /* save starting addr of HDU */
(fptr->Fptr)->hdutype = IMAGE_HDU; /* so that correct fill value is used */
/* ffiblk also increments headstart for all following HDUs */
if (ffiblk(fptr, nblocks, 1, status) > 0) /* insert the blocks */
return(*status);
}
((fptr->Fptr)->maxhdu)++; /* increment known number of HDUs in the file */
for (ii = (fptr->Fptr)->maxhdu; ii > (fptr->Fptr)->curhdu; ii--)
(fptr->Fptr)->headstart[ii + 1] = (fptr->Fptr)->headstart[ii]; /* incre start addr */
if (nexthdu == 0)
(fptr->Fptr)->headstart[1] = nblocks * 2880; /* start of the old Primary array */
(fptr->Fptr)->headstart[nexthdu] = newstart; /* set starting addr of HDU */
/* set default parameters for this new empty HDU */
(fptr->Fptr)->curhdu = nexthdu; /* we are now located at the next HDU */
fptr->HDUposition = nexthdu; /* we are now located at the next HDU */
(fptr->Fptr)->nextkey = (fptr->Fptr)->headstart[nexthdu];
(fptr->Fptr)->headend = (fptr->Fptr)->headstart[nexthdu];
(fptr->Fptr)->datastart = ((fptr->Fptr)->headstart[nexthdu]) + 2880;
(fptr->Fptr)->hdutype = IMAGE_HDU; /* might need to be reset... */
/* write the required header keywords */
ffphprll(fptr, TRUE, bitpix, naxis, naxes, 0, 1, TRUE, status);
/* redefine internal structure for this HDU */
ffrdef(fptr, status);
return(*status);
}
bool Fits2D::writeKeywords(fitsfile *fptr){
/*
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% FITS 2D Keywords template %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1. SIMPLE = T / file does conform to FITS standard
2. BITPIX = 8 / number of bits per pixel
3. NAXIS = 2 / number of data axes
4. NAXIS1 = 1280 / length of data axis 1
5. NAXIS2 = 960 / length of data axis 2
6. EXTEND = T / FITS dataset may contain extensions
7. FILENAME = 'stationOrsay_YYYYMMJJ_HHMMSS_UT.fits' / name of the fits file
8. DATE = 'YYYY-MM-JJT HH:MM:SS.SS' / date of the creation of the fits file
9. DATE-OBS = 'YYYY-MM-JJT HH:MM:SS.SS' / acquisition date of the first frame
10. OBS_MODE = 'SINGLE' / observation method used to get this fits file 'SINGLE' 'SUM' 'AVERAGE' ('MEDIAN')
11. ELAPTIME = 60 / end observation date - start observation date (sec.)
12. EXPOSURE = 0.033 / frame exp (sec). 999999 if different exposure.
13. ONTIME = 0.033 / total integration time, sum of each frame exp (sec)
14. FILTER = "NONE"
15. TELESCOP = "<Code station>" / station <stationName>
16. OBSERVER = "<responsable camera>"
17. INSTRUME = 'FRIPON-CAM'
18. CAMERA = 'BASLER 1300gm'
19. FOCAL = 1.25
20. APERTURE = 2.0
21. SITELONG = 2.1794397 / geodetic WGS84 longitude (DD)
22. SITELAT = 48.7063906 / geodetic WGS84 latitude (DD)
23. SITEELEV = 90 / geodetic WGS84 elevation (m)
24. XPIXEL = 3.75
25. YPIXEL = 3.75
26. GAINDB = 400 / detector gain
27. SATURATE = 4095 / saturation value
28. PROGRAM = 'FreeTure v0.1' / name of the acquisition software
29. CREATOR = 'FRIPON TEAM' / http://fripon.org
30. BZERO = 0
31. BSCALE = 1
32. RADESYS = 'ICRS'
33. TIMESYS = 'UTC'
34. EQUINOX = 2.000000000000E+03 / equinox of equatorial coordinates
35. CTYPE1 = 'RA---ARC' / projection and reference system
36. CTYPE2 = 'DEC--ARC' / projection and reference system
37. TIMEUNIT = 's '
38. CD1_1 = 0.0 / deg/px
39. CD1_2 = 0.17 / deg/px
40. CD2_1 = 0.17 / deg/pix
41. CD2_2 = 0.0 / deg/pix
42. CD3_3 = 30 / fps
43. CRPIX1 = 640
44. CRPIX2 = 480
45. CRVAL1 = / Sidereal time (decimal degree)
46. CRVAL2 = / geodetic WGS84 latitude (DD)
47. K1 =
48. K2 =
*/
int status = 0;
// DELETE DEFAULT COMMENTS.
if(ffdkey(fptr, "COMMENT", &status))
printerror( status );
if(ffdkey(fptr, "COMMENT", &status))
printerror( status );
/// 7. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% FILENAME %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * filename = new char[kFILENAME.length()+1];
strcpy(filename,kFILENAME.c_str());
char * cfilename = new char[cFILENAME.length()+1];
strcpy(cfilename,cFILENAME.c_str());
if(fits_write_key(fptr, TSTRING, "FILENAME", filename, cfilename, &status)){
delete filename;
delete cfilename;
printerror(status, "Error fits_write_key(FILENAME)");
return false;
}
delete cfilename;
delete filename;
/// 8. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% DATE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * date = new char[kDATE.length()+1];
strcpy(date,kDATE.c_str());
char * cdate = new char[cDATE.length()+1];
strcpy(cdate,cDATE.c_str());
if(fits_write_key(fptr,TSTRING,"DATE",date,cdate,&status)){
delete date;
delete cdate;
printerror(status, "Error fits_write_key(DATE)");
return false;
}
delete cdate;
delete date;
/// 9. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% DATE-OBS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * dateobs = new char[kDATEOBS.length()+1];
strcpy(dateobs,kDATEOBS.c_str());
char * cdateobs = new char[cDATEOBS.length()+1];
strcpy(cdateobs,cDATEOBS.c_str());
if(fits_write_key(fptr,TSTRING,"DATE-OBS",dateobs,cdateobs,&status)){
delete dateobs;
delete cdateobs;
printerror(status, "Error fits_write_key(DATE-OBS)");
return false;
}
delete cdateobs;
delete dateobs;
/// 10. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% OBS_MODE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cobsmode = new char[cOBSMODE.length()+1];
strcpy(cobsmode,cOBSMODE.c_str());
char * obsmode = new char[kOBSMODE.length()+1];
strcpy(obsmode,kOBSMODE.c_str());
if(fits_write_key(fptr,TSTRING,"OBS_MODE",obsmode,cobsmode,&status)){
delete cobsmode;
delete obsmode;
printerror(status, "Error fits_write_key(OBS_MODE)");
return false;
}
delete cobsmode;
delete obsmode;
/// 11. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% ELAPTIME %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * celaptime = new char[cELAPTIME.length()+1];
strcpy(celaptime,cELAPTIME.c_str());
if(fits_write_key(fptr,TDOUBLE,"ELAPTIME",&kELAPTIME,celaptime,&status)){
delete celaptime;
printerror(status, "Error fits_write_key(ELAPTIME)");
return false;
}
delete celaptime;
/// 12. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% EXPOSURE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ceposure = new char[cEXPOSURE.length()+1];
strcpy(ceposure,cEXPOSURE.c_str());
if(fits_write_key(fptr,TDOUBLE,"EXPOSURE",&kEXPOSURE,ceposure,&status)){
delete ceposure;
printerror(status, "Error fits_write_key(EXPOSURE)");
return false;
}
delete ceposure;
/// 13. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% ONTIME %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * contime = new char[cONTIME.length()+1];
strcpy(contime,cONTIME.c_str());
if(fits_write_key(fptr,TDOUBLE,"ONTIME",&kONTIME,contime,&status)){
delete contime;
printerror(status, "Error fits_write_key(ONTIME)");
return false;
}
delete contime;
/// 14. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% FILTER %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cfilter = new char[cFILTER.length()+1];
strcpy(cfilter,cFILTER.c_str());
char * f = new char[kFILTER.length()+1];
strcpy(f,kFILTER.c_str());
if(fits_write_key(fptr,TSTRING,"FILTER",f,cfilter,&status)){
delete cfilter;
delete f;
printerror(status, "Error fits_write_key(FILTER)");
return false;
}
delete cfilter;
delete f;
/// 15. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% TELESCOP %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ctelescop = new char[cTELESCOP.length()+1];
strcpy(ctelescop,cTELESCOP.c_str());
char * t = new char[kTELESCOP.length()+1];
strcpy(t,kTELESCOP.c_str());
if(fits_write_key(fptr,TSTRING,"TELESCOP",t,ctelescop,&status)){
delete ctelescop;
delete t;
printerror(status, "Error fits_write_key(TELESCOP)");
return false;
}
delete ctelescop;
delete t;
/// 16. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% OBSERVER %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cobserver = new char[cOBSERVER.length()+1];
strcpy(cobserver,cOBSERVER.c_str());
char * o = new char[kOBSERVER.length()+1];
strcpy(o,kOBSERVER.c_str());
if(fits_write_key(fptr,TSTRING,"OBSERVER",o,cobserver,&status)){
delete cobserver;
delete o;
printerror(status, "Error fits_write_key(OBSERVER)");
return false;
}
delete cobserver;
delete o;
/// 17. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% INSTRUME %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cinstrume = new char[cINSTRUME.length()+1];
strcpy(cinstrume,cINSTRUME.c_str());
char * i = new char[kINSTRUME.length()+1];
strcpy(i,kINSTRUME.c_str());
if(fits_write_key(fptr,TSTRING,"INSTRUME",i,cinstrume,&status)){
delete cinstrume;
delete i;
printerror(status, "Error fits_write_key(OBSERVER)");
return false;
}
delete cinstrume;
delete i;
/// 18. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CAMERA %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccamera = new char[cCAMERA.length()+1];
strcpy(ccamera,cCAMERA.c_str());
char * cam = new char[kCAMERA.length()+1];
strcpy(cam,kCAMERA.c_str());
if(fits_write_key(fptr,TSTRING,"CAMERA",cam,ccamera,&status)){
delete ccamera;
delete cam;
printerror(status, "Error fits_write_key(CAMERA)");
return false;
}
delete ccamera;
delete cam;
/// 19. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% FOCAL %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cfocal = new char[cFOCAL.length()+1];
strcpy(cfocal,cFOCAL.c_str());
if(fits_write_key(fptr,TDOUBLE,"FOCAL",&kFOCAL,cfocal,&status)){
delete cfocal;
printerror(status, "Error fits_write_key(FOCAL)");
return false;
}
delete cfocal;
/// 20. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% APERTURE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * caperture = new char[cAPERTURE.length()+1];
strcpy(caperture,cAPERTURE.c_str());
if(fits_write_key(fptr,TDOUBLE,"APERTURE",&kAPERTURE,caperture,&status)){
delete caperture;
printerror(status, "Error fits_write_key(APERTURE)");
return false;
}
delete caperture;
/// 21. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% SITELONG %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * csitelong = new char[cSITELONG.length()+1];
strcpy(csitelong,cSITELONG.c_str());
if(fits_write_key(fptr,TDOUBLE,"SITELONG",&kSITELONG,csitelong,&status)){
delete csitelong;
printerror(status, "Error fits_write_key(APERTURE)");
return false;
}
delete csitelong;
/// 22. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% SITELAT %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * csitelat = new char[cSITELAT.length()+1];
strcpy(csitelat,cSITELAT.c_str());
if(fits_write_key(fptr,TDOUBLE,"SITELAT",&kSITELAT,csitelat,&status)){
delete csitelat;
printerror(status, "Error fits_write_key(SITELAT)");
return false;
}
delete csitelat;
/// 23. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% SITEELEV %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * csiteelev = new char[cSITEELEV.length()+1];
strcpy(csiteelev,cSITEELEV.c_str());
if(fits_write_key(fptr,TDOUBLE,"SITEELEV",&kSITEELEV,csiteelev,&status)){
delete csiteelev;
printerror(status, "Error fits_write_key(SITEELEV)");
return false;
}
delete csiteelev;
/// 24. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% XPIXEL %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cxpixel = new char[cXPIXEL.length()+1];
strcpy(cxpixel,cXPIXEL.c_str());
if(fits_write_key(fptr,TDOUBLE,"XPIXEL",&kXPIXEL,cxpixel,&status)){
delete cxpixel;
printerror(status, "Error fits_write_key(XPIXEL)");
return false;
}
delete cxpixel;
/// 25. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% YPIXEL %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cypixel = new char[cYPIXEL.length()+1];
strcpy(cypixel,cYPIXEL.c_str());
if(fits_write_key(fptr,TDOUBLE,"YPIXEL",&kYPIXEL,cypixel,&status)){
delete cypixel;
printerror(status, "Error fits_write_key(YPIXEL)");
return false;
}
delete cypixel;
/// 26. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% GAINDB %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cgaindb = new char[cGAINDB.length()+1];
strcpy(cgaindb,cGAINDB.c_str());
if(fits_write_key(fptr,TINT,"GAINDB",&kGAINDB,cgaindb,&status)){
delete cgaindb;
printerror(status, "Error fits_write_key(GAINDB)");
return false;
}
delete cgaindb;
/// 27. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% SATURATE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * csaturate = new char[cSATURATE.length()+1];
strcpy(csaturate,cSATURATE.c_str());
if(fits_write_key(fptr,TDOUBLE,"SATURATE",&kSATURATE,csaturate,&status)){
delete csaturate;
printerror(status, "Error fits_write_key(SATURATE)");
return false;
}
delete csaturate;
/// 28. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% PROGRAM %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cprograme = new char[cPROGRAM.length()+1];
strcpy(cprograme,cPROGRAM.c_str());
char * p = new char[kPROGRAM.length()+1];
strcpy(p,kPROGRAM.c_str());
if(fits_write_key(fptr,TSTRING,"PROGRAM",p,cprograme,&status)){
delete cprograme;
delete p;
printerror(status, "Error fits_write_key(PROGRAM)");
return false;
}
delete cprograme;
delete p;
/// 29. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CREATOR %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccreator = new char[cCREATOR.length()+1];
strcpy(ccreator,cCREATOR.c_str());
char * c = new char[kCREATOR.length()+1];
strcpy(c,kCREATOR.c_str());
if(fits_write_key(fptr,TSTRING,"CREATOR",c,ccreator,&status)){
delete ccreator;
delete c;
printerror(status, "Error fits_write_key(CREATOR)");
return false;
}
delete ccreator;
delete c;
/// 30. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% BZERO %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cbzero = new char[cBZERO.length()+1];
strcpy(cbzero,cBZERO.c_str());
if(fits_write_key(fptr,TDOUBLE,"BZERO",&kBZERO,cbzero,&status)){
delete cbzero;
printerror(status, "Error fits_write_key(BZERO)");
return false;
}
delete cbzero;
/// 31. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% BSCALE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cbscale = new char[cBSCALE.length()+1];
strcpy(cbscale,cBSCALE.c_str());
if(fits_write_key(fptr,TDOUBLE,"BSCALE",&kBSCALE,cbscale,&status)){
delete cbscale;
printerror(status, "Error fits_write_key(BSCALE)");
return false;
}
delete cbscale;
/// 32. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% RADESYS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * radesys = new char[kRADESYS.length()+1];
strcpy(radesys,kRADESYS.c_str());
char * cradesys = new char[cRADESYS.length()+1];
strcpy(cradesys,cRADESYS.c_str());
if(fits_write_key(fptr,TSTRING,"RADESYS",radesys,cradesys,&status)){
delete cradesys;
delete radesys;
printerror(status, "Error fits_write_key(RADESYS)");
return false;
}
delete cradesys;
delete radesys;
/// 33. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% TIMESYS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ctimesys = new char[cTIMESYS.length()+1];
strcpy(ctimesys,cTIMESYS.c_str());
char * timesys = new char[kTIMESYS.length()+1];
strcpy(timesys,kTIMESYS.c_str());
if(fits_write_key(fptr,TSTRING,"TIMESYS",timesys,ctimesys,&status)){
delete ctimesys;
delete timesys;
printerror(status, "Error fits_write_key(TIMESYS)");
return false;
}
delete ctimesys;
delete timesys;
/// 34. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% EQUINOX %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * cequinox = new char[cEQUINOX.length()+1];
strcpy(cequinox,cEQUINOX.c_str());
if(fits_write_key(fptr,TDOUBLE,"EQUINOX",&kEQUINOX,cequinox,&status)){
delete cequinox;
printerror(status, "Error fits_write_key(EQUINOX)");
return false;
}
delete cequinox;
/// 35. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CTYPE1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ctype1 = new char[cCTYPE1.length()+1];
strcpy(ctype1,cCTYPE1.c_str());
char * ktype1 = new char[kCTYPE1.length()+1];
strcpy(ktype1,kCTYPE1.c_str());
if(fits_write_key(fptr,TSTRING,"CTYPE1",ktype1,ctype1,&status)){
delete ctype1;
delete ktype1;
printerror(status, "Error fits_write_key(CTYPE1)");
return false;
}
delete ctype1;
delete ktype1;
/// 36. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CTYPE2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ctype2 = new char[cCTYPE2.length()+1];
strcpy(ctype2,cCTYPE2.c_str());
char * ktype2 = new char[kCTYPE2.length()+1];
strcpy(ktype2,kCTYPE2.c_str());
if(fits_write_key(fptr,TSTRING,"CTYPE2",ktype2,ctype2,&status)){
delete ctype2;
delete ktype2;
printerror(status, "Error fits_write_key(CTYPE2)");
return false;
}
delete ctype2;
delete ktype2;
/// 37. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% TIMEUNIT %%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ctimeunit = new char[cTIMEUNIT.length()+1];
strcpy(ctimeunit,cTIMEUNIT.c_str());
char * ktimeunit = new char[kTIMEUNIT.length()+1];
strcpy(ktimeunit,kTIMEUNIT.c_str());
if(fits_write_key(fptr,TSTRING,"TIMEUNIT",ktimeunit,ctimeunit,&status)){
delete ctimeunit;
delete ktimeunit;
printerror(status, "Error fits_write_key(TIMEUNIT)");
return false;
}
delete ctimeunit;
delete ktimeunit;
/// 38. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CD1_1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccd1_1 = new char[cCD1_1.length()+1];
strcpy(ccd1_1,cCD1_1.c_str());
if(fits_write_key(fptr,TDOUBLE,"CD1_1",&kCD1_1,ccd1_1,&status)){
delete ccd1_1;
printerror(status, "Error fits_write_key(CD1_1)");
return false;
}
delete ccd1_1;
/// 39. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CD1_2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccd1_2 = new char[cCD1_2.length()+1];
strcpy(ccd1_2,cCD1_2.c_str());
if(fits_write_key(fptr,TDOUBLE,"CD1_2",&kCD1_2,ccd1_2,&status)){
delete ccd1_2;
printerror(status, "Error fits_write_key(CD1_2)");
return false;
}
delete ccd1_2;
/// 40. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CD2_1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccd2_1 = new char[cCD2_1.length()+1];
strcpy(ccd2_1,cCD2_1.c_str());
if(fits_write_key(fptr,TDOUBLE,"CD2_1",&kCD2_1,ccd2_1,&status)){
delete ccd2_1;
printerror(status, "Error fits_write_key(CD2_1)");
return false;
}
delete ccd2_1;
/// 41. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CD2_2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccd2_2 = new char[cCD2_2.length()+1];
strcpy(ccd2_2,cCD2_2.c_str());
if(fits_write_key(fptr,TDOUBLE,"CD2_2",&kCD2_2,ccd2_2,&status)){
delete ccd2_2;
printerror(status, "Error fits_write_key(CD2_2)");
return false;
}
delete ccd2_2;
/// 42. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CD3_3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccd3_3 = new char[cCD3_3.length()+1];
strcpy(ccd3_3,cCD3_3.c_str());
if(fits_write_key(fptr,TDOUBLE,"CD3_3",&kCD3_3,ccd3_3,&status)){
delete ccd3_3;
printerror(status, "Error fits_write_key(CD3_3)");
return false;
}
delete ccd3_3;
/// 43. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CRPIX1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccrpix1 = new char[cCRPIX1.length()+1];
strcpy(ccrpix1,cCRPIX1.c_str());
if(fits_write_key(fptr,TINT,"CRPIX1",&kCRPIX1,ccrpix1,&status)){
delete ccrpix1;
printerror(status, "Error fits_write_key(CRPIX1)");
return false;
}
delete ccrpix1;
/// 44. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CRPIX2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccrpix2 = new char[cCRPIX2.length()+1];
strcpy(ccrpix2,cCRPIX2.c_str());
if(fits_write_key(fptr,TINT,"CRPIX2",&kCRPIX2,ccrpix2,&status)){
delete ccrpix2;
printerror(status, "Error fits_write_key(CRPIX2)");
return false;
}
delete ccrpix2;
/// 45. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CRVAL1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccrval1 = new char[cCRVAL1.length()+1];
strcpy(ccrval1,cCRVAL1.c_str());
if(fits_write_key(fptr,TDOUBLE,"CRVAL1",&kCRVAL1,ccrval1,&status)){
delete ccrval1;
printerror(status, "Error fits_write_key(CRVAL1)");
return false;
}
delete ccrval1;
/// 46. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% CRVAL2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ccrval2 = new char[cCRVAL2.length()+1];
strcpy(ccrval2,cCRVAL2.c_str());
if(fits_write_key(fptr,TDOUBLE,"CRVAL2",&kSITELAT,ccrval2,&status)){
delete ccrval2;
printerror(status, "Error fits_write_key(CRVAL2)");
return false;
}
delete ccrval2;
/// 47. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% K1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ck1 = new char[cK1.length()+1];
strcpy(ck1,cK1.c_str());
if(fits_write_key(fptr,TDOUBLE,"K1",&kK1,ck1,&status)){
delete ck1;
printerror(status, "Error fits_write_key(K1)");
return false;
}
delete ck1;
/// 48. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% K2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
char * ck2 = new char[cK2.length()+1];
strcpy(ck2,cK2.c_str());
if(fits_write_key(fptr,TDOUBLE,"K2",&kK2,ck2,&status)){
delete ck2;
printerror(status, "Error fits_write_key(K2)");
return false;
}
delete ck2;
return true;
}