ADI_c: set OutputAccel flag for IfW with MHK (#27)

* ADI_c: set `OutputAccel` flag for IfW with MHK

* ADI cbind: pass the OutputAccel flag through AeroDyn_Inflow.f90

* AD driver: set OutputAccel flag for IfW

* Remove hard coded OutputAccel flag in favor of input file

---------

Co-authored-by: Rafael M Mudafort <[email protected]>

Author: andrew-platt

modules/aerodyn/src/AeroDyn_Driver_Subs.f90

@@ -477,6 +477,7 @@ subroutine Init_ADI_ForDriver(iCase, ADI, dvr, FED, dt, needInitIW, errStat, err
       InitInp%IW_InitInp%RefHt      = dvr%IW_InitInp%RefHt
       InitInp%IW_InitInp%PLExp      = dvr%IW_InitInp%PLExp
       InitInp%IW_InitInp%MHK        = dvr%MHK
+      InitInp%IW_InitInp%OutputAccel= dvr%MHK /= MHK_None
       InitInp%IW_InitInp%FilePassingMethod = 0_IntKi ! read input file instead of passed file data
       ! AeroDyn
       InitInp%AD%Gravity   = 9.80665_ReKi

modules/aerodyn/src/AeroDyn_Inflow.f90

@@ -388,11 +388,11 @@ subroutine ADI_InitInflowWind(Root, i_IW, u_AD, o_AD, IW, dt, InitOutData, errSt
       endif
       InitInData%RootName         = trim(Root)//'.IfW'
       InitInData%MHK              = i_IW%MHK
-      InitInData%OutputAccel      = .TRUE.            ! TODO: Set via input file or c driver flag
       ! OLAF might be used in AD, in which case we need to allow out of bounds for some calcs. To do that
       ! the average values for the entire wind profile must be calculated and stored (we don't know if OLAF
       ! is used until after AD_Init below).
       InitInData%BoxExceedAllow = .true.
+      InitInData%OutputAccel = i_IW%OutputAccel
       CALL InflowWind_Init( InitInData, IW%u, IW%p, &
                      IW%x, IW%xd, IW%z, IW%OtherSt, &
                      IW%y, IW%m, dt,  InitOutData, errStat2, errMsg2 )

modules/aerodyn/src/AeroDyn_Inflow_C_Binding.f90

@@ -386,6 +386,7 @@ subroutine ADI_C_PreInit(                       &
    ! Set whether these are MHK turbines
    InitInp%AD%MHK = MHK_in
    InitInp%IW_InitInp%MHK = MHK_in
+   InitInp%IW_InitInp%OutputAccel      = MHK_in /= MHK_None
 
    ! Offset the origin to account for water depth for MHK turbines
    do iWT=1,Sim%NumTurbines

modules/aerodyn/src/AeroDyn_Inflow_Registry.txt

@@ -41,6 +41,8 @@ typedef    ^        ^                  IntKi                            FilePass
 typedef    ^        ^                  FileInfoType                     PassedFileInfo   -      -     -   "If we don't use the input file, pass everything through this as a FileInfo structure"       -
 typedef    ^        ^                  InflowWind_InputFile             PassedFileData   -      -     -   "If we don't use the input file, pass everything through this as an IfW InputFile structure" -
 typedef    ^        ^                  LOGICAL                          Linearize        -    .FALSE. -   "Flag that tells this module if the glue code wants to linearize." -
+typedef    ^        ^                  LOGICAL                          OutputAccel      -    .FALSE. -   "Flag to output wind acceleration" -
+
 
 
 # ..... InitIn ....................................................................................................................

modules/aerodyn/src/AeroDyn_Inflow_Types.f90

@@ -63,6 +63,7 @@ MODULE AeroDyn_Inflow_Types
     TYPE(FileInfoType)  :: PassedFileInfo      !< If we don't use the input file, pass everything through this as a FileInfo structure [-]
     TYPE(InflowWind_InputFile)  :: PassedFileData      !< If we don't use the input file, pass everything through this as an IfW InputFile structure [-]
     LOGICAL  :: Linearize = .FALSE.      !< Flag that tells this module if the glue code wants to linearize. [-]
+    LOGICAL  :: OutputAccel = .FALSE.      !< Flag to output wind acceleration [-]
   END TYPE ADI_IW_InputData
 ! =======================
 ! =========  ADI_InitInputType  =======
@@ -309,6 +310,7 @@ subroutine ADI_CopyIW_InputData(SrcIW_InputDataData, DstIW_InputDataData, CtrlCo
    call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
    if (ErrStat >= AbortErrLev) return
    DstIW_InputDataData%Linearize = SrcIW_InputDataData%Linearize
+   DstIW_InputDataData%OutputAccel = SrcIW_InputDataData%OutputAccel
 end subroutine
 
 subroutine ADI_DestroyIW_InputData(IW_InputDataData, ErrStat, ErrMsg)
@@ -341,6 +343,7 @@ subroutine ADI_PackIW_InputData(RF, Indata)
    call NWTC_Library_PackFileInfoType(RF, InData%PassedFileInfo) 
    call InflowWind_PackInputFile(RF, InData%PassedFileData) 
    call RegPack(RF, InData%Linearize)
+   call RegPack(RF, InData%OutputAccel)
    if (RegCheckErr(RF, RoutineName)) return
 end subroutine
 
@@ -359,6 +362,7 @@ subroutine ADI_UnPackIW_InputData(RF, OutData)
    call NWTC_Library_UnpackFileInfoType(RF, OutData%PassedFileInfo) ! PassedFileInfo 
    call InflowWind_UnpackInputFile(RF, OutData%PassedFileData) ! PassedFileData 
    call RegUnpack(RF, OutData%Linearize); if (RegCheckErr(RF, RoutineName)) return
+   call RegUnpack(RF, OutData%OutputAccel); if (RegCheckErr(RF, RoutineName)) return
 end subroutine
 
 subroutine ADI_CopyInitInput(SrcInitInputData, DstInitInputData, CtrlCode, ErrStat, ErrMsg)