用origin8.5作图出现一点问题，作图是欧亚极涡的每月所有格点累加的平均值，实在不懂

蓝色蜗牛角

在写论文中画到上面这组图很怪，明显图片是有规律的，但应该是两边蓝色，中间红色，或者两边红色，中间蓝色。这个图画出来完全没有给我修正的机会。其中nc转dat的gs'reinit'
'sdfopen C:\lunwen3\ndvi_1981~2015_mon_0.432_SM_supply.nc'
'set gxout fwrite'
'set fwrite C:\lunwen3\o3.dat'
tn=1
while(tn<=420)
'set t 'tn
'set lon -180 179.568'
'set lat -54.7195 77.024'
'set z 1'
'd ndvi'
tn=tn+1
endwhile
'disable fwrite'

用dat做平均转成txt程序
program main
implicit none
integer,parameter::xl=833,yl=301,zl=1,nt=420
integer::t,i,j,m,a
real::o3(xl,yl)
real::area(420),o3_ave(420)
open(12,file='C:/lunwen3/o3.dat',form='binary')
open(15,file='C:/lunwen3/o31.txt',form='formatted')
o3_ave=0
do t=1,420
a=0
do i=1,833
do j=1,301

read(12)o3(i,j)
end do
end do
do i=393,833
do j=196,301
if(o3(i,j)>0.0001)then
a=a+1
  o3_ave(t)=o3_ave(t)+o3(i,j)
  end if
  end do
end do
write(15,*)o3_ave(t)/a
end do
close(12)
close(15)
end program