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