grads绘制风速风向图

breeze_li · 发表于 2016-9-5 14:51:03

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

用fortran+grads绘制风速风向图该怎么做，流程大概是怎样的？哪里有类似例子

breeze_li · 发表于 2016-9-5 15:05:12

哪位给个流程也行，主要是没接触过，不知道从何下手

lqouc · 发表于 2016-9-5 15:51:28

啥叫风速风向图？给个示例图看看。

breeze_li · 发表于 2016-9-5 15:55:47

lqouc 发表于 2016-9-5 15:51
啥叫风速风向图？给个示例图看看。

我不是学气象的，刚刚接触，所以说的不准确。。就下面这种图
LN@PWW%XD`TS~0~VU1@Y.jpg

lqouc · 发表于 2016-9-5 16:10:38

先看这个帖子
http://bbs.06climate.com/forum.p ... digest%26digest%3D1
把站点数据做出来之后，去grads官网看看 set gxout model的用法就行了

breeze_li · 发表于 2016-9-5 16:35:14

lqouc 发表于 2016-9-5 16:10
先看这个帖子
http://bbs.06climate.com/forum.p ... digest%26digest%3D1
把站点数据做出来之后，去grad ...

这个之前有看过，弄过降水跟温度的图。风速与风向在用fortran生成.grd需要怎么处理，还是直接用
这是我的站点数据的一部分：
53949 109.3167 36.8833 1.3 180
53937 110.0833 34.4833 3.7 120
53927 110.1833 36.0667 9.8 90
这是fortran：

Program sta2grd
Implicit none
!这里是程序的变量声明
Character*8 stid
Real lon,lat,vwind,uwind,tim
integer nlev,flag
!变量声明结束
!程序开始
tim=0.0
nlev=1
flag=1
Open(1,file='wind1.txt',status='old')
open(2,file='wind1.grd',status='replace',form='binary')
read(1,*)
10 Read(1,*,end=100)stid,lon,lat,vwind,uwind
Write(2)stid,lat,lon,tim,nlev,flag,vwind,uwind
Goto 10
100 continue
Close(1)
nlev=0
Write(2)stid,lat,lon,tim,nlev,flag
close(2)
!程序结束
End

这是ctl：
DSET wind1.grd
DTYPE  station
STNMAP wind.map
UNDEF  -9.99E33
TITLE  wind Data Sample
TDEF 1 linear 01may2011 1mo
VARS 2
v 1  99  vwind
u 1  99  uwind
ENDVARS

在生成.map文件时出错：
TNF$@2@@UAU5}G@GRF6ILTF.png

我这个应该怎么改？

lqouc · 发表于 2016-9-5 22:20:37

代码应该都没啥问题，报错估计是你的数据有问题，自己检查下吧。
至于风速风向怎么用我觉得你一定是没有查看我告诉你的命令。这个是要换成uv的，具体看官网说明
http://cola.gmu.edu/grads/gadoc/usingstationdata.html#operating

lizhidong2016 · 发表于 2017-1-25 13:35:23

感谢lz很无私，很好的资源

		自动登录	找回密码
密码			立即注册

[图形美化] grads绘制风速风向图

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