用grads 画图出现 contouring：-100 to 5800 interval 100 求 大神指教

xiaocaoqiqiao

fortran程序
program EX1
real a(144,73,12,65),ave1(144,73),ave7(144,73),ave8(144,73)
open(2,file='d:\hgt500.grd',form='binary')

open(4,file='d:\hgt5001.grd',form='binary')
      do it=1,65
do imo=1,12
do j=1,73
do i=1,144
read(2)a(i,j,imo,it)
enddo;enddo;enddo;enddo

cccccccccccccc  请完成以下的程序
do  j=1,73
do  i=1,144
   ave1(i,j)=0
do  it=1,65
   ave1(i,j)=a(i,j,1,it)+ave1(i,j)

end  do
    ave1(i,j)= ave1(i,j)/65
    ave7(i,j)=a(i,j,1,61)-ave1(i,j)
end  do
end do

do  j=1,73
     ave8(1,j)=0

      do  i=1,144
      
      ave8(1,j)=a(i,j,1,61)+ave8(1,j)
      end  do
     ave8(1,j)=ave8(1,j)/144
      end do
do  j=1,73
do  i=1,144
    ave8(i,j)=a(i,j,1,61)-ave8(1,j)
end  do   
end do

do j=1,73
do i=1,144
write(4)ave1(i,j),ave7(i,j),ave8(i,j)

end do
end do
end
------------------------------------------------------------------------分割线
ctl   
dset d:\hgt5001.grd
undef -9.99E+33
title monthly mean hgt from the NCEP Reanalysis
xdef 144 linear 0.000 2.500
ydef 73 linear -90.000 2.500
zdef 1 levels 500
tdef 3 linear jan1948 1mo
vars 1
hgt 1 -999  monthly mean  hgt
endvars
---------------------------------------------------------------------------分割线
用grads  画图出现  contouring：-100 to 5800 interval  100  求 大神指教为毛图出不来Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟( ]6 \% w7 n) Y! j1 k. B

mofangbao

出现这个怎么了？

river

难道你不认识这几个单词么，查一下也清楚了，没必要来提问的吧。那不就是说 等值线 ：从-100到5800 间隔100画一条 。这有什么啊，正常的画等值线都会有这样的提示

xiaocaoqiqiao

river 发表于 2013-10-28 21:52

                               
登录/注册后可看大图

难道你不认识这几个单词么，查一下也清楚了，没必要来提问的吧。那不就是说 等值线 ：从-100到5800 间隔100 ...

画的是500hpa 一月历年平均位势高度场，指定年月位势距平场，纬偏场 正确的应该是contouring：5000 to 5800 interval 100。 用set  cmin  5000  也不管用，怀疑是数据的问题，500hpa位势高度场求和数字非常大，是不是存在数据溢出问题。个人瞎猜。

xiaocaoqiqiao

mofangbao 发表于 2013-10-28 21:47

                               
登录/注册后可看大图

出现这个怎么了？

画的是500hpa 一月历年平均位势高度场，指定年月位势距平场，纬偏场 正确的应该是contouring：5000 to 5800 interval 100。 用set  cmin  5000  也不管用，怀疑是数据的问题，500hpa位势高度场求和数字非常大，是不是存在数据溢出问题。个人瞎猜。

mofangbao

你写出之前先把数据输出来看看对不对 不要不知道数据算的对不对就开始画图哈

xiaocaoqiqiao

mofangbao 发表于 2013-10-28 22:11

                               
登录/注册后可看大图

你写出之前先把数据输出来看看对不对 不要不知道数据算的对不对就开始画图哈

生成的GRD文件用什么打开看啊 ？我打开都是乱码的

river

xiaocaoqiqiao 发表于 2013-10-28 22:09

                               
登录/注册后可看大图

画的是500hpa 一月历年平均位势高度场，指定年月位势距平场，纬偏场 正确的应该是contouring：5000 to 58 ...

contouring：5000 to 5800 interval 100这只能一定程度上反应你画的对不对，它只是描述一个事实告诉你现在你画的这个图是这么个情况，这句话没有正确不正确一说。

还有你要画的是两种图，一种是平均场，一种是距平场，你自己想一下有什么区别，不可能这三张画的时候都显示从5000到5800的，如果都显示类似这个值左右的话，平均场先不说，那两个距平场肯定是错的

river

xiaocaoqiqiao 发表于 2013-10-28 22:09

                               
登录/注册后可看大图

画的是500hpa 一月历年平均位势高度场，指定年月位势距平场，纬偏场 正确的应该是contouring：5000 to 58 ...

500毫巴的高度场求和以后当然会很大，5000米左右的高度本身就比较大了。不过就算你加个50年100年的也不存在你说的那个数据溢出，现在的电脑和程序什么的又不是几十年前的

xiaocaoqiqiao

river 发表于 2013-10-28 22:23

                               
登录/注册后可看大图

500毫巴的高度场求和以后当然会很大，5000米左右的高度本身就比较大了。不过就算你加个50年100年的也不存 ...

能不能帮我看看fortran程序哪里有错误呢？我都看一晚上了 ，也不知道哪有问题。我是菜鸟