用哈德来环流中心资料画图

ai果汁小鱼

用资料画图，画厄尔尼诺年典型年份1-12月北极海冰变化图
下面是FORTRAN程序，做的区域平均，得出来的数据与事实相反，海冰应该是夏季少冬季多，但是画出来的图却是冬季少，夏季多
program ave
implicit none
integer i,j,it
real sst(361,26,1728),aepjun,aepfeb,aepmar,aepapi,aepmay,aepjune,aepjul,aepaug,aepsep,aepoct,aepnov,aepdec
real epjun(361,26),epfeb(361,26),epmar(361,26),epapi(361,26),epmay(361,26),epjune(361,26),epjul(361,26),epaug(361,26),epsep(361,26),epoct(361,26),epnov(361,26),epdec(361,26)
real aep(12)
real aaepjun,aaepfeb,aaepmar,aaepapi,aaepmay,aaepjune,aaepjul,aaepaug,aaepsep,aaepoct,aaepnov,aaepdec
open(1,file='F:\biye\sst.grd',form='binary')
do it=1,1728
  do i=1,361
    do j=1,26
          read(1)sst(i,j,it)
        enddo
  enddo
enddo
close(1)

do it=1,1728
  do i=1,361
    do j=1,26
if(sst(i,j,it)==-1e+30) then
sst(i,j,it)=0.0
end if
        enddo
  enddo
enddo


aepjun=0.0
aepfeb=0.0
aepmar=0.0
aepapi=0.0
aepmay=0.0
aepjune=0.0
aepjul=0.0
aepaug=0.0
aepsep=0.0
aepoct=0.0
aepnov=0.0
aepdec=0.0
do i=1,361
do j=1,26
epjun(i,j)=(sst(i,j,973)+sst(i,j,997)+sst(i,j,1189)+sst(i,j,1225)+sst(i,j,1273)+sst(i,j,1345)+sst(i,j,1465)+sst(i,j,1537))/8.0
epfeb(i,j)=(sst(i,j,974)+sst(i,j,998)+sst(i,j,1190)+sst(i,j,1226)+sst(i,j,1274)+sst(i,j,1346)+sst(i,j,1466)+sst(i,j,1538))/8.0
epmar(i,j)=(sst(i,j,975)+sst(i,j,999)+sst(i,j,1191)+sst(i,j,1227)+sst(i,j,1275)+sst(i,j,1347)+sst(i,j,1467)+sst(i,j,1539))/8.0
epapi(i,j)=(sst(i,j,976)+sst(i,j,1000)+sst(i,j,1192)+sst(i,j,1228)+sst(i,j,1276)+sst(i,j,1348)+sst(i,j,1468)+sst(i,j,1540))/8.0
epmay(i,j)=(sst(i,j,977)+sst(i,j,1001)+sst(i,j,1193)+sst(i,j,1229)+sst(i,j,1277)+sst(i,j,1349)+sst(i,j,1469)+sst(i,j,1541))/8.0
epjune(i,j)=(sst(i,j,978)+sst(i,j,1002)+sst(i,j,1194)+sst(i,j,1230)+sst(i,j,1278)+sst(i,j,1350)+sst(i,j,1470)+sst(i,j,1542))/8.0
epjul(i,j)=(sst(i,j,979)+sst(i,j,1003)+sst(i,j,1195)+sst(i,j,1231)+sst(i,j,1279)+sst(i,j,1351)+sst(i,j,1471)+sst(i,j,1543))/8.0
epaug(i,j)=(sst(i,j,980)+sst(i,j,1004)+sst(i,j,1196)+sst(i,j,1232)+sst(i,j,1280)+sst(i,j,1352)+sst(i,j,1472)+sst(i,j,1544))/8.0
epsep(i,j)=(sst(i,j,981)+sst(i,j,1005)+sst(i,j,1197)+sst(i,j,1233)+sst(i,j,1281)+sst(i,j,1353)+sst(i,j,1473)+sst(i,j,1545))/8.0
epoct(i,j)=(sst(i,j,982)+sst(i,j,1006)+sst(i,j,1198)+sst(i,j,1234)+sst(i,j,1282)+sst(i,j,1354)+sst(i,j,1474)+sst(i,j,1546))/8.0
epnov(i,j)=(sst(i,j,983)+sst(i,j,1007)+sst(i,j,1199)+sst(i,j,1235)+sst(i,j,1283)+sst(i,j,1355)+sst(i,j,1475)+sst(i,j,1547))/8.0
epdec(i,j)=(sst(i,j,984)+sst(i,j,1008)+sst(i,j,1200)+sst(i,j,1236)+sst(i,j,1284)+sst(i,j,1356)+sst(i,j,1476)+sst(i,j,1548))/8.0


aepjun=aepjun+epjun(i,j)
aepfeb=aepfeb+epfeb(i,j)
aepmar=aepmar+epmar(i,j)
aepapi=aepapi+epapi(i,j)
aepmay=aepmay+epmay(i,j)
aepjune=aepjune+epjune(i,j)
aepjul=aepjul+epjul(i,j)
aepaug=aepaug+epaug(i,j)
aepsep=aepsep+epsep(i,j)
aepoct=aepoct+epoct(i,j)
aepnov=aepnov+epnov(i,j)
aepdec=aepdec+epdec(i,j)

enddo
enddo



aaepjun=aepjun/9386.0
aaepfeb=aepfeb/9386.0
aaepmar=aepmar/9386.0
aaepapi=aepapi/9386.0
aaepmay=aepmay/9386.0
aaepjune=aepjune/9386.0
aaepjul=aepjul/9386.0
aaepaug=aepaug/9386.0
aaepsep=aepsep/9386.0
aaepoct=aepoct/9386.0
aaepnov=aepnov/9386.0
aaepdec=aepdec/9386.0

aep(1)=aaepjun
aep(2)=aaepfeb
aep(3)=aaepmar
aep(4)=aaepapi
aep(5)=aaepmay
aep(6)=aaepjune
aep(7)=aaepjul
aep(8)=aaepaug
aep(9)=aaepsep
aep(10)=aaepoct
aep(11)=aaepnov
aep(12)=aaepdec
do i=1,12
print*,aep(i)
enddo

open(2,file='F:\biye\monep.grd',form='binary')
do i=1,12
write(2) aep(i)
enddo
close(2)



end

huangzhijin

谢谢楼主分享{:eb513:}

kdy379405950

666666666666666666666666666666666666

oucdaqi

{:5_213:}{:5_213:}