heart: GNSSComponent.get_data_mask invert masking scheme (boolean array)

plotting.geodetic.gnss_fits: with n_ensmeble > 1 fix variance reductions input, i.e. take transpose due to changed atleast2d syntax in numpy

beat/heart.py

@@ -1239,8 +1239,10 @@ def to_kite_scene(self, bins=(600, 600)):
  return scene
 
  def get_data_mask(self, corr_config):
+ """
+ Return dataset mask, i.e. boolean numpy array, 0/1 indicating false/true masking,
+ """
  s2idx = self.station_name_index_mapping()
-
  if (
  len(corr_config.station_whitelist) > 0
  and len(corr_config.station_blacklist) > 0
@@ -1266,8 +1268,8 @@ def get_data_mask(self, corr_config):
  "Stations with idxs %s got blacklisted!"
  % utility.list2string(station_blacklist_idxs)
  )
- mask = num.ones_like(self.lats, dtype=num.bool_)
- mask[num.array(station_blacklist_idxs)] = False
+ mask = num.zeros_like(self.lats, dtype=num.bool_)
+ mask[num.array(station_blacklist_idxs)] = True
  return mask
 
  def station_name_index_mapping(self):

beat/plotting/geodetic.py

@@ -379,7 +379,7 @@ def gnss_fits(problem, stage, plot_options):
 
  out_filename = "/tmp/histbounds.txt"
  m.gmt.pshistogram(
- in_rows=num.atleast_2d(all_var_reductions[dataset.component]),
+ in_rows=num.atleast_2d(all_var_reductions[dataset.component]).T,
  W=(imax - imin) / nbins,
  out_filename=out_filename,
  Z=Z,

beat/pytensorf.py

@@ -547,7 +547,7 @@ def perform(self, node, inputs, output):
  num.array(self.lats), num.array(self.lons), pole_lat, pole_lon, omega
  )
 
- if self.data_mask:
+ if len(self.data_mask) > 0:
  velocities[num.array(self.data_mask), :] = 0.0
 
  z[0] = velocities