diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-09-17 16:27:27 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-09-17 16:27:27 +0000 |
| commit | 81feb0023e8ef4585d77dafa9b2ce6430692bc23 (patch) | |
| tree | 690375fc8d26cbe585e274641b08df6dd5dbd2d9 | |
| parent | a25f8ba2499d399d44916f12607167e597076298 (diff) | |
GRHydro_InitData: dump the bondi 1d solution to file.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@166 ac85fae7-cede-4708-beff-ae01c7fa1c26
| -rwxr-xr-x | src/GRHydro_Bondi_new.F90 | 37 |
1 files changed, 22 insertions, 15 deletions
diff --git a/src/GRHydro_Bondi_new.F90 b/src/GRHydro_Bondi_new.F90 index a2eabcd..9bcf7e1 100755 --- a/src/GRHydro_Bondi_new.F90 +++ b/src/GRHydro_Bondi_new.F90 @@ -139,19 +139,26 @@ subroutine GRHydro_Bondi_Iso(CCTK_ARGUMENTS) u_bondi(i) = utmp v_bondi(i) = vtmp enddo + + open (47,file="bondi.asc",form="formatted") + do i=1,N_points + write(47,'(i5,4f22.14)')i,r_bondi(i),rho_bondi(i),& + u_bondi(i),v_bondi(i) + end do + close(47) - write(debug_message,'(a,4f22.14)') "i=1:",r_bondi(1),rho_bondi(1),& - u_bondi(1),v_bondi(1) - call CCTK_INFO(debug_message) - write(debug_message,'(a,4f22.14)') "i=100:",r_bondi(100),rho_bondi(100),& - u_bondi(100),v_bondi(100) - call CCTK_INFO(debug_message) - write(debug_message,'(a,4f22.14)') "i=1000:",r_bondi(1000),rho_bondi(1000),& - u_bondi(1000),v_bondi(1000) - call CCTK_INFO(debug_message) - write(debug_message,'(a,4f22.14)') "i=1500:",r_bondi(1500),rho_bondi(1500),& - u_bondi(1500),v_bondi(1500) - call CCTK_INFO(debug_message) +!!$ write(debug_message,'(a,4f22.14)') "i=1:",r_bondi(1),rho_bondi(1),& +!!$ u_bondi(1),v_bondi(1) +!!$ call CCTK_INFO(debug_message) +!!$ write(debug_message,'(a,4f22.14)') "i=100:",r_bondi(100),rho_bondi(100),& +!!$ u_bondi(100),v_bondi(100) +!!$ call CCTK_INFO(debug_message) +!!$ write(debug_message,'(a,4f22.14)') "i=1000:",r_bondi(1000),rho_bondi(1000),& +!!$ u_bondi(1000),v_bondi(1000) +!!$ call CCTK_INFO(debug_message) +!!$ write(debug_message,'(a,4f22.14)') "i=1500:",r_bondi(1500),rho_bondi(1500),& +!!$ u_bondi(1500),v_bondi(1500) +!!$ call CCTK_INFO(debug_message) !!$ // find the derivative near r=M rnew = 2.25 * M @@ -223,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(abs(z(i,j,k)).lt.tiny.and. & - x(i,j,k).gt.6.7d0.and.x(i,j,k).lt.6.9d0.and. & - y(i,j,k).gt.1.34d0.and.y(i,j,k).lt.1.38d0 ) 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 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),& |