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

[rhaas at tapir.caltech.edu]

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

Modified:
 /trunk/src/
  GRHydro_Bondi_new.F90

Log:
 GRHydro_InitData: dump the bondi 1d solution to file.
 
 From: Bruno Coutinho Mundim <bcmsma at astro.rit.edu>

File Changes:

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

File [modified]: GRHydro_Bondi_new.F90
Delta lines: +22 -15
===================================================================
--- trunk/src/GRHydro_Bondi_new.F90	2012-09-17 16:27:26 UTC (rev 165)
+++ trunk/src/GRHydro_Bondi_new.F90	2012-09-17 16:27:27 UTC (rev 166)
@@ -139,19 +139,26 @@
      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 @@
                 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),&