def read_mixedlayer(self, c_Network):
c_ObsNetDirs = IncDir.c_ObsNetDir
c_ObsNetName = IncDir.c_ObsNetName
idx_Net = c_ObsNetName.index(c_Network)
c_Files = os.listdir(c_ObsNetDirs[idx_Net])
logging.info('Data Directory: %s' % c_ObsNetDirs[idx_Net])
i_count = 0
if 'PBLHmunoz.txt' in c_Files:
c_FileName = 'PBLHmunoz.txt'
logging.info('Reading File: %s' % c_FileName)
f_AuxData = []
d_AuxDate = []
with open(c_ObsNetDirs[idx_Net] + c_FileName, 'r') as f:
for line in (row.split(',') for row in f):
if i_count > 0:
f_AuxData.append(float(line[3]))
d_AuxDate.append(dt.datetime.strptime(line[0], '%d-%m-%Y_%H:%M'))
i_count += 1
f_AuxDate = date2num(d_AuxDate)
f_date_i = date2num(dt.datetime.strptime(IncF.c_Start_Date[0], '%d-%m-%Y'))
f_date_f = date2num(dt.datetime.strptime(IncF.c_Last_Date[0], '%d-%m-%Y'))
f_Stntime = []
d_Stntime = []
f_date_aux = f_date_i
d_date_aux = num2date(f_date_i)
while f_date_aux <= f_date_f + 23 / 24.:
f_Stntime.append(date2num(d_date_aux))
d_Stntime.append(d_date_aux)
d_date_aux = d_date_aux + dt.timedelta(hours=1)
f_date_aux = date2num(d_date_aux)
f_Data = np.empty(len(f_Stntime))
f_Data.fill(IncF.f_FillValue)
f_Data[np.in1d(f_Stntime, f_AuxDate)] = f_AuxData
return f_Data, f_Stntime
评论列表
文章目录