diff options
Diffstat (limited to 'src/AHFinder_output.F')
-rw-r--r-- | src/AHFinder_output.F | 28 |
1 files changed, 11 insertions, 17 deletions
diff --git a/src/AHFinder_output.F b/src/AHFinder_output.F index c9f4f8b..9073c3b 100644 --- a/src/AHFinder_output.F +++ b/src/AHFinder_output.F @@ -131,6 +131,8 @@ c some compilers cannot trim an empty string, so we add at least one char character(len=200) :: filename + call AHFinder_SetReflections(CCTK_ARGUMENTS) + open(1,file=filename,form='formatted', status='replace') write(1,"(A20)") '# GAUSSIAN CURVATURE' @@ -158,13 +160,8 @@ c some compilers cannot trim an empty string, so we add at least one char write(1,"(A9,L1)") '# refy = ',refy write(1,"(A9,L1)") '# refz = ',refz write(1,"(A1)") '#' - write(1,"(A11,I7)") '# ntheta = ',(ntheta+1) - write(1,"(A11,I7)") '# nphi = ',(nphi+1) - - write(1,*) - write(1,"(A11,I4)") '# Time step',cctk_iteration - write(1,"(A6,ES11.3)") '# Time',cctk_time - write(1,"(A6,I4)") '# Call',ahf_ncall + write(1,"(A11,I7)") '# ntheta = ',(ahf_ntheta+1) + write(1,"(A11,I7)") '# nphi = ',(ahf_nphi+1) close(1) @@ -458,21 +455,12 @@ c some compilers cannot trim an empty string, so we add at least one char inquire(file=almf, exist=almf_exists) if (almf_exists) then open(1,file=almf,form='formatted',status='old',position='append') - write(1,*) else open(1,file=almf,form='formatted',status='replace') write(1,"(A21)") '# Radial coefficients' write(1,"(A1)") '#' end if - if (status.and.report) then - out_centerx = xc - out_centery = yc - out_centerz = zc - - write(1,"(A15,3ES14.6)") '# centered on: ', out_centerx, out_centery, out_centerz - endif - write(1,"(A12,I4)") '# Time step ',cctk_iteration write(1,"(A7,ES14.6)") '# Time ',cctk_time write(1,"(A7,I4)") '# Call ',ahf_ncall @@ -503,6 +491,12 @@ c some compilers cannot trim an empty string, so we add at least one char write(1,"(A30,I4)") '# Surface found: Not a horizon' end if end if + + out_centerx = xc + out_centery = yc + out_centerz = zc + + write(1,"(A15,3ES14.6)") '# centered on: ', out_centerx, out_centery, out_centerz write(1,"(A1)") '#' write(1,"(A22)") '# a_lm l m' @@ -525,7 +519,7 @@ c some compilers cannot trim an empty string, so we add at least one char end do end do end if - + write(1,*) else write(1,"(A18,I4)") '# No surface found' end if |