GrADS画图入门.ppt

GrADS 画图入门 * 高拴柱 gaosz@cma.gov.cn 国家气象中心 2008年12月7日 Grid Analysis and Display System Outline 1、grads功能概述(不讲) 2、Grads数据准备： （1）ncep资料的直接读取 （2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制 3、图形显示类型及绘图设置： （1）一维：直方、折线、散点、误差图.... （2）二维：等值线、阴影图、风矢量、流线图、站点天气图 （3）图形设置：坐标、线型、页面、地图投影、地表状况、图形信息开关 4、grads函数的使用： （1）数学函数 （2）统计函数 （3）网格处理函数 5、grads批处理脚本： 6、练习： real ps(72,46,5), z(72,46,7,5),u(72,46,7,5),v(72,46,7,5),q(72,46,5,5) .......read(***)ps,z,u,v,q...... open (8,file=’example.dat’,form=’unformatted’,access=’direct’,recl=72*46*4) irec = 0 do it=1,5 irec=irec+1 write(8,rec=irec) ((ps(j,i,it),j=1,72),i=1,46) do k=1,7 irec=irec+1 write(8,rec=irec+1) ((z(j,i,k,it),j=1,72),i=1,46) enddo ………….u,v,q……………………………………. enddo 2、Grads数据准备： （1）ncep资料的直接读取 （2）mm5_to_Grads,wrf_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 open (9,file=’example.ctl’) write(9, ’(a)’) ‘dset ^example.dat’ write(9, ’(a)’) ‘title 5 days of model output’ write(9, ’(a)’) ‘undef -2.56e33 write(9, ’(a)’) ‘xdef 72 linear 100.0 1.0 write(9, ’(a)’) ‘ydef 46 linear 10.0 1.0 write(9, ’(a)’) ‘zdef 7 levels 1000,850,700,500,300,200,100 write(9, ’(a)’) ‘tdef 5 linear 02jan2008 1dy’ write(9, ’(a)’) ‘vars 5 write(9, ’(a)’) ‘ps 0 99 Surface Pressure’ write(9, ’(a)’) ‘z 7 99 Geopotential Heights’ write(9, ’(a)’) ‘u 7 99 U winds’ write(9, ’(a)’) ‘v 7 99 V winds’ write(9, ’(a)’) ‘q 5 99 Specific Humidity’ write(9, ’(a)’) ‘endvars close(9) 2、Grads数据准备： （1）ncep资料的直接读取 （2）mm5_to_Grads,wrf_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制 data zlev/1000,925,850,700,500,400,300,250,200,150,100/ tim=0.0；ilev=11；nflag=0 open(2,file=filename.dat,form=unformatted,access=direct,recl=73*4) irecord=0 do j=1,irec irecord=irecord+1 write(2,rec=irecord)abc//sta(j)(1:5),lat(j),lon(j),tim,ilev,nflag, + (zlev(i),z(j,i),t(j,i),td(j,i),d(j,i),f(j,i),i=1,11) enddo ilev=0 ire