我用micaps的第二类数据转grads的站点数据怎么转出来没有数据。

笑哈乐呵

这是我的fortran程序
program main
integer,parameter::ns=435
character*8 stid(ns)
integer i,tim,nlev,flag
real lon(ns),lat(ns),hgt,fx(ns),fs(ns),u(ns),v(ns),temp,class,un,hb
open(1,file='d:\tu\70\10071508.000',status='old')
read(1,*)   
read(1,*)
do i=1,ns
     read(1,'(a6,2f8.2,4f5.1,3f5.1)') stid(i),lon(i),lat(i),hb,class,hgt,temp,fs(i),fx(i),un   
enddo
tim=0.0
nlev=1
flag=1
open(2,file='d:\tu\70\1508.dat',status='new',form='binary')
    do i=1,ns
        write(2) stid(i),lon(i),lat(i),tim,nlev,flag
        u(i)=sin((fx(i)-180)*3.14159/180)*abs(fs(i))
     v(i)=cos((fx(i)-180)*3.14159/180)*abs(fs(i))  
           write(2)u(i),v(i)
    enddo
close(1)
nlev=0
flag=0
Write(2) stid(1),lon(1),lat(1),tim,nlev,flag
close(2)
end
我的数据文件是：

diamond 2 10年07月15日08时700百帕高空观测

10 07 15 08  700 435

  1152   14.37   67.27   13    1  302   -3   26  345    7

  1415    5.63   58.88    9    1  302    3    7  170   21

  2185   22.13   65.55   16    1  303    3   25  255   18

  2365   17.45   62.53    6    1  308    3    9  260   13

  2836   26.65   67.37  179    1  301    1   14  240   17

  2963   23.50   60.82  103    1  312    6   14  215   11

  3005   -1.17   60.13   84    1  297   -1  2.2  105    7

  3354   -1.24   53.00  117    1  290   -3    9  195   11

  3808   -5.31   50.22   88    1  292   -2   .0  260   21

ctl：

DSET   d:\tu\70\1508.dat

DTYPE  station  

STNMAP d:\tu\70\700.map

UNDEF  9999.9

TITLE  700hpa steam

TDEF   1 linear 08Z15jul2010 12hr

VARS 2

u    0  99  wind

v    0  99  wind

ENDVARS

'reinit'

'open d:/tu/70/700.ctl'

'set lon 70 140'

'set lat 10 50'

'set gxout  barb'

'd u;v'

'printim  d:/tu/70/6.png white'

我只想u，v的数据，结果什么也没有。

北五味

帮顶~~~~~~

笑哈乐呵

北五味 发表于 2013-4-2 08:48

                               
登录/注册后可看大图

帮顶~~~~~~

谢谢啊！

river

改成这样试试
do i=1,ns
     u(i)=sin((fx(i)-180)*3.14159/180)*abs(fs(i))
     v(i)=cos((fx(i)-180)*3.14159/180)*abs(fs(i))  
           write(2) stid(i),lon(i),lat(i),tim,nlev,flag,u(i),v(i)
    enddo
close(1)
nlev=0
Write(2) stid(1),lon(1),lat(1),tim,nlev,flag
close(2)
end