User: cott Date: 2010/12/04 06:44 PM

Modified: /TOVSolverHot/src/ tov.c

Log: * fix stupid bug introduced by myself in mapping 1D->3D

File Changes:

Directory: /TOVSolverHot/src/ =============================

File [modified]: tov.c Delta lines: +12 -7 =================================================================== --- TOVSolverHot/src/tov.c 2010-12-05 00:09:30 UTC (rev 27) +++ TOVSolverHot/src/tov.c 2010-12-05 00:44:50 UTC (rev 28) @@ -869,10 +869,15 @@ CCTK_WARN(0,"aborting"); }

- temperature_point[star] = xtemp; - entropy_point[star] = xent; - ye_point[star] = xye; - + if(rho_point[star] > nuc_eos_rho_min) { + temperature_point[star] = xtemp; + entropy_point[star] = xent; + ye_point[star] = xye; + } else { + temperature_point[star] = nuc_eos_temperature; + entropy_point[star] = nuc_eos_entropy; + ye_point[star] = nuc_eos_ye; + } } else if(!use_EOS_Omni) { @@ -1056,9 +1061,9 @@ eps[i3D] += eps_point[star_i]; press[i3D] += press_point[star_i]; if(eoskey==4) { - temperature[i3D] += temperature_point[star]; - entropy[i3D] += entropy_point[star]; - Y_e[i3D] += ye_point[star]; + temperature[i3D] += temperature_point[star_i]; + entropy[i3D] += entropy_point[star_i]; + Y_e[i3D] += ye_point[star_i]; }

/* we still have to know if we are inside one star - and which */