Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

HF_breaking #11

Merged
merged 6 commits into from
Apr 8, 2019
Merged

HF_breaking #11

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
0aa92ed
breaking_bugfix: add proper scaling for radian freq. [radHz]
aronroland Apr 3, 2019
3421754
breaking_bugfix: Add some comments ...
aronroland Apr 3, 2019
92e2ec5
update comments ...
aronroland Apr 4, 2019
d0fb168
update comments ...
aronroland Apr 4, 2019
5677a1f
changed version number in header from 7.00 to 6.07
thesser1 Apr 5, 2019
3faaf20
minor formatting fix on header
thesser1 Apr 5, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
40 changes: 24 additions & 16 deletions model/ftn/w3sdb1md.ftn
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -57,27 +57,37 @@
!/ | FORTRAN 90 |
!/ | J. H. Alves |
!/ | H. L. Tolman |
!/ ! A. Roland |
!/ | Last update : 08-Jun-2018 |
!/ +-----------------------------------+
!/
!/ 25-Apr-2007 : Origination of module. ( version 3.11 )
!/ 08-Jun-2018 : Add DEBUGDB1. ( version 6.04 )
!/ 03-Apr-2019 : Rewrite in terms of energy density (A. Roland,version 6.07)
!/ 03-Apr-2019 : Add Thornton & Guza, 1983 (A. Roland,version 6.07)
!/
! 1. Purpose :
!
! Compute depth-induced breaking using Battjes and Janssen bore
! model approach
!
! 2. Method :
!
! Sdb = - CDB * FMEAN * QB * B * B * SPEC
!
! Where CDB = SDBC1 = 0.25 * BJALFA (defaults to BJALFA = 1)
! modified via ww3_grid namelist parameter BJALFA
! B = HM / HRMS
! HM = GAMMA * DEP
! GAMMA = SDBC2 defaults to 0.73 (mean Battjes/Janssen value)
! modified via ww3_grid namelist parameter BJGAM
! 2. Method : Battjes & Janssen (1978),
!
! Sbr = Dtot/Etot*WA = D * WA
! Dtot = 0.25*alpha*Qb*fm*Hmax²
! fm = sigma/2Pi
! BB = Hrms²/Hmax² = 8Etot/Hmax²
! D = Dtot/Etot = BJALFA * sigma / pi * Qb/BB = 2 * BJALFA * fm * Qb/BB
!
! AR: only valid for Hrms .le. Hm, Qb .le. 1, otherwise, in the degenrative regime it is
! due to Qb > 1 that all wave are broken and Hrms .le. Hmax
! MLIM can be used to enforce this conditions, source term will smoothly converge to this limit.
!
! Where CDB = SDBC1 = BJALFA (defaults to BJALFA = 1)
! modified via ww3_grid namelist parameter BJALFA
! HM = GAMMA * DEP
! GAMMA = SDBC2 defaults to 0.73 (mean Battjes/Janssen value)
! modified via ww3_grid namelist parameter BJGAM
!
! And QB is estimated by iterations using the nonlinear expression
!
Expand Down Expand Up @@ -145,7 +155,6 @@
!/
INTEGER, INTENT(IN) :: IX ! Local grid number
REAL, INTENT(IN) :: A(NSPEC)
!AR: This below is not good I know ... we need more inlined methods ...
REAL, INTENT(INOUT) :: EMEAN, FMEAN, WNMEAN, DEPTH
REAL, INTENT(OUT) :: S(NSPEC), D(NSPEC)
REAL, INTENT(IN) :: CG(NK)
Expand Down Expand Up @@ -179,7 +188,7 @@
!
!/T WRITE (NDST,9000) SDBC1, SDBC2, FDONLY
!
! 1. Integral quantities ... only for the case when it is used nonlinear in the solver ...
! 1. Integral quantities. AR: make sure mean quantities are computed, need to move upward
!
ETOT = 0.
FMEAN2 = 0.
Expand Down Expand Up @@ -256,16 +265,16 @@
IF (IWB == 1) THEN
IF ( ( BB .GT. THR) .AND. ( ABS ( BB - QB ) .GT. THR) ) THEN
IF ( BB .LT. 1.0) THEN
CBJ = DBLE(SDBC1) * QB * DBLE(FMEAN) / BB
CBJ = 2 * DBLE(SDBC1) * QB * DBLE(FMEAN) / BB
ELSE
CBJ = DBLE(SDBC1) * DBLE(FMEAN) * BB ! Make sure the energy vanishes ...
CBJ = 2 * DBLE(SDBC1) * DBLE(FMEAN) * BB ! AR: degenerative regime, all waves must be .le. Hmax, we just smoothly let the excessive energy vanish by * BB.
END IF
ELSE
CBJ = 0.d0
ENDIF
D = - CBJ
S = D * A
ELSE IF (IWB == 2) THEN
ELSE IF (IWB == 2) THEN
IF (ETOT .GT. THR) THEN
HRMS = SQRT(8*EMEAN)
FAK = (1+4./SQRT(PI)*(B*BB+1.5*B)*exp(-BB)-ERF(B))
Expand All @@ -275,7 +284,6 @@
ENDIF
D = - CBJ
S = D * A
ELSE IF (IWB == 3) THEN
ENDIF

IF (CBJ .GT. 0.) THEN
Expand Down
2 changes: 1 addition & 1 deletion model/ftn/ww3_grid.ftn
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1888,7 +1888,7 @@
!/DB1 ' Using Hmax/d in Miche style formulation.'
!/DB1 END IF
!/DB1 WRITE (NDSO,*)
!/DB1 SDBC1 = 0.25 * BJALFA
!/DB1 SDBC1 = BJALFA
!/DB1 SDBC2 = BJGAM
!/DB1 FDONLY = BJFLAG
!
Expand Down