爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7695|回复: 7

[脚本编辑] stnmap出不了.map,恳请帮助

[复制链接]

新浪微博达人勋

发表于 2016-6-17 22:08:51 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
我刚开始学习画图,在论坛里努力搜帖子学习,但是还是解决不了这个问题,恳请各位前辈不吝赐教。
这是我按照清风大大的pdf做的sta2grd.f
       Program sta2grd
      Implicit none!这里是程序的变量声明
        !变量生命结束
        !程序开始
      Character*8 stid
      Character*8 lon,lat,var,t24,tim
      integer nlev,flag
      tim =0.0
        nlev=1
        flag=1
      Open(1,file='e:\jwsksj\surface\t24-p\16012208.000',status='old')
      open(2,file='sta.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,var,t24
      write(2)stid,lat,lon,tim,nlev,flag,t24
      !            print*,stid,lon,lat,var,rain

        ! pause         
      Goto 10
      100 continue
      Close(1)
      nlev=0
      write(2)stid,lat,lon,tim,nlev,flag,t24
      close(2)
      !程序结束
      End
pdf里经纬度等是real,结果出的.grd里有乱码,就改成Character*8
出的sta.grd是这样的:
2.jpg
CTL是:
dset e:/jwsksj/surface/t24-p/sta.grd
dtype station
stnmap e:/jwsksj/surface/t24-p/sta.map
undef -999.9
title  2016012208 t24-p
tdef 1 linear 08:00z22Jan2016 24hr
vars 1
t24 0 99 t24-p data
endvars

在使用!stanmap时出现这个问题:
ga->!stnmap -i e:/jwsksj/surface/t24-p/sta.ctl
    Name of binary data set:e:/jwsksj/surface/t24-p/sta.grd
    Number of times in the data set:1
    Number of surface variables:1
    Number of level dependent variables:0


Starting scan of station data binary file.
Binary data file open:e:/jwsksj/surface/t24-p/sta.grd


Processing time step 1
    Invalid station hdr found in station binary file
    Possible causes:Invalid level count in hdr
                            Descriptor file mismatch
                            File not station data
                            Invalid relative time
        levs = 538981684    flag = 0   time = 4.165e-11

1.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-6-18 20:42:10 | 显示全部楼层

回帖奖励 +2 金钱

经纬度等是real,结果出的.grd不能用文本方式的打开。用文本软件打开就是乱码。
密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2016-6-18 21:02:15 | 显示全部楼层
应该写入站号出现错误了: stid(i)=char(i) 试试
                                   
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-6-21 15:03:39 | 显示全部楼层
霹雳电光鼠 发表于 2016-6-18 20:42
经纬度等是real,结果出的.grd不能用文本方式的打开。用文本软件打开就是乱码。

谢谢你,感激不尽!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-6-21 15:04:10 | 显示全部楼层
小小书童 发表于 2016-6-18 21:02
应该写入站号出现错误了: stid(i)=char(i) 试试

谢谢你,我会好好研究的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-7-18 18:45:35 | 显示全部楼层
请问楼主问题解决了么?怎么改的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-7-31 16:49:57 | 显示全部楼层
279995024 发表于 2017-7-18 18:45
请问楼主问题解决了么?怎么改的

问题是解决了,但也不是大问题,第一个问题就像那些高手指教的,不要用文本文件打开,其实也不需要打开,多此一举,还把自己整糊涂了,第二个问题清风大神说是正常现象,不用管。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-1-4 14:14:06 | 显示全部楼层
想问楼主第二个问题是怎么解决的?就是黑框中提示的错误,map文件生不出来也不能不理会吧
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表