新手求助！站点映射文件报错

hello-lyn

stn160.dat (3.05 KB, 下载次数: 0)

2014-2-26 19:05 上传

点击文件名下载附件

r3.txt (49.41 KB, 下载次数: 0)

2014-2-26 19:05 上传

点击文件名下载附件

在学习站点资料处理时，按照书上说的方法做了一次不对，后来又看了版里那个教程之后修改了一下程序，程序运行无错，也生成了相应的grd文件，但在运行grads生成映射文件时总是报错，看了以前有相似问题的帖子修改后还是无果，新手能力有限，实在看不出问题了。。特来求教！

数据是国家气候中心160站61年2月的降水资料 每个站点为：站号、纬度、经度、tim nlev nflag 降水

fortran程序如下：
program read_160
         parameter (n=160,m=61)
         dimension r0(n,m)
   open(3,file='d:\rain-160\r3.txt')  
   read(3,*) ((r0(i,j),i=1,n),j=1,m)
   call stntogrd(r0(n,m))
  end

subroutine stntogrd(x)
real lat(160),lon(160),x(160,61)
character*8 stid(160)
open(4,file='d:\rain-160\stn160.dat')
do k=1,160
read(4,*) stid(k),lat(k),lon(k)
enddo
close(4)
open(5,file='d:\rain-160\stn-160.grd',form='binary')
do k =1,61
tim=0.0
nlev=1
nflag=1
do i=1,160
write(5) stid(i),lat(i),lon(i),tim,nlev,nflag,x(i,k)
enddo
nlev=0
write(5) stid,lat,lon,tim,nlev,nflag
enddo
close(5)
return        
end

ctl：
dset  d:\rain-160\stn-160.grd
dtype   STATION
stnmap d:\rain\stn160-1.map
undef   -999.0
title   stn160
tdef    61 linear feb1951 1yr
vars    1
p    0  99  Rainfall data
endvars

报错：
file:///C:/Users/Administrator/AppData/Roaming/Tencent/Users/821101349/QQ/WinTemp/RichOle/%7BMZ)@]RNPU(0GQ%7B2GJQ0%60JN.jpg


补充说明：之前运行了fortran之后特意生成一个dat文件查看了一下，结果发现写入的最后一列数据就是降水那一列貌似有点问题，怀疑是这里的问题。。。
[img]file:///C:/Users/Administrator/AppData/Roaming/Tencent/Users/821101349/QQ/WinTemp/RichOle/%7BXE6$%60K[$%7B~E1L(0RUQAAOK.jpg[/img]
当然还是希望大家帮忙指教~

hello-lyn

楼上的截图没有出来，接着上图。。

做个霸气的木头

write(5) stid,lat,lon,tim,nlev,nflag
enddo
close(5)
return        
end
这几行的第一行命令貌似不对，应当是write(5) stid(i-1),lat(i-1),lon(i-1),nlev,nflag
是标注每一个时刻结束的最后一条记录，这里i-1=160.

做个霸气的木头

write(5) stid,lat,lon,tim,nlev,nflag
enddo
close(5)
return        
end
这几行的第一行命令貌似不对，应当是write(5) stid(i-1),lat(i-1),lon(i-1),nlev,nflag
是标注每一个时刻结束的最后一条记录，这里i-1=160.

hello-lyn

做个霸气的木头 发表于 2014-2-26 23:10

                               
登录/注册后可看大图

write(5) stid,lat,lon,tim,nlev,nflag
enddo
close(5)

映射文件出来了，又出现了新的问题。。就是运行画图的gs又只能出来地图，感觉这个错误好低级啊、、、看了其它帖子改了还是一样。。请版主指点！
grid.ctl
dset  d:\rain-160\grid.grd
undef  -9.99e33
title   sample grid data
xdef    71 linear 70 1
ydef    41 linear 15 1
zdef    1  linear 850 1
tdef    61 linear feb1951 1yr
vars    1
g    0  99  grid data
endvars


gs：'open d:\rain-160\grid.ctl'
'open d:\rain-160\stn160-1.ctl'
'set grads off'
'enable print d:\rain-160\stn-160.gmf'
'set lon 73 135'
'set lat 15 55'
'set t 61'
'set mpdset hires cnworld'
'define a=oacres(g,p.2,1.5)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set gxout shaded'
'd aa'
'set grid off'
'set gxout contour'
'set clab forced'
'set black -0.2976 0.2976'
'set cint 0.1'
'd aa'
'print'
'disable print'
'reinit'
;
C:\Users\Administrator\Desktop
C:\Users\Administrator\Desktop

hello-lyn

[img]C:\Users\Administrator\Desktop\VRMZKUZ)}]X57[%R6]U%I$H.jpg[/img]
C:\Users\Administrator\Desktop\}M$ZSGULW8RK7ILX(806F2Q.jpg