[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 321)
roland.haas at physics.gatech.edu
roland.haas at physics.gatech.edu
Sun May 13 14:03:31 CDT 2012
User: rhaas
Date: 2012/05/13 02:03 PM
Modified:
/trunk/src/
GRHydro_Con2PrimM_pt.c
Log:
Consistent atmosphere treatment.
patch by Philipp Moesta
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_Con2PrimM_pt.c
Delta lines: +13 -2
===================================================================
--- trunk/src/GRHydro_Con2PrimM_pt.c 2012-05-13 19:01:37 UTC (rev 320)
+++ trunk/src/GRHydro_Con2PrimM_pt.c 2012-05-13 19:03:31 UTC (rev 321)
@@ -45,6 +45,7 @@
#include <complex.h>
#include "cctk.h"
+#include "cctk_Parameters.h"
/* Set this to be 1 if you want debug output */
#define DEBUG_CON2PRIMM (0)
@@ -215,8 +216,10 @@
CCTK_REAL detg = (*det);
CCTK_REAL sqrt_detg = sqrt(detg);
CCTK_REAL inv_sqrt_detg = 1./sqrt_detg;
- CCTK_INT i,j, i_increase ;
-
+ CCTK_INT i,j, i_increase;
+
+ DECLARE_CCTK_PARAMETERS;
+
struct LocGlob lg;
gammaeos = *gamma_eos;
@@ -395,7 +398,15 @@
*vely = g_o_WBsq * ( usy + QdB_o_W*(*By) ) ;
*velz = g_o_WBsq * ( usz + QdB_o_W*(*Bz) ) ;
+ if (*rho <= rho_abs_min*(1.0+GRHydro_atmo_tolerance) ) {
+ *rho = rho_abs_min;
+ *velx = 0.0;
+ *vely = 0.0;
+ *velz = 0.0;
+ *w_lorentz = 1.0;
+ }
+
#if(DEBUG_CON2PRIMM)
fprintf(stdout,"rho = %26.16e \n",*rho );
fprintf(stdout,"epsilon = %26.16e \n",*epsilon );
More information about the Commits
mailing list