-
Notifications
You must be signed in to change notification settings - Fork 14
converter for COWVR for the NASA STPH8 mission ocean surface wind product #1647
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: develop
Are you sure you want to change the base?
Conversation
…nternal/ioda-converters into feature/cowvr_osw4ingest
|
The test passed. Shall the |
haydenlj
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks @BenjaminRuston !
fabiolrdiniz
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks, @BenjaminRuston!
| # use approximate offset of 725846427s between 01Jan1970 and 01Jan1993 | ||
| # missing value assumed minimum (e.g. -9999.0) | ||
| time_offset = obs_time_tai93[:].flatten() | ||
| missing_value = np.min(time_offset) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
A unit test would be helpful, especially in understanding the use of missing_value = np.min(.). Isn't the missing value encoded in the attributes of the file?
| # standard angle conventions (0 degrees at East, increasing counter-clockwise) | ||
|
|
||
| # use minimum as missing value | ||
| missing_value = np.min(windSpeed) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Same here...
| windNorthwardError = np.full_like(windSpeed, missing_value, dtype=np.float64) | ||
|
|
||
| # Convert valid wind directions to radians (meteorological to standard) | ||
| rad = np.deg2rad(windDirection[valid_mask]) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I hope the obs file will indicate the convention being used for the direction (0 degrees = North or East) and "wind blowing from".
Description
The data from COWVR as part of the NASA STPH8 mission has an ocean surface wind product
DOI: 10.5067/COWVR-STPH8-EDR100
from NASA earthdata link early 2025:
https://cmr.earthdata.nasa.gov/virtual-directory/collections/C3237678855-POCLOUD
let's move forward as ingest suite can proceed
a unit test is still required!
To test this converter with Skylab ingest suite please use these dependent branches:
https://github.com/JCSDA-internal/skylab/pull/721
https://github.com/JCSDA-internal/ewok/pull/1136
Issue(s) addressed
Resolves #1646
Dependencies
uses the python bindings to IODA
Impact
provide a converter for NASA data files to IODA
Checklist