Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension


Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions CHANGELOG.md
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,7 @@ and uses [Semantic Versioning](https://semver.org/spec/v2.0.0.html).

### Fixed
* [811](https://github.com/dbekaert/RAiDER/pull/811) - Fixed a north-south flip in the ECMWF/ERA-5 model-level reader. `_load_model_level` reversed only the level axis of the geopotential cube (`z[::-1]`) while `t`, `q`, and `lnsp` had their latitude axis reversed, leaving heights mirrored in latitude relative to the meteorology. Surface height and surface pressure were anti-correlated (r = -0.54 where physics requires ~+1); sea-level ZTD over a 3-degree box spanned 1621-3233 mm instead of 2239-2395 mm. The error is antisymmetric in latitude, so it largely cancels in a domain average while being severe at individual points.
* Update dem to match NISAR DEM as noted here: https://github.com/ACCESS-Cloud-Based-InSAR/dem-stitcher/blob/dev/notebooks/analysis_and_comparison/2_Comparison_with_NISAR_DEM.ipynb

### Changed
* [808](https://github.com/dbekaert/RAiDER/pull/808) - Use `xarray.open_dataset` instead of `load_dataset` for memory efficiency, and guard weather-model file reads with explicit error handling so that datasets are always closed. Dropped the redundant full-cube `shutil.copy` from the ECMWF model-level download, and added regression tests covering the layout contract between `ECMWF._get_from_cds` and `ECMWF._makeDataCubes`.
Expand Down
2 changes: 1 addition & 1 deletion environment.yml
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@ dependencies:
- cdsapi
- cfgrib
- dask
- dem_stitcher>=2.5.8
- dem_stitcher>=3.0.0
- ecmwf-api-client
- geopandas
- h5netcdf
Expand Down
2 changes: 1 addition & 1 deletion tools/RAiDER/dem.py
Original file line number Diff line number Diff line change
Expand Up @@ -62,7 +62,7 @@ def download_dem(
list(bounds),
dem_name='glo_30',
dst_ellipsoidal_height=True,
dst_area_or_point='Area',
dst_area_or_point=None,
)
metadata = cast(RIO.Profile, metadata)
if writeDEM:
Expand Down