snapshot...

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()