站点显示位置与实际站点位置不对应，求助各位大神！！！

大力

画站点图的问题，站点显示位置与实际位置不对应，求助各位大神！！！按照《GrADS站点文件作图详细解决方案》的帖子，进行数据处理和画图。数据处理没有问题，生成映射文件也正常，就是在站点数据显示的时候出现问题。我是把站点的海拔高度信息作为变量进行显示，其显示位置与实际站点位置不对应。如图


我的站点数据是国内的，站点信息如：23175       87.200    68.700    43.000，分别为站号、经度、纬度、海拔高度

Fortran程序如下：
program station2grd
implicit none
integer,parameter :: sumstid=10387
character*8 stid(sumstid)
real lon(sumstid),lat(sumstid),hgt(sumstid)
real tim
integer nlev,flag,i
!hgt,height above sea level
tim=0.0
nlev=1
flag=1
open(1,file='stationlist.txt',status='old')
do i=1,sumstid
  read(1,*)stid(i),lon(i),lat(i),hgt(i)
enddo
close(1)
open(2,file='sta.grd',status='replace',form='binary')
!form='unformatted'  form='binary'
do i=1,sumstid
    write(2)stid(i),lon(i),lat(i),tim,nlev,flag,hgt(i)
enddo
nlev=0
write(2)stid(sumstid-1),lon(sumstid-1),lat(sumstid-1),tim,nlev,flag
close(2)
end


ctl文件为：
dset E:/f90/station/sta.grd
dtype station
index E:/f90/station/sta.map
undef 9999.00
title station data
tdef 1 linear 00Z07jul2017 6hr
vars 1
hgt 0 99 height above sea level
ENDVARS

river

一看就是写反了，你的经纬度写反了好吧。grads站点作图的时候纬度要写在前面

大力

river 发表于 2017-10-10 14:54
一看就是写反了，你的经纬度写反了好吧。grads站点作图的时候纬度要下载前面

感谢感谢！大神又来救我了！
一下子对了，感觉真棒！

river

大力 发表于 2017-10-10 16:07
感谢感谢！大神又来救我了！
一下子对了，感觉真棒！

不是我来救你。本来你自己就能弄好，你么仔细看别人的教程，所以才浪费了好多时间······所以一定要认真，注意细节。加油！

大力

river 发表于 2017-10-10 17:01
不是我来救你。本来你自己就能弄好，你么仔细看别人的教程，所以才浪费了好多时间······所以一定要 ...

是 ！要细心，注意细节！
之前看你的帖子关于斜剖面的绘制，你也给讲解了好多，很受教！当时是另外一个账号登录的

river

大力 发表于 2017-10-10 18:06
是 ！要细心，注意细节！
之前看你的帖子关于斜剖面的绘制，你也给讲解了好多，很受教！当时是另外一个 ...

没事儿换什么账号啊，一直用一个才有很多福利的

大力

river 发表于 2017-10-10 18:49
没事儿换什么账号啊，一直用一个才有很多福利的

噢噢，后面就不换了，一直用这个向大神沟通！

river

大力 发表于 2017-10-10 19:01
噢噢，后面就不换了，一直用这个向大神沟通！

原来的号码密码忘记了？可以加管理员QQ找回的