st_greater_than = 260
st_less_than = None
algo = eustace.surface_temperature.ST_ALGORITHM.IST
elif algorithm == eustace.surface_temperature.ST_ALGORITHM.IST + "_GT_240_LT_260":
st_less_than = 260
st_greater_than = 240
algo = eustace.surface_temperature.ST_ALGORITHM.IST
else:
st_less_than = None
st_greater_than = None
algo = algorithm
y_array = np.array([row[0] for row in db.get_perturbed_values(swath_variables=None,
lat_less_than=args["--lat-lt"],
lat_greater_than=args["--lat-gt"],
tb_11_minus_tb_12_limit=args["--t11-t12-limit"],
st_less_than=st_less_than,
st_greater_than=st_greater_than,
algorithm=algo,
limit=limit)])
LOG.debug("Took: %s" % (str(datetime.datetime.now() - t)))
# Number of samples - total.
LOG.info("Number of samples: %i." %(len(y_array)))
# Make sure that there are no nan in the array.
y_array_is_not_nan = y_array[~np.isnan(y_array)]
# Number of samples.
LOG.info("Number of samples without NaN: %i." %(len(y_array_is_not_nan)))
LOG.debug("Calculating the average.")
x_arrays[variable]=[]
random.seed(1)
LOG.debug("Get the values from the database.")
t = datetime.datetime.now()
with eustace.db.Db(args["<database-filename>"]) as db:
# Where sql used for the title in the plots.
where_sql = db.build_where_sql(lat_less_than=args["--lat-lt"],
lat_greater_than=args["--lat-gt"],
tb_11_minus_tb_12_limit=args["--t11-t12-limit"],
algorithm=args["--algorithm"])
# Get the values.
for row in db.get_perturbed_values(variable_names,
lat_less_than=args["--lat-lt"],
lat_greater_than=args["--lat-gt"],
tb_11_minus_tb_12_limit=args["--t11-t12-limit"],
algorithm=args["--algorithm"],
limit=limit):
y_array.append(row[0])
for i in range(len(variable_names)):
if row[i + 1] == None:
x_arrays[variable_names[i]].append(np.NaN)
else:
x_arrays[variable_names[i]].append(row[i + 1])
LOG.debug("Took: %s" % (str(datetime.datetime.now() - t)))
LOG.info("%i samples" %(len(y_array)))
y_array = np.array(y_array)
y_array_is_not_nan = y_array[~np.isnan(y_array)]
average_all = np.average(y_array_is_not_nan)