Merge pull request #3066 from deslaughter/fastfarm_kdtree

Use k-d tree to speed up wake application to high-res grid in AWAE

Author: deslaughter

glue-codes/fast-farm/src/FAST_Farm_IO_Params.f90

@@ -9618,14 +9618,20 @@ MODULE FAST_Farm_IO_Params
       contains
 
 ! --- Functions not automatically generated
-logical function PointInAABB(x, y, z, x0, y0, z0, x1, y1, z1)
-   real(ReKi), intent(in) :: x,y,z,x0,y0,z0,x1,y1,z1
+logical function PointInAABB(x, y, z, minXYZ, sizeXYZ)
+   real(ReKi), intent(in) :: x,y,z
+   real(ReKi), intent(in) :: minXYZ(3), sizeXYZ(3)
+ 
+   ! Default to false
+   PointInAABB = .false.
 
-      ! default to return false
-   PointInAABB = .false.;
-      !Check if the point is less than max and greater than min
-    if (x >= x0 .and. x <= x1 .and. y >= y0 .and. y <= y1 .and. z >= z0 .and. z <= z1)  PointInAABB = .true.;
+   ! If point is outside box, return false
+   if (x < minXYZ(1) .or. x > minXYZ(1) + sizeXYZ(1)) return
+   if (y < minXYZ(2) .or. y > minXYZ(2) + sizeXYZ(2)) return
+   if (z < minXYZ(3) .or. z > minXYZ(3) + sizeXYZ(3)) return
 
+   ! If we reach here, the point is inside the AABB
+   PointInAABB = .true.
 
 end function PointInAABB
 
@@ -15504,7 +15510,8 @@ SUBROUTINE Farm_SetOutParam(OutList, farm, ErrStat, ErrMsg )
 
       ! Add checks for the WindVel locations based on knowledge of the wind grids and NWindVel
    do i = 1, farm%p%NWindVel
-      if (.not. PointInAABB(farm%p%WindVelX(i), farm%p%WindVelY(i), farm%p%WindVelZ(i), farm%AWAE%p%X0_low, farm%AWAE%p%Y0_low,farm%AWAE%p%Z0_low, farm%AWAE%p%X0_low+(farm%AWAE%p%nX_low-1)*farm%AWAE%p%dX_low, farm%AWAE%p%Y0_low+(farm%AWAE%p%nY_low-1)*farm%AWAE%p%dY_low, farm%AWAE%p%Z0_low+(farm%AWAE%p%nZ_low-1)*farm%AWAE%p%dZ_low) ) then
+      if (.not. PointInAABB(farm%p%WindVelX(i), farm%p%WindVelY(i), farm%p%WindVelZ(i), &
+                            farm%AWAE%p%LowRes%oXYZ, farm%AWAE%p%LowRes%Size)) then
          InvalidOutput( WVAmbX  (i) ) = .true.
          InvalidOutput( WVAmbY  (i) ) = .true.
          InvalidOutput( WVAmbZ  (i) ) = .true.

glue-codes/fast-farm/src/FAST_Farm_Subs.f90

@@ -245,15 +245,15 @@ SUBROUTINE Farm_Initialize( farm, InputFile, ErrStat, ErrMsg )
 
    farm%AWAE%IsInitialized = .true.
 
-   farm%p%X0_Low = AWAE_InitOutput%X0_Low
-   farm%p%Y0_low = AWAE_InitOutput%Y0_low
-   farm%p%Z0_low = AWAE_InitOutput%Z0_low
-   farm%p%nX_Low = AWAE_InitOutput%nX_Low
-   farm%p%nY_low = AWAE_InitOutput%nY_low
-   farm%p%nZ_low = AWAE_InitOutput%nZ_low
-   farm%p%dX_low = AWAE_InitOutput%dX_low
-   farm%p%dY_low = AWAE_InitOutput%dY_low
-   farm%p%dZ_low = AWAE_InitOutput%dZ_low
+   farm%p%X0_Low = AWAE_InitOutput%oXYZ_Low(1)
+   farm%p%Y0_low = AWAE_InitOutput%oXYZ_Low(2)
+   farm%p%Z0_low = AWAE_InitOutput%oXYZ_Low(3)
+   farm%p%nX_Low = AWAE_InitOutput%nXYZ_Low(1)
+   farm%p%nY_low = AWAE_InitOutput%nXYZ_Low(2)
+   farm%p%nZ_low = AWAE_InitOutput%nXYZ_Low(3)
+   farm%p%dX_low = AWAE_InitOutput%dXYZ_Low(1)
+   farm%p%dY_low = AWAE_InitOutput%dXYZ_Low(2)
+   farm%p%dZ_low = AWAE_InitOutput%dXYZ_Low(3)
    farm%p%Module_Ver( ModuleFF_AWAE  ) = AWAE_InitOutput%Ver
 
       !-------------------
@@ -653,10 +653,6 @@ SUBROUTINE Farm_InitFAST( farm, WD_InitInp, AWAE_InitOutput, ErrStat, ErrMsg )
       FWrap_InitInp%tmax          = farm%p%TMax
       FWrap_InitInp%n_high_low    = farm%p%n_high_low + 1   ! Add 1 because the FAST wrapper uses an index that starts at 1
       FWrap_InitInp%dt_high       = farm%p%dt_high
-     
-      FWrap_InitInp%nX_high       = AWAE_InitOutput%nX_high
-      FWrap_InitInp%nY_high       = AWAE_InitOutput%nY_high
-      FWrap_InitInp%nZ_high       = AWAE_InitOutput%nZ_high
 
       if (farm%p%MooringMod > 0) then
          FWrap_Interval = farm%p%dt_mooring    ! when there is a farm-level mooring model, FASTWrapper will be called at the mooring coupling time step
@@ -676,14 +672,17 @@ SUBROUTINE Farm_InitFAST( farm, WD_InitInp, AWAE_InitOutput, ErrStat, ErrMsg )
          FWrap_InitInp%TurbNum       = nt
          FWrap_InitInp%RootName      = trim(farm%p%OutFileRoot)//'.T'//num2lstr(nt)
 
-         
-         FWrap_InitInp%p_ref_high(1) = AWAE_InitOutput%X0_high(nt)
-         FWrap_InitInp%p_ref_high(2) = AWAE_InitOutput%Y0_high(nt)
-         FWrap_InitInp%p_ref_high(3) = AWAE_InitOutput%Z0_high(nt)
+         FWrap_InitInp%nX_high       = AWAE_InitOutput%nXYZ_high(1,nt)
+         FWrap_InitInp%nY_high       = AWAE_InitOutput%nXYZ_high(2,nt)
+         FWrap_InitInp%nZ_high       = AWAE_InitOutput%nXYZ_high(3,nt)
+
+         FWrap_InitInp%p_ref_high(1) = AWAE_InitOutput%oXYZ_high(1,nt)
+         FWrap_InitInp%p_ref_high(2) = AWAE_InitOutput%oXYZ_high(2,nt)
+         FWrap_InitInp%p_ref_high(3) = AWAE_InitOutput%oXYZ_high(3,nt)
 
-         FWrap_InitInp%dX_high       = AWAE_InitOutput%dX_high(nt)
-         FWrap_InitInp%dY_high       = AWAE_InitOutput%dY_high(nt)
-         FWrap_InitInp%dZ_high       = AWAE_InitOutput%dZ_high(nt)
+         FWrap_InitInp%dX_high       = AWAE_InitOutput%dXYZ_high(1,nt)
+         FWrap_InitInp%dY_high       = AWAE_InitOutput%dXYZ_high(2,nt)
+         FWrap_InitInp%dZ_high       = AWAE_InitOutput%dXYZ_high(3,nt)
 
          FWrap_InitInp%Vdist_High   => AWAE_InitOutput%Vdist_High(nt)%data