Changeset 12332 for trunk/Ohana/src/opihi/dvo/catalog.c
- Timestamp:
- Mar 8, 2007, 12:21:18 PM (19 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/opihi/dvo/catalog.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/opihi/dvo/catalog.c
r8427 r12332 197 197 dparse (&D, Ad, tbuffer); 198 198 dparse (&M, Am, tbuffer); 199 catalog.average[nstar].R _PS= R;200 catalog.average[nstar].D _PS= D;201 catalog.average[nstar].M _PS= M;199 catalog.average[nstar].R = R; 200 catalog.average[nstar].D = D; 201 catalog.average[nstar].M = M; 202 202 nstar++; 203 203 if (nstar == NSTARS - 1) { … … 239 239 dparse (&D, Ad, &tbuffer[i*Bytes_Star]); 240 240 dparse (&M, Am, &tbuffer[i*Bytes_Star]); 241 catalog.average[nstar].R _PS= R;242 catalog.average[nstar].D _PS= D;243 catalog.average[nstar].M _PS= M;241 catalog.average[nstar].R = R; 242 catalog.average[nstar].D = D; 243 catalog.average[nstar].M = M; 244 244 nstar++; 245 245 if (nstar == NSTARS - 1) { … … 289 289 dparse (&D, 2, &tbuffer[i*BYTES_STAR]); 290 290 dparse (&M, 3, &tbuffer[i*BYTES_STAR]); 291 catalog.average[nstar].R _PS= R;292 catalog.average[nstar].D _PS= D;293 catalog.average[nstar].M _PS= M;291 catalog.average[nstar].R = R; 292 catalog.average[nstar].D = D; 293 catalog.average[nstar].M = M; 294 294 nstar++; 295 295 if (nstar == NSTARS - 1) { … … 342 342 if (IDclip && (catalog.average[i].code != IDchoice)) 343 343 continue; 344 Zvec.elements[N] = MIN (1.0, MAX (0.01, (catalog.average[i].M _PS- Mz) / Mr));344 Zvec.elements[N] = MIN (1.0, MAX (0.01, (catalog.average[i].M - Mz) / Mr)); 345 345 if (LimExclude && (Zvec.elements[N] > 0.99)) continue; 346 346 if (Zvec.elements[N] < 0.011) continue; 347 while (catalog.average[i].R _PS < Rmin) catalog.average[i].R_PS+= 360.0;348 while (catalog.average[i].R _PS > Rmax) catalog.average[i].R_PS-= 360.0;349 if (fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R _PS, catalog.average[i].D_PS, &graphmode.coords)) N ++;347 while (catalog.average[i].R < Rmin) catalog.average[i].R += 360.0; 348 while (catalog.average[i].R > Rmax) catalog.average[i].R -= 360.0; 349 if (fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R, catalog.average[i].D, &graphmode.coords)) N ++; 350 350 } 351 351 break; 352 352 case (NUMSCALE): 353 353 for (N = i = 0; i < catalog.Naverage; i++) { 354 if (clip && ((catalog.average[i].M _PS > Mz) || (catalog.average[i].M_PS< Mr+Mz)))354 if (clip && ((catalog.average[i].M > Mz) || (catalog.average[i].M < Mr+Mz))) 355 355 continue; 356 356 if (IDclip && (catalog.average[i].code != IDchoice)) … … 360 360 if (Zvec.elements[N] == 0.01) 361 361 continue; 362 while (catalog.average[i].R _PS < Rmin) catalog.average[i].R_PS+= 360.0;363 while (catalog.average[i].R _PS > Rmax) catalog.average[i].R_PS-= 360.0;364 if (fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R _PS, catalog.average[i].D_PS, &graphmode.coords)) N++;362 while (catalog.average[i].R < Rmin) catalog.average[i].R += 360.0; 363 while (catalog.average[i].R > Rmax) catalog.average[i].R -= 360.0; 364 if (fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R, catalog.average[i].D, &graphmode.coords)) N++; 365 365 } 366 366 break; 367 367 case (MISSCALE): 368 368 for (N = i = 0; i < catalog.Naverage; i++) { 369 if (clip && ((catalog.average[i].M _PS > Mz) || (catalog.average[i].M_PS< Mr+Mz)))369 if (clip && ((catalog.average[i].M > Mz) || (catalog.average[i].M < Mr+Mz))) 370 370 continue; 371 371 if (IDclip && (catalog.average[i].code != IDchoice)) … … 375 375 if (Zvec.elements[N] == 0.01) 376 376 continue; 377 while (catalog.average[i].R _PS < Rmin) catalog.average[i].R_PS+= 360.0;378 while (catalog.average[i].R _PS > Rmax) catalog.average[i].R_PS-= 360.0;379 if (fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R _PS, catalog.average[i].D_PS, &graphmode.coords)) N++;377 while (catalog.average[i].R < Rmin) catalog.average[i].R += 360.0; 378 while (catalog.average[i].R > Rmax) catalog.average[i].R -= 360.0; 379 if (fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R, catalog.average[i].D, &graphmode.coords)) N++; 380 380 } 381 381 break;
Note:
See TracChangeset
for help on using the changeset viewer.
