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

Update w3srcemd.F90 for considering the contribution of ice scattering of IS1 #1387

Open
wants to merge 2 commits into
base: develop
Choose a base branch
from
+6 −0
Open

Update w3srcemd.F90 for considering the contribution of ice scattering of IS1 #1387

Changes from 1 commit
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
579c004
Update w3srcemd.F90
Biao-Zhao Mar 7, 2025
f8b6abe
Update w3srcemd.F90 for considering the contribution of ice scatterin…
Biao-Zhao Mar 7, 2025
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
Next Next commit
Loading
Loading status checks…
Update w3srcemd.F90 
I reported a bug in #1386, the current IS1 is not contributing to action density spectrum. Since scattering doesn’t dissipate wave energy, I believe the correct solution is to redistribute the scattered energy to other directions, similar to IS2.
  • Loading branch information
@Biao-Zhao
Biao-Zhao authored Mar 7, 2025

Verified

This commit was created on GitHub.com and signed with GitHub’s verified signature.
GPG key ID: B5690EEEBB952194
Verified
Learn about vigilant mode
commit 579c00485b112d5900156eaf28a4bc1ec933a45b
6 changes: 6 additions & 0 deletions model/src/w3srcemd.F90
View file
Original file line number Diff line number Diff line change
@@ -2124,6 +2124,12 @@ SUBROUTINE W3SRCE ( srce_call, IT, ISEA, JSEA, IX, IY, IMOD, &
!
! Second part of ice term integration: scattering including re-distribution in directions
!

#ifdef W3_IS1
! IS1 doesn't dissipate wave energy, it just redistribute the scattered energy to other directions
SPEC(1+(IK-1)*NTH:NTH+(IK-1)*NTH) = SPEC(1+(IK-1)*NTH:NTH+(IK-1)*NTH) + VSIR(IS)*DTMIN
#endif

#ifdef W3_IS2
IF (IS2PARS(2).GE.0) THEN
IF (IS2PARS(20).GT.0.5) THEN