多时次站点数据作图问题

LNYZ

生成站点映射文件时显示错误



看到许多人提出这个问题是循环有问题，但是我用print输出的数据与我的原数据是一致的，所以想请教一下大神，不胜感激。
这是fortran文件
program main
implicit none
integer::ero,i,j,datasNb,Year
real::datas(6,50),lat(6),lon(6)
open(1,file='E:/SPI/Console1/SPI.txt',status='old')
read(1,*)lat(1),lat(2),lat(3),lat(4),lat(5),lat(6)
read(1,*),lon(1),lon(2),lon(3),lon(4),lon(5),lon(6)
print*,lat(2)
i=3
do
   read(1,*,iostat=ero)Year,datas(1,i),datas(2,i),datas(3,i),datas(4,i),datas(5,i),datas(6,i)
   if(ero/=0)exit
   i=i+1
enddo
print*,datas(2,3)
datasNb=i-1
open(2,file='E:\SPI\Console1\SPI.grd',status='replace',form='binary',access='direct',recl=6*datasNb)
write(2,rec=1) lat(1),lat(2),lat(3),lat(4),lat(5),lat(6),lon(1),lon(2),lon(3),lon(4),lon(5),lon(6),year,((datas(i,j),i=1,6),j=1,datasNb)
close(2)
pause
end program


这是我的描述文件
dset E:\SPI\Console1\SPI.grd
dtype station
stnmap E:\SPI\SPI.map
undef -999.0
title station data
xdef 3 linear 97.54 2.5
ydef 3 linear 32.31 2.5
tdef 42 linear 01JAN1976 1yr
vars 1
spi 0 99 spidata
endvars

river

有很多人分享过多站点，多时次，多层次的站点图做法，请楼主在论坛里搜索，看看相关的帖子。别人提出的循环问题，不是fortran的问题，是说的你的这种循环写出来的站点资料格式是否符合grads对站点数据格式的要求问题。所以建议你先了解一下grads要求的二进制的站点格式资料是什么样的

哟哟呀

应该是生成二进制文件的时候出错啦

me007

他山之石，可以攻玉，楼主论坛找找，应该能够解决问题。我也是菜鸟，没办法帮助您。

yuu777

应该是不符合grads的站点数据格式，需要台站号

老A

站点数据格式写的不对，提示里面说的很清楚

老A

每个站点数据要先写头，再写站号，经纬度层次标记等等后再写值，值要写4字节有符号浮点数，再仔细看看站点数据的格式要求