From 8c5621196e518865679d58476d7df474016f4364 Mon Sep 17 00:00:00 2001 From: andrew-platt Date: Fri, 19 Jul 2024 11:07:18 -0600 Subject: [PATCH] Vis: Fix bug with vtk writing of non-square wave surface When a rectangular SeaState grid was used, the connection list in the vtk file was incorrectly generated due to an indexing error. The list would contain points beyond the total number of points in the point list. For a 5x3 SeaState grid (Nx=3, Ny=2) there are a total of 15 points (0 to 14). The connection list was being incorrectly generated as ``` 0 1 3 1 4 3 1 2 4 2 5 4 5 6 8 6 9 8 6 7 9 7 10 9 10 11 13 11 14 13 11 12 14 12 15 14 15 16 18 16 19 18 16 17 19 17 20 19 ``` The connection list should be ``` 0 1 3 1 4 3 1 2 4 2 5 4 3 4 6 4 7 6 4 5 7 5 8 7 6 7 9 7 10 9 7 8 10 8 11 10 9 10 12 10 13 12 10 11 13 11 14 13 ``` Co-authored-by: luwang00 --- modules/openfast-library/src/FAST_Subs.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/modules/openfast-library/src/FAST_Subs.f90 b/modules/openfast-library/src/FAST_Subs.f90 index 09a26558d0..1c0d0b50bb 100644 --- a/modules/openfast-library/src/FAST_Subs.f90 +++ b/modules/openfast-library/src/FAST_Subs.f90 @@ -8686,7 +8686,7 @@ SUBROUTINE WrVTK_WaveElevVisGrid(t_global, p_FAST, y_FAST, SeaSt) do ix=1,p_FAST%VTK_surface%NWaveElevPts(1)-1 do iy=1,p_FAST%VTK_surface%NWaveElevPts(2)-1 - n = p_FAST%VTK_surface%NWaveElevPts(1)*(ix-1)+iy - 1 ! points start at 0 + n = p_FAST%VTK_surface%NWaveElevPts(2)*(ix-1)+iy - 1 ! points start at 0 WRITE(Un,'(3(i7))') n, n+1, n+p_FAST%VTK_surface%NWaveElevPts(2) WRITE(Un,'(3(i7))') n+1, n+1+p_FAST%VTK_surface%NWaveElevPts(2), n+p_FAST%VTK_surface%NWaveElevPts(2)