Changeset 2096
- Timestamp:
- 05/30/08 07:04:12 (3 months ago)
- Location:
- branches/newmole/source
- Files:
-
- 7 modified
-
conv_base.cpp (modified) (1 diff)
-
molcol.cpp (modified) (1 diff)
-
mole_drive.cpp (modified) (1 diff)
-
mole_reactions.cpp (modified) (1 diff)
-
mole_solve.cpp (modified) (1 diff)
-
mole_species.cpp (modified) (1 diff)
-
radius_increment.cpp (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/newmole/source/conv_base.cpp
r2090 r2096 618 618 iso.xIonSimple[ipISO][nelem]<1e-10) ) 619 619 { 620 const realnum err_tol = 1e-2 ;620 const realnum err_tol = 1e-2f; 621 621 realnum abund = 0.; 622 622 for (long n=0; n<iso.numLevels_local[ipISO][nelem]; ++n) 623 623 { 624 abund += StatesElem[ipISO][nelem][n].Pop;624 abund += (realnum)StatesElem[ipISO][nelem][n].Pop; 625 625 } 626 626 if( conv.lgConvIoniz && iso.pop_ion_ov_neut[ipISO][nelem] > 1e-15 && -
branches/newmole/source/molcol.cpp
r2042 r2096 39 39 for( i=0; i < mole.num_calc; i++ ) 40 40 { 41 mole.list[i]->column += (realnum) mole.list[i]->den*radius.drad_x_fillfac;41 mole.list[i]->column += (realnum)(mole.list[i]->den*radius.drad_x_fillfac); 42 42 } 43 43 } -
branches/newmole/source/mole_drive.cpp
r2055 r2096 133 133 134 134 /* total H2 - all forms */ 135 hmi.H2_total = (realnum) findspecies("H2*")->den + findspecies("H2")->den;135 hmi.H2_total = (realnum)(findspecies("H2*")->den + findspecies("H2")->den); 136 136 /* first guess at ortho and para densities */ 137 137 h2.ortho_density = 0.75*hmi.H2_total; -
branches/newmole/source/mole_reactions.cpp
r2068 r2096 1950 1950 { 1951 1951 rate = p->second; 1952 realnum oldrk = rate->rk;1952 realnum oldrk = (realnum)rate->rk; 1953 1953 ASSERT(rate->fun != NULL); 1954 1954 rate->rk = rate->a*rate->fun(rate); -
branches/newmole/source/mole_solve.cpp
r2076 r2096 279 279 findspecies("H2*")->den = 0.; 280 280 281 hmi.H2_total = findspecies("H2")->den + findspecies("H2*")->den;281 hmi.H2_total = (realnum)(findspecies("H2")->den + findspecies("H2*")->den); 282 282 /* first guess at ortho and para densities */ 283 283 h2.ortho_density = 0.75*hmi.H2_total; -
branches/newmole/source/mole_species.cpp
r2094 r2096 673 673 mole.grain_saturation = 1.0; 674 674 else 675 mole.grain_saturation = (realnum) 4*den_times_area/(mole_cs*adsorbed_density);675 mole.grain_saturation = (realnum)(4.*den_times_area/(mole_cs*adsorbed_density)); 676 676 } 677 677 else -
branches/newmole/source/radius_increment.cpp
r2089 r2096 105 105 * change is to only do this branch if abundance is large enough to be detected by 106 106 * the ionization convergence solvers */ 107 const realnum err_tol = 1e-2 ;107 const realnum err_tol = 1e-2f; 108 108 /* >>chng 05 sep 02, when low-T solver used solution is approximate, 109 109 * and it must not matter (lot-T solver should not be used if it … … 114 114 for (long n=0; n<iso.numLevels_local[ipISO][nelem]; ++n) 115 115 { 116 abund += StatesElem[ipISO][nelem][n].Pop;116 abund += (realnum)StatesElem[ipISO][nelem][n].Pop; 117 117 } 118 118 //fprintf(ioQQQ,"%.2f %ld %ld %11.4g %11.4g\n",fnzone,ipISO,nelem,
