GrADS实习报告要点.pdf

GrADS 实习报告 专业：大气科学 班级：xxxx 姓名：xxxx 注：本人上传实习报告，并不想为大家抄袭提供方便，只是想给大家一 个参考和验证， 我认为每一个用心之人都会做得比我更好，本文档只起抛砖之效。 实习一 利用已给出的二进制资料，即第二类资料，编写简单的 gs 文件， 熟悉所学Grads 基本命令，绘制简单的图形,以掌握grads 绘图基本步 骤。 1）编写简单的 gs 文件，分别绘制一个简单的一维图形（曲线）、 二维图形（等值线），并生成图形文件（三种格式），学习用 gv32.exe 查看并转换gmf 图形文件。 2）熟悉三类坐标的用法，掌握查询命令 Q，特别掌握其坐标转换 命令的用法，并在grads 命令提示符下运行 “Q w2xy 120 20”,看看 其转换结果。 另外，熟悉“Q time/q file”等用法 Gs 文件： **1、一个简单的一维图形（曲线）(u 风场沿100E 的变化曲线） reinit open d:\data.ctl enable print d:\exc1.2.gmf set lon 100 d u 2 print disable print ; 得到的图形： *2、简单的二维图形（v 风场）,并生成图形文件（三种格式） reinit open d:\data.ctl *生成gmf 文件 enable print d:\exc1.1.gmf disable print d v print disable print *生成gif 文件 enable print d:\exc1.1.gif disable print d v 3 print disable print *生成png 文件 enable print d:\exc1.1.png disable print d v print disable print 出图形示例 3、运行 “Q w2xy 120 20”结果： X=27.5 Y=13.15 4 实习二 编写Fortran 程序，将上述第一类资料1）-3）所列资料转成二进 制格式；4）资料单独转成二进制形式；并配以相应的数据描述文件， 即ctl 文件。 Fortran 程序 ******注意修改路径，下面只是一种编写方式 parameter(ii=37,jj=17) real var(ii,jj) open(1,file=u850.dat) open(2,file=v850.dat) open(3,file=u200.dat) open(4,file=v200.dat) open(5,file=h500.dat) open(6,file=tsfc.dat) open(9,file=data.grd,form=unformatted,access=direct, # recl=ii*jj*4) irec=1 DO 200 iy=1,4 5 do 100 m=1,12 ccc write u850 read(1,1000) read(1,2000) ((var(i,j),I=1,ii),j=1,jj) write(9,rec=irec) ((var(i,j),i=1,ii),j=1,jj)