code: def _compute_full_ts(self, data):
          
          # Get results at each desired timestep and spatial point.
          full_ts, dt = self._compute(data)
          # Vertically integrate.
          vert_types = ('vert_int', 'vert_av')
          if self.dtype_out_vert in vert_types and self.var.def_vert:
              dp = self._get_input_data(_DP_VARS[self.dtype_in_vert],
                                        self.start_date, self.end_date)
              full_ts = utils.vertcoord.int_dp_g(full_ts, dp)
              if self.dtype_out_vert == 'vert_av':
                  ps = self._get_input_data(utils.vertcoord.ps,
                                            self.start_date, self.end_date)
                  full_ts *= (GRAV_EARTH / ps)
          return full_ts, dt
docs: Perform calculation and create yearly timeseries at each point.