00001 subroutine IO_output_AHfinder 00002 use def_horizon, only : ahz 00003 use grid_parameter, only : npg, ntg 00004 implicit none 00005 integer :: ipg, itg 00006 ! 00007 ! --- Apparent horizon. 00008 open(13,file='bbhahz.las',status='unknown') 00009 write(13,'(5i5)') ntg, npg 00010 do ipg = 0, npg 00011 do itg = 0, ntg 00012 write(13,'(1p,6e20.12)') ahz(itg,ipg) 00013 end do 00014 end do 00015 close(13) 00016 ! 00017 end subroutine IO_output_AHfinder