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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3126|回复: 2

[图形美化] 参照清风的多时次站点资料处理出现的一个小BUG,求明眼人指点一二。

[复制链接]

新浪微博达人勋

发表于 2016-2-24 13:21:03 | 显示全部楼层 |阅读模式

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

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

x
(不知道为什么选择脚本编辑发不了帖)一切顺利  能出图
因为我并不需要插值图片所以只进行到“!stnmap -i ********.ctl”这一步。

问题在于生成map这一步 GrADS似乎并没有把整个过程走完。详见图1.


但是奇怪的是这并没有影响我作图,见图2.



献上我的数据(附件1,快照见图3)


处理数据的fortran如下:
program main
implicit none
integer,parameter::nt=168,nn=11
real lat(nn),lon(nn),rain(nt,nn),tim
real u(nt,nn),v(nt,nn)
character*8 stan(nn)
integer nlev,nflag,i,k

open(9,file='E:\SCI2\datstation\UVsurface_allstation-2012.txt',status='old')
do i=1,nt
do k=1,nn
read(9,*)stan(k),lat(k),lon(k),u(i,k),v(i,k)
print*,stan(k),lat(k),lon(k),u(i,k),v(i,k)
enddo
print*,"*"
enddo
close(9)
open(12,file='E:\SCI2\datstation\obs_allstation_V_2012.grd',status='replace',form='binary')
do i=1,nt
tim=0.0
nlev=1
nflag=1
do k=1,nn
write(12) stan(k),lat(k),lon(k),tim,nlev,nflag,v(i,k)
print*,stan(k),lat(k),lon(k),tim,nlev,nflag,v(i,k)
enddo
nlev=0
write(12) stan(k-1),lat(k-1),lon(k-1),tim,nlev,nflag
enddo
close(12)
end


ctl(就是那个包涵.map 运行之后生成map的那个ctl)如下:
DSET E:\SCI2\datstation\obs_allstation_V_2012.grd
DTYPE station
STNMAP E:\SCI2\datstation\obs_allstation_V_2012.map
UNDEF -999999.0
TITLE 100m V
tdef 168 linear 00Z20Jun2012 6hr
VARS 1
V 0 99 Rain data
ENDVARS




焦头烂额!!!真相到底是什么呢?




图1

图1

图2

图2

图3

图3

UVsurface_allstation-2012.txt

97.45 KB, 下载次数: 2, 下载积分: 金钱 -5

附件1

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

新浪微博达人勋

发表于 2017-4-29 13:07:20 | 显示全部楼层
      楼主你好!谢谢你的分享。我按照你的程序做了一下,提示有错误,说只有一个时次生成,ctl文件里的时间维度有错。我的降水数据为394*160,394个月,160个测站。请问这是什么原因了?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-29 13:18:06 | 显示全部楼层

回帖奖励 +1 金钱

{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}{:eb315:}
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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