当前位置：文档之家 › GrADS实用绘图--数据准备

GrADS实用绘图--数据准备

  • 格式：pdf
  • 大小：520.35 KB
  • 文档页数：42

zdef 7 levels 1000 850 700 500 300 200 100
• 文件路径推荐用UNIX惯例写
– d:/mydata/test.grd
10
LASG动力论坛版
五. 数据准备(Gr)
– – – – – – – – – yrev zrev sequential 365_day_calendar template byteswapped big_endian little_endian cray_32bit_ieee
17
LASG动力论坛版
五. 数据准备(Gr)
%ix1 %ix3 %iy2 %iy4 %im1 %im2 %imc %id1 %d1 %id2 %ih1 %ih2 %ih3 %in2
— .ctl文件语法9-4
tdef 中起始时间的替代方式
起始的1 位数年代(v1.9) 起始的3 维数年代 (v1.9) 起始的2 位数年 起始的4 位数年 起始的1/2 位数月 起始的2 位数月(必要时前置0) 起始的3 字符月 起始的1/2 位数天(必要时前置0) 1 位数字天 起始的2 位数天 起始的1/2 位数时 起始的2 位数时 起始的3 位数时 起始的2 位数分(必要时前置0)
— .ctl文件语法7-1
pdef关键字
－lcc | lccr －nps | sps －BILIN | file －eta.u | pse －ops
－特殊投影坐标系下的数据 Lambert投影 极射赤面投影 其它直线网格投影 US NMC eta模式 | SSM/I极射赤面投影 CSU RAMS斜极射赤面投影
3
LASG动力论坛版
五. 数据准备(Gr)
本节内容提要
– .ctl文件语法
– – – – – –
GrADS支持的数据格式 Binary数据格式 netCDF数据格式 GRIB数据格式 站点数据格式 进阶向导
4
LASG动力论坛版
五. 数据准备(Gr)
• 3-D空间数据排列方式
– 自西向东 – 从南到北 – 从低层到高层
15
LASG动力论坛版
五. 数据准备(Gr)
文 件 名 中 时 间 格 式 说 明 符
%x1 %x3 %y2 %y4 %m1 %m2 %mc %d1 %d2 %h1 %h2 %h3 %f2 %f3 %n2
— .ctl文件语法9-2
(如: 1910, 1920, 1930) (如: 1980, 1990, 2000) (如: 1, 2, …, 11,12) (如: 01, 03, …, 12) (如: Jan, Feb, Mar, …) (如: 1, 2, …, 31) (如: 01, 02, 10, 31) (如: 1z, 2z, …, 12z) (如: 01, 02, 10, 31) (如: 120, 012) (如: 12, 24, 48, 72, 120) (如: 012, 024, 072, 120) (如: 01, 35, 60)
HDF4(SDS)
20
LASG动力论坛版
五. 数据准备(Gr)
Binary格式数据
— Binary格式1 ♥♥♥
– GrADS最基本, 我们最容易准备 – 常用后缀 .dat .grd .bin ,…… – 需要 .ctl 文件一起使用 – 是无格式文件 (form=‘unformatted’) – 若用Fortran 读写
– 方向: 自西向东, 自南向北 (永远不变) – 单位: 度 纬度/经度 – 正负值: 以东为正/ 以北为正
7
LASG动力论坛版
五. 数据准备(Gr)
• Z 方向的维数说明
— .ctl文件语法3
zdef 7 levels 1000 850 700 500 300 200 100 列举法 给出 坐标值
– 方向: 自下而上 – 单位: hPa / 其他
可用linear 关键字替换
8
LASG动力论坛版
五. 数据准备(Gr)
• T 方向的维数说明
— .ctl文件语法4
1mo 时间 步长 3mn 6hr 1dy 1mo 1yr 分钟 小时 天 月 年
– tdef 12 linear 00:00z21jan2000
6
LASG动力论坛版
五. 数据准备(Gr)
• X, Y 方向的维数说明
– xdef
144 linear 坐标值 线性 增加
— .ctl文件语法2
0 起始 坐标 值 2.5 坐标 间隔 可用levels 关键字替换
Y: gausT62 gausR15 gausR20 gausR40
X(lon) 共有 方向 144个 维数 格点 定义
12
LASG动力论坛版
五. 数据准备(Gr)
pdef关键字
－使用注意
— .ctl文件语法7-2
－u,v分量需与GRIB数据格式中风场变量33和34的说明格式 相同, 以便执行矢量的旋转运算
u 18 V 18 u10 0 v10 0 33,100 34,100 33,105 34,105 U-Wind Components on Pressure Levels V-Wind Components on Pressure Levels 10 Meter U Wind 10 Meter V Wind
— .ctl文件语法6
• options 是可选的关键字, 其选项有
自北向南的数据 自高层到低层的数据 顺序存贮的数据 任何年份均是365天 多个文件 文件存储的字节次序与本机相反 在SUN, SGI等机器上生成的数据 在PC, DEC等机器上生成的数据 在Cray机上生成的数据
11
LASG 数据准备(Gr) 五. 动力论坛版
— .ctl文件语法1
• 格点数据描述文件 .ctl 的写法
dset d:/mydata/test.grd title This is a GrADS test data options sequential yrev undef -9999.0 xdef 144 linear 0 2.5 ydef 73 linear -90 2.5 zdef 7 levels 1000 850 700 500 300 200 100 tdef 12 linear 00:00z01jan2000 1mo vars 5 slp 0 99 Sea Surface Temperature (hPa) t2m 0 99 2m temperature( degC) temp 7 99 upper level temperature (degC) uwnd 7 99 upper level U component of wind (m/s) vwnd 7 99 upper level V component of wind (m/s) endvars
1位数字表示的年代 3位数字表示的年代 2位数字表示的年 4位数字表示的年 1位或2位数字表示的月 2位数字表示的月 3位字符表示的月 1位或2位数字表示的日 2位数字表示的日 1位或2位数字表示的小时 2位数字表示的小时 3位数字表示的小时 2位或3位数字表示的预报时效 3位数字表示的预报时效 2位数字表示的分钟
dset ^pr.%ch.nc chsub 1 600 1851-1900 chsub 601 1800 1901-2000 options template dtype netcdf tdef 1800 linear jan1851 1mo
若两个文件分别位于不同的磁盘上:
dset %ch chsub 1 600 /disk1/pr.1851-1900.nc chsub 601 1800 /disk2/pr.1901-2000.nc
1
LASG动力论坛版
GrADS实用绘图
第三讲 数据准备
＝Funny @ LASG动力论坛＝
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
2
LASG动力论坛版
内容提要
请注意标记♥♥♥的地方
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1 — 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2 — 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
16
LASG动力论坛版
五. 数据准备(Gr)
— .ctl文件语法9-3
dset %ch 和 chsub t1 t2 string 若有如下NetCDF 数据文件, pr.1851-1900.nc包含了50 年 的逐月资料(即第1-600月), pr.1901-2000.nc包含100 年(第 601-1800月), 则:
pdef isize jsize { lcc | lccr } latref lonref iref jref Slat Nlat std_lon dx dy lcc 表示数据中风分量已经旋转成经纬度坐标, lccr 表示风分量仍然是lambert 正型投影,需GrADS 对其进行旋转. isize, jsize X, Y 方向的总点数; latref, lonref 参考点的经, 纬度; iref, jref 参考点的(i,j)位置; Slat, Nlat Lambert 投影m=1 时南, 北两个纬度; std_lon 标准经度; dx, dy 以M 为单位的X 和Y 方向的网格距.
— .ctl文件语法9-1
• 用一个 .ctl 文件描述多个数据文件
dset %y4%m2%d2.grd options template … … tdef 365 linear %id2%im2%iy4 1dy …
.ctl
• %y4%m2%d2 是时间的替代格式符
%y4 | %iy4 4位数的年份 %m2 | %im2 2位数的月份 %d2 | %id2 2位数的日期

合集下载

GrADS基础绘图

GrADS基础绘图


绘图速度较慢，当数据文件较大时，读入很容易内存溢出 远程处理比较麻烦。

grads的优点; （1），最多能读入20个数据文件，数据读取迅速，出图速度快，但是要写描述文件和gs文 件比较麻烦。 （2），在画二维等值线图时，图形非常平滑，这点是我到现在还没大规模用matlab的主要 原因。 （3），它包含了很多气象上所需要的函数，在gs文件里可以直接调用，因此在用它时也省 了不少事。
Grads在win/Linux下的下载、安装、运行

程序：


Windows：grads-1.8sl11-win32e.exe Linux：grads-bin-1.9b4-linuxRH9.tar.gz grads-1.8sl11-linux.tar.gz data.tar.Z

环境变量：

Windows：




Panoply OpenDX (formerly IBM Data Explorer) NCL (NCAR Command Language) ncBrowse (netCDF File Browser) HDF (Hierarchical Data Format) interface GrADS (Grid Analysis and Display System) FERRET NetCDF Toolbox for MATLAB-5 IDV (Integrated Data Viewer) Origin surfer
命令顺序 grads –b|l|p|c open/run set display draw print quit/reinit
显示类型


GrADS绘图实习一数据文件的转换及数据描述文件的建立

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基础绘》课件


GrADS的主要应用领 域
GrADS被广泛应用于气象 学、气候学、海洋学、环 境科学等领域的数据处理 和可视化。
安装及环境配置
1
GrADS的安装步骤
详细介绍如何下载、安装和配置GrADS软件。
2
GrADS环境的配置
探索如何设置和调整GrADS的环境变量以适应个人需求。
3
GrADS的基本设置
了解GrADS的基本设置和首选项，以优化您的工作流程。
GrADS与其他软件集成
GrADS和Python的 数据交互
探索如何结合GrADS和 Python的强大功能，进行数 据交互和自动化处理。
GrADS与Matlab的 数据联动
了解如何在GrADS和Matlab 之间无缝传递数据，实现更 多分析和可视化功能。
GrADS在GIS领域的 应用
发现如何将GrADS与地理信 息系统（GIS）工具结合，实 现空间数据的分析和可视化。
总结
GrADS的优缺点评价
综合评价GrADS的优点和不足，帮助您了解 该工具是否适合自己的数据可视化需求。
GrADS使用心得交流
分享使用GrADS的一些心得体会和实用技巧， 促进经验交流和学习。
GrADS高阶绘图技能
GrADS高级绘图基础
学习如何创建精美且高度可定 制的数据可视化图表，彰显数 据的特点和趋势。
GrADS图像调整及美化
发掘GrADS图像调整和美化工 具，使您的图表更具吸引力和 专业性。
GrADS数据可视化技巧
了解如何创建交互式、动态和 多维数据可视化，提升数据沟 通的效果。
GrADS案例分析
1
GrADS数据处理应用案例
2
分享一些使用GrADS进行数据处理和
GrADS画图入门

GrADS画图入门

设置 set gxout shaded
（缺省）
命令：
display express
ga-> d strp Contouring: 0 to 100 interval 10
设置：
set *** ###
维数参数设置 等值线设置
Set set Set set lon 100 cint 15 120 指定等值线间隔 lat 103 25 指定等值线颜色 ccolor
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
• 3、图形显示类型及绘图设置：
•
2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5_to_Grads,wrf_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
data open(3,file=filename.ctl,form='formatted',status='unknown') zlev/1000,925,850,700,500,400,300,250,200,150,100/ write(3,'(a)')'dset ^'//filename.dat tim=0.0 ；ilev=11；nflag=0 open(2,file=filename.dat,form='unformatted',access='direct',recl=73*4) write(3,'(a)')'dtype station' irecord=0 write(3,'(a)')'stnmap ^sounding.map' do j=1,irec write(3,'(a)')'undef -999.0' irecord=irecord+1 write(3,'(a)')'title sounding' write(2,rec=irecord)'abc'//sta(j)(1:5),lat(j),lon(j),tim,ilev,nflag, write(3,'(a)')'tdef 1 linear 02JAN1987 1DY ' + (zlev(i),z(j,i),t(j,i),td(j,i),d(j,i),f(j,i),i=1,11) write(3,'(a)')'vars 5' enddo write(3,'(a)')'zz 11 99 '//staa//' height1' ilev=0 write(3,'(a)')'tt 11 99 '//staa//' temper1' irecord=irecord+1 write(3,'(a)')'t_td 11 99 '//staa//' t_td1' write(2,rec=irecord)'abcdefgh',lat(1),lon(1),tim,ilev,nflag, write(3,'(a)')'uu 11 99 '//staa//' u wind1' + write(3,'(a)')'vv (zlev(i),i=1,11),(z(1,i),i=1,11),(t(1,i),i=1,11), 11 99 '//staa//' v wind1' + write(3,'(a)')'endvars‘ (td(1,i),i=1,11),(d(1,i),i=1,11),(f(1,i),i=1,11) close(2) close(3)