Skip to content

Bug fixes to MPAS->Albany mesh converter #338

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

Conversation

matthewhoffman
Copy link
Member

This merge fixes some small bugs in the MPAS->Albany mesh converter. With these changes, the grounding line flux difference between MALI and standalone Albany is down to 10 digits of precision.

With these changes, the GLF difference between MALI and standalone Albany is down to
10 digits of precision
* write also the labels (floating/dirichlet) for boundary edges and nodes
* print the mesh and ascii files in double precision.
* Fix a small bug in the MPAS interface. In some very peculiar cases we set the thickness
to 2 or 3 meters. Then we change accordingly the surface elevation in those location to
comply with the floating point condition. However, instead of setting
h = (1-rho_i/rho_o) H, we were setting h = (rho_o/rho_i -1) H. The impact on the GLF
is small, 10e-6 relative error or so, but better to fix it.
@matthewhoffman
Copy link
Member Author

Note: this PR to be replaced by a more extensive update to the interface code in progress by @mperego .

@matthewhoffman
Copy link
Member Author

Replaced by #402

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

Successfully merging this pull request may close these issues.

3 participants