docs: def apb(s, o):
      """
          Absolute Percent Bias
          input:
          s: simulated
          o: observed
          output:
          apb_bias: absolute percent bias
          """
      # s,o = filter_nan(s,o)
      return 100.0*np.sum(np.abs(s-o))/np.sum(o)