NCAS-CMS / NCAS-CMS/cf-python

reading pp files with non-gregorian calendars does not work

Aperta
#170 4 commenti 0 reazioni 0 assegnatari Vedi su GitHub

Nessuno ha ancora preso questa issue.

bug um/pp
Lingua principale
Python
Stelle
150
Fork
23
Merge medio
1g 11h
PR unite (30g)
2

Descrizione

I am attempting to aggregate as set of UM pp files in a directory.
These include a file which appears to have a 30th of February.
I had expected cf-python to warn me and change to a 360_day calendar
(or allow me to specify it in the read statement).

Error message

 File "/home/users/lawrence/miniconda3/envs/macbook/lib/python3.7/site-packages/cf/read_write/um/umread.py", line 678, in __init__
    c = self.time_coordinate(axiscode)
  File "/home/users/lawrence/miniconda3/envs/macbook/lib/python3.7/site-packages/cf/read_write/um/umread.py", line 2201, in time_coordinate
    vtimes = numpy_array([self.vtime(rec) for rec in recs], dtype=float)
  File "/home/users/lawrence/miniconda3/envs/macbook/lib/python3.7/site-packages/cf/read_write/um/umread.py", line 2201, in <listcomp>
    vtimes = numpy_array([self.vtime(rec) for rec in recs], dtype=float)
  File "/home/users/lawrence/miniconda3/envs/macbook/lib/python3.7/site-packages/cf/read_write/um/umread.py", line 2320, in vtime
    time = netCDF4_date2num(cftime.datetime(*LBVTIME),
 "cftime/_cftime.pyx", line 948, in cftime._cftime.datetime.__init__
  File "cftime/_cftime.pyx", line 1492, in cftime._cftime.assert_valid_date
ValueError: invalid day number provided in cftime.datetime(1992, 2, 30, 0, 0, 0, 0, calendar='gregorian')

Environment: sci6.jasmin.ac.uk, private conda install:

Platform: Linux-3.10.0-1127.19.1.el7.x86_64-x86_64-with-centos-7.8.2003-Core
HDF5 library: 1.10.6
netcdf library: 4.7.4
udunits2 library: /home/users/lawrence/miniconda3/envs/macbook/lib/libudunits2.so.0
python: 3.7.9
netCDF4: 1.5.3
cftime: 1.3.0
numpy: 1.19.4
psutil: 5.7.3
scipy: 1.5.3
matplotlib: 3.3.3
ESMF: not available
cfdm: 1.7.11
cfunits: 3.2.4
cfplot: 3.0.38
cf: 3.1.0

Guida per i contributori

Apri la guida per i contributori

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Direzione di ricerca

Inizia con l’aggregazione di cf.read dei file UM pp e segui l’errore attraverso cf/read_write/um/umread.py, in particolare time_coordinate e vtime. Riproduci il caso del 30 febbraio con l’ambiente indicato o con dati di test equivalenti, quindi verifica che la lettura gestisca il calendario non gregoriano come richiesto oppure segnali chiaramente la limitazione.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Valutazione

Stack tecnologico
python
Ambito
data
Tipo di issue
Bug
Difficoltà
4/5
Tempo stimato
3-5 giorni
Stato di attività
Ferma
Chiarezza
Abbastanza chiara
Idoneità per principianti
30/100

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.