sst的dat数据求距平，有缺省值

乌云珠

我下的是HadISST_sst.nc的数据，已经转化为dat的了，要求距平，下面是我的Fortran 程序，感觉应该是缺省值的在if 或循环中的问题，求帮助！求程序呀！
program ex1
integer,parameter:: ix=360, iy=180,it=732 ,ik=12
real sst(ix,iy,it),ssh1(ix,iy,ik),sjp(ix,iy,it)
open(1,file='d:\my4\sst.dat',form='unformatted', access='direct', recl=1)
open(2,file='d:\my4\sjp.dat',form='unformatted',access='direct',recl=1)
irec=0
do t=1,it
  do j=1,iy
   do i=1,ix
  irec=irec+1
read(1,rec=irec) sst(i,j,t)
  enddo
end do
end do
write(*,*) sst(123,:,12)
irec=0
do k=1,12
do t=k,it,12
  do j=1,iy
   do i=1,ix
  irec=irec+1
  if (sst(i,j,t)==-9.99e+08)then
   ssh1(i,j,k)=(-9.99e+08)*61
   else
ssh1(i,j,k)=sst(i,j,t)+ssh1(i,j,k)
  end if
  end do
end do
end do
end do
ssh1=ssh1/61

irec=0
do k=1,12
do t=k,it,12
  do j=1,iy
   do i=1,ix
  irec=irec+1
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,k)
end if
  end do
end do
end do
end do
write(*,*)sjp(123,167,:)
irec=0
do t=1,it
  do j=1,iy
   do i=1,ix
  irec=irec+1
write(2,rec=irec) sjp(i,j,t)

  end do
end do
end do
end

river

找茬游戏吗？您的问题在哪啊？

kongfeng0824

可以贴出来你的错误信息是什么吗？

漂流的屋顶

我同学问的，帮她说一下哈。她的一个站点有的月份是缺省值，有的月份是正常海温。不晓得出了什么怪问题。正常运行出的数据，有的达到了缺省值量级还不是缺省值。是不是数据有错误啊？还是说就是有这样的数据，是需要在处理一下呢？

我就是-猴子

ssh1(i,j,k)=(-9.99e+08)*61
这里为什么要*61  不太明白

bree115

是不是想缺省值在求了距平之后仍然为缺省？

bree115

是不是想缺省值在求了距平之后仍然为缺省？

李小毛123

漂流的屋顶 发表于 2013-4-25 21:34
我同学问的，帮她说一下哈。她的一个站点有的月份是缺省值，有的月份是正常海温。不晓得出了什么怪问题。正 ...

啊啊啊你同学知道什么原因了吗？我也出现了这个问题，同一个站点有些月份缺省有些月份不缺省，最终所有值都达到了缺省值的量级