На случай, если кому-то захочется получить netcdf числа квадратных метров в каждой ячейке сетки длиной в лат. Вероятно, это не самое чистое решение, но в каждой сетке будет создан netcdf (earth_m2.n c) m2 с использованием xarray.
Функция gridsize()
адаптирована из другого вопроса переполнения стека. Затем мы можем создать фиктивный массив и создать поле земли в м2 с использованием долготных расстояний в каждом местоположении.