grads循环读取文件

小浣熊

我想用T639资料画图，因为t639资料包括00时起报的资料的预报时段为0-240小时，且时间间隔不一样。比如2014年2月21日起报的资料为2014022100003、2014022100006,....。现在我想批量画出这个时效的所有时间间隔的图，写了一个脚本，但是提示找不到文件。我把我所有的时间写在了name.txt里。这是我的gs文件：
parameter(m=38)
character time1(m)*13
character time2*13
open(1,file='E:\ZLCL\name.txt',status='unknow')
do i=1,m
read(1,100) time1(i,j)
time2=time1(i,j)                     
'open E:\ZLCL\time2.ctl'
'reint'
'enable print E:\ZLCL\time2.gif gif'
'set lon 60 140'
'set lat 10 60'
'set grads off'
'set grid off'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set mpdset cnworld '
'set map 1 1 6'
'draw map'
'set mpdset cnbase '
'set map 1 1 4'
'draw map'
'set gxout barb'
'set lev 700'
'set gxout shaded'
'set cmin 12'
'd mag(ugrdprs,vgrdprs)'
'd skip(UGRdprs,3);skip(VGRdprs,3)'
'cbarn 0.8 1'
'set gxout contour'
*'set cterp on'
*'set csmooth on'
'set lev 500'
'set clevs 520 524 528 532 536 540 544 548 552 556 560 564 568 572 576 580 584 588 592'
'set ccolor 4'
'set cthick 10'
'set cmin 5'
'd HGTprs/10'
'draw title 2013-10-21-06 500hPa HGT 700hpa uv'
'print'
'printim E:\ZLCL\time2.gif gif white '
'disable print'
end do
100 format(a13)
运行以后提示的错误是‘can't open description file'.grads是不是不能直接打开包含变量名的文件？

river

grads的编程语言和fortran有很多不同的地方，你这种fortran和grads杂交的东西在grads里当然是运行不通的。

小浣熊

river 发表于 2014-2-27 17:15

                               
登录/注册后可看大图

grads的编程语言和fortran有很多不同的地方，你这种fortran和grads杂交的东西在grads里当然是运行不通的。

求指导，新手上路，直接混了，都是边做边试

斥鷃

小浣熊 发表于 2014-2-27 17:43

                               
登录/注册后可看大图

求指导，新手上路，直接混了，都是边做边试

搜《GRADS实用绘图》，倒数第二部分还是哪个是关于脚本和函数使用的，看看吧。

mofangbao

擦。。。楼主应该好好分别看看fortran和grads了。。。这代码

又是那隻貓

牛比啊

又是那隻貓

牛比啊

clare

保持队形。牛比啊

冰枫

楼主这个问题解决了吗？目前也面临着批处理时间间隔不一致的文件问题

dyldcx

学习学习~~