Skip to content

[BugFix] TurbSim calculating grid bottom location #2808

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
andrew-platt merged 4 commits into from
May 22, 2025
Merged

[BugFix] TurbSim calculating grid bottom location #2808

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
9f46214
compute Zbottom from GridHeight, not RotorDiameter
ptrbortolotti May 13, 2025
31918a5
Zbottom is simply HubHt - half GridHeight. skip Ztop
ptrbortolotti May 13, 2025
e676fb5
TurbSim: if grid goes below ground, shift the bottom point
bjonkman May 14, 2025
f1a8ef4
Merge pull request #5 from bjonkman/f/TurbSimGrid
ptrbortolotti May 14, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 3 additions & 1 deletion modules/turbsim/src/TS_FileIO.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4868,7 +4868,9 @@ SUBROUTINE GetDefaultRS( p, OtherSt_RandNum, TmpUstarHub, ErrStat, ErrMsg )


Z(2) = p%grid%HubHt + 0.5*p%grid%RotorDiameter ! top of the grid
Z(1) = Z(2) - p%grid%GridHeight ! bottom of the grid
Z(1) = MAX( Tolerance, Z(2) - p%grid%GridHeight ) ! bottom of the grid
Z(2) = Z(1) + p%grid%GridHeight ! re-calculate just in case Z1 is set to Tolerance

CALL getVelocityProfile(p, p%UHub, p%grid%HubHt, Z, V, ErrStat2, ErrMsg2)
CALL SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, 'GetDefaultRS')

Expand Down
8 changes: 1 addition & 7 deletions modules/turbsim/src/TSsubs.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1349,13 +1349,7 @@ SUBROUTINE CreateGrid( p_grid, p_usr, UHub, AddTower, ErrStat, ErrMsg )

p_grid%Zbottom = p_grid%HubHt + 0.5*p_grid%RotorDiameter ! height of the highest grid points
p_grid%Zbottom = p_grid%Zbottom - p_grid%GridRes_Z * REAL(p_grid%NumGrid_Z - 1, ReKi) ! height of the lowest grid points

IF ( p_grid%Zbottom <= 0.0_ReKi ) THEN
CALL SetErrStat(ErrID_Fatal,'The lowest grid point ('//TRIM(Num2LStr(p_grid%Zbottom))// ' m) must be above the ground. '//&
'Adjust the appropriate values in the input file.',ErrStat,ErrMsg,RoutineName)
RETURN
ENDIF

p_grid%Zbottom = MAX( Tolerance, p_grid%Zbottom) ! make sure it's above the ground

! (2) the tower points:
IF ( AddTower ) THEN
Expand Down