linux 下fortran 建立二进制文件的问题

xiaohundanpipi

我想把站点数据读成二进制的文件，然后用grads做图，根据站里的教程我读入txt然后生成二进制文件，
open(3,file='fallaop.grd',form='binary')
可是在运行fortran的时候linux下显示的错误是：
OPEN 语句中 FORM 指定了非法的值‘binary’

我又把form='binary'改成form='unformatted'，这次fortran程序可以运行，生成了二进制文件，可是到了生成map的时候，ctl如下
dset /home/corre/FALL/fallaop.grd
dtype station
stnmap index.map
undef -999.0
title drought index
tdef 1 linear Jul1951 1mo
vars 1
ind 0 99 drought index
endvars
就出现了这种错误：
Invalid station hdr found in station binary file
  Possible causes:  Invalid level count in hdr
                    Descriptor file mismatch
                    File not station data
                    Invalid relative time
    levs = 0  flag = 1  time = 126.65
上网查了下，我觉得是因为二进制数据格式的问题，应该用binary这个格式，可是linux不认？这是为什么，求帮助。

mofangbao

Linux的Fortran还没测试，binary不知道支不支持。把你输出的部分也贴出来吧

xiaohundanpipi

这是把txt都成二进制文件的fortran程序:

program main
real::lat(160),lon(160),ind(160)
integer::i
character*8 stid(160)

open(1,file='fall_p_ao.txt',status='old')
do i=1,160
  read(1,*) lat(i),lon(i),ind(i)
enddo
close(1)

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

open(3,file='fallaop.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do i=1,160
   write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,ind(i)
enddo
nlev=0
write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
close(3)
end

言深深

linux的没有玩过，不过写有格式文件（binary）的时候建议用形如这样的语句写
(((x(i,j,k),i=1,m),j=1,n),k=1,o)
试试看

xiaohundanpipi

言深深 发表于 2012-10-30 10:59

                               
登录/注册后可看大图

linux的没有玩过，不过写有格式文件（binary）的时候建议用形如这样的语句写
(((x(i,j,k),i=1,m),j=1,n),k ...

还是不行，我编了小程序：
program ii
implicit none
open(1,file='1.grd',form='binary')
close(1)
end
还是出现的那个错误，
(1)处 OPEN 语句中 FORM 指定了非法的值‘binary’
应该是linux下不认binary，unformatted又行不通
谁知道该怎么办？

兰溪之水

open(3,file='fallaop.grd',form='unformatted',access='direct')

银狐

open(10,file='fallaop.grd',form='unformatted',access='direct',status='replace',recl=4)在windows中recl=1，linux中为4.若是在工作站中运行还要注意是否需要高位翻转，convert='big_endian'..

xiaohundanpipi

兰溪之水 发表于 2012-10-31 00:53

                               
登录/注册后可看大图

open(3,file='fallaop.grd',form='unformatted',access='direct')

谢谢兰溪之水，按照你的思路已经作出map文件了

xiaohundanpipi

银狐 发表于 2012-10-31 17:01

                               
登录/注册后可看大图

open(10,file='fallaop.grd',form='unformatted',access='direct',status='replace',recl=4)在windows中rec ...

嗯嗯，刚接触站点资料，还没弄懂高位翻转是什么意思，不过recl=4真的和window中的不同，不久前在这里卡住过，看来linux下还是要尽快学习，尽快熟悉，谢谢啦～

letring

谁知道怎么用fortran读取GrADS的站点降水资料啊