Interpolate ECCO ASTE reanalysis to MALI grid #691
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Rearranges script so that remapping happens for each individual tile. Allows ECCO arrays to remain gridded and improves efficiency of ESMF_RegridWeightGen.
This reverts commit bf4de92db9d8f40500493c15cf33ed145166cb95.
Incorporates ECCO grid file to create custom scrip file for ECCO format
Hardcodes ECCO tile numbers and explicitly defines edge cells along boundary for the two tiles. Currently only tiles 14 and 27 are supported. Not enough information exists in ecco grid files to define boundary cells edges without hardcoding (unconventional MITgcm grid file format)..
Addresses formating issue with xtime variable in initial output file. Still need to fix formatting in 'meshVars' output file.
Conservative remapping averages defined and undefined ocean cells. Use orig3dOceanMask to identify these cells and treat as undefined
Creates a consistent naming scheme so that only one output file is created and the remaining temporary files are removed
Extrapolates iceCellArea to fill ocean cells, and defines icebergFjordMask wherever iceCellArea is above a threshold value
Adds an option to use a geojson file to define a region where icebergFjordMask is permanently 1
alexolinhager
force-pushed
the
interpolate_ecco_to_mali
branch
from
f13ce7c to
daac57e
Compare
Moves saving of xtime time down to where final output file is created in order to preserve its format
Creates a copy of the output file with netcdf3 64-bit offset formatting to be consistent with MALI requirements
Monthly averaged ECCO data is output at the beginning of following month. Need to subtract a month from the time stamp in order to force MALI with the correct monthly ocean conditions.
Previously, a few invalid TS cells were designated as valid in the orig3dOceanMask after remmapping, which caused issues when extrapolating ocean properties. This commit add a check to remove this cells if they occur
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
1 participant
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR introduces
landice/mesh_tools_li/interpolate_ecco_to_mali.py, a script that interpolates ECCO ASTE R1 reanalysis data (https://arcticdata.io/data/10.18739/A2CV4BS5K/) onto a MALI mesh. Currently supported fields are ocean temperature, ocean salinity, and sea ice fraction. ECCO grid files are unconventional for MITgcm simulations and do not contain adequate grid cell corner information for tile edge cells (necessary for creating scrip file and remapping). Therefore, this information needed to be hardcoded intointerpolate_ecco_to_mali.py, using grid corner information from adjacent MITgcm tiles. Currently supported tiles are 14 and 27, which cover the western and northern coasts of Greenland.