fixed bug in UU floor, still an issue with GRMHD monopole

achael · Mar 3, 2022 · 17bbfdc · 17bbfdc
1 parent 04ebb94
commit 17bbfdc
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/PROBLEMS/MONOPOLE_2D/define.h b/PROBLEMS/MONOPOLE_2D/define.h
@@ -61,7 +61,7 @@
 #define CORRECT_POLARAXIS
 #define NCCORRECTPOLAR 2
 
-#define B2RHOFLOORFRAME DRIFTFRAME // ZAMOFRAME 
+#define B2RHOFLOORFRAME ZAMOFRAME 
 #define UURHORATIOMIN 0.
 #define UURHORATIOMAX 50.
 #define B2UURATIOMIN 0.

diff --git a/u2p.c b/u2p.c
@@ -414,7 +414,7 @@ check_floors_mhd(ldouble *pp, int whichvel,void *ggg)
     }
 
     if(verbose) {printf("hd_floors CASE 2 at (%d,%d,%d | %d,%d,%d): %e %e | tijk: %d %d %d\n",geom->ix+TOI,geom->iy+TOJ,geom->iz+TOK,geom->ix,geom->iy,geom->iz,pp[RHO],pp[UU],TI,TJ,TK);}//getchar();}
-      pp[UU]=UURHORATIOMIN*pp[UU]; //increasing uint
+      pp[UU]=UURHORATIOMIN*pp[RHO]; //increasing uint
 
 
     ret=-1;
-Original file line number
+Diff line change
@@ Expand Up / @@ -414,7 +414,7 @@ check_floors_mhd(ldouble *pp, int whichvel,void *ggg) @@
         }
         if(verbose) {printf("hd_floors CASE 2 at (%d,%d,%d | %d,%d,%d): %e %e | tijk: %d %d %d\n",geom->ix+TOI,geom->iy+TOJ,geom->iz+TOK,geom->ix,geom->iy,geom->iz,pp[RHO],pp[UU],TI,TJ,TK);}//getchar();}
-          pp[UU]=UURHORATIOMIN*pp[UU]; //increasing uint
+          pp[UU]=UURHORATIOMIN*pp[RHO]; //increasing uint
         ret=-1;
@@ Expand Down @@