- 积分
- 43
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-4-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
先贴出FORTRAN程序:
real::m(12*176),lat(12*176),lon(12*176)
character*8 stid(12*176)
integer nlev,flag
tim=0.0
nlev=1
flag=1
open(1,file='F:\m.dat')
read(1,*)m
open(2,file='F:\lat.dat')
read(2,*)lat
open(3,file='F:\lon.dat')
read(3,*)lon
close(1)
close(2)
close(3)
do i=1,176*12
if(lat(i)==0.0)then
m(i)=-999.0
endif
enddo
open(4,file='F:\station.txt')
do i=1,176*12
read(4,*)stid(i)
enddo
close(4)
open(5,file='F:\station.grd',form='binary')
do i=1,12*176
write(5,*)stid(i),lat(i),lon(i),tim,nlev,flag,m(i)
enddo
nlev=0
write(5,*)stid(i-1),lat(i-1),lon(i-1),tim,nlev,flag,m(i-1)
close(5)
end
我的CTL:
DSET F:/station.grd
DTYPE station
STNMAP F:/grid.map
UNDEF -999.0
TITLE DUCT Data
TDEF 1 linear 01jan2014 1mo
VARS 1
m 0 99 DUCT MAGTINUTY
ENDVARS
请大家帮帮忙,如果可以进一步帮我,加我QQ吧,我想跟大家多学点,1284569459
|
|