aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authordiener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c>2005-02-03 21:01:58 +0000
committerdiener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c>2005-02-03 21:01:58 +0000
commitb5c183ea222b363d1ad0845cb42d9f88c9338cfa (patch)
tree7fbcbdbd18665e78b36e98811b109dc46a698e56
parent367044234c0f6c70df58464c9cce3252c46bf139 (diff)
Fixed subroutine call parameters so that it compiles.
Still doesn't work, though. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@194 2a26948c-0e4f-0410-aee8-f1d3e353619c
-rwxr-xr-xsrc/EHFinder_IsoSurface_optimized.F9016
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