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