diff --git a/modules/util/gfs_reader.py b/modules/util/gfs_reader.py index b4107744f57dc5dc93931ab955e2c3940b53d76c..3ac5bc78167ed44552a05091c27c181def04fe3c 100644 --- a/modules/util/gfs_reader.py +++ b/modules/util/gfs_reader.py @@ -33,6 +33,7 @@ def get_timestamp(filename): toks = filename.split('.') tstr = toks[1].split('_')[0] dto = datetime.datetime.strptime(tstr, gfs_date_format + '%H').replace(tzinfo=timezone.utc) + dto = dto + datetime.timedelta(hours=12) return dto.timestamp() @@ -43,6 +44,7 @@ def get_time_tuple_utc(timestamp): def get_bounding_gfs_files(timestamp): dt_obj, time_tup = get_time_tuple_utc(timestamp) + dt_obj = dt_obj + datetime.timedelta(hours=12) date_str = dt_obj.strftime(gfs_date_format) dt_obj = datetime.datetime.strptime(date_str, gfs_date_format).replace(tzinfo=timezone.utc) dt_obj_r = dt_obj + datetime.timedelta(days=1) @@ -60,10 +62,10 @@ def get_bounding_gfs_files(timestamp): ftimes = [] for pname in filelist: # TODO: make better with regular expressions (someday) fname = os.path.split(pname)[1] - toks = fname.split('.') - tstr = toks[1].split('_')[0] - dto = datetime.datetime.strptime(tstr, gfs_date_format+'%H').replace(tzinfo=timezone.utc) - ftimes.append(dto.timestamp()) + # toks = fname.split('.') + # tstr = toks[1].split('_')[0] + # dto = datetime.datetime.strptime(tstr, gfs_date_format+'%H').replace(tzinfo=timezone.utc) + ftimes.append(get_timestamp(fname)) tarr = np.array(ftimes) sidxs = tarr.argsort()