add interpolate_mpasOcean_to_mali.py #568
Conversation
Creates new script that interpolates MPAS-Ocean thermal forcing and sub ice shelf melt rates to a MALI mesh for one way coupling. Still a work in progress at this point
|
@matthewhoffman, opening this now so you can get a sense of the script so far. I added "To Do" notes where needed |
Implements the thermal forcing calculation, as well as reorganizes script to object oriented. Further work is labeled with "<< TO DO >>:"
Finishes to-do list of interpolation script. Script capabilities are now all complete, but script still needs to be tested and likely debugged
Two changes to get xtime written correctly: * when the char array is generated, pad it to 64 characters * fix a bug in that the return code to encoding.update was being assigned to xtime
Previous code was giving an error when landIceFreshwaterFlux exists.
minLevelCell does not exist in older meshes, which assumed it was 1.
…al attr The global attribute version could sometimes have the value 'file'.
Verification testingmeltAs a reference for the melt field, I time averaged one year of SORRM monthly output files in one directory: I compare this against the output of the script applied to that same set of files. There is a sign convention error that I will push a fix for, but otherwise results look the same ( So I think this verifies both the time averaging code and the grid remapping code. TFTop layer of thermal forcing in the output: 23rd layer: While the spatial patterns look appropriate for the melt field, the values seem off, with most of the continent having negative values. Also, at shallow layers, there is some strange scalloping in the open ocean regions. I can't think of a straightforward comparison to the MPAS-Ocean data here. There could be issues with either the TF calculations, the vertical remapping, or both. |
matthewhoffman
force-pushed
the
interpolate_mpasOcean_to_mali
branch
from
ad612ad
to
c815154
Compare
… files There will likely be other files in output directories, so ensure we are only working with the relevant ones.
Reorders script to first create separate files containing the interpolated thermal forcing and ice shelf melt for each year. Then reopens each annual file and concatenates into the final output file at the end. Done to save memory
Creates new script that interpolates MPAS-Ocean thermal forcing and sub ice shelf melt rates to a MALI mesh for one way coupling. Still a work in progress at this point