[vip专享]Grads实习报告1

实习一

1.所用资料

“data”文件夹下有ASCII码数据资料文件：月平均风场：u200.dat、u850.dat、 v200.dat、v850.dat月平均高度场：hgt500.dat

时间范围：2002.1-2005.12 共48个月

水平范围：37×17个网格点，分辨率为2.5°×2.5°经纬度范围：60-150°E， 0-40°N

2.fortran程序

program exam01

parameter(nx=37,ny=17,nt=48)

dimension u200(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')do it=1,nt

read(1,*) ((u200(i,j,it),i=1,nx),j=1,ny) read(2,*) ((u850(i,j,it),i=1,nx),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)end dodo it=1,nt

write(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)end doend

3.数据描写文件①uv.ctl

DSET D:\\data\%uv.ctlUNDEF -9.99E33

TITLE Upper Air DataXDEF 37 LINEAR 60.0 2.5YDEF 17 LINEAR 0.0 2.5

ZDEF 2 LEVELS 850 200

TDEF 48 LINEAR JAN2002 1movars 2

u 2 99 u wind(m/s)v 2 99 v wind(m/s)endvars

②hgt.ctl

DEST D:\\data\\hgt.grdUNDEF -9.99E+33

XDEF 37 linear 60.0 2.5YDEF 17 linear 0.0 2.5ZDEF 1 levels 500

TDEF 48 linear JAN2002 1movars 1

hgt 1 99 q1 endvars

4.所绘制图形

①2002年一月850hpa风场图

②2002年7月850hpa风场图

③200hpa uv风场图

实习二

1.所用资料

“data”文件夹下有GrADS数据格式的二进制资料文件model.le.dat以及对应的数据描述文件model.le.ctl

这是一组模式输出的全球5天数值预报结果，包括了多个要素、多层，按经纬度网格存放的数据，具体见数据描述文件model.le.ctl

2.①一维 line.gs

'open D:\\data\\model.le.ctl''set lon 120''set lat 0 90''set lev 200''set t 1'

'set gxout line ''set ccolor 2''set cmark 3''set cstyle 1''set cthick 6''set grid on 3 3''d z'

'enable print D:\\data.line.gmf white';

②二维 contour.gs

'open D:\\data\\model.le.ctl''set lat -90 90''set lon 0 360''set lev 500''set t 1'

'set gxout contour''set cint 8''set cterp on''set csmooth on'

'set rgb 20 160 240 40''set ccolor 20''set clab %.if''set clskip 4'

'set clopts 3 0.1 0.2''d z/9.8'

'set clevs 588''set ccols 2''set cthick 7''d z/9.8'

'enable print D:\\data\\contour.gmf white';