diff options
| -rwxr-xr-x | src/EHFinder_IsoSurface_optimized.F90 | 16 |
1 files changed, 7 insertions, 9 deletions
diff --git a/src/EHFinder_IsoSurface_optimized.F90 b/src/EHFinder_IsoSurface_optimized.F90 index fcadde3..1733e20 100755 --- a/src/EHFinder_IsoSurface_optimized.F90 +++ b/src/EHFinder_IsoSurface_optimized.F90 @@ -29,17 +29,15 @@ subroutine EHFinder_IsoSurface(CCTK_ARGUMENTS) call CCTK_INFO ('Entered IsoSurface') -! print*,nx,ny,nz -! print*,size(f) -! print*,Xf0,Xf1,Xf2,f_length -! print*,size(x) -! allocate(x_coord(5),y_coord(5),z_coord(5),t_index(20)) + nx = cctk_lsh(1); ny = cctk_lsh(2); nz = cctk_lsh(3) + print*,nx,ny,nz + print*,size(f) + print*,Xf0,Xf1,Xf2,f_length + print*,size(x) + allocate(x_coord(5),y_coord(5),z_coord(5),t_index(20)) do l = 1, eh_number_level_sets -! call triangularization ( n_coords, n_triangles, nx, ny, nz, x_coord, & -! y_coord, z_coord, t_index, f(:,:,:,l), x, y, z ) -! call triangularization ( n_coords, n_triangles, nx, ny, nz, t_index, f(:,:,:,l), x, y, z ) - call triangularization(l,contor,nx,ny,nz,t_index,f(:,:,:,l),x,y,z) + call triangularization(l,contor,nx,ny,nz,x_coord,y_coord,z_coord,t_index,f(:,:,:,l),x(:,1,1),y(1,:,1),z(1,1,:)) print*,'The surface consists of ',n_triangles/3,' triangles and ',n_coords,' coords' end do end subroutine EHFinder_IsoSurface |