GrADS显示图形问题

windys

我用fortran将一个500hPa高度场资料，（资料以ASCII码形式存放，文件名h500.dat，范围：60~150E，0~40N.时段：1982.1~1985.12共48个月。水平分辨率：2.5*2.5，格点数：37*17）转换成二进制数据文件h500.grd,再写了一个ctl 文件，然后用GrADS画出其中一个时次的高度场图，可是出来的图很奇怪，图线都挤成一堆了，要不就是很稀疏。。。不明白什么原因，求大神指教！

下面是fortran程序：
program turn
  implicit none
  integer nx,ny,nt,i,j,it
  parameter(nx=37,ny=17,nt=48)
  real h(nx,ny,nt)
  open(1,file='F:\statistics\practice1\h500.dat')
  open(2,file='F:\statistics\practice1\h500.grd',form='binary')
  do 10 it=1,nt
  read(1,*)((h(i,j,it),i=1,nx),j=1,ny)
10 continue
  do 100 it=1,nt
  write(2)((h(i,j,it),i=1,nx),j=1,ny)
100 continue
  end


然后是h500.grd的ctl文件
dset F:\statistics\practice1\h500.grd
undef -9.99E+33
title NCEP/NCAR REANALYSIS PROJECT
xdef   37 linear   60.0  2.5
ydef   17 linear   0.0  2.5
zdef    1 levels   500
tdef 48 linear   JAN1982   1mo
vars  1
h  1 99 h500
endvars


最后是绘图的gs文件：

'reinit'
'open F:\statistics\practice1\h500.ctl'
'set t 1'
'd h'
;

出来就是这个鬼样。。。图线都挤在左下角啦，左上角和右上角各有一根。。。



求指教!是我的gs文件设置有问题吗？

男紫汗

等值线间隔设个4试试。。。ascii格式这么读就可以吗？

windys

不是，这原本是个ASCII码格式的资料没错，但我已经用fortran把它转换成二进制的了，但是用grads显示就出现那种问题了，等值线间隔设成4也不行。。。

windys

男紫汗 发表于 2017-3-29 20:37
等值线间隔设个4试试。。。ascii格式这么读就可以吗？

不是，这原本是个ASCII码格式的资料没错，但我已经用fortran把它转换成二进制的了，但是用grads显示就出现那种问题了，等值线间隔设成4也不行。。。

男紫汗

windys 发表于 2017-3-30 19:49
不是，这原本是个ASCII码格式的资料没错，但我已经用fortran把它转换成二进制的了，但是用grads显示就出 ...

fortran里print出来看看数据对不对