Changeset 2419 for branches/newmole
- Timestamp:
- 10/23/08 16:31:27 (3 months ago)
- Files:
-
- 1 modified
-
branches/newmole/source/dense.cpp (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/newmole/source/dense.cpp
r1851 r2419 47 47 2e-5*dense.gas_phase[nelem] ) 48 48 { 49 fprintf(ioQQQ,"Problem ions %g moles %g error %g of %g nelem %ld(%s)\n",50 sum,dense.xMolecules[nelem],51 sum+dense.xMolecules[nelem]-dense.gas_phase[nelem],52 dense.gas_phase[nelem],nelem,53 elementnames.chElementSym[nelem]);54 49 fprintf(ioQQQ,"DEBUG non-conv nelem\t%li\tsum\t%e\ttot gas\t%e\trel err\t%.3e\n", 55 50 nelem, 56 sum ,51 sum+dense.xMolecules[nelem], 57 52 dense.gas_phase[nelem], 58 (sum -dense.gas_phase[nelem])/dense.gas_phase[nelem]);53 (sum+dense.xMolecules[nelem]-dense.gas_phase[nelem])/dense.gas_phase[nelem]); 59 54 fprintf(ioQQQ,"Conservation should be better than this\n"); 60 55 lgOK=false;
