[HydroDyn] Error when reading QTF file .12d for multiple wave headings

[getChaos22]

Hi all,

When I use QTF file .12d with multiple wave headings, OpenFAST will give a error message which is:

Part of my .12d file is like:

I noticed that on line 20, the real part and imaginary part (column 8 and 9) are truly different from those on line 8, whose frequency and directions are same with line 20.

However, even though I changed column 6 to 9 of line 20 to the same with line 8, the error still exist. Only when I make line 20 totally same with line 8 (as is shown in the following screenshot), the error of line 20 will disappear, and the error of line 21 came out.

The QTF files I used were not directly obtained by Wamit. Instead, I employed Nemoh to get the QTF results and used BEMRosetta to output Wamit format QTF files. So I'm not sure whether the format of the QTF files I obtained is completely right. Could you give me some information on how the Wamit QTF files (12d and 12s) are like for multiple wave headings?

Thank you.

Chao Li

[getChaos22]

Hi all,

Even though I still haven't figured out how a correct QTF file with multiple wave directions is like, I do find a place that is not so make sense in the source code of WAMIT2.f90 (Line 4252 ~ Line 4257):

I tried to modify this place in a way that I think maybe more reasonable:

By doing this, the program can run normally without error.

I would like to know your ideas. Is this the right way to solve this issue?

Thanks

Chao Li

[luwang00]

Hi @getChaos22, what you have here does look like a bug. In your original example, I think what happened was that when reading in line 8, OpenFAST tried to automatically populate the entry associated with Line 20. However, it copied the wrong value as shown in your code change. The wrong value triggered an error of miss-matching values when OpenFAST later tried to read line 20. We should get this fixed. Thanks for finding this. Hi @andrew-platt, can you take a quick look at this as well? Thanks!

[andrew-platt]

Hi @getChaos22, that looks like the correct solution to the issue! I suspect I had a copy/paste error there when I wrote that part of the code 10 years ago. Do you want to make a pull request for this fix?

[getChaos22]

Hi @andrew-platt
Sure, which branch do you want me to make the pull request?

[andrew-platt]

@getChaos22, thanks for being willing to do that. Start with the dev branch and make a pull request there. Thanks!

[getChaos22]

No problem @andrew-platt, please check the pull request and see whether it is ready to merge.