數(shù)據(jù)文件的轉(zhuǎn)換及數(shù)據(jù)描述文件的建立試驗二

1、南京信息工程大學(xué)實驗（實習(xí)）報告實驗（實習(xí)）名稱數(shù)據(jù)文件的轉(zhuǎn)換及數(shù)據(jù)描述文件的建立日期2016.3.18指導(dǎo)老師馬利系專業(yè)年級_班次姓名學(xué)號實習(xí)二數(shù)據(jù)文件的轉(zhuǎn)換及數(shù)據(jù)描述文件的建立1 .實習(xí)資料：現(xiàn)有“data”文件夾ASCII碼數(shù)據(jù)資料文件：月平均風(fēng)場（u200.dat、u850.datv200.datv850.dat）;月平均高度場（hgt500.dat）。資料的水平網(wǎng)格范圍是MXN個網(wǎng)格點（M=37N=17），分辨率為2.50X2.5,范圍為：自西向東經(jīng)度60-1500E,由南至北緯度0-40No時段為2002.12005.12共48個月。2 .實習(xí)要求：1）安裝GrADS運行軟件（1.

2、8版或2.0版），熟悉GrADS!行環(huán)境；2）根據(jù)所提供的數(shù)據(jù)文件，編寫出將ASCII碼數(shù)據(jù)資料文件（月平均風(fēng)場和月平均高度場）分別轉(zhuǎn)換成二進(jìn)制無格式直接存?。℅rAD繳據(jù)格式）文件（uv.grd和hgt.grd文件）的Fortran程序；3）寫出相應(yīng)的數(shù)據(jù)描述文件（uv.ctl和hgt.ctl文件），并利用GrADSS本操作命令（open*.ctl;d*）顯示該地區(qū)850hPa和200hPauv風(fēng)場及500hPa位勢高度場圖形，最后保存圖形（printim路徑*png）。3 .實習(xí)目的：通過編寫該程序，進(jìn)一步熟悉數(shù)據(jù)GrADS的數(shù)據(jù)格式，熟練編寫將ASCII碼數(shù)據(jù)轉(zhuǎn)換為二進(jìn)制格式的Fortr

3、an程序以及建立相應(yīng)的CTL文件。4 .實習(xí)步驟：4.1 安裝GrADS運行軟件（1.8版或2.0版）；4.2 利用“寫字板”程序打開*.dat數(shù)據(jù)文件，熟悉該數(shù)據(jù)資料；4.3 按要求編寫Fortran程序，將所給的ASCII碼數(shù)據(jù)資料文件轉(zhuǎn)換成二進(jìn)制無格式直接存取文件，結(jié)果保存為uv.grd和hgt.grd;4.4 通過“寫字板”或“記事本”程序編寫相應(yīng)的數(shù)據(jù)描述文件：uv.ctl和hgt.ctl;4.5 畫出2002年1月850hPa風(fēng)場圖，與圖1比較驗證數(shù)據(jù)正確性；4.6 利用GrADSS本操作命令（open*.ctl,d*）顯示2002年7月850hPa和200hPauv風(fēng)場圖形及50

4、0hPa位勢高度場圖形，最后保存圖形（printim路徑*.png）；4.7 完成實習(xí)報告1）說明所用資料2）給出所編寫的Fortran程序3）給出所編寫的數(shù)據(jù)描述文件4）給出所繪制的圖形10圖12002年1月850hpa風(fēng)場圖5 .實驗結(jié)果所用資料“data”文件夾下有ASCII碼數(shù)據(jù)資料文件:月平均風(fēng)場：u200.dat、u850.dat、v200.dat、v850.dat月平均高度場：hgt500.dat時間范圍：2002.1-2005.12共48個月水平范圍：37X17個網(wǎng)格點，分辨率為2.5X2.5經(jīng)緯度范圍：60-150E,0-40N2 .fortran程序programexam0

5、1parameter(nx=37,ny=17,nt=48)dimensionu200(nx,ny,nt),u850(nx,ny,nt),v200(nx,ny,nt),v850(nx,ny,nt)open(1,file=D:u200.dat)open(2,file=D:u850.dat)open(3,file=D:v200.dat)open(4,file=D:v850.dat)open(12,file=D:uv.grd,form=binary)doit=1,ntread(1,*)(u200(i,j,it),i=1,nx),j=1,ny)read(2,*)(u850(i,j,it),i=1,nx)

6、,j=1,ny)read(3,*)(v200(i,j,it),i=1,nx),j=1,ny)read(4,*)(v850(i,j,it),i=1,nx),j=1,ny)enddodoit=1,ntwrite(12)(u850(i,j,it),i=1,nx),j=1,ny)write(12)(u200(i,j,it),i=1,nx),j=1,ny)write(12)(v850(i,j,it),i=1,nx),j=1,ny)write(12)(v200(i,j,it),i=1,nx),j=1,ny)enddoend3 .數(shù)據(jù)描寫文件uv.ctlDSETD:datauv.ctlUNDEF-9.99E33TITLEUpperAirDataXDEF37LINEAR60.02.5YDEF17LINEAR0.02.5ZDEF2LEVELS850200TDEF48LINEARJAN20021movars2u299uwind(m/s)v299vwind(m/s)endvarshgt.ctlDESTD:datahgt.grdUNDEF-9.99E+33XDEF37linear60.02.5YDEF17linear0.02.5ZDEF1levels500