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.

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

Add wave mesh generation for global_ocean mesh test cases #818

Open
wants to merge 26 commits into
base: main
Choose a base branch
from

Conversation

sbrus89
Copy link
Collaborator

@sbrus89 sbrus89 commented Apr 19, 2024

This PR is a collaboration with @erinethomas that adds the workflow used to create coastal-refined wave mesh that conforms to the coastlines for a given ocean mesh. It includes several steps to:

  • build the base mesh
  • cull the base mesh
  • create the wave mesh scrip file
  • generate remapping files between the wave and ocean mesh
  • rotate the mesh
  • generate the files for the unresolved obstacle source term (uost) in WAVEWATCH III

This work depends on the MPAS-Tools PR: MPAS-Dev/MPAS-Tools#579

Checklist

  • User's Guide has been updated
  • Developer's Guide has been updated
  • API documentation in the Developer's Guide (api.rst) has any new or modified class, method and/or functions listed
  • Documentation has been built locally and changes look as expected
  • Document (in a comment titled Testing in this PR) any testing that was used to verify the changes
  • New tests have been added to a test suite

@sbrus89 sbrus89 added the ocean label Apr 19, 2024
@sbrus89 sbrus89 requested a review from xylar April 19, 2024 22:02
@xylar
Copy link
Collaborator

xylar commented Jul 7, 2024

@sbrus89, I had the impression you would ping me when this is ready to review. Just wanted to check that I hadn't missed that ping.

sbrus89 added 5 commits July 15, 2024 11:49
 - This means Antarctic coastline is between ocean and ice
   instead of ocean and grounding line
 - Also add the ability to (optionally) specify paths to existing
   base and culled ocean meshes in the config file
data = np.loadtxt('angled.d')
angled = data[:, 2]

filename_local_in = 'obstructions_local.waves_mesh_culled.in'

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

file can't be found. The file is placed in output directory ... suggest changing the file name to be 'output/obstructions_local.waves_mesh_culled.in'

self.write_alpha_beta(filename_local_out, header, nodes, lon_local, lat_local, sizes, alpha_interp, beta_interp)
self.plot_alpha_beta_spectra(Theta, Freq, alpha_spec, beta_spec, alpha_interp, beta_interp, angled, nodes, 'local')

filename_shadow_in = 'obstructions_shadow.waves_mesh_culled.in'

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

again, file can't be found. The file is placed in output directory ... suggest changing the file name to be 'output/obstructions_shadow.waves_mesh_culled.in'

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants