Reading NSIDC Sea Ice Edge NetCDF Files #1718
-
|
Dear METplus team, We tried to open a NSIDC sea ice edge netCDF file using plot_data_plane, but got the following error messages:
This file is in a predefined polar stereographic grid. It's a curvilinear grid that CDO doesn't recognize and cannot be remapped into latlon grid using "cdo remapbil". This data set is specified as a validation source for the EVS project at EMC, so we would like to find a way to read it into MET/METplus. I put two sample files at ftp.rap.ucar.edu:/incoming/irap/met_help/chen_data. Please advise. Thanks, |
Beta Was this translation helpful? Give feedback.
Replies: 2 comments 2 replies
-
|
Hi Gwen: The longitude of projection origin is not an attribute of the crs variable which is the specific error message you are getting. However, I am not sure what the best way to get this data into MET would be and I need to double check with a few people here. I will get back to you on this. Christina |
Beta Was this translation helpful? Give feedback.
-
|
@GwenChen-NOAA, I work with @CPKalb on providing support for METplus. FYI, I wasn't able to run the two NetCDF files you provided through ncdump or ncview. I assume they got corrupted somewhere along the way. But I did find and test this file available publicly here: I'm pretty sure you've uncovered a bug in MET in the parsing of Polar Stereographic grids from CF-compliant NetCDF files. The CF-Convention documentation defines Polar Stereographic grids here. And this is parsed by the get_grid_mapping_polar_stereographic() function in MET. I see 2 issues in the code:
This will require a bugfix to patch up the logic. |
Beta Was this translation helpful? Give feedback.
-
|
Thank you both! And apologize for the corrupted files. I'm not sure what happened to the files, but the one you used is from the same source where I downloaded the data. When will this bugfix be available in MET? Thanks, |
Beta Was this translation helpful? Give feedback.
-
|
@GwenChen-NOAA sorry for the late response here. As we discussed in the last NOAA METplus telecon, we've created a bugfix issue (dtcenter/MET#2218) in MET to correct this. @hsoh-u is working on it now, and we plan to include the fix in the next development version for MET, MET-11.0.0-beta3. Looks like we do also plan to include a fix in the main_v10.1 branch so that it could be included in a MET-10.1.3 bugfix release. We don't have any immediate plans for the date of that bugfix release. But please feel free to comment directly on the dtcenter/MET#2218 issue to let us know if/when it's needed. I'll go ahead and mark this discussion as being answered and recommend that any future discussion for it be done via comments on dtcenter/MET#2218. Thanks! |
Beta Was this translation helpful? Give feedback.
@GwenChen-NOAA, I work with @CPKalb on providing support for METplus.
FYI, I wasn't able to run the two NetCDF files you provided through ncdump or ncview. I assume they got corrupted somewhere along the way. But I did find and test this file available publicly here:
I'm pretty sure you've uncovered a bug in MET in the parsing of Polar Stereographic grids from CF-compliant NetCDF files.
The CF-Convention documentation defines Polar Stereographic grids here. And this is parsed by the get_grid_mapping_polar_stereographic() function in MET.
I see 2 issues in the code:
longitude_of_project…