def prepareSites(cl3files, dt_start, dt_end, args):
#=====================================================================
# Work out how many station models need to be created
#=====================================================================
numModels = 0
params = []
for f in range(0, np.size(cl3files)):
filename = os.path.basename(cl3files[f])
siteID = filename[0:4]
sdata = gsf.parseSite(args.station_file, siteID.upper())
changes = gsf.determineESMChanges(dt_start, dt_stop, sdata)
sitepos = gapr.getStationPos(args.apr_file, siteID)
numModels = numModels + np.size(changes['ind']) + 1
info = {}
info['filename'] = cl3files[f]
info['basename'] = filename
info['site'] = siteID
info['numModels'] = np.size(changes['ind']) + 1
info['changes'] = changes
info['sitepos'] = sitepos
params.append(info)
return params, numModels
def prepareSites(cl3files, dt_start, dt_end, args):
# =====================================================================
# Work out how many station models need to be created
# =====================================================================
numModels = 0
params = []
for f in range(0, np.size(cl3files)):
filename = os.path.basename(cl3files[f])
siteID = filename[0:4]
sdata = gsf.parseSite(args.station_file, siteID.upper())
changes = gsf.determineESMChanges(dt_start, dt_stop, sdata)
sitepos = gapr.getStationPos(args.apr_file, siteID)
numModels = numModels + np.size(changes["ind"]) + 1
info = {}
info["filename"] = cl3files[f]
info["basename"] = filename
info["site"] = siteID
info["numModels"] = np.size(changes["ind"]) + 1
info["changes"] = changes
info["sitepos"] = sitepos
params.append(info)
return params, numModels
dt_stop = dt.datetime(int(args.eyyyy),01,01) + dt.timedelta(days=int(args.edoy)-1)
else:
print("")
print("Warning:")
print("\tusing:",args.resfile,"to work out the time period to determine how many satellites were operating.")
print("")
dt_start = gt.unix2dt(site_residuals[0,0])
res_start = int(dt_start.strftime("%Y") + dt_start.strftime("%j")-1)
dt_stop = gt.unix2dt(site_residuals[-1,0])
res_stop = int(dt_stop.strftime("%Y") + dt_stop.strftime("%j")-1)
print("\tResiduals run from:",res_start,"to:",res_stop)
filename = os.path.basename(args.resfile)
siteID = filename[0:4]
sdata = gsf.parseSite(args.station_file,siteID.upper())
changes = gsf.determineESMChanges(dt_start,dt_stop,sdata)
sitepos = gapr.getStationPos(args.apr_file,siteID)
numModels = np.size(changes['ind']) + 1
info = {}
info['filename'] = args.resfile
info['basename'] = filename
info['site'] = siteID
info['numModels'] = np.size(changes['ind']) + 1
info['changes'] = changes
info['sitepos'] = sitepos
params.append(info)
for s in range(0,info['numModels']):
siteIDList.append(info['site']+"_model_"+str(s+1))
print(
"\tusing:", args.resfile,
"to work out the time period to determine how many satellites were operating."
)
print("")
dt_start = gt.unix2dt(site_residuals[0, 0])
res_start = int(
dt_start.strftime("%Y") + dt_start.strftime("%j") - 1)
dt_stop = gt.unix2dt(site_residuals[-1, 0])
res_stop = int(dt_stop.strftime("%Y") + dt_stop.strftime("%j") - 1)
print("\tResiduals run from:", res_start, "to:", res_stop)
filename = os.path.basename(args.resfile)
siteID = filename[0:4]
sdata = gsf.parseSite(args.station_file, siteID.upper())
changes = gsf.determineESMChanges(dt_start, dt_stop, sdata)
sitepos = gapr.getStationPos(args.apr_file, siteID)
numModels = np.size(changes['ind']) + 1
info = {}
info['filename'] = args.resfile
info['basename'] = filename
info['site'] = siteID
info['numModels'] = np.size(changes['ind']) + 1
info['changes'] = changes
info['sitepos'] = sitepos
params.append(info)
for s in range(0, info['numModels']):
siteIDList.append(info['site'] + "_model_" + str(s + 1))
