Update HD template for OpenFAST#1008

OpenFAST/openfast#1008
bjonkman · Jun 22, 2022 · 6eae813 · 6eae813
1 parent a38a5bc
commit 6eae813
Show file tree

Hide file tree

Showing 3 changed files with 18 additions and 5 deletions.
diff --git a/ConvertFASTversions/InputConversions/newInputs_HD.m b/ConvertFASTversions/InputConversions/newInputs_HD.m
@@ -18,6 +18,18 @@
 
  % v3 Specification
 
+ if isfield(HDP,'AxCoefs')
+     nc = length(HDP.AxCoefs.Headers);
+     if ~any(strcmpi(HDP.AxCoefs.Headers,'AxFDMod'))
+         HDP.AxCoefs.Headers{nc+1} = 'AxFDMod';
+         HDP.AxCoefs.Headers{nc+2} = 'AxVnCOff';
+         HDP.AxCoefs.Headers{nc+3} = 'AxFDLoFSc';
+
+         HDP.AxCoefs.Table(:,nc+1) = 0;   %AxFDMod
+         HDP.AxCoefs.Table(:,nc+2) = 0.0; %AxVnCOff
+         HDP.AxCoefs.Table(:,nc+3) = 1.0; %AxFDLoFSc      
+     end
+ end
 
     %----------------------------------------------------------------------
     % Sort output channels:

diff --git a/ConvertFASTversions/TemplateFiles/HydroDyn.dat b/ConvertFASTversions/TemplateFiles/HydroDyn.dat
@@ -56,12 +56,13 @@ False            Echo           - Echo the input file data (flag)
              0             0             0             0             0             0
 ---------------------- STRIP THEORY OPTIONS --------------------------------------
              1   WaveDisp       - Method of computing Wave Kinematics {0: use undisplaced position, 1: use displaced position) } (switch)
+             0   AMMod          - Method of computing distributed added-mass force. (0: Only and always on nodes below SWL at the undisplaced position. 2: Up to the instantaneous free surface) [overwrite to 0 when WaveMod = 0 or 6 or when WaveStMod = 0 in SeaState]
 ---------------------- AXIAL COEFFICIENTS --------------------------------------
              2   NAxCoef        - Number of axial coefficients (-)
-AxCoefID  AxCd     AxCa     AxCp
-   (-)    (-)      (-)      (-)
-    1     0.00     0.00     1.00
-    2     9.60     0.00     1.00
+AxCoefID  AxCd     AxCa     AxCp    AxFDMod   AxVnCOff  AxFDLoFSc
+   (-)    (-)      (-)      (-)      (-)       (-)        (-)
+    1     0.00     0.00     1.00      0        0.0        1.0
+    2     9.60     0.00     1.00      0        0.0        1.0
 ---------------------- MEMBER JOINTS -------------------------------------------
             44   NJoints        - Number of joints (-)   [must be exactly 0 or at least 2]
 JointID   Jointxi     Jointyi     Jointzi  JointAxID   JointOvrlp   [JointOvrlp= 0: do nothing at joint, 1: eliminate overlaps by calculating super member]

diff --git a/MATLAB2FAST/Matlab2HD.m b/MATLAB2FAST/Matlab2HD.m
@@ -115,7 +115,7 @@ function Matlab2HD(HDPar,TemplateFile,OutputFilename,hdrLines)
                 disp(  ['WARNING: heave coefficients table not found in the HD data structure.'] );
                 printTable = true;
             else
-                frmt = ' %4i %8.2f %8.2f %8.2f';
+                frmt = ' %4i %8.2f %8.2f %8.2f %6i    %8.2f   %8.2f';
                 WriteFASTTable(line, fidIN, fidOUT, HDPar.AxCoefs, newline, frmt);
                 continue; %let's continue reading the template file            
             end