Changeset 1838 for branches/newmole/source/mole_co_atom.cpp

Timestamp:

03/09/08 16:59:05 (10 months ago)

Author:

rjrw

Message:

Minor changes to resulting from trying to chase down differences on
pdf_leiden_f1 --

ion_solver: allow molecule effects in first iteration,
iso_level: defensive programming for z[][]
mole_eval_balance: move some definitions for sp
mole_co_atom: add SDIV (aka fixit...) to prevent FPE in some circumstances

Files:

• branches/newmole/source/mole_co_atom.cpp (modified) (1 diff)

branches/newmole/source/mole_co_atom.cpp

@@ -333 +333 @@
                         (*Rotate)[j].Emis->ColOvTot = 0.0;
                 else
-                        (*Rotate)[j].Emis->ColOvTot = (realnum)(CollRate[j][j+1]/(CollRate[j][j+1]+(*Rotate)[j].Emis->pump) );
+                        (*Rotate)[j].Emis->ColOvTot = (realnum)(CollRate[j][j+1]/SDIV(CollRate[j][j+1]+(*Rotate)[j].Emis->pump) );
 
                 /* two cases - collisionally excited (usual case) or