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 an iceMask field to the interpolation script #467

Merged

Conversation

trhille
Copy link
Collaborator

@trhille trhille commented Jun 22, 2022

Add an iceMask field to the interpolation script, which can be used
for trimming down fields to the original extent of ice after
interpolating. This is useful when fields must be extrapolated before
using the interpolation script in order to avoid including zeros, NaNs,
or large missing values in the interpolation at the ice edge.

Add an iceMask field to the interpolation script, which can be used
for trimming down fields to the original extent of ice after
interpolating. This is useful when fields must be extrapolated before
using the interpolation script in order to avoid including zeros, NaNs,
or large missing values in the interpolation at the ice edge.
@@ -649,6 +649,7 @@ def vertical_interp_MPAS_grid(mpas_grid_input_layers, destVertCoord, input_layer
if filetype=='cism':

fieldInfo['thickness'] = {'InputName':'thk', 'scalefactor':1.0, 'offset':0.0, 'gridType':'x1', 'vertDim':False}
fieldInfo['iceMask'] = {'InputName':'iceMask', 'scalefactor':1.0, 'offset':0.0, 'gridType':'x1', 'vertDim':False}
Copy link
Member

Choose a reason for hiding this comment

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

Did you mean to keep this outside of the thicknessOnly check?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Yes, but I don't have a strong preference for where it should go.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Although thinking about it more, I think it is placed appropriately here. The value of iceMask will give some indication of how bad the interpolation ramp is in the thickness field, while you can't know that from looking at the interpolated thickness field alone.

Copy link
Member

Choose a reason for hiding this comment

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

Ok, that makes sense, thanks.

@matthewhoffman matthewhoffman merged commit 554a338 into MPAS-Dev:master Jun 22, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants