请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 11475|回复: 13

[分享资料] grads绘图中遇到的问题

[复制链接]

新浪微博达人勋

发表于 2012-6-30 19:55:47 | 显示全部楼层 |阅读模式

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

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

x
数据文件为      文件1.   40 个站的 站点信息、经度、纬度           文件2        两列数据 1列为最高气温,1列为最低气温
Fortran, Ctl, gs 文件如下, grads中提示的错误信息如下图所示,不胜感激~


Fortran 程序
      PARAMETER(num=40,ks=2)
c    num=stn,ks=nt it=2(1=eof;2=reof)
     DIMENSIONxlon(num),xlat(num),vis(num,ks)
      CHARACTERSTA(num)*8,FL*35,F5*5,F2*2,FL1*35
       FL='c:\abc\40stn.txt'
     OPEN (11,FILE=FL,STATUS='OLD')
       Fl='win-ave.txt';FL1='win-ave.grd'
     OPEN (1,FILE=FL,status='old')
     open(3,file=fl1,form='binary')
     DO  i=1,num
     READ(11,*) ist,xlat(i),xlon(i)
     WRITE(sta(i),'(I8.8)')ist
      read(1,*)(vis(i,j),j=1,ks)
     end do
     close(11)
   
      close(1)
     
     tim=0.0
      nlev=1
      nflag=1
      do i=1,num
      write(3)sta(i),xlat(i),xlon(i),tim,nlev,nflag,(vis(i,k),k=1,ks)
      enddo
      
       nlev=0
      write(3)sta(num),xlat(num),xlon(num),tim,nlev,nflag
      close(3)
     STOP 999
     END
Ctl
dset C:\abc\win-ave.grd
dtype station
stnmap C:\abc\zj.map
undef 3276.6
*undef 1.0
title average
tdef 1 linear 01Jan1979 1yr
vars 2
V1 0 99 eigenvalue 1
V2 0 99 eigenvalue 2
endvars
Gs
'REINIT'
'enable print C:\abc\win-ave.gmf'
i=1
while(i<=2)
'set mpdset cnworld'
'set map 1'
'OPEN C:\abc\win-ave.ctl';a0='10 7';a1='1.57'
'open C:\abc\prc5.ctl'
'set parea 1. 10.5 1.0 8.'
'set lon 108 119';'set lat 20 27'
'set grid off';'set grads off'
'set dignum 0'
'SET DIGSIZ 0.1'
'set xlopts 1 5 0'
'set ylopts 1 5 0'
'set clab forced'
*'set clevs -28.0 -24.0 -20.0 -16.0 -12.0-8.0 -4.0 4.0 8.0 12.0 16.0 20.0 24.0 28.0'
'set ccolor 1'
'set clopts 1 5 0.14';'set cthick 5';'setcint 1'
'Dsmth9(MASKOUT(OACRES(prc.2(T=1),v'i',20,16,12,8,4,2),prc.2(T=1)-0.5))'
*'xylab2 2 1 0.16'
if(i=1);fs='a';endif
if(i=2);fs='b';endif
'set string 1 c 6'
'set strsiz 0.35';'draw string 1.7 7.5`1'fs
'print'
'reset'
'c'
'close 2'
'close 1'
i=i+1
endwhile
'disable print'

QQ截图20120630121746.gif
QQ截图20120630121808.gif


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

新浪微博达人勋

发表于 2012-6-30 22:32:38 来自手机 | 显示全部楼层
论坛里有好多站点数据处理的帖子,你先看看吧~根据报的错一点点改……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-7-1 05:19:30 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2012-7-1 09:56:03 | 显示全部楼层
ctl~文件有问题,好好看基础知识吧通过自己找现有知识完全能解决,你现在的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-1 10:15:05 来自手机 | 显示全部楼层
り。3角戀※ 发表于 2012-7-1 05:19  请问no hardcopy metafile open 是怎么回事啊

这不是报错,因为gs头一句是reinit造成的,没有问题。我是手机上的,不方便帮你看,你gs有一个明显的错误,打开文件的顺序反了,应该先打开用于插值的格点数据,再打开站点数据,,其它还有错的话,你可以按照这样方法检查:1、先检查你编程生成的站点数据有没有错,画画站点图,图上有非0数据说明没错,有错的话,按照我写过的一个贴子里的错误提醒修改,2、第一步完成,再有问题,你就把你gs里面修饰图形的部分都去掉,一行行的检查错误!再有问题的话,上来问
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-3 21:30:36 | 显示全部楼层
学习了,好多高人在这里。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-8-26 16:49:55 | 显示全部楼层
高手!{:5_231:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2016-5-16 16:14:05 | 显示全部楼层
{:5_235:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-1-5 10:09:53 | 显示全部楼层
写的很好
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-1-5 13:20:03 | 显示全部楼层
随便看看,学习学习
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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