Changeset 1198

Show
Ignore:
Timestamp:
06/16/07 20:38:18 (1 year ago)
Author:
gary
Message:

bugfix merge of 1197 - did not print incident continuum correctly in both emergent and intrinsic continuum

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • branches/c07.02_branch/source/prt_lines_continuum.cpp

    r1130 r1198  
    226226        continuum.cn4861 /= (float)radius.dVeff; 
    227227        linadd(continuum.cn4861,4860,"Inci",'i'); 
     228        /* >>chng 07 jun 13, at this point nsum is the number of lines in the stack 
     229         * so nsum-1 is the entry created by the above call to linadd 
     230         * we want to set the emergent incident continuum to the intrinsic 
     231         * incident continuum - there are reports of the incident continuum  
     232         * striking the cloud and the emergent - incident distinction does 
     233         * not apply */ 
     234        if( LineSave.ipass > 0 ) 
     235        { 
     236                LineSv[LineSave.nsum-1].sumlin[1] = LineSv[LineSave.nsum-1].sumlin[0]; 
     237        } 
    228238 
    229239        /* incident continuum nu*f_nu near Ly-alpha, at illuminated face of cloud */ 
    230240        continuum.cn1216 /= (float)radius.dVeff; 
    231241        linadd(continuum.cn1216,1215,"Inci",'i'); 
     242        /* see comment concerning parallel code immediately above */ 
     243        if( LineSave.ipass > 0 ) 
     244        { 
     245                LineSv[LineSave.nsum-1].sumlin[1] = LineSv[LineSave.nsum-1].sumlin[0]; 
     246        } 
    232247 
    233248        continuum.cn1216 *= (float)radius.dVeff;