160站降水生成map文件时发生这样的错误

好猫

原文件是在气候中心上下载的txt形式，月平均的降水，总共63年。
做完fortran程序转为dat格式：
program preci
parameter(sta=160,tt=63)
integer d(sta,tt),a(tt,sta),nlev,flag
real lon(sta),lat(sta),tim
character*8 stid(sta)
open(11,file='F:\graduate\work\data\station.txt')
open(22,file='F:\graduate\work\data\1.txt')
open(33,file='F:\graduate\sta.dat',status='replace',form='binary')

do i=1,sta
        read(11,*) stid(i),lon(i),lat(i)
enddo

do i=1,160
stid(i)=char(i)
end do

read(22,*) ((d(i,j),i=1,sta),j=1,tt)

a=transpose(d) !对原降水序列进行转置，使行为时次（外循环），列为站名（内循环，空间）

do i=1,tt
tim=0.0
nlev=1
flag=1
        do j=1,sta
                write(33)stid(j),lon(j),lat(j),tim,nlev,flag,a(i,j)
        enddo
nlev=0
write(33)stid(sta),lon(sta),lat(sta),tim,nlev,flag   
!print*, stid(sta),lon(sta),lat(sta),tim,nlev,flag
enddo
close(33)
close(22)
close(11)
end


然后是ctl：
DSET   e:/160/sta.dat
DTYPE  station   
STNMAP e:/160/data_sample.map
UNDEF  -999.0  
TITLE  Station Data Sample  
TDEF   756 linear 01jan1951 1mo
VARS 1  
rain    0  99  rainfall data of 160station
endvars


最后运行ctl时发生的错误是，查了一下，此处的EOF意为：end of file

还请各位多多帮忙了！我是新手，按着清风写的套路做的，翻来覆去也没看出问题！

river

lat在前，lon在后。改一下试试吧

好猫

river 发表于 2013-11-28 22:21

                               
登录/注册后可看大图

lat在前，lon在后。改一下试试吧

我确实把经纬度搞错了，不过改完还是一样的问题。
在fortran中试着在屏幕上输出也没有问题，请问您觉得可能是什么原因？

linikiku

把 这句 write(33)stid(sta),lon(sta),lat(sta),tim,nlev,flag
改成 write(33) stid(j-1),lat(j-1),lon(j-1),tim,nlev,flag

试一试。

river

好猫 发表于 2013-11-28 22:31

                               
登录/注册后可看大图

我确实把经纬度搞错了，不过改完还是一样的问题。
在fortran中试着在屏幕上输出也没有问题，请问您觉得可 ...

tim=0.0
nlev=1
flag=1
这几句应该放在循环里面的吧，就是放在时间和站点循环的最里面

好猫

linikiku 发表于 2013-11-29 00:34

                               
登录/注册后可看大图

把 这句 write(33)stid(sta),lon(sta),lat(sta),tim,nlev,flag
改成 write(33) stid(j-1),lat(j-1),lon(j ...

找到问题了！原来是我只处理的63年每年一月的数据，所以ctl文件的时次应该是63而不是756！
还是谢谢了！

好猫

river 发表于 2013-11-29 08:22

                               
登录/注册后可看大图

tim=0.0
nlev=1
flag=1

找到问题了！原来是我只处理的63年每年一月的数据，所以ctl文件的时次应该是63而不是756！
还是谢谢了！以后真的要细心！

菜菜黑

为什么要进行转制呢？

apple973

请问map文件如何生成