站点文件出错

希尔博拉斯

小弟最近在做160站站点文件的降水，按照清风大大的教学贴一步步做的，可是画出的图一直有问题，找了好久找不到原因。。。只有左上角。。。求帮助。。下面贴上程序。。

希尔博拉斯

program main
implicit none
character*8::stid
real::lon,lat,rain,tim
integer::nlev,flag

tim=0.0
nlev=1
flag=1

open(11,file='e:\qingzang\chi\stid.txt',status='old')
open(12,file='e:\qingzang\chi\bu1.txt',status='old')
open(13,file='e:\qingzang\chi\sta.grd',status='replace',form='binary')
open(14,file='e:\qingzang\chi\st1.txt',status='replace')

10 read(11,*,end=100)stid,lon,lat
   read(12,*,end=100)rain
   
   write(13)stid,lat,lon,tim,nlev,flag,rain
   write(14,*)stid,lat,lon,tim,nlev,flag,rain
  
   
goto 10
100 continue
close(11)
close(12)

nlev=0
flag=0
write(13)stid,lat,lon,tim,nlev,flag
write(14,*)stid,lat,lon,tim,nlev,flag
close(13)
close(14)
end

这个是fortran。。

希尔博拉斯

dset e:/qingzang/chi/sta.grd
dtype station
stnmap e:/qingzang/chi/sta.map
undef -999.0
title rain data sample
tdef 1 linear 01jan0000 1mo
vars 1
rain 0 99 rain data
endvars


这个是ctl

希尔博拉斯

stnmap文件也能正常生成（虽然不知道有没有问题），可是图就是不在正常位置啊

一个不小心

没看懂你那个FORTRAN在做什么，能解释下么？

希尔博拉斯

一个不小心 发表于 2014-5-16 20:45
没看懂你那个FORTRAN在做什么，能解释下么？

把站点数据的txt写到dat里面去

希尔博拉斯

我突然发现国家气候中心原来是纬度在前，经度在后的。。。

river

160站是中国的吧，那你那图都不知道跑哪去了，经纬度有问题吧

希尔博拉斯

river 发表于 2014-5-16 21:21
160站是中国的吧，那你那图都不知道跑哪去了，经纬度有问题吧

后来发现的确是的。。中国气候中心上经度在纬度后面。当时没注意到。。