Changeset 15208 for trunk/Ohana/src/opihi/dvo/gstar.c
- Timestamp:
- Oct 3, 2007, 4:55:24 PM (19 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/opihi/dvo/gstar.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/opihi/dvo/gstar.c
r15038 r15208 11 11 double Mcat, Mrel; 12 12 float *RA, *DEC; 13 int i, j, k, m, N, *N1, Nsecfilt, NPTS, QUIET ;13 int i, j, k, m, N, *N1, Nsecfilt, NPTS, QUIET, FULL_OUTPUT; 14 14 int Nstars, found, GetMeasures, Nlo, Nhi; 15 15 int SaveVectors; … … 42 42 } 43 43 44 FULL_OUTPUT = FALSE; 45 if ((N = get_argument (argc, argv, "-full"))) { 46 FULL_OUTPUT = TRUE; 47 remove_argument (N, &argc, argv); 48 } 49 44 50 GetMeasures = FALSE; 45 51 if ((N = get_argument (argc, argv, "-m"))) { 52 GetMeasures = TRUE; 53 remove_argument (N, &argc, argv); 54 } 55 if ((N = get_argument (argc, argv, "-meas"))) { 46 56 GetMeasures = TRUE; 47 57 remove_argument (N, &argc, argv); … … 139 149 if (!QUIET) { 140 150 gprint (GP_LOG, "star: %d\n", k); 141 gprint (GP_LOG, "%11.7f %11.7f %3d of %3d %4.1f %5d\n", catalog.average[k].R, catalog.average[k].D, 142 catalog.average[k].Nm, catalog.average[k].Nn + catalog.average[k].Nm, 143 0.01*catalog.average[k].Xp, catalog.average[k].code); 151 gprint (GP_LOG, "%11.7f ", catalog.average[k].R); 152 gprint (GP_LOG, "%11.7f ", catalog.average[k].D); 153 gprint (GP_LOG, "%3d ", catalog.average[k].Nm); 154 gprint (GP_LOG, "%4.1f ", 0.01*catalog.average[k].Xp); 155 gprint (GP_LOG, "%5d", catalog.average[k].code); 156 157 if (FULL_OUTPUT) { 158 gprint (GP_LOG, "%f", catalog.average[k].dR); 159 gprint (GP_LOG, "%f", catalog.average[k].dD); 160 gprint (GP_LOG, "%f", catalog.average[k].uR); 161 gprint (GP_LOG, "%f", catalog.average[k].uD); 162 gprint (GP_LOG, "%f", catalog.average[k].duR); 163 gprint (GP_LOG, "%f", catalog.average[k].duD); 164 gprint (GP_LOG, "%f", catalog.average[k].P); 165 gprint (GP_LOG, "%f", catalog.average[k].dP); 166 gprint (GP_LOG, "%x", catalog.average[k].objID_hi); 167 gprint (GP_LOG, "%x", catalog.average[k].objID_lo); 168 } 169 170 gprint (GP_LOG, "\n"); 144 171 145 172 /* filter names */ … … 169 196 if (GetMeasures && !QUIET) { 170 197 date = ohana_sec_to_date (catalog.measure[m].t); 171 gprint (GP_LOG, "%6.3f %6.3f %5.3f %20s %5.2f %5.2f %2d %3x %3d %-20s\n", 172 Mcat, Mrel, catalog.measure[m].dM, 173 date, catalog.measure[m].dR, catalog.measure[m].dD, 174 catalog.measure[m].dophot, catalog.measure[m].dbFlags, 175 catalog.measure[m].photcode, GetPhotcodeNamebyCode (catalog.measure[m].photcode)); 198 gprint (GP_LOG, "%6.3f ", Mcat); 199 gprint (GP_LOG, "%6.3f ", Mrel); 200 gprint (GP_LOG, "%5.3f ", catalog.measure[m].dM); 201 gprint (GP_LOG, "%20s ", date); 202 gprint (GP_LOG, "%7.4f ", catalog.measure[m].dR); 203 gprint (GP_LOG, "%7.4f", catalog.measure[m].dD); 204 gprint (GP_LOG, "%2d ", catalog.measure[m].dophot); 205 gprint (GP_LOG, "%3x ", catalog.measure[m].dbFlags); 206 gprint (GP_LOG, "%3d ", catalog.measure[m].photcode); 207 gprint (GP_LOG, "%-20s ", GetPhotcodeNamebyCode (catalog.measure[m].photcode)); 208 209 if (FULL_OUTPUT) { 210 gprint (GP_LOG, "%f", catalog.measure[m].Mcal); 211 gprint (GP_LOG, "%f", catalog.measure[m].Mgal); 212 gprint (GP_LOG, "%f", pow(10.0, 0.4*catalog.measure[m].dt)); 213 gprint (GP_LOG, "%f", 1.0 + catalog.measure[m].airmass); 214 gprint (GP_LOG, "%f", catalog.measure[m].az); 215 gprint (GP_LOG, "%f", catalog.measure[m].Xccd); 216 gprint (GP_LOG, "%f", catalog.measure[m].Yccd); 217 gprint (GP_LOG, "%d", catalog.measure[m].dXccd); 218 gprint (GP_LOG, "%d", catalog.measure[m].dYccd); 219 gprint (GP_LOG, "%f", catalog.measure[m].Sky); 220 gprint (GP_LOG, "%f", catalog.measure[m].dSky); 221 gprint (GP_LOG, "%d", catalog.measure[m].averef); 222 gprint (GP_LOG, "%d", catalog.measure[m].detID); 223 gprint (GP_LOG, "%d", catalog.measure[m].imageID); 224 gprint (GP_LOG, "%f", catalog.measure[m].qPSF); 225 gprint (GP_LOG, "%f", catalog.measure[m].psfProb); 226 gprint (GP_LOG, "%f", catalog.measure[m].crNsigma); 227 gprint (GP_LOG, "%f", catalog.measure[m].extNsigma); 228 gprint (GP_LOG, "%f", 0.01*catalog.measure[m].FWx); 229 gprint (GP_LOG, "%f", 0.01*catalog.measure[m].FWy); 230 gprint (GP_LOG, "%f", (360.0/(float)0xffff)*catalog.measure[m].theta); 231 232 gprint (GP_LOG, "%x", catalog.measure[m].photFlags); 233 gprint (GP_LOG, "%d", catalog.measure[m].stargal); 234 } 235 gprint (GP_LOG, "\n"); 236 176 237 free (date); 177 238 }
Note:
See TracChangeset
for help on using the changeset viewer.
