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.

Sign up for GitHub

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

Jump to bottom

cfa can't read PP file #165

Closed
JonathanGregory opened this issue Dec 8, 2020 · 4 comments · Fixed by #166
Closed

cfa can't read PP file #165

JonathanGregory opened this issue Dec 8, 2020 · 4 comments · Fixed by #166
Labels
bug Something isn't working CFA Relating to CFA datasets
Milestone
3.8.0

Comments

@JonathanGregory
Copy link

cfa is unable to read the PP file attached in example.zip. It says

racc-login-0-1$ cfa jaaac.003003m.example.pp
cfa ERROR reading /storage/basic/baobab/jonathan/general/jaaac.003003m.example.pp: TODO
racc-login-0-1$ which cfa
/home/users/pg902583/jonathan/scripts/cfa

Please could you help? Thanks. Jonathan
@JonathanGregory JonathanGregory added the bug Something isn't working label Dec 8, 2020
@davidhassell
Copy link
Collaborator

Thanks, Jonathan. Having improved the error message somewhat, I now get:

>> cfa -vm jaaac.003003m.example.pp 
cfa ERROR reading jaaac.003003m.example.pp: Incompatible units: <Units: days since 7481-1-1 360_day>, <Units: days since 7481-9-1 gregorian>
>>

I don't know what this means, but I am investigating ...

@JonathanGregory
Copy link
Author

Thanks! That's an improvement. This is a multiannual seasonal mean. Perhaps that is causing the problem.

@davidhassell
Copy link
Collaborator

OK - I've got to the bottom of it. I now get:

>> cfa -vm jaaac.003003m.example.pp 
Field: toa_incoming_shortwave_flux (ncvar%UM_m01s01i207_vn405)
--------------------------------------------------------------
Data            : toa_incoming_shortwave_flux(latitude(37), longitude(48)) W m-2
Cell methods    : time(1): mean within years time(1): mean over years
Dimension coords: time(1) = [7481-07-16 00:00:00] 360_day
                : latitude(37) = [90.0, ..., -90.0] degrees_north
                : longitude(48) = [0.0, ..., 352.5] degrees_east
>>

The problem was an interaction between improved error trapping in cf.Data and the recent changes in calendar expectation in cftime. The fix will be in v3.8.0, out next week.

@JonathanGregory
Copy link
Author

JonathanGregory commented Dec 8, 2020 •
edited
Loading

Lovely, thanks. Could it be made available sooner via /home/users/pg902583/jonathan?

@davidhassell davidhassell mentioned this issue Dec 8, 2020
Merged
@davidhassell davidhassell added this to the 3.8.0 milestone Apr 26, 2021
@davidhassell davidhassell added the CFA Relating to CFA datasets label Jun 23, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
bug Something isn't working CFA Relating to CFA datasets
Projects
None yet
Milestone
3.8.0
Development

Successfully merging a pull request may close this issue.

specify calendar in cftime.datetime calls davidhassell/cf-python
2 participants
@JonathanGregory @davidhassell