hadley资料计算海冰范围

啦啦啦啦啦啦

我想问一下，hadley中心的海冰百分数(sea ice area fraction)资料，怎么计算海冰范围(sea ice extent,10^6 km^2)呀?


我的资料是这样的：


我的计算方法是这样：
pi = 4.0*atan(1.0)
nx = 89;lat
ny = 359;lon
nt = 1799;time

do i = 0,nt
    sien(i) = 0.0
   do j = 0,nx
      do k = 0,ny
         if (.not.ismissing(sic(i,j,k))) then
            if (sic(i,j,k) .gt. 0.15) then
                sien(i) = sien(i) + 12351.77853*cos((89.5-j)*pi/180)           
            end if
        end if
      end do
   end do
end do


也是从站上别的那里学到的计算方法，但是结果显示，在北半球，3月份的海冰范围最小，9月份反而最大。
不知道是程序出错了还是哪边有问题，求求各位大佬指点一下，拜托拜托。。

602自习室

请问楼主问题解决了吗？急需要计算海冰范围，求助

萤火矢车菊

计算海冰范围，收钱将海冰覆盖率大于15%的格点通过where语句设置成数值为1，其他格点设置为缺省值或者为0，再通过wgt_areasum2进行区域求和