Changeset 1197

Timestamp:

06/16/07 20:35:27 (1 year ago)

Author:

gary

Message:

prt_lines_continuum.cpp - bugfix - incident continuum near Lya and Hbeta were zero in the emergent continuum. there is actually no distinction between emergent and intrinsic incident continua. this fix sets both to the correct incident value

Files:

• trunk/source/prt_lines_continuum.cpp (modified) (1 diff)

trunk/source/prt_lines_continuum.cpp

@@ -229 +229 @@
-        linadd(continuum.cn4861,4860,"Inci",'i',
-                "incident continuum nu*f_nu at H-beta, at illuminated face of cloud ");
+        /* >>chng 07 jun 13, at this point nsum is the number of lines in the stack
+         * so nsum-1 is the entry created by the above call to linadd
+         * we want to set the emergent incident continuum to the intrinsic
+         * incident continuum - there are reports of the incident continuum 
+         * striking the cloud and the emergent - incident distinction does
+         * not apply */
+        if( LineSave.ipass > 0 )
+        {
+                LineSv[LineSave.nsum-1].sumlin[1] = LineSv[LineSave.nsum-1].sumlin[0];
+        }
-
-        continuum.cn1216 /= (float)radius.dVeff;
-        linadd(continuum.cn1216,1215,"Inci",'i',
-                "incident continuum nu*f_nu near Ly-alpha, at illuminated face of cloud");
+        /* see comment concerning parallel code immediately above */
+        if( LineSave.ipass > 0 )
+        {
+                LineSv[LineSave.nsum-1].sumlin[1] = LineSv[LineSave.nsum-1].sumlin[0];
+        }
-
-        continuum.cn1216 *= (float)radius.dVeff;