- Timestamp:
- Mar 4, 2007, 3:45:18 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dvo-mods-2007-02/Ohana/src/relastro/src/bcatalog.c
r12048 r12220 19 19 if (catalog[0].average[i].Nm < 2) continue; 20 20 21 /* XXX this limitation is absurd22 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 28 21 /* start with all stars good */ 29 22 subcatalog[0].average[Naverage] = catalog[0].average[i]; … … 34 27 35 28 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? 41 30 subcatalog[0].average[Naverage].code &= ~ID_STAR_FEW; 42 31 subcatalog[0].average[Naverage].code &= ~ID_STAR_POOR; … … 59 48 60 49 /* select measurements by quality */ 61 if ( catalog[0].measure[offset].dophot != 1) continue;50 if (DophotSelect && (catalog[0].measure[offset].dophot != DophotValue)) continue; 62 51 63 52 /* select measurements by mag limit */ … … 66 55 67 56 /* 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; 69 58 70 59 /* select measurements by mag limit */ … … 111 100 subcatalog[0].Nsecfilt = catalog[0].Nsecfilt; 112 101 102 // XXX really do this or not??? 113 103 if (catalog[0].Naverage > 0) free (catalog[0].average); 114 104 if (catalog[0].Nmeasure > 0) free (catalog[0].measure);
Note:
See TracChangeset
for help on using the changeset viewer.
