Fixed a bug on line 1265 of MoorDyn_Line.f90 that lead to bending mom…

…ent for rod EndB to be computed incorrectly. This fix makes the fortran version of MoorDyn v2 consistent with the C++ version.
OpenFAST · Sep 18, 2023 · 56c63e8 · 56c63e8
1 parent 357d827
commit 56c63e8
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/modules/moordyn/src/MoorDyn_Rod.f90 b/modules/moordyn/src/MoorDyn_Rod.f90
@@ -827,7 +827,7 @@ SUBROUTINE Rod_DoRHS(Rod, m, p)
             Rod%Mext = Rod%Mext + (/ Mtemp*sinBeta, -Mtemp*cosBeta, 0.0_DbKi /) 
 
             ! axial drag
-            Rod%Dq(:,I) = Rod%Dq(:,I) + VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq
+            Rod%Dq(:,I) = Rod%Dq(:,I) + 0.5 * VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq
 
             ! >>> what about rotational drag?? <<<   eqn will be  Pi* Rod%d**4/16.0 omega_rel?^2...  *0.5 * Cd...
 
@@ -857,7 +857,7 @@ SUBROUTINE Rod_DoRHS(Rod, m, p)
             Rod%Mext = Rod%Mext + (/ Mtemp*sinBeta, -Mtemp*cosBeta, 0.0_DbKi /) 
 
             ! axial drag
-            Rod%Dq(:,I) = Rod%Dq(:,I) + VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq
+            Rod%Dq(:,I) = Rod%Dq(:,I) + 0.5 * VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq
 
             ! long-wave diffraction force
             Rod%Aq(:,I) = Rod%Aq(:,I) + VOF * p%rhoW* Rod%CaEnd * (2.0/3.0*Pi*Rod%d**3 /8.0) * aq