From 81feb0023e8ef4585d77dafa9b2ce6430692bc23 Mon Sep 17 00:00:00 2001 From: rhaas Date: Mon, 17 Sep 2012 16:27:27 +0000 Subject: GRHydro_InitData: dump the bondi 1d solution to file. From: Bruno Coutinho Mundim git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@166 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_Bondi_new.F90 | 37 ++++++++++++++++++++++--------------- 1 file 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),& -- cgit v1.2.3