def read_periodfile(c_Model, c_Var):
# Routine to read variable from hourly (HH) file
c_Files = ops.find(IncF.c_ModelDir, 'PP_' + c_Model + '-2D*')
if c_Files is None:
return None
c_FileName = c_Files[0]
c_ModelFile = c_FileName.split(os.sep)[-1]
logging.info('Reading File: %s', c_ModelFile)
new_data = ncfile.read_ncfile(c_FileName, c_Var)
if new_data is not None:
f_Data_aux, f_lat, f_lon, d_Time_aux, c_Units = new_data
else:
return None
# CHECKING FOR TIME FREQUENCIES BETWEEN TIME STEPS
c_TimeFreq, Dt = delta_time_freq(d_Time_aux)
f_Time_aux = map(float, date2num(d_Time_aux))
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_Data, d_Time = time_crop(f_date_i, f_date_f, Dt, f_Time_aux, f_Data_aux)
return f_Data, f_lat, f_lon, d_Time, c_ModelFile, c_TimeFreq, c_Units
# Funciones auxiliares
评论列表
文章目录