爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5814|回复: 14

[脚本编辑] grads出图错误

[复制链接]

新浪微博达人勋

发表于 2018-6-25 12:40:24 | 显示全部楼层 |阅读模式

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

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

x
捕获.PNG
grads里先将站点数据转为二进制数据,之后运行时出现了这个错误,不知道是什么原因,求指导啊
这是ctl
dset F:\VIS\VIS-PIC\VIS-1980-1989\vis-1980-1989.grd
dtype station
stnmap F:\VIS\VIS-PIC\VIS-1980-1989\AH.map
undef -9999
title rain data sample
tdef 1 linear 01may1980 1mo
vars 1
day 0 99 rain data
endvars

这是gs
'reinit'
'open F:\VIS\VIS-PIC\VIS-1980-1989\AH.ctl'  
'open F:\VIS\VIS-PIC\VIS-1980-1989\first.ctl'
'set mpdset ahs'
'set lon 114.5 123';'set lat 27 35.5'
'define a=oacres(mask.1,day.2,50)'
'define a1=maskout(a,mask.1-1.0)'
'define aa=smth9(a1)'
'set grads off'
'set grid off'
'set xlopts 1 6 0.18';'set ylopts 1 6 0.18'
'set xlint 2';'set ylint 2'


'set rgb 20 255 244 197'
'set rgb 21 255 227 101'
'set rgb 22 255 198 42'
'set rgb 23 255 167 42'
'set rgb 24 255 144 42'
'set rgb 25 255 100 42'
'set rgb 26 255 74 0'
'set rgb 27 227 32 23'
'set rgb 28 194 23 15'

'set rbcols 20 21 22 23 24 25 26 27 28 '

'set clevs 15 20 25 30 35 40 45 50 '
'set gxout shaded'
'd aa'
'cbar 1.0 1'
'printim F:\VIS\VIS-PIC\VIS-1980-1989\AH-1980-1989.png white'

这是转换数据的fortran
Character*8 stid
Real lat,lon,day,tim
integer nlev,flag
!变量声明结束
!程序开始
tim=0.0
nlev=1
flag=1
PRINT*,111
Open(1,file='F:\VIS\VIS-PIC\VIS-1980-1989\vis-1980-1989.txt',status='old')
PRINT*,111
open(2,file='F:\VIS\VIS-PIC\VIS-1980-1989\vis-1980-1989.grd',status='replace',form='binary')
10 Read(1,*,end=100)stid,lat,lon,day
   print*,stid,lat,lon,day
Write(2)stid,lat,lon,tim,nlev,flag,day
!write(2)stid,lat,lon,day
Goto 10
100 continue
Close(1)
nlev=0
Write(2)stid,lat,lon,tim,nlev,flag
close(2)
!程序结束
End

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

新浪微博达人勋

发表于 2018-6-25 12:42:43 | 显示全部楼层
把你的ctl重新另存一下成gs,要Unicode编码的试试,可能是编码有问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-25 20:25:18 | 显示全部楼层
river 发表于 2018-6-25 12:42
把你的ctl重新另存一下成gs,要Unicode编码的试试,可能是编码有问题

你好,有别的方法吗,不是很懂你的方法
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-25 20:25:27 | 显示全部楼层
river 发表于 2018-6-25 12:42
把你的ctl重新另存一下成gs,要Unicode编码的试试,可能是编码有问题

你好,有别的方法吗,不是很懂你的方法
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-25 20:25:34 | 显示全部楼层
river 发表于 2018-6-25 12:42
把你的ctl重新另存一下成gs,要Unicode编码的试试,可能是编码有问题

你好,有别的方法吗,不是很懂你的方法
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-6-25 21:00:33 | 显示全部楼层
believefuture 发表于 2018-6-25 20:25
你好,有别的方法吗,不是很懂你的方法

那就新建一个txt文件,把后缀名改成gs ,然后用记事本打开。把你原来的脚本复制,然后粘贴在新建的那个文件里,不是直接粘贴,选择那个只粘贴文本那个试试,然后保存。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-25 21:21:24 | 显示全部楼层
river 发表于 2018-6-25 21:00
那就新建一个txt文件,把后缀名改成gs ,然后用记事本打开。把你原来的脚本复制,然后粘贴在新建的那个文 ...

你好,我准备将原来的脚本复制在你说的新建的gs文件里,可是没有只粘贴文本
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-25 21:24:18 | 显示全部楼层
believefuture 发表于 2018-6-25 21:21
你好,我准备将原来的脚本复制在你说的新建的gs文件里,可是没有只粘贴文本

只有粘贴选项
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-6-26 09:12:51 | 显示全部楼层

哦,那你直接粘贴吧,然后再运行试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-26 09:58:56 | 显示全部楼层
river 发表于 2018-6-26 09:12
哦,那你直接粘贴吧,然后再运行试试

直接粘贴还是不行,奇怪的是我用个gs去运行我其他的数据,图却能出来。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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