GEUS-Glaciology-and-Climate / mass_balance

Greenland ice sheet mass balance from 1840 through next week
https://doi.org/10.5194/essd-13-5001-2021
GNU Lesser General Public License v3.0
6 stars 2 forks source link

SMB uncertainty is >> 9 % (Table 3) #19

Closed mankoff closed 2 years ago

mankoff commented 2 years ago

This study SMB uncertainty is reported in Table 3 as "9 %: Average of 15 % SMB uncertainties above, assuming uncorrelated.". However, that 9 % was leftover from when SMB uncertainties were originally assumed to be 5 %, not 15 %.

import uncertainties
from uncertainties import unumpy
err = unumpy.uarray([1,1,1], [5,5,5])
print('{:.4f}'.format(err.sum()))
>  3.0000+/-8.6603

When using 15 % uncertainties, the average error increases to 26 %.

err = unumpy.uarray([1,1,1], [15,15,15])
print('{:.4f}'.format(err.sum()))
>  3.0000+/-25.9808
mankoff commented 2 years ago

Incorrect. The code was not updated, but the table value of 9 % is correct. Need to be better at writing literate documents where the numbers come from the code.

The 9 % uncertainty is from the mean of the three RCMS with 15 % uncertainty:

err = unumpy.uarray([1,1,1], [15,15,15])
print('{:.4f}'.format(err.mean()))
> 1.0000+/-8.6603