OpenFAST · andrew-platt · May 16, 2025 · May 16, 2025 · May 16, 2025 · Copilot
diff --git a/docs/source/dev/index.rst b/docs/source/dev/index.rst
@@ -61,7 +61,7 @@ with the NREL OpenFAST team to define the scope of the work and coordinate
 development efforts. This is particularly important since many groups
 work on OpenFAST simultaneously. By engaging early, all developers can
 stay up to date and minimize conflicts during the code merge.
-The prefered method of communication is `GitHub Issues <https://github.com/openfast/openfast/issues>`_.
+The preferred method of communication is `GitHub Issues <https://github.com/openfast/openfast/issues>`_.
 An initial post should contain all relevant information about the planned
 development work, the areas of the software that will be impacted,
 and any model validation materials. See :ref:`development_plan`

diff --git a/docs/source/install/index.rst b/docs/source/install/index.rst
@@ -478,7 +478,7 @@ The CMake options specific to OpenFAST and their default settings are:
     ORCA_DLL_LOAD                  - Enable OrcaFlex library load (Default: ON)
     USE_DLL_INTERFACE              - Enable runtime loading of dynamic libraries (Default: ON)
     USE_LOCAL_STATIC_LAPACK        - Enable downloading and building static LAPACK and BLAS libs (Default: OFF)
-    VARIABLE_TRACKING              - Enables variable tracking for better runtime debugging output. May increase compile time. Valid only for GNU. (Defualt: ON)
+    VARIABLE_TRACKING              - Enables variable tracking for better runtime debugging output. May increase compile time. Valid only for GNU. (Default: ON)
 
 
 Additional system-specific options may exist for a given system, but those

diff --git a/docs/source/user/aerodyn-olaf/OLAFTheory.rst b/docs/source/user/aerodyn-olaf/OLAFTheory.rst
@@ -94,7 +94,7 @@ a line of varying circulation. The line follows the motion of the blade and is
 referred to as “bound” circulation. The bound circulation does not follow the
 same dynamic equation as the free vorticity of the wake. Instead, the intensity
 is linked to airfoil lift via the Kutta-Joukowski theorem. Spanwise variation of
-the bound circulation results in vorticity being emitted into the the wake. This
+the bound circulation results in vorticity being emitted into the wake. This
 is referred to as “trailed vorticity”. Time changes of the bound circulation are
 also emitted in the wake, referred to as “shed” vorticity. The subsequent
 paragraphs describe the representation of the bound vorticity.

diff --git a/docs/source/user/aerodyn/input.rst b/docs/source/user/aerodyn/input.rst
@@ -187,7 +187,7 @@ Determines the kind of BEM algorithm to use.
 
 .. warning::
 
-    ``BEM_Mod`` currently governs the coordinate system used for "ill-defined" outputs (outputs that don't have a specified coordinate system) such as the ones that ends with "x" and "y". Other ill-defined outputs are the typical BEM quantities such as "AxInd", "TnInd", "Phi", etc. These are defined in a different coordinate system depending on `BEM_Mod`. For consistency accross differents `Wake_Mod` (even when `Wake_Mod/=1`), we use `BEM_Mod` to determine the coordinate system of the ill-defined outputs. 
+    ``BEM_Mod`` currently governs the coordinate system used for "ill-defined" outputs (outputs that don't have a specified coordinate system) such as the ones that ends with "x" and "y". Other ill-defined outputs are the typical BEM quantities such as "AxInd", "TnInd", "Phi", etc. These are defined in a different coordinate system depending on `BEM_Mod`. For consistency across differents `Wake_Mod` (even when `Wake_Mod/=1`), we use `BEM_Mod` to determine the coordinate system of the ill-defined outputs. 
 
 The following inputs in this section are only used when ``Wake_Mod = 1``.
 

diff --git a/docs/source/user/aerodyn/theory.rst b/docs/source/user/aerodyn/theory.rst
@@ -4,7 +4,7 @@
 AeroDyn Theory
 ==============
 
-This theory manual is work in progress, please refer to the AeroDyn 14 manual for more details :cite:`ad-AeroDyn:manual`. Many changes have occured since AeroDyn 14 (e.g. BEM formulation, coordinate system used in the BEM equations, dynamic stall, dynamic BEM), but these changes are not yet documented here.
+This theory manual is work in progress, please refer to the AeroDyn 14 manual for more details :cite:`ad-AeroDyn:manual`. Many changes have occurred since AeroDyn 14 (e.g. BEM formulation, coordinate system used in the BEM equations, dynamic stall, dynamic BEM), but these changes are not yet documented here.
 
 
 

diff --git a/docs/source/user/beamdyn/theory.rst b/docs/source/user/beamdyn/theory.rst
@@ -82,7 +82,7 @@ Local blade coordinate system
 
 The local blade coordinate system is used for some input and output
 quantities, for example, the cross-sectional mass and stiffness matrices
-and the the sectional force and moment resultants. This coordinate
+and the sectional force and moment resultants. This coordinate
 system is different from the blade reference coordinate system in that
 its :math:`Z_l` axis is always tangent to the blade axis as the blade
 deflects. Note that a subscript :math:`l` denotes the local blade

diff --git a/docs/source/user/cppapi/files/cDriver.i b/docs/source/user/cppapi/files/cDriver.i
@@ -27,7 +27,7 @@ n_substeps: 1            # Number of substeps per timestep of the glue-code
 
 n_checkpoint: 160        # Restart files will be written every so many time steps
 
-set_exp_law_wind: false  # Set velocity at the the turbine using an exponential law profile.
+set_exp_law_wind: false  # Set velocity at the turbine using an exponential law profile.
 
 Turbine0:
 

diff --git a/docs/source/user/elastodyn/theory.rst b/docs/source/user/elastodyn/theory.rst
@@ -5,7 +5,7 @@ ElastoDyn Theory
 ================
 
 Note this document is work in progress and is greatly incomplete. 
-This documentation was started to document some code changes to the the tail furl and rotor furl part of ElastoDyn. 
+This documentation was started to document some code changes to the tail furl and rotor furl part of ElastoDyn. 
 Please refer to the different ressources provided in :numref:`ed_intro` for additional documents.
 
 

diff --git a/docs/source/user/fast.farm/FFarmTheory.rst b/docs/source/user/fast.farm/FFarmTheory.rst
@@ -1278,7 +1278,7 @@ where :math:`D` is a reference diameter, and :math:`\bar{U}` is the mean
 velocity taken as the filtered velocity at the turbine location normal to the
 rotor disk.  The coordinates :math:`x,y,z` and :math:`r,\theta` are taken in the
 meandering frame of reference.  The parameters :math:`k_\text{def}^\text{WAT}`
-and :math:`k_\text{grad}^\text{WAT}` are tuning paramters of the model
+and :math:`k_\text{grad}^\text{WAT}` are tuning parameters of the model
 respectively multiplying the quasi-steady wake deficit and the gradient of the
 wake deficit. These are based on an eddy-viscosity filter with five calibrated
 parameters to give a more realistic dependence on downstream position.  The

diff --git a/docs/source/user/fast.farm/OutputFiles.rst b/docs/source/user/fast.farm/OutputFiles.rst
@@ -138,5 +138,5 @@ module levels, as specified within the OpenFAST input files) are
 described in the OpenFAST documentation and include summary (*.sum*)
 files, time-series results (ASCI *.out* or binary *.outb*) files,
 visualization (*.vtk*) files, etc. FAST.Farm simulations will generate
-these same files, but with the the path/rootname changed to *<RootName
+these same files, but with the path/rootname changed to *<RootName
 of WT_FASTInFile>.T<n*\ :sub:`t`\ *>*.
diff --git a/docs/source/user/subdyn/theory.rst b/docs/source/user/subdyn/theory.rst
@@ -1879,7 +1879,7 @@ This applies e.g.: to :math:`F_{R,e}, F_{L,e}, F_{R,g}, F_{L,g}`, where the foll
 
 
 The dependency of the load vectors on :math:`U_{TP}` introduces some complications for the state space representation, where for instance the :math:`B` and  :math:`F_X` matrices should be modified to account for the dependency in :math:`U_{TP}` in Eq. :eq:`ABFx`.
-The equation remains valid even if :math:`F_{L,e}` and :math:`F_{L,g}` contains a dependency in :math:`U_{TP}`, but the matrix :math:`B` shouldn't be used for the linearization (numerical differentiation is then prefered for simplicity).
+The equation remains valid even if :math:`F_{L,e}` and :math:`F_{L,g}` contains a dependency in :math:`U_{TP}`, but the matrix :math:`B` shouldn't be used for the linearization (numerical differentiation is then preferred for simplicity).
 Similar considerations apply for Eq. :eq:`bigY2`.
 
 

diff --git a/glue-codes/fast-farm/src/FAST_Farm.f90 b/glue-codes/fast-farm/src/FAST_Farm.f90
@@ -49,8 +49,8 @@ PROGRAM FAST_Farm
 REAL(ReKi)                            :: PrevClockTime      !< Previous clock time in seconds past midnight
 INTEGER                               :: SimStrtTime (8)    !< An array containing the elements of the start time (after initialization).
 INTEGER                               :: ProgStrtTime (8)   !< An array containing the elements of the program start time (before initialization).
-REAL(ReKi)                            :: SimStrtCPU         !< User CPU time for simulation (without intialization)
-REAL(ReKi)                            :: ProgStrtCPU        !< User CPU time for program (with intialization)
+REAL(ReKi)                            :: SimStrtCPU         !< User CPU time for simulation (without initialization)
+REAL(ReKi)                            :: ProgStrtCPU        !< User CPU time for program (with initialization)
 
 ! these should probably go in the FAST.Farm registry:
 type(All_FastFarm_Data)               :: farm  

diff --git a/glue-codes/fast-farm/src/FAST_Farm_Subs.f90 b/glue-codes/fast-farm/src/FAST_Farm_Subs.f90
@@ -341,7 +341,7 @@ SUBROUTINE WAT_init( p, WAT_IfW, AWAE_InitInput, ErrStat, ErrMsg )
    HAWC_InitInput%WindFileName(3) = trim(BoxFileRoot)//trim(FileEnding(3))
 
    ! HAWC spatial grid
-   if (p%WAT == Mod_WAT_PreDef) then       ! from libary of WAT files, set the NxNyNz and DxDyDz terms
+   if (p%WAT == Mod_WAT_PreDef) then       ! from library of WAT files, set the NxNyNz and DxDyDz terms
       call MannLibDims(BoxFileRoot, p%RotorDiamRef, p%WAT_NxNyNz, p%WAT_DxDyDz, ErrStat2, ErrMsg2);  if (Failed()) return
       write(sDummy, '(3(I8,1X))') p%WAT_NxNyNz
       call WrScr('  WAT: NxNyNz set to: '//trim(sDummy)//' (inferred from filename)')

diff --git a/glue-codes/openfast-cpp/src/OpenFAST.H b/glue-codes/openfast-cpp/src/OpenFAST.H
@@ -81,7 +81,7 @@ struct turbineDataType {
     int nBRfsiPtsTwr;
     //! The mean azimuth at which the loads are blended between AeroDyn and CFD
     double azBlendMean;
-    //! The delta azimuth over which the the loads are blended between AeroDyn and CFD
+    //! The delta azimuth over which the loads are blended between AeroDyn and CFD
     double azBlendDelta;
     //! Mean velocity at reference height
     double velMean;
@@ -708,7 +708,7 @@ private:
 
   //! Allocate memory for data structures for all turbines on this processor
   void allocateMemory_preInit();
-  //! Allocate more memory for each turbine after intialization/restart
+  //! Allocate more memory for each turbine after initialization/restart
   void allocateMemory_postInit(int iTurbLoc);
 
   //! Get the nacelle drag coefficient of local turbine number 'iTurbLoc'

diff --git a/modules/aerodyn/src/AeroDyn_Driver.f90 b/modules/aerodyn/src/AeroDyn_Driver.f90
@@ -26,8 +26,8 @@ program AeroDyn_Driver
    ! Program variables
    REAL(ReKi)                       :: PrevClockTime ! Clock time at start of simulation in seconds [(s)]
    REAL(ReKi)                       :: UsrTime1      ! User CPU time for simulation initialization [(s)]
-   REAL(ReKi)                       :: UsrTime2      ! User CPU time for simulation (without intialization) [(s)]
-   INTEGER(IntKi) , DIMENSION(1:8)  :: StrtTime      ! Start time of simulation (including intialization) [-]
+   REAL(ReKi)                       :: UsrTime2      ! User CPU time for simulation (without initialization) [(s)]
+   INTEGER(IntKi) , DIMENSION(1:8)  :: StrtTime      ! Start time of simulation (including initialization) [-]
    INTEGER(IntKi) , DIMENSION(1:8)  :: SimStrtTime   ! Start time of simulation (after initialization) [-]
    REAL(DbKi)                       :: t_global         ! global-loop time marker
    REAL(DbKi)                       :: t_final         ! global-loop time marker

diff --git a/modules/aerodyn/src/AeroDyn_Inflow_C_Binding.f90 b/modules/aerodyn/src/AeroDyn_Inflow_C_Binding.f90
@@ -422,10 +422,10 @@ SUBROUTINE ADI_C_Init( ADinputFilePassed, ADinputFileString_C, ADinputFileString
    ! Input file info
    integer(c_int),            intent(in   )  :: ADinputFilePassed                      !< Write VTK outputs [0: none, 1: init only, 2: animation]
    type(c_ptr),               intent(in   )  :: ADinputFileString_C                    !< Input file as a single string with lines deliniated by C_NULL_CHAR
-   integer(c_int),            intent(in   )  :: ADinputFileStringLength_C              !< lenght of the input file string
+   integer(c_int),            intent(in   )  :: ADinputFileStringLength_C              !< length of the input file string
    integer(c_int),            intent(in   )  :: IfWinputFilePassed                     !< Write VTK outputs [0: none, 1: init only, 2: animation]
    type(c_ptr),               intent(in   )  :: IfWinputFileString_C                   !< Input file as a single string with lines deliniated by C_NULL_CHAR
-   integer(c_int),            intent(in   )  :: IfWinputFileStringLength_C             !< lenght of the input file string
+   integer(c_int),            intent(in   )  :: IfWinputFileStringLength_C             !< length of the input file string
    character(kind=c_char),    intent(in   )  :: OutRootName_C(IntfStrLen)              !< Root name to use for echo files and other
    character(kind=c_char),    intent(in   )  :: OutVTKDir_C(IntfStrLen)                !< Directory to put all vtk output
    ! Environmental
@@ -2110,7 +2110,7 @@ end subroutine ADI_C_GetDiskAvgVel
 !===================================================================================================================================
 
 !> This routine is operating on module level data.  Error handling here in case checks added
-!! NOTE: the OriginInit is not included in the data passed in and must be added to the the position info here
+!! NOTE: the OriginInit is not included in the data passed in and must be added to the position info here
 subroutine Set_MotionMesh(iWT, ErrStat3, ErrMsg3)
    integer(IntKi),            intent(in   )  :: iWT         !< current rotor/turbine
    integer(IntKi),            intent(  out)  :: ErrStat3
@@ -2138,7 +2138,7 @@ end subroutine Set_MotionMesh
 
 !> Map the motion of the intermediate input mesh over to the input meshes
 !! This routine is operating on module level data, hence few inputs
-!! NOTE: the OriginInit is not included in the data passed in and must be added to the the position info here
+!! NOTE: the OriginInit is not included in the data passed in and must be added to the position info here
 subroutine AD_SetInputMotion( iWT, u_local,        &
          HubPos_C, HubOri_C, HubVel_C, HubAcc_C,   &
          NacPos_C, NacOri_C, NacVel_C, NacAcc_C,   &

diff --git a/modules/aerodyn/src/AirfoilInfo.f90 b/modules/aerodyn/src/AirfoilInfo.f90
@@ -1391,8 +1391,8 @@ SUBROUTINE ComputeUA360_AttachedFlow(p, ColUAf, cn_cl, iLower, iUpper)
       !------------------------------------------------
       call Compute_iLoweriUpper(p, iLower, iUpper)
 
-      p%UA_BL%alphaLower = p%alpha(iLower) ! note we are overwritting values here to make them consistent in the linear equation
-      p%UA_BL%alphaUpper = p%alpha(iUpper) ! note we are overwritting values here to make them consistent in the linear equation
+      p%UA_BL%alphaLower = p%alpha(iLower) ! note we are overwriting values here to make them consistent in the linear equation
+      p%UA_BL%alphaUpper = p%alpha(iUpper) ! note we are overwriting values here to make them consistent in the linear equation
 
       p%UA_BL%c_alphaLower = cn_cl(iLower) ! for vortex calculations (x5, HGMV model)
       p%UA_BL%c_alphaUpper = cn_cl(iUpper) ! for vortex calculations (x5, HGMV model)

diff --git a/modules/aerodyn/src/DBEMT.f90 b/modules/aerodyn/src/DBEMT.f90
@@ -844,7 +844,7 @@ SUBROUTINE DBEMT_ABM4( i, j, t, n, u, utimes, p, x, OtherState, m, ErrStat, ErrM
       x_in = x%element(i,j)
 
 
-         ! predict: (note that we are overwritting x%element(i,j)%vind and x%element(i,j)%vind_1 here):
+         ! predict: (note that we are overwriting x%element(i,j)%vind and x%element(i,j)%vind_1 here):
       CALL DBEMT_AB4( i, j, t, n, u, utimes, p, x, OtherState, m, ErrStat2, ErrMsg2 )
          CALL SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
          IF ( ErrStat >= AbortErrLev ) RETURN

diff --git a/modules/aerodyn/src/FVW.f90 b/modules/aerodyn/src/FVW.f90
@@ -480,7 +480,7 @@ subroutine FVW_FinalWrite(u, p, x, z, m, ErrStat, ErrMsg)
    ErrMsg  = ""
    ! Place any last minute operations or calculations here:
    if (p%WrVTK>0 .and. m%VTKstep<FINAL_STEP) then
-      call WrScr('OLAF: writting final VTK outputs')
+      call WrScr('OLAF: writing final VTK outputs')
       t=-1.0_ReKi
       if (p%WrVTK==1) then
          if (m%VTKstep<m%iStep+1) then
@@ -1529,7 +1529,7 @@ subroutine WriteVTKOutputs(t, force, VTKstep, u, p, x, z, m, ErrStat, ErrMsg)
             ! Global coordinate system, ALL VTK will be exported in global
             call WrVTK_FVW(p, x, z, m, trim(p%VTK_OutFileBase)//'FVW_Glb', VTKstep, 9, bladeFrame=.FALSE.)
          endif
-         m%VTKstep=VTKstep ! We save the step at which writing occured
+         m%VTKstep=VTKstep ! We save the step at which writing occurred
       endif
    endif
    ! --- Write VTK grids
@@ -1549,7 +1549,7 @@ subroutine WriteVTKOutputs(t, force, VTKstep, u, p, x, z, m, ErrStat, ErrMsg)
             call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
             m%GridOutputs(iGrid)%tLastOutput = t
             call WrVTK_FVW_Grid(p, m, iGrid, trim(p%VTK_OutFileBase)//'FVW_Grid', VTKstep, 9)
-            m%VTKstep=VTKstep ! We save the step at which writing occured
+            m%VTKstep=VTKstep ! We save the step at which writing occurred
          endif
       enddo
    endif

diff --git a/modules/aerodyn/src/FVW_Subs.f90 b/modules/aerodyn/src/FVW_Subs.f90
@@ -146,7 +146,7 @@ subroutine ReadAndInterpGamma(CirculationFileName, s_CP_LL, L, Gamma_CP_LL, ErrS
       read(iUnit,*, iostat=istat) sPrescr(i), GammaPrescr(i)
       if (istat/=0) then
          errStat2=ErrID_Fatal
-         errMsg2='Error occured while reading line '//num2lstr(i+1)//' of circulation file: '//trim(CirculationFileName)
+         errMsg2='Error occurred while reading line '//num2lstr(i+1)//' of circulation file: '//trim(CirculationFileName)
          if(Failed()) return
       endif
    enddo

diff --git a/modules/aerodyn/src/UnsteadyAero.f90 b/modules/aerodyn/src/UnsteadyAero.f90
@@ -3179,7 +3179,7 @@ SUBROUTINE UA_ABM4( i, j, t, n, u, utimes, p, x, OtherState, AFInfo, m, ErrStat,
       ! save copy of x(t):
       x_in     = x%element(i,j)
 
-         ! predict: (note that we are overwritting x%element(i,j) here):
+         ! predict: (note that we are overwriting x%element(i,j) here):
       CALL UA_AB4( i, j, t, n, u, utimes, p, x, OtherState, AFInfo, m, ErrStat2, ErrMsg2 )
          CALL SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
          IF ( ErrStat >= AbortErrLev ) RETURN

diff --git a/modules/beamdyn/src/BeamDyn.f90 b/modules/beamdyn/src/BeamDyn.f90
@@ -2282,7 +2282,7 @@ SUBROUTINE BD_QuadraturePointDataAt0( p )
 
    DO nelem = 1,p%elem_total
        DO idx_qp = 1,p%nqp
-            !> ### Calculate the the initial displacement fields in an element
+            !> ### Calculate the initial displacement fields in an element
             !! Initial displacement field \n
             !!    \f$   \underline{u_0}\left( \xi \right) =
             !!                \sum_{k=1}^{p+1} h^k\left( \xi \right) \underline{\hat{u}_0}^k
@@ -2370,7 +2370,7 @@ SUBROUTINE BD_DisplacementQP( nelem, p, x, m )
    INTEGER(IntKi)                :: idx_qp            !< index to the current quadrature point
    INTEGER(IntKi)                :: elem_start        !< Node point of first node in current element
 
-   !> ### Calculate the the displacement fields in an element
+   !> ### Calculate the displacement fields in an element
    !! Using equations (27) and (28) \n
    !!    \f$   \underline{u}\left( \xi \right) =
    !!                \sum_{i=1}^{p+1} h^i\left( \xi \right) \underline{\hat{u}}^i