- Timestamp:
- Mar 7, 2007, 11:54:39 AM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dvo-mods-2007-02/Ohana/src/opihi/dvo/calextract.c
r12046 r12295 6 6 7 7 int i, j, m, N, Nr, mode[2]; 8 int Nsec , NSTAR;8 int Nsecfilt, NSTAR; 9 9 char *RegionName, *RegionList; 10 10 double M1, M2, dM2, color; … … 26 26 /* load photcode information */ 27 27 if (!InitPhotcodes ()) return (FALSE); 28 Nsec = GetPhotcodeNsecfilt ();28 Nsecfilt = GetPhotcodeNsecfilt (); 29 29 30 30 /* command line arguments */ … … 78 78 m = catalog.average[i].offset; 79 79 80 if (code[0][0].c1 && code[0][0].c2 && !PhotColor (&catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], code[0][0].c1, code[0][0].c2, &color)) continue;80 if (code[0][0].c1 && code[0][0].c2 && !PhotColor (&catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], code[0][0].c1, code[0][0].c2, &color)) continue; 81 81 82 82 /* find data for filter 2 (PHOT_REF) */ … … 91 91 92 92 /* find data for filter 1 */ 93 M1 = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], AVE_MAG);93 M1 = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], AVE_MAG); 94 94 if (M1 == NO_MAG) continue; 95 95 … … 99 99 vec[NR ][0].elements[N] = catalog.average[i].R; 100 100 vec[ND ][0].elements[N] = catalog.average[i].D; 101 vec[Nd1][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], AVE_dMAG);101 vec[Nd1][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], AVE_dMAG); 102 102 vec[Nd2][0].elements[N] = dM2; 103 vec[Nx ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], AVE_Xm);104 vec[Nc ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], AVE_NCODE);105 vec[Np ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], AVE_NPHOT);106 vec[Nt ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsec ], &catalog.measure[m], AVE_NCRIT);103 vec[Nx ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], AVE_Xm); 104 vec[Nc ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], AVE_NCODE); 105 vec[Np ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], AVE_NPHOT); 106 vec[Nt ][0].elements[N] = ExtractAverages (code[0], mode[0], &catalog.average[i], &catalog.secfilt[i*Nsecfilt], &catalog.measure[m], AVE_NCRIT); 107 107 N ++; 108 108 if (N == NSTAR) {
Note:
See TracChangeset
for help on using the changeset viewer.
