grads初学者，急急急，把txt数据用fortran转化成grd文件出错

weiqiyun0049

fortran程序参照气象家园教程的：

program transform

    implicit none

    Character*8 stid
    real lat
    real lon
    real tim
    integer nlev
    integer flag
    real var

    tim=0.0
    nlev=1
    flag=0

    open(1,file='transform.txt',status='old')
    open(2,file='result.grd',status='replace',form='binary')

    read(1,*)
    read(1,*)

10  read(1,*,end=100)stid,lon,lat,var
    write(2)stid,lat,lon,tim,nlev,flag,var

    goto 10
100 continue

    close(1)
    nlev=0
    write(2)stid,lat,lon,tim,nlev,flag

    close(2)

    end program transform

weiqiyun0049

四叶草 发表于 2015-8-31 17:37
那楼主继续把报错信息贴出来呗，难道还是一样的信息？

没事了已经解决了谢谢

weiqiyun0049

数据格式是这样的
stid  lon  lat  var
    0    0    0  154
    1    1    0  154
    2    2    0  154
    3    3    0  154
    4    4    0  154
    5    5    0  153
    6    6    0  153
    7    7    0  153
    8    8    0  153
    9    9    0  154
   10   10    0  154
   11   11    0  154
   12   12    0  154

四叶草

    open(1,file='transform.txt',status='old')
    open(1,file='result.grd',status='replace',form='binary')
楼主，好好看看FORTRAN的基本知识。

weiqiyun0049

四叶草 发表于 2015-8-31 17:24
open(1,file='transform.txt',status='old')
    open(1,file='result.grd',status='replace',form='b ...

擦。。。笔误

weiqiyun0049

然而改了也并没有卵用啊，还是一样的报错信息，只不过unit=1变成unit=2而已

weiqiyun0049

已解决，原来是无格式的grd，write的时候不能用（2，*）

四叶草

那楼主继续把报错信息贴出来呗，难道还是一样的信息？

just_for_you！

{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}

str@nger