Skip to content
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

filter null ocean colour data #122

Closed
wants to merge 2 commits into from

Conversation

twsearle
Copy link
Collaborator

@twsearle twsearle commented Oct 11, 2024

Description

75% of ocean colour observations are null values in ODB files. We can improve performance and comparibility with the OPS by filtering these observations from the ODB file before JOPA runs. Ideally, we would do this during creation of the ODB file, however I am not sure how to do this, and this may not be in scope for the ODB creation executables.

I attempted to add this to the base query file, however ODB is unable to handle files with a different number of valid varno's. My change was removed as it caused various crashes on different platforms:
https://github.com/JCSDA-internal/ioda/pull/1345

As an alternative, I am temporarily adding the query file to orca-jedi on the assumption that we only need Chlorophyll data (true for PS46) so that we filter down to non-null values.

In the future, I hope that upstream changes by @adammaycock or @DJDavies2 might render this file unnecessary.

Impact

Enables comparison of OPS and JOPA ocean colour outputs.

@twsearle twsearle requested a review from s-good October 11, 2024 09:33
@twsearle twsearle self-assigned this Oct 11, 2024
@twsearle twsearle force-pushed the feature/filter-null-oceancolour branch from 73f737e to 64e7d4f Compare October 11, 2024 10:33
Copy link
Collaborator

@s-good s-good left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm not sure how this file will get used by orca-jedi. Please could you clarify?

I'm wondering if we can't put the query back into ioda and remove the second varno if we don't need it?

@@ -0,0 +1,15 @@
where:
varno: [288,]
query: initial_obsvalue is not null
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is the indentation correct for this line? In https://github.com/JCSDA-internal/ioda/pull/1345/files the query is indented the same as varno.

@twsearle
Copy link
Collaborator Author

I'm not sure how this file will get used by orca-jedi. Please could you clarify?

I would be using this in the AMM7 sith change, so that we can compare the outputs across OPS and JOPA.

I'm wondering if we can't put the query back into ioda and remove the second varno if we don't need it?

I think I can't put it in ioda as it doesn't have a test. I hesitate to put it in as the change was only just removed from ioda. Another option would be to put it in jjdocs/jjaux. I honestly don't know what the preferred protocol is for this or what others are doing. As I say I think ideally, we wouldn't have these large ODB files with mostly missing data so I would like to resolve it upstream.

@s-good
Copy link
Collaborator

s-good commented Oct 11, 2024

Would it work if it were inserted into the ODB->ODB2 conversion a bit like is done in https://github.com/MetOffice/jjdocs/blob/develop/docs/share/jopa/ocean/apps/ocean_ops_codb_odb_to_odb2/opt/rose-app-o-ts.conf#L2?

If not, I'm happy to approve this.

@twsearle
Copy link
Collaborator Author

Would it work if it were inserted into the ODB->ODB2 conversion a bit like is done in https://github.com/MetOffice/jjdocs/blob/develop/docs/share/jopa/ocean/apps/ocean_ops_codb_odb_to_odb2/opt/rose-app-o-ts.conf#L2?

If not, I'm happy to approve this.

I will give it a go! That at least should be a bit better in performance. Ideally we don't want it in the ODB in the first place, but this does seem like a better place to put it thanks.

@twsearle
Copy link
Collaborator Author

Looks like using the ODB->ODB2 step works closing this one.

@twsearle twsearle closed this Oct 11, 2024
@twsearle twsearle deleted the feature/filter-null-oceancolour branch October 11, 2024 15:25
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants