《GrADS绘图与编程》课程实习
GrADS绘图与编程(值得收藏)
二、文本格式数据
GrADS格式
现有如下资料： 名称： U850，V850，U200，V200，H500 和TSFC 范围： 60~150E，0~40N 分辨率：2.5×2.5 时间： 1982年1月~1985年12月的逐月资料。
如何将上述资料写成GrADS下的数据格式？
lon) y (或lat) z (或lev) t (或time)
从西向东的水平坐标 从南到北的水平坐标 从地面到高空的垂直坐标 时序坐标。
如 ： set lon 40 160 ；set lat 0 60 定义了水平变化范围 set lat 30；set lon 0 180 定义了沿30N的纬向变化范围 set time jan81 dec94 定义了从81年1月到94年12月的时段
z
(X,Y) (X,Y)
v
t
• T=1 • T=2 • T=3
U V H T
U200
分(z)
时(v)
5月20日(t)
秒(x,y)
18
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=„v850.dat‟) Open(5,file=„h500.dat‟) Open(6,file=„tsfc.dat‟)
7
简
介
• GrADS 是Grid Analysis and Display System的缩写,它由美国马里兰大学气象 系Brian E. Doty开发。利用该系统可实现 包括格点数据和站点数据的彩色气象图 形分析和显示。具有操作简便、功能丰 富、图形美观、显示快速的特点。在国 内外气象界得到广泛应用。 • 版本：V1.5, V1.6.8, V1.7.7，V1.7.9等
实验三 GrADS软件的应用(2)
实验三 GrADS软件的应用（2）一.实验目的1.掌握一维线形图的绘制。
2.掌握fnl数据的转换方法。
3.熟悉GrADS的生成ctl文件的常用命令。
4.熟悉GrADS的生成idx文件的常用命令。
5.掌握GrADS画图的方法。
6.掌握查看画好的图片的方法。
二.实验内容1.实验代码解释'open ' 打开数据文件'set lat 40' 设定经度变化为东经40度'set lon -180 0' 设定纬度变化为南纬180度到0度'set t 1' 设定时次固定为数据集中第一个时次'set gxout line' 设定gxout线'set ccolor 2' 设置等值线颜色为红色'set cmark 3' 设置线上的标记符号为实心圆'set cstyle 1' 设置等值线形为实线'set cthick 6' 设置等值线线宽为6'set grid on 3 3' 设定网格线为红色实现2.绘制线性图3.转换fnl文件4.运行fnl文件生成图像三.实验步骤1.绘制线性图把事先做好的.gs,.ctl和.dat文件复制到D:\grads19下面，打开grads.exe,将.gs文件拖动到dos界面，回车就会在rcl窗口显示所画图像，对应在D:\grads19也会生成一个.png文件；2..转换fnl文件（1）生成描述文件.ctl：打开grads.exe，对应输入gs->grib2ctl grib_file>grib_file.ctl,这样就可以生成描述文件.ctl（2）利用grads自带的grabmap.exe生成索引文件gs-> !gribmap -v -i D:/grads19/grib_file.ctl > a.txt,这样就会在D:/grads19路径下生成.idx 的索引文件，并会把整个映射过程写入a.txt文件中，这样就可以开始用grads读取fnl文件进行画图3.运行fnl文件生成图像将fnl_20090808_00_00_c文件对应的.gs文件拖到里面，即可生成它的.gmf文件。
实验二 GrADS软件的应用
实验二 GrADS软件的应用一.实验目的1.掌握GrADS(Grid Analysis and Display System)软件的安装。
2.掌握GrADS软件的启动方法。
3.熟悉GrADS的常用命令。
4.了解GrADS软件包内的常用文件。
5.掌握GrADS画图的方法。
6.掌握查看画好的图片的方法。
二.实验内容1.软件安装2.启动软件方法开始运行cmd 然后输入C：\grads19 \win32 \grads.exe或打开C盘-grads19-win32-grads.exe3.常用命令Reinit 重新初始化GrADS设置，同时关闭所有打开的文件，并释放所有定义的对象open 打开数据文件set t 9 设定时次固定为数据集中第九个时次set lon 90 140 设定经度变化从东经90度至140度set lat 10 60 设定纬度变化从北纬10度至北纬60度set mpdset cnworld 设置地图数据集set map color style thickness 用定制的颜色，线性和线宽绘背景地图draw map 低级绘图指令，直接进行所指定图形元素操作set lev 200 设定高度维数固定为200h Pa 等压面enable print 打开磁盘文件，用于存放随后print命令转换生成的当前屏幕图形的图元数据set grads off 不打印出GrADS标记set grid off' 不画网格set cthick 6 设置等值线线宽set clopts 1 6 0.12 设置等值线标记的颜色set xlopts 1 6 0.18 设置等值线x轴标记的颜色set ylopts 1 6 0.18 设置等值线y轴标记的颜色set gxout contour 二维等值线绘图set cint 4 设置等值线间隔d smth9(u) 温度Print 将当前图形窗口重的图形转化为图元数据，存放在先前enable 命令指定的磁盘文件中，每执行一次prin t 即向该文件中附加一幅窗口图形4.GrADS软件包内的常用文件.dat 二进制无格式记录的原始气象数据文件，其数据既可以是格点数据，也可以是站点数据。
GRAS实习作业
实习一1）编写简单的gs文件，分别绘制一个简单的一维图形（曲线）、二维图形（等值线），并生成图形文件（三种格式），学习用gv32.exe查看并转换gmf图形文件。
<1>'reinit''open d:\grads\data.ctl''enable print d:\grads\homework\ex01_1.gmf''set gxout line''set x 1 37''set y 1''set z 1''set t 1''d u''print''disable print'*'printim d:\grads\homework\ex01_1.png';<2>'reinit''open d:\grads\data.ctl''enable print d:\grads\homework\ex01_2.gmf''set gxout contour''set x 1 37''set y 1 17''set z 1''set t 1''d u''print''disable print'*'printim d:\grads\homework\ex01_2.png';法，并在grads命令提示符下运行“Q w2xy 120 20”,看看其转换结果。
实习三1）掌握line\bar\contour\shaded形式的绘图；学会fwrite的用法。
注：line与contour参考实习一<1>'reinit''open d:\grads\data.ctl''enable print d:\grads\homework\ex03_1.gmf''set gxout bar''set x 1 37''set y 1''set z 1''set t 1''d u''print''disable print'*'printim d:\grads\homework\ex03_1.png';<2>'reinit''open d:\grads\data.ctl''enable print d:\grads\homework\ex03_2.gmf' 'set gxout shaded''set x 1 37''set y 1 17''set z 1''set t 1''d u''print''disable print'*'printim d:\grads\homework\ex03_2.png';2）掌握阴影图和等值线一起绘图时的顺序；'reinit''open d:\grads\data.ctl''enable print d:\grads\homework\ex03_2.gmf''set gxout shaded''set x 1 37''set y 1 17''set z 1''set t 1''d u''set gxout contour''d u''print''disable print'*'printim d:\grads\homework\ex03_2.png';3）自定义绘制一张line\bar图（同样资料叠加）。
Grads实习报告2
实习三1、所用资料“data”文件夹下二进制月平均风场文件：uv.grd 数据描述文件：uv.ctl时间范围：2002.1-2005.12 共48个月水平范围：37×17个网格点，分辨率为2.5°×2.5°经纬度范围：60-150°E，0-40°N垂直层次：850hPa、200hPa物理量：u，v2、gs文件①500hpa高度场shixi3.1.gs'open D:\data\hgt.ctl''printim D:\data\hejiayue20111368018.gif''set grid off''set time jul2002''set lon 60 150''set lat 0 40''set lev 500''set t 7''set gxout contour''set cint 4''set ccolor rainbow''set clskip 1''draw title 500hpa gaoduchang hejiayue20111368018'；②850hpa风场shixi3.2.gs'reinit''open D:\data\uv.ctl''set parea 1 8 0.5 2.5''printim D:\data\hejiayue20111368018.gif''set grid off''set time jul2002''set lon 60 150''set lat 0 40''set lev 850''set t 7''set gxout vector''d u;v''set parea 1 8 3 5''set gxout barb''d u;v''set parea 1 8 5.5 7.5''set gxout stream''d u;v''draw title 850hpa fengchang hejiayue20111368018'；3、①500hpa高度场图②850hpa风场图实习四1、所用资料①地形高度资料：dxgd.grd\dxgd.ctl 范围：2.5-360 E，-90-90N（全球）层次：1层，地面；时间：2002年1月变量名：h分辨率：2.5*2.5，格点数：144*73 ②风场资料：uv.grd\uv.ctl范围：60-150E，0-40N层次：风场，850 hPa、200hPa；时段：2002.1~2005.12共48个月变量名：u、v分辨率：2.5*2.5，格点数：37*17 2、gs文件shixi4.gs'reinit''open D:\data\dxgd.ctl''open D:\data\uv.ctl''printim D:\data\shixi4.gif white''set grid off''set gxout shaded''set t 1''set lon 60 150''set lat 0 40''set clevs 1500''set ccols 0 1''d h(z=1)''set gxout stream''set lev 850''set time jul2003''set t 7''set lon 60 150''set lat 0 40''draw title 850hpa liuchangtu hejiayue20111368018' 'd u.2;v.2';3、850hpa 流场图4、实习小结通过这次实习，我进一步掌握了gs文件的编写方法，懂得了如何绘制相应时刻和高度的风场图、高度场图及流场图，懂得了怎样绘制三种不同的风场图并将他们放在同一张图中，懂得了怎样给图片添加标题，还知道了怎样绘制阴影区及如何运行gs文件。
GrADS绘图学习技巧与实例(阿木)
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！（笔者：阿木）1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。
软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推荐学习者安装2.0版本，选择默认安装路径就可以。
2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。
3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20) sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,nz,nt)real sstopen(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件do it=1,ntdo iz=1,nzread(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！read后只有文件号，数据是无格式的enddoenddodo it=1,ntdo iz=1,nzwrite(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddoclose(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。
绘图实习报告
绘图实习报告绘图实习报告六篇随着社会不断地进步，报告对我们来说并不陌生，不同种类的报告具有不同的用途。
那么什么样的报告才是有效的呢？以下是小编帮大家整理的绘图实习报告7篇，欢迎大家分享。
绘图实习报告篇1选择了装修公司作为实习地点，说实在的，选择在装修公司实习我感觉很对不起我所学的计算机辅助设计与制造。
但是我真的很不想将来在工厂里度过我的后半生。
不过还好我专业所涉及到的软件知识一样有用武之地。
虽然不是从事机械设计，但是我很喜欢这份工作。
在装修公司里除了可以了解建筑方面的工作，还有机会下工地了解施工方面的知识。
公司只有三十个人，虽然不大，但听说接下的工程还不少呢，每个人都有一台电脑，和我原先想的差不多。
很幸运的是，我刚来就有地方坐，还有台电脑用，正好我可以在这里制图。
现在在公司里有问题可以百度和请教师傅，学起来一定比在家里一个人摸索要快多了。
这周主要是熟悉环境和熟悉CAD在建筑上的表现方法。
主要收获与体会在步入社会的第一步，我要学会稳扎稳打，充分利用学校里学习CAD的经验，把自己最大的优势体现出来。
当然我们在学校中学习的也只是蜻蜓点水一般，很多东西和学校中学习的理论知识还是有一定的区别，所以我需要摆正态度，虚心求教。
绘图实习报告篇2转眼间，实习期即将结束。
停下脚步，回首此程，收获还是很多，也让我成熟了不少了。
现总结此过程的点点滴滴，查漏补缺。
20xx年3月23日，我来到了云南昆明实习，所在的部门是公司的工程部，现所承接的工程项目是云南省博物馆新馆建设项目中的外墙建设项目，也就做外墙的施工。
虽然说我的专业是建筑设计，但施工技术也是一个设计师及建筑师不可缺少的的知识，当时就是抱着这个想法，同时也觉得博物馆项目是也一个很难得的学习机会，来到这个遥远而陌生的成市来挑战自己。
一开始来，还是不适应这里的气候的，海拔高，又干燥，还遇到刚好降温的那阵子，让我一上来就感冒了，虽然感冒，但出来实习不比在家里，不能像在家里那样耍娇气，该工作还是得坚持去工作。
GrADS实习报告
南京信息工程大学 实验（实习）报告实验（实习）名称 基于NCAR/NCEP 再分析资料的降水和气温的气候特征分析日期 2016.10 得分 指导教师 孙晓娟 系 大气科学 专业大气科学班次 长望实验班 姓名 夏新露一、 实习目的（1） 了解NCEP/NCAR 再分析数据记得数据结构。
（2） 通过NCEP/NCAR 再分析数据集，掌握FORTRAN 有格式输入输出语句（Read 、Print 、Write ）的读写四维数据的方法。
（3） 掌握FORTRAN 多维数组的定义和使用方法。
（4） 掌握隐式do 循环的使用方法。
（5） 巩固FORTRAN 顺序结构概念和顺序结构程序设计方法。
（6） 掌握GrADS 对“*.nc ”文件的读取方法。
（7） 掌握GrADS 中fwrite 写出二进制数据的方法。
（8） 掌握GrADS 中While 循环的使用方法。
（9） 掌握GrADS 中变量的定义和使用方法。
（10） 掌握四维“*.grd ”数据格式建立数据描述文件（*.ctl ）的方法。
（11） 巩固GrADS 的基本操作命令的使用方法。
（12） 掌握GrADS 对contour 和shaded 图形类型的设置方法以及两种绘图类型的图形要素设置方法。
（13） 掌握GrADS 中基础绘图命令的设置方法。
（14） 掌握GrADS 中printim 生成文件（*.gmf 、*gif 、*.png ）的方法。
二、 实习内容问题描述利用1948-2010年NCAR/NCEP 月平均气温和降水再分析资料（nc 格式文件），分析1月气温和降水气候特征。
要求利用Fortran 提取1948-2010年1月数据并进行计算，以“*.grd ”格式保存1948-2010年1月气温、降水气候场数据，用GrADS 生成1948-2010年1月气温、降水气候二维等值线和二维填色图两种格式叠加的图像，书写标题，并在高温、低温中心标记“H ”和“L ”。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《GrADS绘图与编程》课程实习学生姓名x x学号xxxxxxxxxxx院系大气科学学院专业大气科学任课教师x x x二Ｏ一四年六月十六日GrADS实习1、资料介绍：范围：60~150E，0~40N.要素及层次：850、200hPa经（v）、纬(u)向风场；500h高度场(h)；地表海表温度场（TSFC）时段：1982.1~1985.12共48个月。
变量名：u v h tsfc实习一1）编写简单的gs文件，分别绘制一个简单的一维图形（曲线）、二维图形（等值线），并生成图形文件（三种格式），学习用gv32.exe 查看并转换gmf图形文件。
①二维图形（等值线）：'reinit''open d:\grads\mine\data.ctl''set lon 60 150''set lat 0 40''set t 48''set lev 850''set gxout contour''set grads off'*'enable print d:\grads\mine\pictures\1contour.gmf''d u'*'print'*'disable print''printim d:\grads\mine\pictures\1contour.png png white' *'printim d:\grads\mine\pictures\1contour.gif gif white';②一维图形（曲线）'reinit''open d:\grads\mine\data.ctl''set lon 120''set lat 25 ''set time jan1982 dec1985''set z 1''set gxout line'*'set ccolor 2''set cstyle 1''set cmark 3''set cthick 6'*'set grid off''set grads off'*'enable print d:\grads\mine\pictures\1line.gmf''d u'*'print'*'disable print''printim d:\grads\mine\pictures\1line.png png white'*'printim d:\grads\mine\pictures\1line.gif gif white';2）熟悉三类坐标的用法，掌握查询命令Q，特别掌握其坐标转换命令的用法，并在grads命令提示符下运行“Q w2xy 120 20”,看看其转换结果。
另外，熟悉“Q time/q file”等用法实习二编写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=1DO 200 iy=1,4do 100 m=1,12! write u850read(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! write u200read(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! write v850read(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! write v200read(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! write h500read(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! write tsfcread(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+1100 continue200 continue1000 format(2i7)2000 FORMAT(37F6.2)3000 format(37f8.1)4000 format(37f7.2)End②ctl文件:dset d:\grads\mine\dat\data.grdundef -9.99E+33title NCEP/NCAR REANALYSIS PROJECT xdef 37 linear 60.000 2.500ydef 17 linear 0.000 2.500zdef 2 levels 850 200tdef 48 linear JAN1982 1movars 4U 2 99 u wind (m/s)V 2 99 v wind (m/s)H 1 99 H500TSFC 1 99 TSFC DATAendvars实习三1）掌握line\ bar \contour\shaded形式的绘图；学会fwrite的用法。
①bar'reinit''set grads off''open d:\grads\mine\data.ctl''set lon 120''set lat 25''set time jan1982 dec1985''set z 1''set gxout bar''set bargap 20''d u''printim d:\grads\mine\pictures\3bar.png png white'；②shaded'reinit''open d:\grads\mine\data.ctl''set lon 60 150''set lat 0 40''set t 48''set lev 850''set gxout shaded'*'set cint 2''set clevs 0'*'set ccols 14''set csmooth on''set grads off''d u''printim d:\grads\mine\pictures\3shaded.png png white'；2）掌握阴影图和等值线一起绘图时的顺序；3）自定义绘制一张line\bar图（同样资料叠加）。
'reinit''open d:\grads\mine\data.ctl''set lon 120''set lat 25''set time jan1982 dec1985''set z 1''set gxout bar''set bargap 20''set grads off''d u''set gxout line''set cthick 6''set ccolor 1''d u''printim d:\grads\mine\pictures\3barline.png png white';4）自定义绘制一张contour\shaded图（同样资料叠加）'reinit''open d:\grads\mine\data.ctl''set lon 60 150''set lat 0 40''set t 48''set lev 850''set gxout shaded''d u''set gxout contour'''set cterp on''d u''printim d:\grads\mine\pictures\1countourshaded.png png white';5）从原始资料中写出20-30N，80-100E范围内资料，并配以数据描述文件，简单绘图验证正确与否。
①gs文件：'reinit''open d:\grads\mine\data.ctl''set fwrite d:\grads\mine\pictures\fwrite.grd''set gxout fwrite''set lon 80 100''set lat 20 30''set t 1 48''set lev 200''d u''disable fwrite';②数据描述文件dset d:\grads\mine\pictures\fwrite.grdundef -9.99E+33title NCEP/NCAR REANALYSIS PROJECT xdef 37 linear 80.000 0.540ydef 17 linear 20.000 0.588zdef 2 levels 850 200tdef 48 linear JAN1982 1movars 2U 2 99 u wind (m/s)V 2 99 v wind (m/s)endvars实习四1）求出500hPa高度场的气候场、距平场，并给出图形。