From 21c7316a83e5e5bdfda385eb7a939e2fd9b23268 Mon Sep 17 00:00:00 2001 From: Fabrice Ardhuin Date: Thu, 5 Mar 2020 15:57:46 -0800 Subject: [PATCH] Bug correction: MSSSUM was not computed in some cases --- model/ftn/w3src4md.ftn | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/model/ftn/w3src4md.ftn b/model/ftn/w3src4md.ftn index 865ed05cb..75894dc09 100755 --- a/model/ftn/w3src4md.ftn +++ b/model/ftn/w3src4md.ftn @@ -1725,7 +1725,7 @@ ! ! 1.b MSS parameters used for Modulation factors for B or lambda ! - IF (SSDSC(8).GT.0.OR.SSDSC(11).GT.0) THEN + IF (SSDSC(8).GT.0.OR.SSDSC(11).GT.0.OR.SSDSC(18).GT.0) THEN DO IK=1,NK IMSSMAX (IK) = 1 MSSP = 0. @@ -2093,7 +2093,6 @@ ! C=SIG(IK)/K(IK) BTH0(IK)=sum(BTH(1:NTH)*DTH) - IF (SSDSC(18).GT.0) THEN ! Applies modulation factor on Lambda DO ITH=1,NTH IF (SSDSC(11).GT.0) THEN @@ -2105,6 +2104,7 @@ (ECOS(ITH)*COS(DLWMEAN)+ESIN(ITH)*SIN(DLWMEAN))**2) ENDIF LMODULATION(ITH)= FACSTRAINL**SSDSC(19) +!WRITE(998,*) IK,ITH,SSDSC(19),SSDSC(11),FACSTRAINL,LMODULATION(ITH) END DO ELSE LMODULATION(:)= 1. @@ -2123,8 +2123,8 @@ SRHS(IS0+1:IS0+NTH)= BRM12(IK)/GRAV**2*BRLAMBDA(IS0+1:IS0+NTH)*C**5 ! diagonal DDIAG(IS0+1:IS0+NTH) = SRHS(IS0+1:IS0+NTH)*SSDSBR/MAX(1.e-20,BTH(1:NTH))/MAX(1e-20,A(IS0+1:IS0+NTH)) ! - END DO -! Breaking probability (Is actually the breaking rate) + END DO +! Breaking probability (Is actually the breaking rate) PB = BRLAMBDA *C ! END SELECT