Skip to content
This repository was archived by the owner on Nov 27, 2023. It is now read-only.

Commit c3649de

Browse files
ebranlardebranlard
committed
IO: update to read moordyn v2 (Fixes #54)
1 parent 4bc73f7 commit c3649de

File tree

4 files changed

+132
-16
lines changed

4 files changed

+132
-16
lines changed

‎pyFAST/input_output/fast_input_file.py

Lines changed: 62 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -351,12 +351,19 @@ def _read(self):
351351

352352
# --- Tables that can be detected based on the "Label" (second entry on line)
353353
# NOTE: MJointID1, used by SubDyn and HydroDyn
354-
NUMTAB_FROM_LAB_DETECT = ['NumAlf' , 'F_X' , 'MemberCd1' , 'MJointID1' , 'NOutLoc' , 'NOutCnt' , 'PropD' ,'Diam' ,'Type' ,'LineType' ]
355-
NUMTAB_FROM_LAB_DIM_VAR = ['NumAlf' , 'NKInpSt' , 'NCoefMembers' , 'NMembers' , 'NMOutputs' , 'NMOutputs' , 'NPropSets' ,'NTypes' ,'NConnects' ,'NLines' ]
356-
NUMTAB_FROM_LAB_VARNAME = ['AFCoeff' , 'TMDspProp' , 'MemberProp' , 'Members' , 'MemberOuts' , 'MemberOuts' , 'SectionProp' ,'LineTypes' ,'ConnectionProp' ,'LineProp' ]
357-
NUMTAB_FROM_LAB_NHEADER = [2 , 2 , 2 , 2 , 2 , 2 , 2 , 2 , 2 , 2 ]
358-
NUMTAB_FROM_LAB_NOFFSET = [0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ]
359-
NUMTAB_FROM_LAB_TYPE = ['num' , 'num' , 'num' , 'mix' , 'num' , 'sdout' , 'num' ,'mix' ,'mix' ,'mix' ]
354+
NUMTAB_FROM_LAB_DETECT = ['NumAlf' , 'F_X' , 'MemberCd1' , 'MJointID1' , 'NOutLoc' , 'NOutCnt' , 'PropD' ]
355+
NUMTAB_FROM_LAB_DIM_VAR = ['NumAlf' , 'NKInpSt' , 'NCoefMembers' , 'NMembers' , 'NMOutputs' , 'NMOutputs' , 'NPropSets' ]
356+
NUMTAB_FROM_LAB_VARNAME = ['AFCoeff' , 'TMDspProp' , 'MemberProp' , 'Members' , 'MemberOuts' , 'MemberOuts' , 'SectionProp' ]
357+
NUMTAB_FROM_LAB_NHEADER = [2 , 2 , 2 , 2 , 2 , 2 , 2 ]
358+
NUMTAB_FROM_LAB_NOFFSET = [0 , 0 , 0 , 0 , 0 , 0 , 0 ]
359+
NUMTAB_FROM_LAB_TYPE = ['num' , 'num' , 'num' , 'mix' , 'num' , 'sdout' , 'num' ]
360+
# MoorDyn Version 1 and 2 (with AUTO for LAB_DIM_VAR)
361+
NUMTAB_FROM_LAB_DETECT += ['Diam' ,'Type' ,'LineType' , 'Attachment']
362+
NUMTAB_FROM_LAB_DIM_VAR += ['NTypes:AUTO','NConnects' ,'NLines:AUTO' , 'AUTO']
363+
NUMTAB_FROM_LAB_VARNAME += ['LineTypes' ,'ConnectionProp' ,'LineProp' , 'Points']
364+
NUMTAB_FROM_LAB_NHEADER += [ 2 , 2 , 2 , 2 ]
365+
NUMTAB_FROM_LAB_NOFFSET += [ 0 , 0 , 0 , 0 ]
366+
NUMTAB_FROM_LAB_TYPE += ['mix' ,'mix' ,'mix' , 'mix']
360367
# SubDyn
361368
NUMTAB_FROM_LAB_DETECT += ['GuyanDampSize' , 'YoungE' , 'YoungE' , 'EA' , 'MatDens' ]
362369
NUMTAB_FROM_LAB_DIM_VAR += [6 , 'NPropSets', 'NXPropSets', 'NCablePropSets' , 'NRigidPropSets']
@@ -480,6 +487,19 @@ def _read(self):
480487
i+=1;
481488
self.readBeamDynProps(lines,i)
482489
return
490+
elif line.upper().find('OUTPUTS')>0:
491+
if 'Points' in self.keys() and 'dtM' in self.keys():
492+
OutList,i = parseFASTOutList(lines,i+1)
493+
d = getDict()
494+
d['label'] = 'Outlist'
495+
d['descr'] = ''
496+
d['tabType'] = TABTYPE_FIL # TODO
497+
d['value'] = OutList
498+
self.addComment('------------------------ OUTPUTS --------------------------------------------')
499+
self.data.append(d)
500+
self.addComment('END')
501+
self.addComment('------------------------- need this line --------------------------------------')
502+
return
483503

484504
# --- Parsing of standard lines: value(s) key comment
485505
line = lines[i]
@@ -597,7 +617,6 @@ def _read(self):
597617
self.data.append(dd)
598618

599619
d['label'] = NUMTAB_FROM_LAB_VARNAME[ii]
600-
d['tabDimVar'] = NUMTAB_FROM_LAB_DIM_VAR[ii]
601620
if d['label'].lower()=='afcoeff' :
602621
d['tabType'] = TABTYPE_NUM_WITH_HEADERCOM
603622
else:
@@ -607,10 +626,28 @@ def _read(self):
607626
d['tabType'] = TABTYPE_NUM_SUBDYNOUT
608627
else:
609628
d['tabType'] = TABTYPE_MIX_WITH_HEADER
610-
if isinstance(d['tabDimVar'],int):
629+
# Finding table dimension (number of lines)
630+
tabDimVar = NUMTAB_FROM_LAB_DIM_VAR[ii]
631+
if isinstance(tabDimVar, int): # dimension hardcoded
632+
d['tabDimVar'] = tabDimVar
611633
nTabLines = d['tabDimVar']
612634
else:
613-
nTabLines = self[d['tabDimVar']+labOffset]
635+
# We either use a variable name or "AUTO" to find the number of rows
636+
tabDimVars = tabDimVar.split(':')
637+
for tabDimVar in tabDimVars:
638+
d['tabDimVar'] = tabDimVar
639+
if tabDimVar=='AUTO':
640+
# Determine table dimension automatically
641+
nTabLines = findNumberOfTableLines(lines[i+nHeaders:], break_chars=['---','!','#'])
642+
break
643+
else:
644+
try:
645+
nTabLines = self[tabDimVar+labOffset]
646+
break
647+
except KeyError:
648+
#print('Cannot determine table dimension using {}'.format(tabDimVar))
649+
# Hopefully this table has AUTO as well
650+
pass
614651

615652
d['label'] += labOffset
616653
#print('Reading table {} Dimension {} (based on {})'.format(d['label'],nTabLines,d['tabDimVar']));
@@ -755,10 +792,13 @@ def mat_tostring(M,fmt='24.16e'):
755792
s+='\n'.join('\t'.join('{:15.8e}'.format(x) for x in y) for y in d['value'])
756793
elif d['tabType']==TABTYPE_FIL:
757794
#f.write('{} {} {}\n'.format(d['value'][0],d['tabDetect'],d['descr']))
795+
label = d['label']
796+
if 'kbot' in self.keys(): # Moordyn has no 'OutList' label..
797+
label=''
758798
if len(d['value'])==1:
759-
s+='{} {} {}'.format(d['value'][0],d['label'],d['descr']) # TODO?
799+
s+='{} {} {}'.format(d['value'][0], label, d['descr']) # TODO?
760800
else:
761-
s+='{} {} {}\n'.format(d['value'][0],d['label'],d['descr']) # TODO?
801+
s+='{} {} {}\n'.format(d['value'][0], label, d['descr']) # TODO?
762802
s+='\n'.join(fil for fil in d['value'][1:])
763803
elif d['tabType']==TABTYPE_NUM_BEAMDYN:
764804
# TODO use dedicated sub-class
@@ -1183,6 +1223,17 @@ def detectUnits(s,nRef):
11831223
return Units
11841224

11851225

1226+
def findNumberOfTableLines(lines, break_chars):
1227+
""" Loop through lines until a one of the "break character is found"""
1228+
for i, l in enumerate(lines):
1229+
for bc in break_chars:
1230+
if l.startswith(bc):
1231+
return i
1232+
# Not found
1233+
print('[FAIL] end of table not found')
1234+
return len(lines)
1235+
1236+
11861237
def parseFASTNumTable(filename,lines,n,iStart,nHeaders=2,tableType='num',nOffset=0, varNumLines=''):
11871238
"""
11881239
First lines of data starts at: nHeaders+nOffset

‎pyFAST/input_output/tests/example_files/FASTIn_MD.dat renamed to ‎pyFAST/input_output/tests/example_files/FASTIn_MD-v1.dat

File renamed without changes.

‎pyFAST/input_output/tests/example_files/FASTIn_MD-v2.dat

Lines changed: 39 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,39 @@
1+
--------------------- MoorDyn Input File ------------------------------------
2+
Mooring system for OC4-DeepCwind Semi
3+
FALSE Echo - echo the input file data (flag)
4+
----------------------- LINE TYPES ------------------------------------------
5+
Name Diam MassDen EA BA/-zeta EI Cd Ca CdAx CaAx
6+
(-) (m) (kg/m) (N) (N-s/-) (-) (-) (-) (-) (-)
7+
main 0.0766 113.35 7.536E8 -1.0 0 2.0 0.8 0.4 0.25
8+
---------------------- POINTS --------------------------------
9+
ID Attachment X Y Z M V CdA CA
10+
(-) (-) (m) (m) (m) (kg) (m^3) (m^2) (-)
11+
1 Fixed 418.8 725.383 -200.0 0 0 0 0
12+
2 Fixed -837.6 0.0 -200.0 0 0 0 0
13+
3 Fixed 418.8 -725.383 -200.0 0 0 0 0
14+
4 Vessel 20.434 35.393 -14.0 0 0 0 0
15+
5 Vessel -40.868 0.0 -14.0 0 0 0 0
16+
6 Vessel 20.434 -35.393 -14.0 0 0 0 0
17+
---------------------- LINES --------------------------------------
18+
ID LineType AttachA AttachB UnstrLen NumSegs Outputs
19+
(-) (-) (-) (-) (m) (-) (-)
20+
1 main 1 4 835.35 20 -
21+
2 main 2 5 835.35 20 -
22+
3 main 3 6 835.35 20 -
23+
---------------------- SOLVER OPTIONS ---------------------------------------
24+
0.001 dtM - time step to use in mooring integration (s)
25+
3.0e6 kbot - bottom stiffness (Pa/m)
26+
3.0e5 cbot - bottom damping (Pa-s/m)
27+
2.0 dtIC - time interval for analyzing convergence during IC gen (s)
28+
60.0 TmaxIC - max time for ic gen (s)
29+
4.0 CdScaleIC - factor by which to scale drag coefficients during dynamic relaxation (-)
30+
0.01 threshIC - threshold for IC convergence (-)
31+
------------------------ OUTPUTS --------------------------------------------
32+
FairTen1
33+
FairTen2
34+
FairTen3
35+
AnchTen1
36+
AnchTen2
37+
AnchTen3
38+
END
39+
------------------------- need this line --------------------------------------

‎pyFAST/input_output/tests/test_fast_input.py

Lines changed: 31 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -2,7 +2,7 @@
22
import os
33
import numpy as np
44

5-
from .helpers_for_test import MyDir, reading_test
5+
from pyFAST.input_output.tests.helpers_for_test import MyDir, reading_test
66

77
import pyFAST
88
from pyFAST.input_output import FASTInputFile
@@ -67,10 +67,6 @@ def test_FASTIn(self):
6767
F.test_ascii(bCompareWritesOnly=True,bDelete=True)
6868
self.assertEqual(F['PitManRat(1)'],2)
6969

70-
F=FASTInputFile(os.path.join(MyDir,'FASTIn_MD.dat'))
71-
F.test_ascii(bCompareWritesOnly=True,bDelete=True)
72-
self.assertEqual(float(F['LineTypes'][0,1]),0.02)
73-
7470
def test_FASTWnd(self):
7571
F=FASTWndFile(os.path.join(MyDir,'FASTWnd.wnd'))
7672
F.test_ascii(bCompareWritesOnly=True,bDelete=True)
@@ -87,7 +83,37 @@ def test_FASTInGraph(self):
8783
#graph = F.toGraph()
8884
# self.assertEqual(len(graph.Nodes), 2)
8985
# self.assertEqual(len(graph.Elements), 1)
86+
def test_FASTInMoorDyn(self):
87+
# MoorDyn version 1
88+
F=FASTInputFile(os.path.join(MyDir,'FASTIn_MD-v1.dat'))
89+
F.test_ascii(bCompareWritesOnly=True,bDelete=True)
90+
self.assertEqual(float(F['LineTypes'][0,1]),0.02)
91+
92+
# MoorDyn version 2
93+
F=FASTInputFile(os.path.join(MyDir,'FASTIn_MD-v2.dat'))
94+
#F.write(os.path.join(MyDir,'FASTIn_MD-v2.dat---OUT'))
95+
self.assertTrue('Points' in F.keys())
96+
self.assertTrue('LineTypes' in F.keys())
97+
self.assertTrue('LineProp' in F.keys())
98+
self.assertEqual(F['LineProp'].shape , (3,7))
99+
self.assertEqual(F['LineTypes'].shape , (1,10))
100+
self.assertEqual(F['Points'].shape , (6,9))
101+
self.assertEqual(len(F['Outlist']) , 6)
102+
self.assertEqual(F['Outlist'][0] , 'FairTen1')
103+
self.assertEqual(F['LineProp'][0,0] , '1')
104+
self.assertEqual(F['LineProp'][0,1] , 'main')
105+
self.assertEqual(F['LineProp'][0,6] , '-')
106+
107+
def test_FASTInAirfoil(self):
108+
F=FASTInputFile(os.path.join(MyDir,'FASTIn_AD15_arfl.dat'))
109+
F.test_ascii(bCompareWritesOnly=True,bDelete=True)
110+
self.assertTrue('InterpOrd' in F.keys())
111+
self.assertTrue('AFCoeff' in F.keys())
112+
self.assertEqual(F['AFCoeff'].shape, (30,4))
90113

91114
if __name__ == '__main__':
115+
from welib.tools.clean_exceptions import *
92116
#Test().test_FASTIn()
117+
#Test().test_FASTInAirfoil()
118+
#Test().test_FASTInMoorDyn()
93119
unittest.main()

0 commit comments

Comments
 (0)