- 积分
- 2229
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-28
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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
焦头烂额!!!真相到底是什么呢?
|
|