IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Mar 4, 2007, 3:45:18 PM (19 years ago)
Author:
eugene
Message:

finishing up the basic code; almost compiles

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/dvo-mods-2007-02/Ohana/src/relastro/src/bcatalog.c

    r12048 r12220  
    1919    if (catalog[0].average[i].Nm < 2) continue;
    2020
    21     /* XXX this limitation is absurd
    22     if (catalog[0].average[i].R < fullregion[0].RA[0]) continue;
    23     if (catalog[0].average[i].R > fullregion[0].RA[1]) continue;
    24     if (catalog[0].average[i].D < fullregion[0].DEC[0]) continue;
    25     if (catalog[0].average[i].D > fullregion[0].DEC[1]) continue;
    26     */
    27 
    2821    /* start with all stars good */
    2922    subcatalog[0].average[Naverage] = catalog[0].average[i];
     
    3427
    3528    if (RESET) {
    36       float *p;
    37       p = (PhotPrimary) ? &subcatalog[0].average[Naverage].M : &subcatalog[0].secfilt[PhotNsec*Naverage+PhotSec].M;
    38       *p = NO_MAG;
    39       p = (PhotPrimary) ? &subcatalog[0].average[Naverage].dM : &subcatalog[0].secfilt[PhotNsec*Naverage+PhotSec].dM;
    40       *p = NO_MAG;
     29      // XXX reset the ra,dec coords?
    4130      subcatalog[0].average[Naverage].code &= ~ID_STAR_FEW;
    4231      subcatalog[0].average[Naverage].code &= ~ID_STAR_POOR;
     
    5948
    6049      /* select measurements by quality */
    61       if (catalog[0].measure[offset].dophot != 1) continue;
     50      if (DophotSelect && (catalog[0].measure[offset].dophot != DophotValue)) continue;
    6251
    6352      /* select measurements by mag limit */
     
    6655
    6756      /* select measurements by measurement error */
    68       if (catalog[0].measure[offset].dM > SIGMA_LIM) continue;
     57      if ((SIGMA_LIM > 0) && (catalog[0].measure[offset].dM > SIGMA_LIM)) continue;
    6958
    7059      /* select measurements by mag limit */
     
    111100  subcatalog[0].Nsecfilt = catalog[0].Nsecfilt;
    112101
     102  // XXX really do this or not???
    113103  if (catalog[0].Naverage > 0) free (catalog[0].average);
    114104  if (catalog[0].Nmeasure > 0) free (catalog[0].measure);
Note: See TracChangeset for help on using the changeset viewer.