Skip to content

Commit 30adaca

Browse files
author
stephanie
committed
update gapfill
1 parent 507d852 commit 30adaca

1 file changed

Lines changed: 16 additions & 13 deletions

File tree

odmtools/odmservices/edit_service.py

Lines changed: 16 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -149,22 +149,25 @@ def fill_gap(self, gap, period):
149149
points = []
150150
series= self.memDB.series
151151
timegap = np.timedelta64(gap, self.time_units[period])
152-
152+
#if gaps is not of type dataframe- put it in a dataframe
153+
#if not isinstance(gaps, pd.DataFrame
153154
for g in gaps.iterrows():
154-
row= g[0]
155-
e = row["datetime"]
156-
s= row["prevdate"]
157-
print (s)
158-
print (e)
155+
row = g[1]
156+
e = row.datetime
157+
s = row.dateprev
159158

159+
#prime the loop
160+
s = s + timegap
160161
# for each gap time period in the larger gap ( until datetime = prev value)
161-
while s != e:
162-
s= s + timegap
163-
points.append(('-9999', None, s, series.begin_date_time_utc, s, None, None, u'nc', None, None, series.site_id, series.variable_id, series.method_id, series.source_id, series.quality_control_level_id))
164-
#add a row points.append() [('-9999', None, datetime.datetime(2007, 9, 28, 0, 10), '-7', datetime.datetime(2007, 9, 28, 7, 10), None, None, u'nc', None, None, 1, 4, 2, 1, 0)
165-
#('-9999', None, DATE, series.begin_date_time_utc, UTCDATE, None, None, u'nc', None, None, series.site_id, series.variable_id, series.method_id, series.source_id, series.quality_control_level_id
166-
167-
162+
while s < e:
163+
utc_offset = (series.begin_date_time-series.begin_date_time_utc).total_seconds()/3600
164+
points.append((-9999, None, s, utc_offset, s, None, None, u'nc', None, None, series.site_id, series.variable_id, series.method_id, series.source_id, series.quality_control_level_id))
165+
#('-9999', None, DATE, series.begin_date_time_utc, UTCDATE, None, None, u'nc', None, None,
166+
# series.site_id, series.variable_id, series.method_id, series.source_id,
167+
# series.quality_control_level_id
168+
169+
s = s + timegap
170+
print points
168171
self.add_points(points)
169172

170173
time_units = {

0 commit comments

Comments
 (0)