code: def calcperc_sym(b, perc=(0.1,99.9)):
      
      clim = np.max(np.abs(calcperc(b, perc)))
      #clim = (-clim, clim)
      return -clim, clim
docs: Get symmetrical percentile values
      Useful for determining clim centered on 0 for difference maps