站点数据处理ctl报错

hello-lyn

数据是9个站35年月平均相对湿度资料，由于不是一般的站点存放格式，自己做了一些处理，把站点信息（站号及经纬度）和相对湿度数据分开存放了（fortran中的rh.txt和stid.txt）
fortran读取资料：

Program sta2grd
parameter(nx=9,tt=420)
character*8 stid(9),lon(9),lat(9)
real rh(9,420),tim
integer nlev,flag
open(3,file='d:\0512\RH.txt')
open(4,file='d:\0512\stid.txt')
open(7,file='d:\0512\RH.grd',form='binary')
do i=1,nx
do j=1,tt
   read(3,*)rh(i,j)
   print*,rh(i,j)
enddo
enddo
do i=1,nx
read(4,*)stid(i),lat(i),lon(i)
enddo
do i=1,tt
tim=0.0
nlev=1
flag=1
do j=1,nx
write(7)stid(j),lat(j),lon(j),tim,nlev,flag,rh(I,J)
enddo
nlev=0
write(7)stid(i-1),lat(i-1),lon(i-1),tim,nlev,flag
enddo
close(7)
end

程序运行ok生成了相应的grd

ctl：
dset d:\0512\rh.grd
dtype station
stnmap d:\0512\rh.map
undef 32766
title relative humidity
tdef 420 linear JAN1979 1mo
VARS 1
rh 0 99 rhSTN data
endvars
生成站点映射文件时
打开ctl就显示open error missing xdef/ydef/zdef/tdef record。。。
看了之前一些帖子也没有什么头绪，由于刚开始做站点数据处理，不是很明白，又把这个问题放上来求助！

sun_shine_Xia

做站点映射map文件的时候不是应该用！stnmap么？不是open啊

lqouc

lon，lat要用real，修改一下再试试。

lqouc

sun_shine_Xia 发表于 2014-5-21 11:31
做站点映射map文件的时候不是应该用！stnmap么？不是open啊

可能他说的是报错的信息，不是命令。

river

你怎么操作出现你说的那个错误提示的？问题说不清楚，让别人猜还行！

sun_shine_Xia

lqouc 发表于 2014-5-21 12:27
可能他说的是报错的信息，不是命令。

搞不懂。。只有等他试试再看是啥问题了

lqouc

sun_shine_Xia 发表于 2014-5-21 13:12
搞不懂。。只有等他试试再看是啥问题了

嘿嘿，其实我一直想知道你头像的菊花旁边到底写的什么字，真真看不清。

sun_shine_Xia

lqouc 发表于 2014-5-21 13:14
嘿嘿，其实我一直想知道你头像的菊花旁边到底写的什么字，真真看不清。

为啥会对那朵菊花旁边的字感兴趣呢。。。花才是重点

sun_shine_Xia

sun_shine_Xia 发表于 2014-5-21 13:20
为啥会对那朵菊花旁边的字感兴趣呢。。。花才是重点

其实我也忘了。。。好久以前的图片，好像是“我希望躺在向日葵上，即使劳累也可以向着阳光”！

hello-lyn

sun_shine_Xia 发表于 2014-5-21 11:31
做站点映射map文件的时候不是应该用！stnmap么？不是open啊

是用的! stnmap命令，但显示就是open error的错误。