diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-09-17 16:27:29 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-09-17 16:27:29 +0000 |
| commit | 902036305f1ac06a5a20786c7de02b82a340d2c6 (patch) | |
| tree | 8cd80d14dc0fea6c0d02df7f1456b1df3b6c91ac | |
| parent | 81feb0023e8ef4585d77dafa9b2ce6430692bc23 (diff) | |
GRHydro_InitData: pump up the number of 1d grid points for bondi solution.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@167 ac85fae7-cede-4708-beff-ae01c7fa1c26
| -rwxr-xr-x | src/GRHydro_Bondi_new.F90 | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/GRHydro_Bondi_new.F90 b/src/GRHydro_Bondi_new.F90 index 9bcf7e1..7bb5548 100755 --- a/src/GRHydro_Bondi_new.F90 +++ b/src/GRHydro_Bondi_new.F90 @@ -52,7 +52,7 @@ subroutine GRHydro_Bondi_Iso(CCTK_ARGUMENTS) 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 @@ subroutine GRHydro_Bondi_Iso(CCTK_ARGUMENTS) 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),& |