《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 ”。
GrADS绘图实习
实习一数据文件的转换及数据描述文件的建立一、实习资料:现在“data”文件夹下有ASCII码数据资料文件:月平均风场 (u200.dat、u850.dat、v200.dat、v850.dat)、月平均高度场 (hgt500.dat),2002.1-2005.12共48个月,资料的水平网格范围是M×N个网格点(M=37,N=17),分辨率为2.5°×2.5°,自西向东经度60-150°E,由南至北纬度0-40°N。
二、实习要求:编写出将ASCII码数据资料文件:月平均风场和月平均高度场分别转换成二进制无格式直接存取(GrADS数据格式)文件(uv.grd和hgt.grd文件)的Fortran程序,给出相应的数据描述文件(uv.ctl和hgt.ctl文件),并利用GrADS基本操作命令显示该地区850hPa和200hpa uv风场,最后保存图形。
三、实习目的:通过编写该程序,进一步熟悉数据GrADS的数据格式,熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。
四、实习步骤:4.1 利用“写字板”程序打开*.dat数据文件,熟悉该数据资料;4.2 按要求编写Fortran程序,将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件,结果保存为mhy.grd;4.3 通过“写字板”或“记事本”程序编写相应的数据描述文件,保存为mhy.ctl;4.4 画出2002年1月850hpa风场图,与图1比较验证数据正确性;4.5 利用GrADS基本操作命令(open,display,print等)显示2002年7月850hpa和200hpa uv 风场图形;4.5 完成实习报告1)说明所用资料资料名称:月平均风场(u200.dat、u850.dat、v200.dat、v850.dat);月平均高度场(hgt500.dat)范围:60~150E,0~40N.层次:风场,850 hPa、200hPa;高度场,500hPa;时段:2002.1~2005.12共48个月。
GrADS实习报告
南京信息工程大学滨江学院GrADS实习报告院系:大气与遥感系专业:大气科学年级:14级班级:1学号:***********姓名:**实习一第一章气象程序设计与绘图软件安装与运行1.1实验实习目的1.1.1了解并掌握FORTRAN 90软件开发环境及基本操作1.1.2了解并掌握GrADS2.0的软件环境及基本操作方法(1)GrADS的安装(2)GrADS的启动与退出(3)GrADS文件处理(4)GrADS的使用流程1.2实验实习内容11.2.1问题的描述现有北京1951—1980年1月月平均气温文件beijingT.dat,编写程序计算北京30年的平均气温,并输出结果。
北京1951—1980年1月平均气温为1.0,—5.3,—2.0,—5.7,—2.1,—0.6,—1.7,—3.6,—3.0,0.1,—2.6,—1.4,—3.9,—4.7,—6.0,—1.7,—3.4,—3.1,—3.8,—2.0,—1.7,—3.6,—2.7,—2.4,—0.9,—2.7,—1.6,—3.91.2.2算法设计通过对问题的分析,设计求解算法,并绘制流图1.2.3程序编写program maininteger,parameter:: n=30realtave,sum,temp(n)open(1,file='c:users\bjxy\desktop\grads\chap1\data\beijingT.TXT',FORM='formatted' )do i=1,nread(1,*)temp(i)end doclose(1)open(2,file='c:users\bjxy\desktop\grads\chap1\data\tave.dat',FORM='formatted') sum=0.0do i=1,nsum=sum+temp(i)end dotave=sum/nwrite(2,*)taveclose(2)End1.2.4实验实习要求(1)创建新工作空间test1,工作空间文件夹创建在D盘上(2)创建新工作空间test1内创建新项目xml,项目文件夹创建在工作空间文件夹内(3)在项目xml内创建源程序文件test1.f90,文件创建在项目文件夹内,编辑输入给定的源程序文本,在前三行“”处输入班级、姓名、日期信息(以后要求相同)(4)在项目xml内创建辅助文档文件miaoshu1.txt,文件创建夹内,在文件中编辑输入问题描述文本(5)在项目xml内创建辅助文件suanfal.doc,文件创建在项目文件夹内,在文件中编辑绘制图1.1所示的流程图(6)编译源程序test1.f90构建可执行程序xml.exe,运行可执行程序xml.exe (7)强输入数据和输出结果以注释形式编辑输入到源程序末尾1.2.5 实验实习步骤(1)Compaq Visual FORTRAN 6.6系统安装后,在桌面创建一个Developer Studio 图标,在“开始”→“程序”菜单中创建Compaq Visual FORTRAN子菜单项,通过桌面图标或“开始”菜单可快速启动Microsoft Developer Studio。
绘图实习报告5篇_1
绘图实习报告5篇绘图实习报告篇1时间过的很快,一转眼就到了CAD实训周结束的时候了,记得在实训的第一天,老师给我们将了这周实训的任务安排,讲述了本周实训的主要内容,实训目的以及意义所在,然后交代了一些细节方面的问题,强调应当注意的一些地方,以及考试成绩打分等。
在这次实训中,让我体会最深的是理论联系实际,实践是检验真理的唯一标准。
理论知识固然重要,可是无实践的理论就是空谈。
真正做到理论与实践的相结合,将理论真正用到实践中去,才能更好的将自己的才华展现出来。
我以前总以为看书看的明白,也理解就得了,经过这次的实训,我现在终于明白,没有实践所学的东西就不属于你的。
俗话说:“尽信书则不如无书”我们要读好书,而不是读死书。
现在的社会需要的是动手能力强的人,而不是理论好的人,对于我们高职生来说这更重要了,我们学校的目标就是培养一流的职业技能的人。
虽然时间很短很紧,但是我一直认认真真去绘制每一个图,思考每一个细节,作图步骤,哪怕是一个很小的问题。
的确,在作图的过程中我遇到了不少的难题,但都在同学的帮助下,一个个的被我解决,自己难免会感到有点高兴,从而增加了对CAD的浓厚兴趣。
虽然在生活中存在很大的不足,但很多方面对我的印象还是很深刻。
比方说,一位来自贵港的小伙子(其年龄只有17岁)一次离奇不归,大家焦急不安,分别发挥各自能力通过不同的手段去寻找不归的他,回来后大家还曾为他平安归来喝酒庆祝。
这让我看到人性的光辉,凡此种种,数不胜举。
在实训中,我们首先掌握了书写技能,如金额大小写,日期,收付款人等等的填写。
接着填写相关的数据资料。
再继续审核和填写原始凭证。
然后根据各项经济业务的明细账编制记账凭证。
根据有关记账凭及所附原始凭证逐日逐笔的登记现金日记账、银行存款日记账以及其他有关明细账;往后我们会编制科目汇总表,进行试算平衡。
根据科目汇总表登记总分类账。
并与有关明细账、日记账相核对;最后要编制会计报表;将有关记账凭证、账页和报表加封皮并装订成册并写实训总结。
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 ”。
保存并分析之。
问题分析已知:1948-2010年NCAR/NCEP 月降水和气温再分析资料。
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,大气的数据要根据数据的层次确定几层。
GRADS学习心得
GRADS学习心得grads是当今气象界广泛使用的一种数据处理和显示软件系统。
该软件系统通过其集成环境可以对气象数据进行读取,加工,图形显示和打印输出。
他在进行数据处理时,所有数据在grads中均被视为纬度,经度,层次和时间的4维场,而具有操作简单,功能强大,显示快速,出图类型多样化,图形美观等特点。
正因为如此,他才成为国内外气象界通用的标准图形环境之一。
grads作为一门应用性的课程,内容包含了若干命令,在学时有限的情况下,我们如何在短时间内掌握其中最核心的内容,做到尽快入门,是有一定的难度的。
老师帮我们理清思路,凝练重点,促进知识的巩固。
这学期,我们在徐老师的指导下花了大约八周的时间联系和学习这门课程。
我开始下载数据,自己用软件画画。
在老师的严格要求下,我们一次又一次地互相帮助。
在绘画成功的那一刻,我们真的感到了巨大的成就感。
首先,从头回忆这个故事。
一开始刚刚上课,我们了解到了一大堆关于grads可以应用的各种文件。
例如(1).dat二进制无格式记录的原始气象数据文件,其既是以是格点数据,也可以是站点数据。
它们是从其它气象数据(如站点气象报、格点气象报、模式格点输出结果)转换生成的。
对格点数据而言,其格式为二进制无格式直接或顺序记录格式。
总之对格点数据,整个数据集是一个大的五维数据场,包括三维物理空间、一维物理变量、一维时间变量,存放时以二维数组片的形式按水平、垂直、物理变量、时间序列的顺序排放,(2).ctl原始数据描述文件。
该文件为纯ascii文件,用以描述原始数据集的基本信息,包括数据集文件名、数据类型、数据结构、变量描述等等,详见§4在grads环境中至少得首先打开(open)一个数据描述文件,以便后续的操作有数据对象。
(3).gsgrads控制文件,用grads命令run执行之。
也是一个纯文本文件。
还了解到十进制的原始数据文件(*.txt)、二进制的数据文件(*.dat)、数据描述文件(*.ctl)、批处理文件(*.gs)和图形文件(*.gmf)。
Grads实习
Grads实习实习二实习要求:1)利用所提供的数据资料,绘制一张第2时刻温度垂直剖面图(纬度固定于30ºN,经度范围80ºW-0º,层次由第1层变化至第7层),要求将温度正值区域用阴影显示,负值区域用等值线显示,并加粗0值等值线,同时Y轴标注为“Z/hPa”,标题标注为“Temperature/ºC”。
2)采用流线方式绘制第3天西北半球(-180w-0,0-90N)850hPa风场图,并将图片输出到850uv.gmf图片文件中。
3)恰当利用所提供的数据资料,绘制一张北半球地图底图,即图形中只显示地图信息,不含其它形式的数据信息。
'reinit''open D:\grads\grads1\test2\model.le.ctl''set lat 30''set lon -80 0''set z 1 7''set t 2''set grads of''set gxout shaded''set cmin 0''d t-273.15''cbar''set gxout contour''set cmax 0''d t-273.15''set clevs 0''set ccols 3''set cthick 7''d t-273.15''draw ylab Z/hPa''draw title Temperature/`3.`1C''printim D:\grads\grads1\test2\1.gif gif';'reinit''open D:\grads\grads1\test2\model.le.ctl' 'set lat 0 90''set lon -180 0''set lev 850''set t 3''set gxout stream''d u;v''printim D:\grads\grads1\test2\aa.gif';;‘open d:\grads\grads1\test2\model.le.ctl’‘set lat 0 90’‘set grads off’‘set cmax 0’‘set mpdset hires’‘set map auto’‘set mpdraw on’‘set poli on’‘d z’‘printim d:\grads\grads1\test2\ssss.gif gif’;实习三实习要求:1)根据所提供的数据资料,利用fwrite命令将该范围内(纬度:0-90ºN,经度:180ºW-0º,层次为500hPa,第2时刻)高度场资料写入名为z.grd的二进制数据文件中。
GrADS绘图与编程(值得收藏)
由GrADS.exe文件调用 产生格点资料映射文件(*.idx)
看Grib码资料,转为TXT格式文件
显示*.gmf格式的图形文件 产生台站资料的映射文件(*.map)
Grib码资料的解读程序
将*.gmf格式的图形文件转换为*.wmf 格式
10
启动和退出GrADS
100 continue 200 continue 1000 format(2i7) 2000 format(37f6.2) 3000 format(37f8.1) 4000 format(37f7.2) end
三、数据描述文件
Data.ctl 文件内容:
Dset data.grd Undef –9.99e+33 Title Ncep/Ncar reanalysis project Xdef 37 linear 60 2.5 Ydef 17 linear 0 2.5 Zdef 2 levels 850 200 Tdef 48 linear jan1982 1mo Vars 4 U 2 99 u wind (m/s) V 2 99 v wind(m/s) H 1 99 H500 T 1 99 Tsfc data endvars
11 8.5
11
所需预备知识
DOS和Windows基础 FORTRAN编程 WORD基础
学习要求
强调上机编程实践
12
第二章 GrADS 数据格式
• 目标 (1)能将文本格式数据转为GrADS格式
(2)能写出数据描述文件(*.ctl)
重要概念
一、维数环境的概念
1、含义:GrADS视每一个变量(VAR)场为一个四维数据集, 即包括三维空间(x,y,z)和一维时间(t)。 2、作用:说明和指定随后的分析或图形操作时参加的原 始数据集的维数范围。 3、定义方法: Set lat|lon|lev|time val1 <val2> 地球坐标 Set x|y|z|t val1 <val2> 格点坐标 13
编程教育实习报告记录
一、实习背景随着我国信息技术的飞速发展,编程教育逐渐成为我国教育改革的重要方向。
为了更好地适应社会需求,提升自身的编程技能和教学能力,我于2021年7月至9月在某知名编程培训机构进行了为期两个月的实习。
此次实习让我对编程教育有了更加深刻的认识,以下是我实习期间的详细记录。
二、实习内容1. 教学观摩在实习初期,我主要进行了教学观摩。
通过观摩资深教师的授课过程,我学习了他们的教学方法和技巧,了解了不同年龄段学生的需求。
在教学观摩过程中,我发现以下几点:(1)教学内容要贴近实际,让学生在学习过程中感受到编程的实用价值。
(2)教学过程要生动有趣,激发学生的学习兴趣。
(3)注重培养学生的逻辑思维能力,提高编程能力。
2. 课程设计在实习中期,我参与了课程设计。
在课程设计过程中,我学习了如何根据学生的年龄、兴趣爱好和认知水平,设计适合他们的编程课程。
以下是我设计的一堂针对小学生入门的Python课程:课程名称:Python入门——认识变量和运算符课程目标:1. 了解Python编程语言的基本概念。
2. 掌握变量的定义和赋值。
3. 熟悉Python的运算符及其使用。
课程内容:1. Python简介2. 变量的定义和赋值3. 运算符及其使用4. 实践操作3. 课堂教学在实习后期,我担任了助教,协助教师进行课堂教学。
以下是我在课堂教学过程中的一些心得体会:(1)关注学生的个体差异,因材施教。
针对不同学生的接受能力,调整教学进度和难度。
(2)善于运用多媒体教学手段,提高课堂效果。
通过图片、动画等形式,让学生更容易理解抽象的编程概念。
(3)鼓励学生提问,培养他们的自主学习能力。
在课堂上,鼓励学生提出问题,教师及时解答,让学生在解决问题的过程中提高编程技能。
(4)关注学生的学习状态,及时发现并解决问题。
在课堂教学中,关注学生的情绪和表现,及时发现并解决他们在学习过程中遇到的问题。
三、实习收获1. 提升了编程技能。
通过实习,我掌握了Python编程语言的基本语法和常用库,提高了自己的编程能力。
grad实习
(2)绘制1948年7月,90E的纬度-高度剖面 图(90S-90N,1000hPa-100hPa)
实习资料: 1948年1月-2009年12月高度场资料
实习报告内容 目的要求 主要内容 资料方法 结果包括: (1) 500hPa1月气候平均高度场 (2) 1948年7月,90E的纬度-高度剖面图
实习内容:
第一部分 掌握GrADS的安装、启动、退出 熟悉GrADS 安装文件夹内容 掌握GrAD文件的建立 实习内容:将给定的十进制原始数据转换为 GrADS要求的二进制形式,并写出相应的描 述文件。 实习资料:
十进制存放格式数据资料文件u.txt、v.txt和 h.txt,空间范围(60-150°E,0-40°N);层 次(850hPa、 500hPa 200hPa);时段(2009.12009..12)逐月;分辨率:2.5*2.5。
ctl文件内容、 1、7月500hPa高度场图 1、7月850hPa矢量风场图
第三部分 GrADS绘图
实习内容:
根据2001年1月-2005年12月逐月资料,编写gs程序完 成以下绘图内容:
(1)求出500hPa高度场距平,绘制2001年12月份和2002年7月份 的位势高度场距平图,同时负值区标出阴影。 (2)求出2001-05年逐月850hPa等压面100E-120E的经向平均u 值,输出数据文件,根据新生成的数据文件及其描述文件绘制 纬 度-时间剖面图。 (3)计算u(850hPa)-u(200hPa),2001-05年逐月的南亚地区 ()的区域平均值,绘制单线图并在峰、谷处标记。 (4)极射赤面北极投影绘制2001年1月200hpa U场,并标出 90E,180E, 90W,180W
2001-05年逐月南亚区域平均u(850hPa)-u(200hPa)时间序列图 2001年1月200hpa U场
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《GrADS绘图与编程》课程实习学生姓名x x学号xxxxxxxxxxx院系大气科学学院专业大气科学任课教师x x x二O一四年六月十六日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高度场的气候场、距平场,并给出图形。