update waveSpread calculation in irregular waves (WEC-Sim#1290)

* fix on pDis function call * preliminary fix for issue WEC-Sim#1288
kmruehl · Aug 9, 2024 · 0520f46 · 0520f46
1 parent d9bd629
commit 0520f46
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/source/functions/simulink/model/irregExcF.m b/source/functions/simulink/model/irregExcF.m
@@ -8,7 +8,7 @@
     B1= sin(bsxfun(@plus,A1,phaseRand(:,ii)));
     B11 = sin(bsxfun(@plus,w*time,phaseRand(:,ii)));
     C0 = bsxfun(@times,A*spread(ii),dw);
-    C1 = sqrt(bsxfun(@times,A*spread(ii),dw));
+    C1 = sqrt(bsxfun(@times,A*spread(ii).^2,dw));
     D0 =bsxfun(@times,squeeze(fExtMD(ii,:,:)),C0);
     D1 =bsxfun(@times,squeeze(fExtRE(ii,:,:)),C1);
     D11 = bsxfun(@times,squeeze(fExtIM(ii,:,:)),C1);

diff --git a/source/functions/simulink/model/irregnLYaw.m b/source/functions/simulink/model/irregnLYaw.m
@@ -75,7 +75,7 @@
     B1= sin(bsxfun(@plus,A1,phaseRandint));
     B11 = sin(bsxfun(@plus,w*time,phaseRandint));
     C0 = bsxfun(@times,A*WaveSpreadint,dw);
-    C1 = sqrt(bsxfun(@times,A*WaveSpreadint,dw));
+    C1 = sqrt(bsxfun(@times,A*WaveSpreadint.^2,dw));
     D0 =bsxfun(@times,fExtMDint,C0);
     D1 =bsxfun(@times,fExtREint,C1);
     D11 = bsxfun(@times,fExtIMint,C1);