cross merge rc-5.0.1 into dev

.gitignore

@@ -64,3 +64,8 @@ varcache
 # Python cache files
 openfast_io/dist/
 openfast_io/openfast_io/_version.py
+
+# docker with VScode config
+*.code-workspace
+docker-compose.yml
+.devcontainer/

docs/source/user/aerodyn-olaf/InputFiles.rst

@@ -312,6 +312,59 @@ of a box of shape 5x20x30 and dimension 1200x300x295.  The grid contains both th
 The two other grids are vertical and horizontal planes containing only the velocity.
 
 
+Non-equidistant grid points
+~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+By default, grid points along each axis (X, Y, Z) are equidistant, defined by
+the ``Start``, ``End``, and ``n`` columns of the grid output table.
+
+Alternatively, a ``Start`` cell can be given as a quoted filename instead of a
+number. In that case, the corresponding axis uses an explicit, user-defined
+list of coordinates read from that file, and the ``End``/``n`` columns for
+that axis are ignored.
+
+This choice is made independently for each axis, so a single grid can mix
+equidistant and list-defined axes (e.g., equidistant in X and Z, list-defined
+in Y), and a filename can be given for ``XStart``, ``YStart``, ``ZStart``, or
+any combination of the three.
+
+The referenced file must:
+
+- contain exactly one coordinate value per line,
+- list values in strictly ascending order (no duplicates),
+- contain no comments or blank lines,
+- be located in the same directory as the main OLAF input file.
+
+Example, requesting a non-equidistant Y-axis::
+
+    GridName  GridType  TStart  TEnd     DTOut     XStart    XEnd   nX      YStart     YEnd    nY   ZStart   ZEnd   nZ
+    (-)         (-)      (s)     (s)      (s)        (m)      (m)   (-)       (m)       (m)   (-)    (m)      (m)  (-)
+    "Yline"      1     default default  default      90       90     1   "Ypoints.dat"   -     -     100      100   1
+
+with ``Ypoints.dat`` (located next to the OLAF input file) containing::
+
+    -125
+    -100
+    -75
+    -50
+    -37.5
+    -25
+    -12.5
+     0
+     12.5
+     25
+     37.5
+     50
+     75
+     100
+     125
+
+.. note::
+   Vorticity output (``GridType=2``) requires equidistant spacing
+   and is not available for a grid that uses a list-defined axis.
+   Use ``GridType=1`` (velocity only) for such grids.
+
+
 Advanced Options
 ~~~~~~~~~~~~~~~~
 

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

@@ -1047,6 +1047,7 @@ SUBROUTINE Farm_ValidateInput( p, WD_InitInp, AWAE_InitInp, ErrStat, ErrMsg )
    CHARACTER(*),   PARAMETER     :: RoutineName = 'Farm_ValidateInput'
    INTEGER(IntKi)                :: n_disDT_dt
    character(60)                 :: tmpStr
+   logical                       :: file_exists                               ! Flag for file existence check
 
    ErrStat = ErrID_None
    ErrMsg  = ""
@@ -1066,8 +1067,11 @@ SUBROUTINE Farm_ValidateInput( p, WD_InitInp, AWAE_InitInp, ErrStat, ErrMsg )
    ! TODO : Verify that the DLL file exists
 
    ! --- SHARED MOORING SYSTEM ---
-   ! TODO : Verify that p%MD_FileName file exists
-   if ((p%DT_mooring <= 0.0_ReKi) .or. (p%DT_mooring > p%DT_high)) CALL SetErrStat(ErrID_Fatal,'DT_mooring must be greater than zero and no greater than dt_high.',ErrStat,ErrMsg,RoutineName)
+   if (p%MooringMod == 3) then
+      inquire(file=trim(p%MD_FileName), exist=file_exists)
+      if (.not. file_exists) call SetErrStat(ErrID_Fatal,'Cannot find MoorDyn input file '//trim(p%MD_FileName),ErrStat,ErrMsg,RoutineName)
+      if ((p%DT_mooring <= 0.0_ReKi) .or. (p%DT_mooring > p%DT_high)) CALL SetErrStat(ErrID_Fatal,'DT_mooring must be greater than zero and no greater than dt_high.',ErrStat,ErrMsg,RoutineName)
+   end if
 
    ! --- AMBIENT WIND: INFLOWWIND MODULE --- [used only for Mod_AmbWind=2 or 3] ---
    ! FIXME: this really should be checked with the turbine specific size diameter -- maybe relocate this check to AWAE or in FF after initializing all turbines?

modules/aerodyn/src/AeroAcoustics.f90

@@ -165,6 +165,7 @@ subroutine SetParameters( InitInp, InputFileData, p, AFInfo, ErrStat, ErrMsg )
     character(*), parameter :: RoutineName = 'SetParameters'
     REAL(ReKi)              :: val1,val10,f2,f4, dist1, dist10
     REAL(ReKi)              :: BladeSpanUsedForNoise
+    REAL(ReKi)              :: LastElemPct
 
     ! Initialize variables for this routine
     ErrStat  = ErrID_None
@@ -271,15 +272,30 @@ subroutine SetParameters( InitInp, InputFileData, p, AFInfo, ErrStat, ErrMsg )
     p%BlSpn   = InitInp%BlSpn
     p%BlChord = InitInp%BlChord
 
-    p%startnode = max(1, p%NumBlNds - 1)
+    ! Calculate last element size as percentage of blade span
+    IF (p%NumBlNds > 1) THEN
+        LastElemPct = 100.0 * (p%BlSpn(p%NumBlNds,1) - p%BlSpn(p%NumBlNds-1,1)) / p%BlSpn(p%NumBlNds,1)
+    ELSE
+        LastElemPct = 100.0  ! Single node means element spans entire blade
+    ENDIF
+
+    IF (InputFileData%AA_Bl_Prcntge .lt. LastElemPct) THEN
+        CALL SetErrStat(ErrID_Warn, 'BldPrcnt is smaller than the last blade element size. '// &
+            'OpenFAST will move on assuming the last blade element, which is the minimum blade span used for noise calculations. '// &
+            'Either increase BldPrcnt in your aeroacoustic input file '// &
+            'or refine the aerodynamic mesh in your blade AeroDyn input file.', ErrStat2, ErrMsg2, RoutineName )
+        if(Failed()) return
+    endif
+
+    p%startnode = min(p%NumBlNds, 2)
     BladeSpanUsedForNoise = p%BlSpn(p%NumBlNds,1)*(1.0 - InputFileData%AA_Bl_Prcntge/100.0)
     do j=p%NumBlNds-1,2,-1
         IF ( p%BlSpn(j,1) .lt. BladeSpanUsedForNoise )THEN
             p%startnode=j
             exit ! exit the loop
         endif
     enddo
-    p%startnode = max(min(p%NumBlNds,2),p%startnode)
+
 
    p%BlElemSpn = 0;
    DO I = 1,p%numBlades