Changeset 20190
- Timestamp:
- Oct 16, 2008, 9:18:27 AM (18 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/relphot/src/load_catalogs.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/relphot/src/load_catalogs.c
r15743 r20190 3 3 Catalog *load_catalogs (SkyList *skylist, int *Ncatalog) { 4 4 5 int i, N star;5 int i, Nmeas, Nstar, Nmeas_total, Nstar_total; 6 6 Catalog *catalog, tcatalog; 7 7 … … 9 9 10 10 ALLOCATE (catalog, Catalog, skylist[0].Nregions); 11 12 Nmeas_total = Nstar_total = 0; 11 13 12 14 // load data from each region file, only use bright stars … … 25 27 } 26 28 if (VERBOSE && !tcatalog.Naves_disk) fprintf (stderr, "no data in %s, skipping\n", tcatalog.filename); 29 Nstar_total += tcatalog.Naverage; 30 Nmeas_total += tcatalog.Nmeasure; 27 31 28 32 // select only the brighter stars … … 33 37 34 38 // XXX keep this test? 35 Nstar = 0;39 Nstar = Nmeas = 0; 36 40 for (i = 0; i < skylist[0].Nregions; i++) { 37 41 Nstar += catalog[i].Naverage; 42 Nmeas += catalog[i].Nmeasure; 38 43 } 39 44 if (Nstar < 2) { 40 45 fprintf (stderr, "warning: insufficient stars %d\n", Nstar); 41 46 } 47 48 fprintf (stderr, "using %d of %d stars (%d of %d measurements)\n", Nstar, Nstar_total, Nmeas, Nmeas_total); 42 49 43 50 // XXX consider only returning the populated catalogs
Note:
See TracChangeset
for help on using the changeset viewer.
