文档介绍：实****二
编写Fortran程序,将上述第一类资料1)-3)所列资料转成二进制格式;4)资料单独转成二进制形式;并配以相应的数据描述文件,即ctl文件。
注意:验证转换后的资料是否正确,可以再配以ctl文件以后,对资料进行作图,应该与提供的第二类资料图形一致。
在fortran编译器中输入以下代码:
program main
integer i,j,k
parameter(ii=37,jj=17,nz=2,nt=48)
dimension a(ii,jj),b(ii,jj),c(ii,jj),d(ii,jj),e(ii,jj),f(ii,jj)
open(1,file='d:\llp\dat\')
open(2,file='d:\llp\dat\')
open(3,file='d:\llp\dat\')
open(4,file='d:\llp\dat\')
open(5,file='d:\llp\dat\')
open(6,file='d:\llp\dat\')
open(12,file='d:\llp\',form='binary')
do k=1 ,48
read(1,1000)
read(1,2000) ((a(i,j),I=1,ii),j=1,jj)
read(3,1000)
read(3,2000) ((c(i,j),I=1,ii),j=1,jj)
read(2,1000)
read(2,2000) ((b(i,j),I=1,ii),j=1,jj)
read(4,1000)
read(4,2000) ((d(i,j),I=1,ii),j=1,jj)
read(5,1000)
read(5,3000) ((e(i,j),I=1,ii),j=1,jj)
read(6,1000)
read(6,4000) ((f(i,j),I=1,ii),j=1,jj)
1000 format(2i7)
2000 FORMAT()
3000 format()
4000 format()
write(12) ((a(i,j),i=1,ii),j=1,jj)
write(12) ((b(i,j),i=1,ii),j=1,jj)
write(12) ((c(i,j),i=1,ii),j=1,jj)
write(12) ((d(i,j),i=1,ii),j=1,jj)
write(12) ((e(i,j),i=1,ii),j=1,jj)
write(12) ((f(i,j),i=1,ii),j=1,jj)
enddo
end;,与原来ctl文件相配,画图验证即可。
实****三
求出高度场 850hPa、82-85年1月、沿100-120E的平均值,并画出纬度-时间剖面。
编写gs文件如下:
'reinit'
'open e:\shiyan3\'
'set gxout fwrite'
'set fwrite e:\shiyan8\'
'set lon 110'
'set lat 0 40'
'set lev 850'
i=1
while(i<=48)
'set t 'i''
'define ha=ave(h,lon=100,lon=120)'
'd ha'
i=i+12
endwhile
'disable fwrite'
;
运行e:\shiyan3\
得到一个grd文件,
为得到的grd文件配ctl文件:
dset e:\shiyan3\
undef -+33
title NCEP/NCAR REANALYSIS PROJECT
xdef 1 linear
ydef 17 linear
zdef 1 levels 850
tdef 4 linear JAN1982 1yr
vars 1
ha 0 99 height
endvars
再编写gs文件如下:
'reinit'
'open e:\shiyan3\'
'enable print e:\shiyan3\'
'set lon 100'
'set lat 0 40'
'set lev 850'
'set t 1 4'
'