IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Sep 5, 2012, 4:25:56 PM (14 years ago)
Author:
eugene
Message:

merging change from eam_branches/ipp-20120805: adding psfQF, psfQFperf to both average and measure; deprecate average.Xp

Location:
trunk/Ohana
Files:
18 edited
1 copied

Legend:

Unmodified
Added
Removed

  • trunk/Ohana

  • trunk/Ohana/src/libdvo/include/dvo.h

    r34291 r34405  
    623623int dvo_image_save_raw (FITS_DB *db, int VERBOSE);
    624624int dvo_image_addrows (FITS_DB *db, Image *new, off_t Nnew);
     625int dvo_image_createID (Header *header);
    625626void dvo_image_create (FITS_DB *db, double ZeroPoint);
    626627

  • trunk/Ohana/src/libdvo/src/dbExtractAverages.c

    r34260 r34405  
    147147    case AVE_TRANGE:
    148148      value.Flt = TimeValue (average[0].Trange, 0, TimeFormat);
    149       break;
    150 
    151     case AVE_Xp:
    152       value.Flt = 0.01*average[0].Xp;
    153149      break;
    154150

  • trunk/Ohana/src/libdvo/src/dbExtractMeasures.c

    r34260 r34405  
    303303    case MEAS_SECFILT_FLAGS: /* OK */
    304304      equiv = GetPhotcodeEquivbyCode (measure[0].photcode);
     305      if (!equiv) break;
    305306      Nsec = GetPhotcodeNsec (equiv->code);
    306307      if (Nsec == -1) break;
    … …  
    479480      break;
    480481    case MEAS_PSF_QF: /* OK */
    481       value.Flt = measure[0].psfQual;
     482      value.Flt = measure[0].psfQF;
    482483      break;
    483484    case MEAS_PSF_QF_PERFECT: /* OK */

  • trunk/Ohana/src/libdvo/src/dbFields.c

    r34260 r34405  
    335335  if (!strcasecmp (fieldName, "EXTID_LO"))  ESCAPE (AVE_EXTID_LO,  MAG_NONE, OPIHI_INT);
    336336
    337   // if (!strcasecmp (fieldName, "Xp"))      ESCAPE (AVE_Xp,        MAG_NONE, OPIHI_FLT);
    338   // if (!strcasecmp (fieldName, "TYPE"))    ESCAPE (AVE_TYPE,      MAG_NONE, OPIHI_INT);
    339   // for words that don't parse, try a photcode
    340 
    341337  // check for code:mode in photcode name
    342338  code = ParsePhotcodeField (fieldName, &mode, MAG_AVE);

  • trunk/Ohana/src/libdvo/src/dvo_catalog.c

    r34260 r34405  
    8989  average->Trange          = 0;
    9090
    91   average->Xp              = 0;
     91  average->stargal         = 0.0;
    9292  average->Npos            = 0;
    9393
    … …  
    174174 measure->extID     = 0;
    175175
    176  measure->psfQual   = NAN;
     176 measure->psfQF     = NAN;
    177177 measure->psfChisq  = NAN;
    178178 measure->psfNdof   = 0;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_DEV_1.c

    r34260 r34405  
    2727    out[i].detID      = in[i].detID;
    2828    out[i].imageID    = in[i].imageID;
    29     out[i].psfQual    = in[i].psfQual;
     29    out[i].psfQF      = in[i].psfQF;
    3030    out[i].psfChisq   = in[i].psfChisq;
    3131    out[i].crNsigma   = in[i].crNsigma;
    … …  
    6565    out[i].FluxKron   = NAN;
    6666    out[i].dFluxKron  = NAN;
     67    out[i].psfQFperf  = NAN;
    6768  }
    6869  return (out);
    … …  
    9394    out[i].detID      = in[i].detID;
    9495    out[i].imageID    = in[i].imageID;
    95     out[i].psfQual    = in[i].psfQual;
     96    out[i].psfQF      = in[i].psfQF;
    9697    out[i].psfChisq   = in[i].psfChisq;
    9798    out[i].crNsigma   = in[i].crNsigma;
    … …  
    126127    out[i].R             = in[i].R;     
    127128    out[i].D             = in[i].D;     
    128     out[i].Xp            = in[i].Xp;     
    129129    out[i].dR            = in[i].dR;
    130130    out[i].dD            = in[i].dD;
    … …  
    159159    out[i].photFlagsUpper = 0;
    160160    out[i].photFlagsLower = 0;
     161
     162    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     163    out[i].psfQF         = NAN;
     164    out[i].psfQFperf     = NAN;
     165    out[i].stargal       = NAN;
    161166  }
    162167  return (out);
    … …  
    174179    out[i].R             = in[i].R;     
    175180    out[i].D             = in[i].D;     
    176     out[i].Xp            = in[i].Xp;     
    177181    out[i].dR            = in[i].dR;
    178182    out[i].dD            = in[i].dD;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_DEV_2.c

    r34260 r34405  
    2929    out[i].detID      = in[i].detID;
    3030    out[i].imageID    = in[i].imageID;
    31     out[i].psfQual    = in[i].psfQual;
     31    out[i].psfQF      = in[i].psfQF;
    3232    out[i].psfChisq   = in[i].psfChisq;
    3333    out[i].crNsigma   = in[i].crNsigma;
    … …  
    6262    out[i].FluxKron   = NAN;
    6363    out[i].dFluxKron  = NAN;
    64   }
     64    out[i].psfQFperf  = NAN;
     65 }
    6566  return (out);
    6667}
    … …  
    9293    out[i].detID      = in[i].detID;
    9394    out[i].imageID    = in[i].imageID;
    94     out[i].psfQual    = in[i].psfQual;
     95    out[i].psfQF      = in[i].psfQF;
    9596    out[i].psfChisq   = in[i].psfChisq;
    9697    out[i].crNsigma   = in[i].crNsigma;
    … …  
    130131    out[i].P             = in[i].P;
    131132    out[i].dP            = in[i].dP;
    132     out[i].Xp            = in[i].Xp;     
    133133    out[i].Nmeasure      = in[i].Nmeasure;     
    134134    out[i].Nmissing      = in[i].Nmissing;     
    … …  
    153153    out[i].photFlagsUpper = 0;
    154154    out[i].photFlagsLower = 0;
     155
     156    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     157    out[i].psfQF         = NAN;
     158    out[i].psfQFperf     = NAN;
     159    out[i].stargal       = NAN;
    155160  }
    156161  return (out);
    … …  
    176181    out[i].P             = in[i].P;
    177182    out[i].dP            = in[i].dP;
    178     out[i].Xp            = in[i].Xp;     
    179183    out[i].Nmeasure      = in[i].Nmeasure;     
    180184    out[i].Nmissing      = in[i].Nmissing;     

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_REF.c

    r34260 r34405  
    3333    out[i].catID      = in[i].catID;
    3434    out[i].extID      = 0;
    35     out[i].psfQual    = 1.0; // we do not save a value, but some operations (relphot/average) filter on this value..
     35    out[i].psfQF      = 1.0; // we do not save a value, but some operations (relphot/average) filter on this value..
    3636    out[i].psfChisq   = NAN;
    3737    out[i].psfNdof    = 0;
    … …  
    6060    out[i].FluxKron   = NAN;
    6161    out[i].dFluxKron  = NAN;
     62    out[i].psfQFperf  = NAN;
    6263  }
    6364  return (out);
    … …  
    105106    out[i].P             = NAN;
    106107    out[i].dP            = NAN;
    107     out[i].Xp            = NAN_S_SHORT;
    108108    out[i].ChiSqAve      = NAN;
    109109    out[i].ChiSqPM       = NAN;
    … …  
    126126    out[i].photFlagsUpper = 0;
    127127    out[i].photFlagsLower = 0;
     128
     129    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     130    out[i].psfQF         = NAN;
     131    out[i].psfQFperf     = NAN;
     132    out[i].stargal       = NAN;
    128133  }
    129134  return (out);

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V1.c

    r34260 r34405  
    3333    out[i].catID      = in[i].catID;
    3434    out[i].extID      = in[i].extID;
    35     out[i].psfQual    = in[i].psfQual;
     35    out[i].psfQF      = in[i].psfQF;
    3636    out[i].psfChisq   = in[i].psfChisq;
    3737    out[i].psfNdof    = in[i].psfNdof;
    … …  
    6060    out[i].FluxKron   = NAN;
    6161    out[i].dFluxKron  = NAN;
     62    out[i].psfQFperf  = NAN;
    6263  }
    6364  return (out);
    … …  
    9495    out[i].catID      = in[i].catID;
    9596    out[i].extID      = in[i].extID;
    96     out[i].psfQual    = in[i].psfQual;
     97    out[i].psfQF      = in[i].psfQF;
    9798    out[i].psfChisq   = in[i].psfChisq;
    9899    out[i].psfNdof    = in[i].psfNdof;
    … …  
    137138    out[i].P             = in[i].P;
    138139    out[i].dP            = in[i].dP;
    139     out[i].Xp            = in[i].Xp;     
    140140    out[i].ChiSqAve      = in[i].ChiSq;     
    141141    out[i].ChiSqPM       = NAN;
    … …  
    158158    out[i].photFlagsUpper = 0;
    159159    out[i].photFlagsLower = 0;
     160
     161    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     162    out[i].psfQF         = NAN;
     163    out[i].psfQFperf     = NAN;
     164    out[i].stargal       = NAN;
    160165  }
    161166  return (out);
    … …  
    181186    out[i].P             = in[i].P;
    182187    out[i].dP            = in[i].dP;
    183     out[i].Xp            = in[i].Xp;     
    184188    out[i].ChiSq         = in[i].ChiSqAve;     
    185189    out[i].Npos          = in[i].Npos;     

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V2.c

    r34260 r34405  
    3333    out[i].catID      = in[i].catID;
    3434    out[i].extID      = in[i].extID;
    35     out[i].psfQual    = in[i].psfQual;
     35    out[i].psfQF      = in[i].psfQF;
    3636    out[i].psfChisq   = in[i].psfChisq;
    3737    out[i].psfNdof    = in[i].psfNdof;
    … …  
    6161    out[i].FluxKron   = NAN;
    6262    out[i].dFluxKron  = NAN;
     63    out[i].psfQFperf  = NAN;
    6364  }
    6465  return (out);
    … …  
    9596    out[i].catID      = in[i].catID;
    9697    out[i].extID      = in[i].extID;
    97     out[i].psfQual    = in[i].psfQual;
     98    out[i].psfQF      = in[i].psfQF;
    9899    out[i].psfChisq   = in[i].psfChisq;
    99100    out[i].psfNdof    = in[i].psfNdof;
    … …  
    138139    out[i].P             = in[i].P;
    139140    out[i].dP            = in[i].dP;
    140     out[i].Xp            = in[i].Xp;     
    141141    out[i].ChiSqAve      = in[i].ChiSqAve;     
    142142    out[i].ChiSqPM       = in[i].ChiSqPM;     
    … …  
    159159    out[i].photFlagsUpper = 0;
    160160    out[i].photFlagsLower = 0;
     161
     162    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     163    out[i].psfQF         = NAN;
     164    out[i].psfQFperf     = NAN;
     165    out[i].stargal       = NAN;
    161166  }
    162167  return (out);
    … …  
    182187    out[i].P             = in[i].P;
    183188    out[i].dP            = in[i].dP;
    184     out[i].Xp            = in[i].Xp;     
    185189    out[i].ChiSqAve      = in[i].ChiSqAve;     
    186190    out[i].ChiSqPM       = in[i].ChiSqPM;     

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V3.c

    r34260 r34405  
    3333    out[i].catID      = in[i].catID;
    3434    out[i].extID      = in[i].extID;
    35     out[i].psfQual    = in[i].psfQual;
     35    out[i].psfQF      = in[i].psfQF;
    3636    out[i].psfChisq   = in[i].psfChisq;
    3737    out[i].psfNdof    = in[i].psfNdof;
    … …  
    6161    out[i].FluxKron   = NAN;
    6262    out[i].dFluxKron  = NAN;
     63    out[i].psfQFperf  = NAN;
    6364  }
    6465  return (out);
    … …  
    9596    out[i].catID      = in[i].catID;
    9697    out[i].extID      = in[i].extID;
    97     out[i].psfQual    = in[i].psfQual;
     98    out[i].psfQF      = in[i].psfQF;
    9899    out[i].psfChisq   = in[i].psfChisq;
    99100    out[i].psfNdof    = in[i].psfNdof;
    … …  
    138139    out[i].P             = in[i].P;
    139140    out[i].dP            = in[i].dP;
    140     out[i].Xp            = in[i].Xp;     
    141141    out[i].ChiSqAve      = in[i].ChiSqAve;     
    142142    out[i].ChiSqPM       = in[i].ChiSqPM;     
    … …  
    157157    out[i].catID         = in[i].catID;
    158158    out[i].extID         = in[i].extID;
     159
     160    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     161    out[i].psfQF         = NAN;
     162    out[i].psfQFperf     = NAN;
     163    out[i].stargal       = NAN;
    159164  }
    160165  return (out);
    … …  
    180185    out[i].P             = in[i].P;
    181186    out[i].dP            = in[i].dP;
    182     out[i].Xp            = in[i].Xp;     
    183187    out[i].ChiSqAve      = in[i].ChiSqAve;     
    184188    out[i].ChiSqPM       = in[i].ChiSqPM;     

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V4.c

    r34260 r34405  
    3939    out[i].catID      = in[i].catID;
    4040    out[i].extID      = in[i].extID;
    41     out[i].psfQual    = in[i].psfQual;
     41    out[i].psfQF      = in[i].psfQF;
     42    out[i].psfQFperf  = in[i].psfQFperf;
    4243    out[i].psfChisq   = in[i].psfChisq;
    4344    out[i].psfNdof    = in[i].psfNdof;
    … …  
    99100    out[i].catID      = in[i].catID;
    100101    out[i].extID      = in[i].extID;
    101     out[i].psfQual    = in[i].psfQual;
     102    out[i].psfQF      = in[i].psfQF;
     103    out[i].psfQFperf  = in[i].psfQFperf;
    102104    out[i].psfChisq   = in[i].psfChisq;
    103105    out[i].psfNdof    = in[i].psfNdof;
    … …  
    142144    out[i].P             = in[i].P;
    143145    out[i].dP            = in[i].dP;
    144     out[i].Xp            = in[i].Xp;     
     146    out[i].psfQF         = in[i].psfQF;
     147    out[i].psfQFperf     = in[i].psfQFperf;
     148    out[i].stargal       = in[i].stargal;     
    145149    out[i].ChiSqAve      = in[i].ChiSqAve;     
    146150    out[i].ChiSqPM       = in[i].ChiSqPM;     
    … …  
    184188    out[i].P             = in[i].P;
    185189    out[i].dP            = in[i].dP;
    186     out[i].Xp            = in[i].Xp;     
     190    out[i].psfQF         = in[i].psfQF;
     191    out[i].psfQFperf     = in[i].psfQFperf;
     192    out[i].stargal       = in[i].stargal;     
    187193    out[i].ChiSqAve      = in[i].ChiSqAve;     
    188194    out[i].ChiSqPM       = in[i].ChiSqPM;     

  • trunk/Ohana/src/libdvo/src/dvo_convert_elixir.c

    r34260 r34405  
    3535    out[i].Sky       = 0;
    3636    out[i].dSky      = 0;
    37     out[i].psfQual   = 0;
     37    out[i].psfQF     = 0;
    3838    out[i].psfChisq  = 0;
    3939    out[i].crNsigma  = 0;
    … …  
    7272    out[i].FluxKron   = NAN;
    7373    out[i].dFluxKron  = NAN;
     74    out[i].psfQFperf  = NAN;
    7475  }
    7576  return (out);
    … …  
    127128    out[i].R             = in[i].R;     
    128129    out[i].D             = in[i].D;     
    129     out[i].Xp            = in[i].Xp;     
    130130
    131131    // changed for PANSTARRS_DEV_0 (moved from Average to Measure)
    … …  
    170170    out[i].photFlagsUpper = 0;
    171171    out[i].photFlagsLower = 0;
     172
     173    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     174    out[i].psfQF         = NAN;
     175    out[i].psfQFperf     = NAN;
     176    out[i].stargal       = NAN;
    172177  }
    173178  return (out);
    … …  
    185190    out[i].R       = in[i].R;     
    186191    out[i].D       = in[i].D;     
    187     out[i].Xp      = in[i].Xp;     
    188192
    189193    // changed for PANSTARRS_DEV_0 (moved from Average to Measure)

  • trunk/Ohana/src/libdvo/src/dvo_convert_loneos.c

    r34260 r34405  
    3939    out[i].Sky       = 0;
    4040    out[i].dSky      = 0;
    41     out[i].psfQual   = 0;
     41    out[i].psfQF     = 0;
    4242    out[i].psfChisq  = 0;
    4343    out[i].crNsigma  = 0;
    … …  
    7474    out[i].FluxKron   = NAN;
    7575    out[i].dFluxKron  = NAN;
     76    out[i].psfQFperf  = NAN;
    7677  }
    7778  return (out);
    … …  
    118119    out[i].R             = in[i].R;     
    119120    out[i].D             = in[i].D;     
    120     out[i].Xp            = in[i].Xp;     
    121121
    122122    // added for ELIXIR
    … …  
    162162    out[i].photFlagsUpper = 0;
    163163    out[i].photFlagsLower = 0;
     164
     165    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     166    out[i].psfQF         = NAN;
     167    out[i].psfQFperf     = NAN;
     168    out[i].stargal       = NAN;
    164169  }
    165170  return (out);
    … …  
    177182    out[i].R       = in[i].R;     
    178183    out[i].D       = in[i].D;     
    179     out[i].Xp      = in[i].Xp;     
    180184
    181185    // changed for PANSTARRS_DEV_0 (moved from Average to Measure)

  • trunk/Ohana/src/libdvo/src/dvo_convert_panstarrs_DEV_0.c

    r34260 r34405  
    3636
    3737    // changed or added for PS1_DEV_1 (2008.02.26)
    38     out[i].psfQual    = (in[i].psfQual == NAN_S_SHORT) ? NAN : in[i].psfQual;
     38    out[i].psfQF      = (in[i].psfQF == NAN_S_SHORT) ? NAN : in[i].psfQF;
    3939    out[i].dbFlags    = in[i].flags;
    4040    out[i].detID      = in[i].detID_lo;
    … …  
    6565    out[i].Mkron      = NAN;
    6666    out[i].dMkron     = NAN;
    67     out[i].FluxPSF     = NAN;
    68     out[i].dFluxPSF    = NAN;
    69     out[i].FluxKron    = NAN;
    70     out[i].dFluxKron   = NAN;
     67    out[i].FluxPSF    = NAN;
     68    out[i].dFluxPSF   = NAN;
     69    out[i].FluxKron   = NAN;
     70    out[i].dFluxKron  = NAN;
     71    out[i].psfQFperf  = NAN;
    7172  }
    7273  return (out);
    … …  
    105106    // changed or added for PS1_DEV_1 (2008.02.26)
    106107    out[i].flags      = in[i].dbFlags;
    107     out[i].psfQual    = isnan(in[i].psfQual) ? NAN_S_SHORT : in[i].psfQual;
     108    out[i].psfQF      = isnan(in[i].psfQF) ? NAN_S_SHORT : in[i].psfQF;
    108109    out[i].detID_hi   = 0;
    109110    out[i].detID_lo   = in[i].detID;
    … …  
    133134    out[i].R             = in[i].R;     
    134135    out[i].D             = in[i].D;     
    135     out[i].Xp            = in[i].Xp;     
    136136    out[i].dR            = in[i].dR;
    137137    out[i].dD            = in[i].dD;
    … …  
    166166    out[i].photFlagsUpper = 0;
    167167    out[i].photFlagsLower = 0;
     168
     169    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     170    out[i].psfQF         = NAN;
     171    out[i].psfQFperf     = NAN;
     172    out[i].stargal       = NAN;
    168173  }
    169174  return (out);
    … …  
    182187    out[i].R        = in[i].R;     
    183188    out[i].D        = in[i].D;     
    184     out[i].Xp       = in[i].Xp;     
    185189    out[i].dR       = in[i].dR;
    186190    out[i].dD       = in[i].dD;

  • trunk/Ohana/src/libdvo/src/dvo_convert_panstarrs_DEV_1.c

    r34260 r34405  
    3636
    3737    // changed or added for PS1_DEV_1 (2008.02.26)
    38     out[i].psfQual    = (in[i].psfQual == NAN_S_SHORT) ? NAN : in[i].psfQual;
     38    out[i].psfQF      = (in[i].psfQF == NAN_S_SHORT) ? NAN : in[i].psfQF;
    3939    out[i].dbFlags    = in[i].flags;
    4040    out[i].detID      = in[i].detID_lo;
    … …  
    6565    out[i].Mkron      = NAN;
    6666    out[i].dMkron     = NAN;
    67     out[i].FluxPSF     = NAN;
    68     out[i].dFluxPSF    = NAN;
    69     out[i].FluxKron    = NAN;
    70     out[i].dFluxKron   = NAN;
     67    out[i].FluxPSF    = NAN;
     68    out[i].dFluxPSF   = NAN;
     69    out[i].FluxKron   = NAN;
     70    out[i].dFluxKron  = NAN;
     71    out[i].psfQFperf  = NAN;
    7172  }
    7273  return (out);
    … …  
    105106    // changed or added for PS1_DEV_1 (2008.02.26)
    106107    out[i].flags      = in[i].dbFlags;
    107     out[i].psfQual    = isnan(in[i].psfQual) ? NAN_S_SHORT : in[i].psfQual;
     108    out[i].psfQF      = isnan(in[i].psfQF) ? NAN_S_SHORT : in[i].psfQF;
    108109    out[i].detID_hi   = 0;
    109110    out[i].detID_lo   = in[i].detID;
    … …  
    133134    out[i].R             = in[i].R;     
    134135    out[i].D             = in[i].D;     
    135     out[i].Xp            = in[i].Xp;     
    136136    out[i].dR            = in[i].dR;
    137137    out[i].dD            = in[i].dD;
    … …  
    166166    out[i].photFlagsUpper = 0;
    167167    out[i].photFlagsLower = 0;
     168
     169    // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway)
     170    out[i].psfQF         = NAN;
     171    out[i].psfQFperf     = NAN;
     172    out[i].stargal       = NAN;
    168173  }
    169174  return (out);
    … …  
    182187    out[i].R        = in[i].R;     
    183188    out[i].D        = in[i].D;     
    184     out[i].Xp       = in[i].Xp;     
    185189    out[i].dR       = in[i].dR;
    186190    out[i].dD       = in[i].dD;

  • trunk/Ohana/src/libdvo/src/dvo_image.c

    r34260 r34405  
    200200  gfits_modify (&db[0].header, "NIMAGES", "%d", 1, 0);
    201201  gfits_modify (&db[0].header, "ZERO_PT", "%lf", 1, ZeroPoint);
     202 
     203  dvo_image_createID (&db[0].header);
    202204
    203205  if (db[0].format == DVO_FORMAT_INTERNAL)        gfits_modify (&db[0].header, "FORMAT", "%s", 1, "INTERNAL");
    … …  
    217219  return;
    218220}
     221
     222// given a dvo image db, add an ID to the header
     223int dvo_image_createID (Header *header) {
     224 
     225  char dbID[33];
     226
     227  if (!header->buffer) return FALSE;
     228
     229  int status = gfits_scan (header, "DVO_DBID", "%s", 1, dbID);     
     230  if (status) {
     231    // do not add a DVO_DBID to a table which already has one
     232    return TRUE;
     233  }
     234
     235  // I have a header, I want to add ONE line.  double check there is enough space
     236  char *p = gfits_header_field (header, "END", 1);
     237  if (p == NULL) {
     238    fprintf (stderr, "header is missing END\n");
     239    return FALSE;
     240  }
     241  if (p - header->buffer + 80 == header->datasize) {
     242    fprintf (stderr, "no free space in block, can't insert keyword\n");
     243    return FALSE;
     244  }
     245
     246  int start_size = header->datasize;
     247
     248  long A = time(NULL);
     249  srand48(A);
     250 
     251  int i;
     252  for (i = 0; i < 32; i++) {
     253    sprintf (&dbID[i], "%1x", (int)(16.0*drand48()));
     254  }
     255
     256  gfits_modify (header, "DVO_DBID", "%s", 1, dbID);     
     257  if (start_size != header->datasize) {
     258    fprintf (stderr, "error: header buffer grew? (%d to %d bytes)\n", (int) start_size, (int) header->datasize);
     259    return FALSE;
     260  }
     261  return TRUE;
     262}
     263
Note: See TracChangeset for help on using the changeset viewer.