【求助】关于Fortran读取站点资料为grads识别格式

cjj

我将62个站点38年的月降水量整理到txt文件中，但是用Fortran读取为二进制格式后，用grads生成.map文件出现问题，一直找不到解决方案，恳请论坛大佬帮我解答一下！
ctl文件如下：
dset D:\data\rain\1979-2016.dat
dtype station
stnmap D:\data\rain\1979-2016.map
undef -999.0
title monthly rainfall
tdef 456 linear JAN1979 1mo
vars 1
rain 0 99 rainfall data
endvars


站点数据格式：

                               
登录/注册后可看大图

Fortran读取为.dat二进制格式：
Program sta2grd
    Implicit none
    Character(len=8) ::  stid
    real lat,lon,rain
    integer year,month
    integer nlev,flag,yrold,mnold

    flag=1
    yrold=1979
    mnold=1
    Open(11,file='D:\data\rain\1979-2016.txt',status='old')
    Open(12,file='D:\data\rain\1979-2016.dat',status='replace',form='binary')
10  Read(11,*) year,month,stid,lat,lon,rain
    if(yrold.NE.year.OR.mnold.NE.month) then
       nlev=0
       write(12) stid,lat,lon,nlev,flag
    end if
    nlev=1
    write(12) stid,lat,lon,nlev,flag,rain
    yrold=year
    mnold=month
    goto 10   
    close(11)
    close(12)
    End program


Fortran运行过程也出现了问题：

                               
登录/注册后可看大图

安达充

写成站点数据的时候头信息写错了，漏了一个相对时刻。
write(12) stid,lat,lon,nlev,flag应该改为write(12) stid,lat,lon,t,nlev,flag

（不会贴图，楼主自己找一下Grads站点数据的文件格式吧，GrADS使用手册里就有）

cjj

安达充 发表于 2020-3-24 11:37
写成站点数据的时候头信息写错了，漏了一个相对时刻。
write(12) stid,lat,lon,nlev,flag应该改为write(12 ...

非常感谢您的回复！我按您的意见修改了一下Fortran代码，也把goto换成了do循环就不会出现运行时错误了。但是在用grads生成.map文件时出现了错误。修改后代码：Program sta2grd
    Implicit none
    integer,parameter :: num=28262
    Character(len=8) ::  stid
    real lat,lon,rain,t
    integer year,month
    integer nlev,flag,yrold,mnold,i

    flag=1
    yrold=1979
    mnold=1
    Open(11,file='D:\data\rain\1979-2016.txt',status='old')
    Open(12,file='D:\data\rain\1979-2016.dat',status='replace',form='binary')
    do i=1,num
    Read(11,*) year,month,stid,lat,lon,rain
    if(yrold.NE.year.OR.mnold.NE.month) then
       nlev=0
       write(12) stid,lat,lon,nlev,flag
    end if
    t=0.0
    nlev=1
    flag=1
    write(12) stid,lat,lon,t,nlev,flag,rain
    yrold=year
    mnold=month  
    end do
    close(11)
    close(12)
    End program

安达充

cjj 发表于 2020-3-24 12:08
非常感谢您的回复！我按您的意见修改了一下Fortran代码，也把goto换成了do循环就不会出现运行时错误了。 ...

你有两个地方都有输出，都是要有t的，你看GrADS里的报错信息就知道你写入的数据和GrADS要求的是不一致的。