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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/11/08 07:59 PM

Modified:
 /trunk/src/
  GRHydro_BondiM_new.F90, GRHydro_Bondi_new.F90

Log:
 GRHydro_Init_Data: only write bondi data file on root process
 
 From: Roland Haas <roland.haas at physics.gatech.edu>

File Changes:

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

File [modified]: GRHydro_BondiM_new.F90
Delta lines: +8 -6
===================================================================
--- trunk/src/GRHydro_BondiM_new.F90	2012-11-09 01:59:19 UTC (rev 183)
+++ trunk/src/GRHydro_BondiM_new.F90	2012-11-09 01:59:21 UTC (rev 184)
@@ -149,12 +149,14 @@
      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)
+  if (CCTK_MyProc(cctkGH) == 0) then 
+    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)
+  end if
 
 !!$  write(debug_message,'(a,4f22.14)') "i=1:",r_bondi(1),rho_bondi(1),&
 !!$                                            u_bondi(1),v_bondi(1)

File [modified]: GRHydro_Bondi_new.F90
Delta lines: +8 -6
===================================================================
--- trunk/src/GRHydro_Bondi_new.F90	2012-11-09 01:59:19 UTC (rev 183)
+++ trunk/src/GRHydro_Bondi_new.F90	2012-11-09 01:59:21 UTC (rev 184)
@@ -142,12 +142,14 @@
      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)
+  if(CCTK_MyProc(cctkGH) .eq. 0) then
+     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)
+  end if
 
 !!$  write(debug_message,'(a,4f22.14)') "i=1:",r_bondi(1),rho_bondi(1),&
 !!$                                            u_bondi(1),v_bondi(1)