[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