IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 35005


Ignore:
Timestamp:
Jan 21, 2013, 1:32:12 PM (13 years ago)
Author:
eugene
Message:

add -append option, add -extname option, fix Ndof, Npix (where allowed)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/eam_branches/ipp-20121219/Ohana/src/addstar/src/mkcmf.c

    r35002 r35005  
    4040  Coords coords;
    4141
     42  int APPEND = FALSE;
     43  if ((N = get_argument (argc, argv, "-append"))) {
     44    remove_argument (N, &argc, argv);
     45    APPEND = TRUE;
     46  }
     47
    4248  static char *photcode = "SIMTEST.r.Chip";
    4349  if ((N = get_argument (argc, argv, "-photcode"))) {
    … …  
    6672    remove_argument (N, &argc, argv);
    6773    time = strcreate (argv[N]);
     74    remove_argument (N, &argc, argv);
     75  }
     76
     77  char extname[80], exthead[80];
     78  strcpy (extname, "Chip.psf");
     79  strcpy (exthead, "Chip.hdr");
     80  if ((N = get_argument (argc, argv, "-extname"))) {
     81    remove_argument (N, &argc, argv);
     82    snprintf (extname, 80, "%s.psf", argv[N]);
     83    snprintf (exthead, 80, "%s.hdr", argv[N]);
    6884    remove_argument (N, &argc, argv);
    6985  }
    … …  
    275291
    276292  PutCoords (&coords, &header);
    277   gfits_modify (&header, "EXTNAME",   "%s", 1, "Chip.hdr");
     293  gfits_modify (&header, "EXTNAME",   "%s", 1, exthead);
    278294
    279295  ftable.header = &theader;
    … …  
    310326  }
    311327
    312   gfits_modify (ftable.header, "EXTHEAD",   "%s", 1, "Chip.hdr");
    313   gfits_modify (ftable.header, "EXTNAME",   "%s", 1, "Chip.psf");
     328  gfits_modify (ftable.header, "EXTHEAD",   "%s", 1, exthead);
     329  gfits_modify (ftable.header, "EXTNAME",   "%s", 1, extname);
    314330
    315331  // open for append (if you want a new file, need to blow the old object)
    316   fits = fopen (argv[2], "a+");
     332  if (APPEND) {
     333    fits = fopen (argv[2], "a+");
     334  } else {
     335    fits = fopen (argv[2], "w");
     336  }
    317337  if (fits == NULL) {
    318338    fprintf (stderr, "ERROR: can't open output file %s\n", argv[2]);
    … …  
    320340  }
    321341
    322   /* position to begining of file to write header */
    323   off_t Nend = ftello (f);
    324 
    325342  /* if we are appending, fix up header */
    326   if (Nend > 0) {
     343  if (APPEND) {
    327344    static char simple[] = "XTENSION= 'IMAGE  '            / Image extension";
    328345    int Ns, No;
    … …  
    542559    stars[i].dSky      = DSKY;
    543560    stars[i].psfChisq  = PSFCHI;
     561    stars[i].psfNdof   = 1;
     562    stars[i].psfNpix   = 2;
    544563    stars[i].crNsigma  = CRN;
    545564    stars[i].extNsigma = EXTN;
    … …  
    596615    stars[i].dSky      = DSKY;
    597616    stars[i].psfChisq  = PSFCHI;
     617    stars[i].psfNdof   = 1;
     618    stars[i].psfNpix   = 2;
    598619    stars[i].crNsigma  = CRN;
    599620    stars[i].extNsigma = EXTN;
    … …  
    650671    stars[i].dSky      = DSKY;
    651672    stars[i].psfChisq  = PSFCHI;
     673    stars[i].psfNdof   = 1;
     674    stars[i].psfNpix   = 2;
    652675    stars[i].crNsigma  = CRN;
    653676    stars[i].extNsigma = EXTN;
    … …  
    716739    stars[i].dSky      = DSKY;
    717740    stars[i].psfChisq  = PSFCHI;
     741    stars[i].psfNdof   = 1;
     742    stars[i].psfNpix   = 2;
    718743    stars[i].crNsigma  = CRN;
    719744    stars[i].extNsigma = EXTN;
Note: See TracChangeset for help on using the changeset viewer.