爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5066|回复: 19

为什么图是这样的?

[复制链接]

新浪微博达人勋

发表于 2013-11-10 21:08:40 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 男紫汗 于 2013-11-11 11:05 编辑

箭头长度好长。。。。求救各位大神

gs:
'reinit'
'set mapdraw off'
'open u3.ctl'
'open v3.ctl'
'set dfile 2'
'set lon 0 357.5'
'set lat -90 90'
'd u.1;v.2'
;
u3.ctl:
DSET   ^u3.dat
OPTIONS sequential
UNDEF  -2.56E33
TITLE 5 Days of Sample Model Output
XDEF 144 LINEAR  0.0 2.5
YDEF 73 LINEAR  -90.0 2.5
ZDEF 1 LEVELS 1000
TDEF 42 LINEAR 02JAN1960 1yr
VARS 1
u      1   99   Temperature
ENDVARS
Fortran程序:program ssimplicit none
integer i,j
real(8)::u1(73,144),u2(73,144)=0,v1(73,144),v2(73,144)=0,u3(73,144)=0,v3(73,144)=0
OPEN(11,FILE='vwind.dat')
do i=1,73
do j=1,144
  read(11,"(E21.14)") v1(i,j)

  enddo
  enddo
  close(11)
open(20,file='uwind.dat')
do i=1,73
do j=1,144
  read(20,"(E21.14)") u1(i,j)

  enddo
enddo
close(20)
print*,u1(1,1)
print*,u1(73,144)

do j=1,144
do i=1,73
if(2.le.j.and.j.le.143.and.2.le.i.and.i.le.72) then
  u2(i,j)=u1(i,j)/2.+(u1(i,j-1)+u1(i,j+1)+u1(i-1,j)+u1(i+1,j))/8.
   v2(i,j)=v1(i,j)/2.+(v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+v1(i+1,j))/8.
else  
u2(i,j)=u1(i,j)
v2(i,j)=v1(i,j)
   endif
  enddo
enddo
u3=u1-u2
v3=v1-v2

open(12,file='u3.dat',form='unformatted',access='SEQUENTIAL')
open(13,file='v3.dat',form='unformatted',access='SEQUENTIAL')
open(14,file='u2.dat',form='unformatted',access='SEQUENTIAL')
open(15,file='v2.dat',form='unformatted',access='SEQUENTIAL')
open(16,file='u1.dat',form='unformatted',access='SEQUENTIAL')
open(17,file='v1.dat',form='unformatted',access='SEQUENTIAL')
do i=1,144
do j=1,73
  write(12) u3(j,i)
  write(13) v3(j,i)
  write(14) u2(j,i)
  write(15) v2(j,i)
  write(16) u1(j,i)
  write(17) v1(j,i)
   enddo
   enddo
end
这个有问题吗?


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

新浪微博达人勋

发表于 2013-11-10 21:14:18 | 显示全部楼层
上gs文件吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-10 21:36:13 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 19710
发表于 2013-11-11 07:53:25 | 显示全部楼层
缺省值弄错
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-11 08:45:37 | 显示全部楼层
你的数据里面出现了NAN~~~看看缺省值吧~~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-11 09:57:11 | 显示全部楼层
男紫汗 发表于 2013-11-10 21:36
感觉是ctl的问题

UNDEF  -2.56E33
缺测值有问题,用fortran读出来资料看一下缺测值是多少,改一下ctl
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-11 10:42:46 | 显示全部楼层
andrewsoong 发表于 2013-11-11 08:45
你的数据里面出现了NAN~~~看看缺省值吧~~~

数据我输到屏幕上,没出现NAN啊。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-11 10:44:01 | 显示全部楼层
兰溪之水 发表于 2013-11-11 07:53
缺省值弄错

能不能具体一点啊~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-11 10:55:13 | 显示全部楼层
river 发表于 2013-11-11 09:57
UNDEF  -2.56E33
缺测值有问题,用fortran读出来资料看一下缺测值是多少,改一下ctl

全球风速场资料,怎会有缺测值
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-11 11:19:09 | 显示全部楼层
男紫汗 发表于 2013-11-11 10:42
数据我输到屏幕上,没出现NAN啊。。。

那你是怎么设置缺省值的???
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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