In the last post, I forecast next year’s September sea ice extent using not only extent data, but latitude of the ice edge. To compute that from extent, I applied a formula given by Eisenman in the supplemental info of his publication.

A reader pointed out that there seemed to be a discrepancy between his numerical values and mine — and he’s right. I’ve double-checked the approximation formula and my program, but I can’t resolve the discrepancy.

Therefore, I’d be grateful if some of you would reproduce the calculation of latitude based on extent data, using Eisenman’s approximation formula. The formula as given is this:

The data I used for sea ice extent is here (it’s monthly averages from NSIDC):

If you could report the *range* of latitude values you get from the extent data in that file, using the formula as given by Eisenman, that would confirm or deny that I’ve implemented his formula as stated.

Thanks!

P.S. The figure I used for the radius of the earth is 6371 km.