Changeset 4949 for trunk/psphot/src/psphotOutput.c
- Timestamp:
- Sep 5, 2005, 8:44:00 PM (21 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotOutput.c (modified) (18 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotOutput.c
r4946 r4949 52 52 53 53 sprintf (name, "%s.psf.dat", filename); 54 p sModelWritePSFs (sources, name);54 pmModelWritePSFs (sources, name); 55 55 56 56 sprintf (name, "%s.flt.dat", filename); 57 p sModelWriteFLTs (sources, name);57 pmModelWriteFLTs (sources, name); 58 58 59 59 sprintf (name, "%s.nul.dat", filename); 60 p sModelWriteNULLs (sources, name);60 pmModelWriteNULLs (sources, name); 61 61 62 62 sprintf (name, "%s.mnt.dat", filename); 63 p sMomentsWriteText (sources, name);63 pmMomentsWriteText (sources, name); 64 64 65 65 psFree (name); … … 71 71 72 72 int i, type, status; 73 p sModel *model;73 pmModel *model; 74 74 psF32 *PAR, *dPAR; 75 75 float sky, dmag, apMag, fitMag; … … 85 85 // write sources with models first 86 86 for (i = 0; i < sources->n; i++) { 87 p sSource *source = (psSource *) sources->data[i];87 pmSource *source = (pmSource *) sources->data[i]; 88 88 89 89 // use the correct model … … 162 162 163 163 int i, status, flags; 164 p sModel *model;164 pmModel *model; 165 165 psF32 *PAR, *dPAR; 166 166 float sky, dmag, apMag, fitMag; … … 176 176 // write sources with models first 177 177 for (i = 0; i < sources->n; i++) { 178 p sSource *source = (psSource *) sources->data[i];178 pmSource *source = (pmSource *) sources->data[i]; 179 179 flags = 0; 180 180 … … 260 260 261 261 int i, type, status; 262 p sModel *model;262 pmModel *model; 263 263 psF32 *PAR, *dPAR; 264 264 float sky, dmag, apMag, fitMag, lsky, Theta; … … 282 282 // write sources with models first 283 283 for (i = 0; i < sources->n; i++) { 284 p sSource *source = (psSource *) sources->data[i];284 pmSource *source = (pmSource *) sources->data[i]; 285 285 286 286 // use the correct model … … 395 395 396 396 // write the moments to an output file 397 bool p sMomentsWriteText (psArray *sources, char *filename) {397 bool pmMomentsWriteText (psArray *sources, char *filename) { 398 398 399 399 int i; 400 400 FILE *f; 401 p sSource *source;401 pmSource *source; 402 402 403 403 f = fopen (filename, "w"); 404 404 if (f == NULL) { 405 psLogMsg ("p sMomentsWriteText", 3, "can't open output file for moments%s\n", filename);406 return false; 407 } 408 409 for (i = 0; i < sources->n; i++) { 410 source = (p sSource *) sources->data[i];405 psLogMsg ("pmMomentsWriteText", 3, "can't open output file for moments%s\n", filename); 406 return false; 407 } 408 409 for (i = 0; i < sources->n; i++) { 410 source = (pmSource *) sources->data[i]; 411 411 if (source->moments == NULL) continue; 412 412 fprintf (f, "%5d %5d %7.1f %7.1f %7.1f %6.3f %6.3f %8.1f %7.1f %7.1f %7.1f %4d %2d\n", … … 423 423 424 424 // dump the sources to an output file 425 bool p sModelWritePSFs (psArray *sources, char *filename) {425 bool pmModelWritePSFs (psArray *sources, char *filename) { 426 426 427 427 double dP, flux; … … 430 430 psVector *params; 431 431 psVector *dparams; 432 p sModel *model;432 pmModel *model; 433 433 434 434 f = fopen (filename, "w"); 435 435 if (f == NULL) { 436 psLogMsg ("p sModelWritePSFs", 3, "can't open output file for moments%s\n", filename);436 psLogMsg ("pmModelWritePSFs", 3, "can't open output file for moments%s\n", filename); 437 437 return false; 438 438 } … … 440 440 // write sources with models first 441 441 for (i = 0; i < sources->n; i++) { 442 p sSource *source = (psSource *) sources->data[i];443 model = (p sModel *) source->modelPSF;442 pmSource *source = (pmSource *) sources->data[i]; 443 model = (pmModel *) source->modelPSF; 444 444 if (model == NULL) continue; 445 445 … … 459 459 dP = sqrt (dP); 460 460 461 p sModelFlux modelFluxFunc = psModelFlux_GetFunction (model->type);461 pmModelFlux modelFluxFunc = pmModelFlux_GetFunction (model->type); 462 462 flux = modelFluxFunc (params); 463 463 … … 489 489 490 490 // dump the sources to an output file 491 bool p sModelWriteFLTs (psArray *sources, char *filename) {491 bool pmModelWriteFLTs (psArray *sources, char *filename) { 492 492 493 493 double dP; … … 496 496 psVector *params; 497 497 psVector *dparams; 498 p sModel *model;498 pmModel *model; 499 499 500 500 f = fopen (filename, "w"); 501 501 if (f == NULL) { 502 psLogMsg ("p sModelWriteFLTs", 3, "can't open output file for moments%s\n", filename);502 psLogMsg ("pmModelWriteFLTs", 3, "can't open output file for moments%s\n", filename); 503 503 return false; 504 504 } … … 506 506 // write sources with models first 507 507 for (i = 0; i < sources->n; i++) { 508 p sSource *source = (psSource *) sources->data[i];509 model = (p sModel *) source->modelFLT;508 pmSource *source = (pmSource *) sources->data[i]; 509 model = (pmModel *) source->modelFLT; 510 510 if (model == NULL) continue; 511 511 if (source->type == PS_SOURCE_GALAXY) goto valid; … … 551 551 552 552 // dump the sources to an output file 553 bool p sModelWriteNULLs (psArray *sources, char *filename)553 bool pmModelWriteNULLs (psArray *sources, char *filename) 554 554 { 555 555 … … 563 563 } 564 564 565 p sMoments *empty = pmMomentsAlloc ();565 pmMoments *empty = pmMomentsAlloc (); 566 566 567 567 // write sources with models first 568 568 for (i = 0; i < sources->n; i++) { 569 p sSource *source = (psSource *) sources->data[i];569 pmSource *source = (pmSource *) sources->data[i]; 570 570 571 571 // skip these sources (in PSF or FLT) … … 598 598 599 599 // translations between psphot object types and dophot object types 600 int pmSourceDophotType (p sSource *source) {600 int pmSourceDophotType (pmSource *source) { 601 601 602 602 switch (source->type) {
Note:
See TracChangeset
for help on using the changeset viewer.
