如何把两个数据文件画到一张图中？用了两个d但是等值线出不来

周扬澄

Fortran部分程序：

do j=1,n
  do i=1,m
  dp(i,j)=elpj(i,j)-lapj(i,j)
  end do
end do

do j=1,n
  do i=1,m
  tt(i,j)=(elpj(i,j)-lapj(i,j))/sqrt((3*sel(i,j)**2+4*sla(i,j)**2)/7*(1/5+1/4))
  end do
end do

open(2,file='d:/dp.grd',form='binary')
write(2)((dp(i,j),i=1,m),j=1,n)

write(2)((tt(i,j),i=1,m),j=1,n)
close(2)

end

ctl文件
dset D:/dp.grd
undef  1.0E+33
title  dp data
xdef  48  levels   5          12.5            20          27.5            35          42.5            50          57.5            65          72.5            80          87.5            95         102.5           110         117.5           125         132.5           140         147.5           155         162.5           170         177.5           185         192.5           200         207.5           215         222.5           230         237.5           245         252.5           260         267.5           275         282.5           290         297.5           305         312.5           320         327.5           335         342.5           350        356.25
ydef  24  levels   -87.5          -80        -72.5          -65        -57.5          -50        -42.5          -35        -27.5          -20        -12.5           -5          2.5           10         17.5           25         32.5           40         47.5           55         62.5           70         77.5           85
zdef   1  levels   1000
tdef  30  linear  0z01Feb1978  1yr
vars   2
dp   0   99  different
tt   0   99  tt
endvars

gs文件
'reinit'
'enable print D:\dp.gmf'
'open D:\dp.ctl'
'set grads off'
'set lon 0 360'
'set lat -90 90'
'set t 1'
'set gxout shaded'
'set clab on'
'set string 1 c 2 0'
'set strsiz 0.3 0.6'
'set csmooth on'
'd dp'
'set gxout contour'
'set clevs 2.37'
'd tt'
'draw title slp dp '
'print'
'disable print'
；
另外，我尝试了把数据放在两个文件，写两个ctl把变量d dp.1和d dp.2还是不行，希望有好心人可以教教我

yiyiyi

请问解决了吗

coolcode

请问解决了么

zlh

楼主请问这个问题解决了吗，我刚好也遇到这个问题