GrADS绘图和编程值得收藏.ppt

GrADS绘图与编程;图例1 1月份500hPa位势高度场;图例2 500hPa高度场（1、4、7、10月）;图例3 亚洲季风区7月850hPa流场图;图例4 SST和Sea Level分布图;图例5（130E，25N）处的感热和潜热的时间演变;第一章 概述 ;简 介;安 装;GrADS 软件包的组成;启动和退出GrADS; 所需预备知识DOS和Windows基础 FORTRAN编程 WORD基础 学习要求强调上机编程实践;第二章 GrADS 数据格式; x (或lon) 从西向东的水平坐标y (或lat) 从南到北的水平坐标z (或lev) 从地面到高空的垂直坐标t (或time) 时序坐标。  ;;说出下列维数环境的含义;二、文本格式数据 GrADS格式;顺序：(x,y) z v t; 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=‘v850.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 Do 100 m=1,12 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) Irec=irec+1;Read(3,1000) Read(3,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(2,1000) Read(2,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(4,1000) Read(4,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(5,1000) Read(5,3000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(6,1000) Read(6,4000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 ;100 continue 200 continue 1000 format(2i7) 2000 format(37f6.2) 3000 format(37f8.1) 4000 format(37f7.2) end;三、数据描述文件;;;;;;地形高度资料数据描述文件;第三章 GrADS绘图和参数设置 ;以BAR形式画出沿35N的地形高度剖面图（基点取为3000米）;画出地形高度等值线图（间隔500M，等值线标注为500M字样，大于3000M的区域打阴影）;画出1、7月份的H500气候场;画出84年1月、7月的H500距平场（负值区加阴影）;画出UV850的1月份气候流场图（在气旋和反气旋中心位置标注字符C和A，其中地形高度大于1500M的区域画阴影）;???出82年3月、83年4月、85年7月合成的H500场;画出H500单点相关图相关基点为（100~120E，20~30N）平均;逐月显示82年1月~85年12月的UV850风矢量图;画出（60~100E，5~10N）平均的U850-U200的逐月变化曲线在大于30和小于-15的峰谷值点位置画●;在北半球极地投影下标经度值