关于！stnmap出错，不能生成map文件

黄鸿惠

fortran小白最近学习使用站点数据插值，参考了清风的帖子，一步步做的。用清风提供的raindata.txt文件进行学习时都没有问题的。但是用自己的数据处理时，在用！stnmap -i xxxx.ctl 生成map文件时卡住了【GrADS站点资料作图详细解决方案http://bbs.06climate.com/forum.p ... &fromuid=107364】


我用的是diamond3的6小时降水量资料 大概是这样的
读取的部分是从第15行开始的 依次是站点号 经度 纬度 海拔高度 6小时降水量

                               
登录/注册后可看大图

******************************************************************
以下是我读取数据的fortran程序：
Program sta2grd
Implicit none
!这里是程序的变量声明
Character*8 stid      !stid是站点  是字符型，长度为8
character*6 rain
Real lon,lat,alti,tim  !经纬度，alti海拔高度，降水，时次这些变量名声明它们是单精度实型
integer nlev,flag     !总层次，有无地面资料
!变量声明结束

!程序开始
tim=0.0
nlev=1
flag=1
!=======================打开文件====读取数据=======================!
Open(1,file='e:/bian/jiayuan/grads/chazhi/rumen/07061416.000',status='old')   
Open(2,file='e:/bian/jiayuan/grads/chazhi/rumen/07061416.grd',status='replace',form='binary')
read(1,*)   
read(1,*)   
read(1,*)
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)   
read(1,*)     
10 Read(1,*,end=100)stid,lon,lat,alti,rain  
!=========用于测试输出所读取的数据
!Print*,stid,lon,lat,alti,rain   
!pause   
!==========写入==========
write(2)stid,lat,lon,tim,nlev,flag,rain
Goto 10
100 continue
Close(1)
nlev=0  
write(2)stid,lat,lon,tim,nlev,flag
close(2)
!程序结束
End


********************************************
下面是写的ctl
dset e:\bian\jiayuan\grads\chazhi\rumen\07061416.grd
dtype station
stnmap e:\bian\jiayuan\grads\chazhi\rumen\07061416.map
title rain data sample
undef -32766
tdef  1 linear 16z14jun2007 1hr
vars 1
rain 0 99 Rain data
endvars


*************************************************
接下来是在grads中运行！stnmap -i xxxx.ctl 生成map文件  显示如下报错，想问各位大神能看出是那里出问题了吗？
【我试过将ctl文件里面的斜杠‘/’  '\'的表达都试了一遍 都不行 都出现同样的报错  感觉上是数据读取写入的时候有问题 但也不会改 故发出求助】

                               
登录/注册后可看大图

黄鸿惠

不知道为啥贴不了图  我这里补充一下报错的内容

h:\error.png

lqouc

rain应为real类型

黄鸿惠

lqouc 发表于 2022-1-5 09:46
rain应为real类型

之前这么写是因为看diamond3的说明文件说是字符串就把rain定义为character 没想到是这里出问题 后来把rain定义成real 果然问题解决 实在是太感谢您了！