Grads实习报告文档格式.docx
《Grads实习报告文档格式.docx》由会员分享，可在线阅读，更多相关《Grads实习报告文档格式.docx（16页珍藏版）》请在冰点文库上搜索。
\llp\data2.grd'
form='
binary'
dok=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）
1000format（2i7）
2000FORMAT（37F6.2）
3000format（37f8.1）
4000format（37f7.2）
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；
得到data2.grd，与原来ctl文件相配，画图验证即可。
实习三
求出高度场850hPa、82-85年1月、沿100-120E的平均值，并画出纬度-时间剖面。
编写gs文件如下：
shiyan3_1.gs
'
reinit'
'
opene:
\shiyan3\data.ctl'
setgxoutfwrite'
setfwritee:
\shiyan8\h.grd'
setlon110'
setlat040'
setlev850'
i=1
while（i<
=48）
sett'
i'
defineha=ave（h,lon=100,lon=120）'
dha'
i=i+12
endwhile
disablefwrite'
;
运行e:
\shiyan3\shiyan3.gs
得到一个grd文件，即h.grd
为得到的grd文件配ctl文件：
shiyan3.ctl
dsete:
\shiyan3\h.grd
undef-9.99E+33
titleNCEP/NCARREANALYSISPROJECT
xdef1linear100.02.5
ydef17linear0.02.5
zdef1levels850
tdef4linearJAN19821yr
vars1
ha099height
endvars
再编写gs文件如下：
shiyan3_2.gs
\shiyan3\shiyan3.ctl'
enableprinte:
\shiyan3\shiyan3.gmf'
setlon100'
sett14'
setgxoutcontour'
print'
disableprint'
运行shiyan3_2.gs
得到图，shiyan3.gmf
实习四
画出地形图，3000到5000之间划阴影，其余画等值线。
首先把orog.dat的资料转成二进制在fortran中输入以下代码得到orog.grd
programex_04
integerii,jj,i,j
parameter（ii=37,jj=17）
realvar（ii,jj）
open（1,file='
c:
\OROG.dat'
read（1,3000）（（var（i,j）,i=1,ii）,j=1,jj）
open（2,file='
\OROG.grd'
write
（2）（（var（i,j）,i=1,ii）,j=1,jj）
3000format（37f8.2）
end配orog.ctl：
\shiyan4\orog.grd
xdef37linear60.0002.500
ydef17linear0.0002.500
tdef48linearJAN19821mo
orog00（0或99均可）
_
之后运行以下命令：
e:
\shiyan4\shiyan4.gs
\shiyan4\orog.ctl'
\shiyan4\shiyan4.gmf'
setgxoutshaded'
setcmin3000'
setcmax5000'
dorog'
实习五
画出35N的地形剖面，并画出3000米线。
本题为两题：
利用上题得到的orog.ctl及orog.grd
画出3000米线:
在grads中键入以下命令：
opene:
\orog.ctl
setclevs3000
dorog
得到shiyan5_1.gmf
35N的地形剖面:
setlat35
得到shiyan5_2.gmf
gs文件：
\shiyan5\shiyan5_1.gs
\orog.ctl'
\shiyan\shiyan5_1.gmf'
setclevs3000'
e:
\shiyan5\shiyan5_2.gs
setlat35'
实习六
从原始资料中写出20-30N，80-10E范围内资料，并配以数据描述文件，简单绘图验证正确与否。
关于实验6，在grads19中键入以下命令：
reinit
openc:
\data.ctl
setfwrited:
\shiyan6\shiyan6.grd
setgxoutfwrite
while（i<
setti
setlat2030
setlon80100
j=1
while（j<
=2）
setzj
du
dv
j=j+1
endwhile
setz1
dh
dtsfc
i=i+1
disablefwrite
通过以上过程可以得到一个从原始资料中写出的20-30N，80-10E范围内资料，然后再配以数据描述文件shiyan6.ctl，
然后在grads19中绘图验证是否正确。
dsetd:
xdef9linear80.0002.500
ydef5linear20.0002.500
zdef2levels850200
vars4
U299uwind（m/s）
V299vwind（m/s）
H199H500
TSFC199TSFCDATA
gs文件
\data.ctl'
\shiyan6\shiyan.grd'
setlat2030'
setlon80100'
setz'
j'
du'
dv'
setz1'
dh'
dtsfc'
shiyan6.gmf为使用shiyan6.grd及shiyan6.ctl所画图
shiyan6_1.gmf为使用data.ctl及data.grd画图之后与原图叠加而成。
实习七
画出u850的沿100E的时间-纬度剖面。
\shiyan7\data.ctl'
\shiyan7.gmf'
du'
实习八
循环显示所有层次、所有时次的v场。
\shiyan8\data.ctl'
\shiyan8\shiyan8.gmf'
i=1
dv'
c'
i=i+1
求出500hpa高度场的距平场，并画出82年12月和83年7月的距平场，要求<
0的值标出阴影。
enableprintc:
\shiyan9.gmf'
setparea11018'
setgradsoff'
setgridoff'
setxlopts150.16'
setylopts150.16'
setmap1110'
sett112'
have=ave（h,t=0,t=48,12）'
modifyhaveseasonal'
sett148'
hano=h-have'
settimedec1982'
setcmax0'
dhano'
setclevs0'
setcthick10'
drawtitle500hpaanomalousheightfieldsforDec1982'
drawxlabLON'
drawylabLAT'
print'
pulldummy'
settimejul1983'
drawtitle500hpaheightanomalousheightfieldsforJul1983'
循环输出五天高度场图，并将所得图形存于指定文件中。
\shiya8.gmf'
=5）
dh'
显示第一时次高度场图，并在低压中心标注字符D：
setlon60150'
sett1'
qw2xy7020'
x1=subwrd（result,3）'
y1=subwrd（result,6）'
drawstring'
x1'
y1'
D'
printimc:
\shi.pngwhite'
升级会员 