Index: source/prt_header.cpp =================================================================== --- source/prt_header.cpp (revision 9455) +++ source/prt_header.cpp (working copy) @@ -243,13 +243,19 @@ /* magnitudes */ if( radius.lgRadiusKnown ) { - solar = log10(continuum.TotalLumin) + radius.pirsq - 33.5828; - absbol = 4.75 - 2.5*solar; + solar = absbol= -38.; + if( continuum.TotalLumin> SMALLFLOAT ) + { + solar = log10(continuum.TotalLumin) + radius.pirsq - 33.5828; + absbol = 4.75 - 2.5*solar; + } /* absv = 4.79 - 2.5 * (LOG10(MAX(1e-30,continuum.fluxv)) + pirsq - 18.742 - * 1 0.016) * allen 76, page 197 */ - absv = -2.5*(log10(MAX2(1e-30,continuum.fluxv)) + radius.pirsq - 20.654202); + absv = -38.; + if( continuum.fluxv>SMALLFLOAT ) + absv = -2.5*(log10(MAX2(1e-30,continuum.fluxv)) + radius.pirsq - 20.654202); /* >>chng 97 mar 18, following branch so zero returned when no radiation at all */ if( continuum.fbeta > 0. ) Index: source =================================================================== --- source (revision 9455) +++ source (working copy) Property changes on: source ___________________________________________________________________ Modified: svn:mergeinfo Merged /trunk/source:r9455 Index: . =================================================================== --- . (revision 9455) +++ . (working copy) Property changes on: . ___________________________________________________________________ Modified: svn:mergeinfo Merged /trunk:r9455