Added TurbSim to OpenFAST

CMakeLists.txt

@@ -79,6 +79,7 @@ set(FAST_MODULES_LOCAL
   extptfm
   openfoam
   supercontroller
+  turbsim
   fast-library
   )
 

modules-local/turbsim/CMakeLists.txt

@@ -0,0 +1,37 @@
+#
+# Copyright 2016 National Renewable Energy Laboratory
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#     http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+#
+
+set(MODULE_SOURCES
+  src/BlankModVKM.f90
+  src/CohStructures.f90
+  src/Profiles.f90
+  src/RandNum.f90
+  src/Root_Searching.f90
+  src/TS_FileIO.f90
+  src/TSsubs.f90
+  src/TurbSim.f90
+  src/TurbSim_Types.f90
+  src/VelocitySpectra.f90
+  )
+
+add_library(turbsimlib ${MODULE_SOURCES})
+target_link_libraries(turbsimlib nwtclibs)
+
+install(TARGETS turbsimlib
+  EXPORT "${CMAKE_PROJECT_NAME}Libraries"
+  RUNTIME DESTINATION bin
+  LIBRARY DESTINATION lib
+  ARCHIVE DESTINATION lib)

modules-local/turbsim/README.md

@@ -0,0 +1,15 @@
+# TurbSim
+A stochastic, full-field, turbulence simulator, primarialy for use with [InflowWind](https://nwtc.nrel.gov/InflowWind "InflowWind")-based simulation tools 
+
+**Authors**: Neil Kelley and [Bonnie Jonkman](mailto:[email protected]), NREL
+
+TurbSim is a stochastic, full-field, turbulent-wind simulator. It uses a statistical model (as opposed to a physics-based model) to 
+numerically simulate time series of three-component wind-speed vectors at points in a two-dimensional vertical rectangular 
+grid that is fixed in space. 
+
+Spectra of velocity components and spatial coherence are defined in the frequency domain, and 
+an inverse Fourier transform produces time series. The underlying theory behind this method of 
+simulating time series assumes a stationary process.
+
+For more information, please refer to documentation on the [TurbSim web site](https://nwtc.nrel.gov/TurbSim).
+

modules-local/turbsim/src/BlankModVKM.f90

@@ -0,0 +1,75 @@
+MODULE ModifiedvKrm_mod
+
+   USE NWTC_Library
+
+   IMPLICIT NONE
+
+CONTAINS
+!=======================================================================
+SUBROUTINE Mod_vKrm ( Ht, Ucmp, Spec )
+
+   ! This subroutine defines the "Improved" von Karman PSD model.
+   ! The use of this subroutine requires that all variables have the units of meters and seconds.
+   IMPLICIT                   NONE
+
+      !Passed variables
+
+   REAL(ReKi), INTENT(IN)    :: Ht                   ! height
+   REAL(ReKi), INTENT(IN)    :: UCmp                 ! wind speed
+   REAL(ReKi), INTENT(  OUT) :: Spec   (:,:)
+
+
+RETURN
+END SUBROUTINE Mod_vKrm
+!=======================================================================
+SUBROUTINE ScaleMODVKM(Ht,UCmp, LambdaU, LambdaV, LambdaW)
+
+! THIS SUBROUTINE DEFINES HUB SCALE PARMS FOR Modified von KARMAN PSD MODEL 
+   IMPLICIT NONE
+
+   REAL(ReKi), INTENT(IN)  :: Ht     ! height
+   REAL(ReKi), INTENT(IN)  :: UCmp   ! wind speed
+
+   REAL(ReKi)              :: LambdaU
+   REAL(ReKi)              :: LambdaV
+   REAL(ReKi)              :: LambdaW
+
+RETURN
+END SUBROUTINE ScaleMODVKM
+!=======================================================================
+FUNCTION FindZ0(z, sigma, U, f)
+
+    ! This function is used in the Modified von Karman model to 
+    ! determine the necessary surface roughness length for a given sigma.
+   IMPLICIT                NONE
+
+   REAL(ReKi)            :: FindZ0        ! Estimated surface roughness length
+   REAL(ReKi),INTENT(IN) :: z             ! Hub height
+   REAL(ReKi),INTENT(IN) :: sigma         ! Target std deviation
+   REAL(ReKi),INTENT(IN) :: U             ! Hub height wind speed
+   REAL(ReKi),INTENT(IN) :: f             ! Coriolis parameter
+
+   FindZ0 = 1.0   ! a default value
+
+RETURN
+END FUNCTION FindZ0
+!=======================================================================
+FUNCTION CalcDiff(z0Guess, z, sigma, U, f)
+
+      ! This function calculates the difference between the specified
+      ! sigma and the calculated one.
+   IMPLICIT                NONE
+
+   REAL(ReKi)             :: CalcDiff  ! Output - will be nearly zero if surface roughness is correct
+   REAL(ReKi), INTENT(IN) :: z0Guess   ! estimated surface roughness
+   REAL(ReKi), INTENT(IN) :: z         ! Hub height (m)
+   REAL(ReKi), INTENT(IN) :: sigma     ! Target standard deviation (m/s)
+   REAL(ReKi), INTENT(IN) :: U         ! Mean hub-height wind speed (m/s)
+   REAL(ReKi), INTENT(IN) :: f         ! Coriolis parameter 
+
+   CalcDiff = 0.0
+
+RETURN
+END FUNCTION CalcDiff
+!=======================================================================
+END MODULE ModifiedvKrm_mod