Interpolate ECCO ASTE reanalysis to MALI grid#691
Open
alexolinhager wants to merge 20 commits intoMPAS-Dev:masterfrom
Open
Interpolate ECCO ASTE reanalysis to MALI grid#691alexolinhager wants to merge 20 commits intoMPAS-Dev:masterfrom
alexolinhager wants to merge 20 commits intoMPAS-Dev:masterfrom
Conversation
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
Changes default mapping method to bilinear. Better for remapping coarse resolution ECCO data onto finer MALI mesh
Contributor
There was a problem hiding this comment.
Pull request overview
This PR adds a new script to interpolate ECCO ASTE R1 reanalysis data (ocean temperature, salinity, sea ice fraction) onto a MALI mesh. It handles the unconventional ECCO grid format by hardcoding corner information for tiles 14 and 27 (covering western and northern Greenland coasts).
Changes:
- New
interpolate_ecco_to_mali.pyscript that merges MITgcm tile files into an unstructured grid, creates SCRIP files, and remaps ECCO data onto MALI meshes using ESMF regridding.
💡 Add Copilot custom instructions for smarter, more guided reviews. Learn how to get started.
Addressing mostly minor reviews from copilt review. Most substantial change is a switch to salt/temp as required variables, while SIArea remains an optional variable. This is to circumvent potential issues setting up orig3dOceanMask, but shouildn't impact performance as salt/temp would almost always be remapped anyway
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
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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. Relevant ECCO files are currently stored at:/global/cfs/cdirs/fanssie/users/ahager/gisMeshSetup/forcingFiles/ASTEon Perlmutter, although these should eventually be moved to a permanent location. 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.Script can be tested with (on permutter):
python /global/cfs/cdirs/fanssie/users/ahager/MPAS-Tools/landice/mesh_tools_li/interpolate_ecco_to_mali.py --eccoDir /global/cfs/cdirs/fanssie/users/ahager/gisMeshSetup/forcingFiles/ASTE --maliMesh /global/cfs/cdirs/fanssie/MALI_input_files/UummannaqDisko_1to10km_r01/optimization/FO/UummannaqDisko_1to10km_r01_20251030_FO.nc --eccoVars SALT THETA SIarea --meshVars --geojson /global/cfs/cdirs/fanssie/users/ahager/gisMeshSetup/forcingFiles/UummannaqDisko/geojson/Ilulissat_Basin.geojson