Skip to content
Snippets Groups Projects
Commit e3c679d2 authored by Anne Philipp's avatar Anne Philipp
Browse files

BugFix for Ticket #286

parent 697b8d00
No related branches found
No related tags found
No related merge requests found
......@@ -19,5 +19,4 @@ ETA 1
FORMAT GRIB2
PREFIX EA
CWC 1
RRINT 1
ECTRANS 1
......@@ -504,13 +504,23 @@ class MarsRetrieval(object):
date = datetime.strptime(attrs['date'], '%Y%m%d')
newattrs['year'] = date.year
newattrs['month'] = date.month
newattrs['day'] = date.day
newattrs['day'] = date.day
# need to correct the time apperance for CDS surface field retrievals
if attrs['type'] == 'FC': # for EA5 only flux fields are retrieved as FC type
# need to convert fc start times 06/18 to usual AN times
# since the surface fields can only be access through their validity time
start, end, step = map(int,attrs['step'].split('/')[::2])
newattrs['time'] = [ "{0:0=2d}".format(s) for s in range(0,24,step) ]
elif '/' in attrs['time']: # we expect a list of times separated by /
newattrs['time'] = attrs['time'].split('/')
elif isinstance(attrs['time'], str): # we expect a single time in here
newattrs['time'] = [ attrs['time'] ]
newattrs['product_type'] = 'reanalysis'
newattrs['area'] = attrs['area'].split('/')
newattrs['grid'] = list(map(float,attrs['grid'].split('/')))
newattrs['param'] = attrs['param'].split('/')
newattrs['time'] = list(map(str,range(0,24,3)))
newattrs['param'] = attrs['param'].split('/')
newattrs['format'] = 'grib'
return newattrs
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment