grads对站点数据创建地图映射文件

樱花阿喵

这是我写入站点文件的程序
open(2,file='e:\practise\6\b.grd',form='binary')
        do i=1,160
    tim=0.0
    nlev=1
    nflag=1
    write(2) sta(i),jd(i),wd(i),tim,nlev,nflag,b(i)
    end do
    nlev=0
    write(2) sta(160),wd(160),jd(160),tim,nlev,nflag
    close(2)
        end
这是我的ctl文件
dset e:\practise\6\b.grd
dtype station
stnmap e:\practise\6\b.map
title 线性倾向性
undef -99999.9
tdef 1 linear jan1982 1yr
vars 1
a 0 99
endvars
图片是grads的错误，有没有大神指点啊

爱到花开

两点可能：1.第一次write的经纬度写反了，grads判断没有纬度100多的站点，所以说你这不是站点数据；2、站点sta(i)不是字符串格式，这个不确定

樱花阿喵

sta(i)确实不是字符串，因为老师给的资料站点就是以整数形式读出来的，我同学把站点以整数形式写进去也可以出图，我把经纬度调整过来还是不行

river

GrADS对站点数据有明确的格式要求，站号必须是字符型*8 ，纬度在前，经度在后。你这两个要求都没达到，所以报错很正常。
你的程序变量申明那一部分呢？

樱花阿喵

你好，这是经过我修改后的整个程序，已经可以生成地图映射文件了，但d  b之后只有地图底图，没有数据,还是很苦恼，还想请您再帮帮我可以吗  
program main
        parameter(m=160,n=25)
      integer t
      character*8 sta(m)
      integer nlev,nflag,tim,st(m)
    real jd(m),wd(m),sumx(160),f0(m,n)
        real sumtt,avex(160),avet,b(160),f(160)
        sumx=0
        sumt=0
    sumtt=0
        sumxx=0
        f=0
!reading data
      open(1,file='e:\practise\6\160zhan-rainfall-summer.txt')
        do I=1,m
        read(1,*) sta(i),jd(i),wd(i),(f0(i,j),j=1,n)
      enddo
      close(1)

    do i=1,160
        do j=1,25
        sumx(i)=sumx(i)+f0(i,j)
    end do
    avex(i)=sumx(i)/25
    end do
    !print*, avex(160)

        do j=1,25
        t=1981+j
        sumtt=sumtt+t*t
        end do

        do i=1,160
        do j=1,25
        t=1981+j
        f(i)=f(i)+f0(i,j)*t
    end do
      b(i)=(f(i)-25*avex(i)*1994)/(sumtt-25*1994*1944)
    end do
   
open(20,file='e:\practise\6\b.grd',form='binary')
    tim=0.0
    nlev=1
    nflag=1
    do 100 i=1,160
        write(20) sta(i),wd(i),jd(i),tim,nlev,nflag,b(i)
100 continue      
    nlev=0
    write(20) sta(i-1),wd(i-1),jd(i-1),tim,nlev,nflag
    close(20)
        end

樱花阿喵

b表示的是25年的降水线性倾向值

樱花阿喵

问题已经解决了，原来是经纬度读反了，但还是谢谢大神对新人的关注