绘制站点分布图（比较菜鸟，但是能用）

这只个名字

第一次发帖子，今天第一次自己用Fortran，grads画出了江淮流域站点分布图，给大家分享一下程序，比较菜鸟，需要的大家可以参考一下
Fortran：
program station
      parameter (mt=428)    ! st means nums of station
      character stid(mt)*8
      integer nflag,nlev
      real lat(mt),lon(mt),tim,t2(mt)
open(30,file='C:\biancheng\lybjw.txt')
open(50,file='C:\biancheng\sta.txt')

open(40,file='C:\biancheng\lybjw.dat',form='binary')
      do i=1,mt
      read(30,*)stid(i),lat(i),lon(i)  
      write(*,*)stid(i),lat(i),lon(i)
      end do   
      close(30)
      read(50,*) (t2(i),i=1,mt)
      close(50)
      tim=0.0
      nflag=1
      nlev=1
      do i=1,mt
      write(40)stid(i),lat(i),lon(i),tim,nlev,nflag,t2(i)
      end do
      nlev=0
      write(40)stid(mt-1),lat(mt-1),lon(mt-1),tim,nlev,nflag
      close(40)
      end

------------------------------------------------------------------------------
CTL：
dset C:\biancheng\lybjw.dat
dtype station
stnmap C:\biancheng\grads\st428.map
undef 999.9
title temp duration
tdef 1 linear 14Jul2002 1hr
vars 1
t2 0 99 station number
endvars

-----------------------------------------
gs:
'reinit'
'open C:\biancheng\grads\station.ctl'
'enable print C:\biancheng\grads\sta428.gmf'
'set lon 100.0 130.0'
'set lat 25.0 38.0'
'set grads off'
'set grid off'
'set mpdset cn cnriver'
'set gxout stnmark'
'set rbcols 3 3 3'
'set cmark 3'
'set digsiz 0.1'
'set stid off '
'd t2'
'print'
'disable print'
;

CRYJKW

感谢楼主分享，我研究研究

梦呼呼

我现在这个水准都达不到，向你学习

snail_t

入门级别操作，加油

aptx-4869

感谢

fy861123

入门级操作，加油