计算球坐标下1×1°每个网格面积

折磨死

大家好，我需要做出的是北极海冰面积在30年内的变化曲线，资料是水平分辨率1×1°的逐月海冰密集度的nc文件。
定义海冰面积为海冰密集度>15%的网格面积之和，因此我需要先确定该格点的海冰密集度是否>15%，然后计算该格点的面积，然后把所有海冰密集度>15%的格点面积求和。所以我的问题是如何根据经纬度计算出水平分辨率1×1°的每个网格的面积呢？

wjy_ecnu

等间距经纬度可以近似用cos latitude来表示面积的比例

折磨死

好哒，谢谢你，我之前还想着积分不会编呢

wsx

wjy_ecnu 发表于 2020-3-30 18:33
等间距经纬度可以近似用cos latitude来表示面积的比例

请问可以详细说明下嘛？

wsx

折磨死 发表于 2020-3-30 19:22
好哒，谢谢你，我之前还想着积分不会编呢

请问可以分享下计算过程吗？

折磨死

DO K=1,360  
READ(10) ((sic(i,j,k),i=1,nx),j=1,ny)
ENDDO
DO K=1,360
area(K)=0.0
DO i=1,nx   
DO j=1,ny        
IF(sic(i,j,k)>=0.15.and.sic(i,j,k)<=1.0)THEN                
;print*,i,j,sic(i,j,k)                                                                              
      area(K)=area(K)+(6371**2)*(sin((FLOAT(J)+0.5)*3.1416/180)-sin((FLOAT(J)-0.5)*3.14/180))*((i+1)*3.14/180-i*3.14/180)

我当初是这么算的