IfW: add VTK output of slice in XY to driver #1643
Conversation
This allows visualizing a slice through hub-height (or other height) when compairing with a simulation.
|
|
||
| enddo | ||
|
|
||
| close (unit) |
There was a problem hiding this comment.
Should close (unit) be inside the time step do loop?
There was a problem hiding this comment.
yes it should! I had removed the file loop, then re-added it.
| character(*), intent(in) :: FileRootName !< RootName for output files | ||
| real(ReKi), intent(in) :: XYslice_height | ||
| integer(IntKi), intent(in) :: unit !< Error status of the operation | ||
| integer(IntKi), intent(out) :: ErrStat !< Error status of the operation |
There was a problem hiding this comment.
Are ErrStat and ErrMsg initialized in this function? Maybe I missed it.
There was a problem hiding this comment.
They aren't explicitly initialized, but get set either by OpenFOutFile or by the check on the slice height. For completeness I can add them.
I also just noticed that there is no upper bounds check on the height, which could lead to an invalid index on iz. I'll fix that also.
|
If it's not too much work, maybe XZ and YZ could be added too (like fast farm). On a side note, plotting slices of TurbSim and Mann boxes can easily be done using the python-toolbox: |
|
Adding the XY and YZ slices should be pretty straightforward. I'll give that a shot. Also, it might make sense to allow exporting multiple slices instead of just one. |
andrew-platt
force-pushed
the
f/ifw_driver_XYslice
branch
from
efa2dc6
to
82eec99
Compare
For generating InflowWind_Driver_Types.f90
andrew-platt
force-pushed
the
f/ifw_driver_XYslice
branch
from
da194da
to
cecb3e1
Compare
This PR is ready to merge.
Feature or improvement description
This allows visualizing a slice through hub-height (or other height) when comparing with a simulation. The user can specify the height at which to export the slice. At the moment it is a very slow process due to the looping and because it exports over the entire time range (this could be improved later).
This figure shows the XY plane at the hub height for a Mann wind field (HAWC format). The vertical plane is exported from the driver using the
WrVTKoption. Both are taken at 23.625 seconds.Related issue, if one exists
none
Impacted areas of the software
InflowWind driver only
Additional supporting information
During testing of the wake added turbulence feature, I ran into issues in verifying that the background turbulence was being added correctly. To check this, I needed some way to visualize the XY plane through the hub from the Mann box turbulence.
Test results, if applicable
No test cases are affected. I did not add a test case for this since it is a very slow code to execute.