-
Notifications
You must be signed in to change notification settings - Fork 3
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.
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
WW3-MPAS-SI Coupling for E3SMv3 #1
base: e3sm
Are you sure you want to change the base?
Conversation
…ring. Set as default in switch file E3SMv3
@sbrus89 - your review is requested. |
model/ftn/w3fld1md.ftn
Outdated
DO T=1,NTH | ||
INSPC(K,T)=BT(K)*INSPC(K,T)/TPI/(WN2(K)**3.0)/AVG | ||
ENDDO | ||
IF (AVG .GT. 0.0) THEN |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There is a slightly different coding of this fix in the WW3 main.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
thanks for noticing @dabail10 - I will double check
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
thans Dave - fixed.
Re-opening PR with updated WW3 code. |
Pull Request Summary
This PR implements several changes to the WW3 model for wave-sea ice coupling in E3SM.
Description
This PR makes several changes to the WW3 code base. The major changes include:
adding a new 'switch' file with the sea ice settings for wave-sea ice coupling in E3SM (IC4, IS0).
implementing the new wave attenuation scheme for waves in sea ice based on Meylan et al. (https://doi.org/10.1016/j.ocemod.2021.101779). In the E3SM branch of WW3: this setting is "IC4 Method 8".
This attenuation scheme takes into consideration damping in sea ice due to scattering, so scattering term must be turned off (IS0 used in new switch file) to avoid 'double counting' effects due to scattering in sea ice.
Lastly, this PR includes an update to the numerics calculation for IC4 based on calculations from Cecilia Bitz. It is numerically more accurate to determine the attenuation of waves in sea ice THEN determine the breaking of ice due to the damped waves. This is implemented through a new logical namelist parameter under the SIC4 namelist group. The new parameter is called "IC4_NUMERICS" and if this is set to "TRUE" the new numerics fix will be implemented (note this is only implemented for IC4)
Note: This numerics fix is expected to be merged into the official NOAA/EMC branch of WW3 and is being discussed here: Improve numerical solution for waves in sea ice NOAA-EMC/WW3#738