利用nc 数据计算某一时段夏季高度场的合成 图形有点问题

我就是-猴子

利用nc数据计算某一时段夏季高度场的合成，发现合成后的图形与ncc网站的图形结果不太一致 ，量级上也相差10倍，请高手帮我看看我的程序啊？
谢谢
program ex1
integer,parameter:: ix=144, iy=73,it=53
real sst(ix,iy,it),ssh1(ix,iy),sjp(ix,iy,it),sjp1(ix,iy,it)
open(1,file='e:/renyj/for/hec/hgt.grd',form='binary')
open(2,file='e:/renyj/for/hec/hgt.jp.grd',form='binary')
   do t=1,it
     read(1) ((sst(i,j,t),i=1,ix),j=1,iy)
!  pause
!  write(*,*)((sst(i,j,t),i=1,ix),j=1,iy)
   enddo
   
!!!!!计算1981-2010年30年的平均!!!!!!!!!!!!!!
do t=21,50
   do j=1,iy
    do i=1,ix
   if (sst(i,j,t)==-9.99e+08)then
     ssh1(i,j)=(-9.99e+08)*30
    else
      ssh1(i,j)=sst(i,j,t)+ssh1(i,j)
   end if
   end do
  end do
end do
ssh1=ssh1/30
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!求1961-2013年 逐年的 距平!!!!!!
do t=1,it
   do j=1,iy
    do i=1,ix
if (sst(i,j,t)==-9.99e+08) then
   sjp(i,j,t)=-9.99e+08
else
   sjp(i,j,t)=sst(i,j,t)-ssh1(i,j)
end if
   end do
end do
end do

!!!!!!!!!!!!!!求制定年份的合成!!!!!!!!!!!!!!
do t=1,16
   do j=1,iy
     do i=1,ix
if (sst(i,j,t)==-9.99e+08) then
   sjp1(i,j,t)=-9.99e+08
else
   sjp1(i,j,t)=sjp1(i,j,t)+sjp(i,j,t)
end if
     enddo
enddo
enddo
sjp1=sjp1/16
do t=1,1
   write(2) ((sjp1(i,j,t),i=1,ix),j=1,iy)
  end do
end