[Commits] [svn:einsteintoolkit] incoming/TOVSolverHot/src/ (Rev. 28)

[cott at tapir.caltech.edu]

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 */