GrADS实习报告范文

试验（实习）汇报6专业大气科学年级14班级1姓名钱亚楠学号2346079试验目的1.掌握气象要素互相影响研究的基本措施。2.深入巩固FORTRAN次序构造和循环构造的程序设计措施和变量，函数，数组的使用措施，重点掌握外部文献对数据进行输入输出操作的措施。3.掌握FORTRAN中WLF函数库的使用措施。4.理解站点资料的数据构造。5.掌握站点资料数据描述文献和站点映射文献的建立措施。6.掌握站点资料生成格点资料的措施。7.学习并掌握GrADS函数的使用措施，加强对oacres(),maskout(),smth()等函数的理解。8.掌握GrADS变量的定义和使用措施。9.巩固GrADS数据处理流程，绘图要素设置，基础绘图命令的使用措施，深入掌握描述语言的应用措施。二.问题描述（1）问题分析三．试验环节四．试验关键及措施（1）计算1951-1月蒙古高压强度与中国气温同期有关Ptcor01.f90USEWFLinteger,parameter::n=60,start=1951character*10char1integeriinteger(4)station,nstationreala(n),b(n),c(160,n),d(n),e(n),coopen(1,file='D:\p.dat',form='formatted')open(2,file='D:\t1601.txt',form='formatted')open(3,file='D:\mh-t-1.dat',form='formatted')open(4,file='D:\mh-t-1.grd',form='binary')read(1,*)(a(i),i=1,n)close(1)read(2,*)((c(i,j),i=1,160),j=1,n)close(2)doi=1,160doj=1,60b(j)=c(i,j)enddocallCOR2(n,a,b,co,colev)write(3,'(f12.5)')cowrite(4)coenddoclose(3)close(4)End站点数据转换成格点数据Cor-stntogrd.f90USEWFLinteger,parameter::n=60,start=1951character*10char1integeriinteger(4)station,nstationreala(n),b(n),c(160,n),d(n),e(n),coopen(1,file='D:\p.dat',form='formatted')open(2,file='D:\t1601.txt',form='formatted')open(3,file='D:\mh-t-1.dat',form='formatted')open(4,file='D:\mh-t-1.grd',form='binary')read(1,*)(a(i),i=1,n)close(1)read(2,*)((c(i,j),i=1,160),j=1,n)close(2)doi=1,160doj=1,60b(j)=c(i,j)enddocallCOR2(n,a,b,co,colev)write(3,'(f12.5)')cowrite(4)coenddoclose(3)close(4)End生成160站的格点文献grid.f90programmainparameter(nx=71,ny=41)reallat(ny),lon(nx)reals(nx,ny)open(1,file='D:\data\grid.grd',form='binary')lat(1)=15.0lon(1)=70.0doj=1,ny-1lat(j+1)=lat(j)+1.0enddodoi=1,nx-1lon(i+1)=lon(i)+1.0enddodoi=1,nxdoj=1,nys(i,j)=1enddoenddowrite(1)sclose(1)end编写“mh-t-gr.grd”的数据描述文献mh-t-gr.ctldsetD:\data\mh-t-gr.grddtypestationstnmapD:\data\china.mapundef-999.0titleCorrelationoftandmhtdef1linearJAN19511monvars1r099correlationofmhpandtEndvars编写“grid.grd”的数据描述文献Grid.ctldsetD:\data\grid.grdundef-999.0titleSampleGRIDDataxdef71linear701ydef41linear151zdef1linear10001tdef1linearJAN19511monvars1gd099griddatapreparedforoacresfunctionEndvars绘制1951-1月蒙古高压强度与中国气温同期有关图ptcor1.gs'reinit''openD:\grid.ctl''openD:\mh-t-gr.ctl''setmap111''setlon72.5137.5''setlat17.555''sett1''setmpdsethirescnworld''setgridoff''setgradsoff''definea=oacres(gd,r.2,1.5)''definea1=maskout(a,gd-0.5)''defineaa=smth9(a1)''setgxoutshaded''setxlopts160.15''setylopts160.15''setclopts060.12''setclevs-0.418-0.333-0.2560.2560.3330.418''setccols4370734''daa''setgxoutcontour''setclabforced''setcthick8''daa''runD:\southsea.gs''setmpdsetcnriver''setmap218''drawmap''printimD:\ptcor.pngwhite';1951-1月蒙古高压强度与中国气温同期有关图试验总结学习并掌握了同期有关系数的求解措施，运用函数求1951-1月蒙古高压强度与中国气温同期有关的160个有