ax1.set_ylabel('km$^3$/year')
km3yearToSv = 3.16887646e-05
ax2 = ax1.twinx()
ax2.set_ylim(mn * km3yearToSv, mx * km3yearToSv)
ax2.set_ylabel('Sv')
Plot single data with two Y axes (two units) in matplotlib