[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 167)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/09/17 11:27 AM

Modified:
 /trunk/src/
  GRHydro_Bondi_new.F90

Log:
 GRHydro_InitData: pump up the number of 1d grid points for bondi solution.
 
 From: Bruno Coutinho Mundim <bcmsma at astro.rit.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Bondi_new.F90
Delta lines: +4 -4
===================================================================
--- trunk/src/GRHydro_Bondi_new.F90	2012-09-17 16:27:27 UTC (rev 166)
+++ trunk/src/GRHydro_Bondi_new.F90	2012-09-17 16:27:29 UTC (rev 167)
@@ -52,7 +52,7 @@
   
   CCTK_INT :: i, j, k, nx, ny, nz, imin, jb,N_points
   CCTK_REAL :: ONEmTINY, tiny
-  PARAMETER (N_points=2000,ONEmTINY=0.999999d0,tiny=1.0d-12)
+  PARAMETER (N_points=100000,ONEmTINY=0.999999d0,tiny=1.0d-12)
   CCTK_REAL :: M, Msq, Mdot, rs, gam, rmin_bondi, rmax_bondi, cs_sq,cs,vs_sq,vs,rhos,gtemp,hs, Kval, Qdot
   CCTK_REAL :: logrmin,dlogr,rhotmp,utmp,vtmp,rspher
   CCTK_REAL :: r_bondi(N_points), logr_bondi(N_points), rho_bondi(N_points), u_bondi(N_points), v_bondi(N_points)
@@ -230,9 +230,9 @@
                 eps(i,j,k),press(i,j,k),w_lorentz(i,j,k))
            
 !           if(riso.gt.1.014d0.and.riso.lt.1.015) then 
-           if(z(i,j,k).ge.-2.0d-1.and.z(i,j,k).le.2.0d-1.and. &
-              x(i,j,k).ge.6.8d0.and.x(i,j,k).le.7.0d0.and. &
-              y(i,j,k).ge.1.2d0.and.y(i,j,k).le.1.4d0 ) then 
+           if(z(i,j,k).ge.-1.0d0.and.z(i,j,k).le.1.0d0.and. &
+              x(i,j,k).ge.6.5d0.and.x(i,j,k).le.7.5d0.and. &
+              y(i,j,k).ge.1.0d0.and.y(i,j,k).le.1.5d0 ) then 
              write(debug_message,'(a,15f22.14)') 'Point to check:', &
                x(i,j,k),y(i,j,k),z(i,j,k),riso,gxx(i,j,k),dens(i,j,k),&
                tau(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),rho(i,j,k),eps(i,j,k),&