aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorrhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26>2012-11-09 01:59:21 +0000
committerrhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26>2012-11-09 01:59:21 +0000
commit96387687ac1bdb56c345fa736e666f78f6ccb407 (patch)
tree1988fbc42d7f2ac7c2cea4f415ada31ef1d4e6e0
parent410956b99bc1dcc6b86447662ad9a2cefc8c32da (diff)
GRHydro_Init_Data: only write bondi data file on root process
From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@184 ac85fae7-cede-4708-beff-ae01c7fa1c26
-rw-r--r--src/GRHydro_BondiM_new.F9014
-rw-r--r--src/GRHydro_Bondi_new.F9014
2 files changed, 16 insertions, 12 deletions
diff --git a/src/GRHydro_BondiM_new.F90 b/src/GRHydro_BondiM_new.F90
index 346bc19..d03a1c2 100644
--- a/src/GRHydro_BondiM_new.F90
+++ b/src/GRHydro_BondiM_new.F90
@@ -149,12 +149,14 @@ subroutine GRHydro_BondiM_Iso(CCTK_ARGUMENTS)
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)
diff --git a/src/GRHydro_Bondi_new.F90 b/src/GRHydro_Bondi_new.F90
index b31ba7c..1fbdca6 100644
--- a/src/GRHydro_Bondi_new.F90
+++ b/src/GRHydro_Bondi_new.F90
@@ -142,12 +142,14 @@ subroutine GRHydro_Bondi_Iso(CCTK_ARGUMENTS)
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)