![](https://img-blog.csdnimg.cn/20201014180756724.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
GrADS
文章平均质量分 62
chht
这个作者很懒,什么都没留下…
展开
-
grads精致绘图说
[程序见GrADS Script Library] 1. basemap.gs: basemap L | O | U 在低分辨率海岸廓线范围内用颜色覆盖陆地/海洋。适用于各种投影方式,需lpoly.asc, lpoly_hires.asc, lpoly_US.asc, opoly.asc, opoly_hires.asc文件。 其中:L(l):覆盖陆地,O(o转载 2010-01-25 22:01:00 · 6557 阅读 · 0 评论 -
分省地图使用
GrADS的全国各省地图,只需要把文件解压后复制到grads/dat/文件夹下,使用方法同“mres,hires...”2008年3月对地图进行了重新处理,修正了地图有缺口的问题。 地图数据取自国家测绘局国家基础地理信息中心1:400万基础地理信息国界数据(http://sms.webmap.cn/default.asp) 使用方法为: set mpdset map1 map2转载 2010-01-25 22:13:00 · 1438 阅读 · 0 评论 -
grads 中读写数据fwrite
reinitsdfopen d:/gfdl21/va_A2.19960101-20001231.ncset fwrite F:/homework/synoptic/model/vwnd/v1996.dat*set lat -30 30*set lon 119.5 288.5set x 1 144set y 1 90tt=1yr=0while(ttset t tt*-转载 2010-01-25 22:13:00 · 6378 阅读 · 0 评论 -
关于grads中国地区地图绘制
1、在进入绘图前,首先 输入:set mpdset lowres*表示低精度格式set lon 73 135set lat 18 55set poli off*表示取消行政边界2、上述工作完成后,输入SET gxout shadedset ccolor 1set ccols 10 1SET cmin 2.7set cint 3转载 2010-01-25 22:15:00 · 7669 阅读 · 0 评论 -
grads 批量处理nc文件
若想要提取从1951-2006年56年nc文件中的某些数据,一个一个处理非常麻烦,这里介绍种较为简易的方法。例如想提取6-8月的位势高度资料。reinitt5=1951*作文件名循环while(t5 set gxout fwrite set fwrite D:/sichuan/hgt1/%t5%.dat sdfopen e:/ncep1/hgt/hgt转载 2010-01-25 22:17:00 · 7485 阅读 · 1 评论 -
画底图mask.gs程序的使用
1.打开grads,输入命令d:/mask.gs,产生mask图片,使用点选功能将所要绘制的底图点击出来,然后点击quit,自动产生maskarea.dat二进制文件;2.对maskarea.dat用ctl文件描述maskarea.ctl: DSET D:/maskarea.datUNDEF -999.0TITLE CPC Soil MoistureXDEF 127 LINEAR 72转载 2010-01-25 22:18:00 · 1306 阅读 · 1 评论 -
关于站点数据stnmap命令的使用
特别提示:!stnmap -i d:/eof/td.ctl中路径用反斜杠/,而不是/!1、首先使用fortran程序,将数据按照规定格式读入,包含站点名称、经纬度、观测数据值program bintd parameter (st=548) ! st means nums of station character stid(st)*8 integer nflag,nlev rea转载 2010-01-25 22:19:00 · 3641 阅读 · 0 评论 -
竖放6张图片的grads文件
reinitoutfile=./outfile_6*----------------input--------------------------infile1=./infile.ctlopen infile1*===============(左上1)==================set parea 0.9 4.4 7.3 10.8 set grads offs转载 2010-01-25 22:48:00 · 2964 阅读 · 1 评论 -
浅谈GRADS气象绘图软件的使用
摘 要:通过编制具体的应用程序,介绍了气象绘图软件GRADS的使用方法。 关键词:GRADS 软件 使用 GRADS是当今气象界广泛使用的一种数据处理和显示软件系统。该软件系统通过其集成环境,可以对气象数据进行读取、加工、图形显示和打印输出。它在进行数据处理时,所有数据在GRADS中均被视为纬度、经度、层次和时间的4维场,而数据可以是格点资料,也可以是站点资料;数据格式可以是二进制,转载 2010-01-25 22:52:00 · 4746 阅读 · 0 评论 -
用grads读取txt中经纬度画站点图
用grads读取txt中经纬度,用于画站点图 txt中的站点信息格式如下: stn_num lon lat data1 data2 .... (每个变量用空格格开一行为一个记录) ****************draw grid************************************** open ******.ctl *打开一个格点资料用于画底图的转载 2010-01-25 23:05:00 · 5460 阅读 · 1 评论 -
grads里批处理的使用
1。数据命名是anth+4位年+2位月+‘.dat’,如anth200012.dat这些数据我原来画都是使用相同的ctl和gs文件,但是在画每一个文件时要每次修改它的输入数据名称和输出的图像名称,具体修的地方见如下ctl和gs文件;********************************************************************ctl文件dse转载 2010-01-25 23:13:00 · 4450 阅读 · 0 评论 -
grads计算散度注意事项及其它
reinitsdfopen e:/data/ncep/data/uwnd.ncsdfopen e:/data/ncep/data/vwnd.ncset gxout fwriteset fwrite e:/data/ncep/data/divmav.grdset lev 1000k=1967while(kset lon 10 240转载 2010-01-25 23:14:00 · 2833 阅读 · 0 评论 -
读取全国160个站资料用于GrADS绘站点图(源程序)
! ######################################################################################! Description: This programe is designed for reading the temperature &! rainfall data from 160 national meteo转载 2010-01-25 23:49:00 · 3327 阅读 · 0 评论 -
maskout()和tmave()做合成分析
maskout()的作用: 逐步剔除不需要的点 tmave()的作用: 使用mask序列, 只求所需时间点的平均 关键点: 将script的脚本数组变成GrADS的变量 范例: 50年的资料序列, 经过初步的分析, 已经分成三组, 分别以1,2,3标识, 并保存在数据文件group.txt中: [1,] 2 [2,] 2 [3,] 3 ...... [18,]转载 2010-01-25 23:07:00 · 2573 阅读 · 0 评论 -
basemap /province-basemap/cnbasemap命令使用
经过测试,终于发现basemap、provincemap-basemap、cnbasemap不能使用的原因,问题主要集中于路径的设置,这是由于程序编写者往往设置的是自己电脑中默认路径!1、basemap使用注意在原始gs文件中添加使用路径,修正使用basemap不能使用的情况2、province-basemap注意将gs文件中23行添加*号,修正调用边界文件而没有的情况3、cnbase转载 2010-01-25 22:12:00 · 3101 阅读 · 0 评论 -
Grads使用点滴
1、多组数据: q file 1----显示文件1的信息;q files------显示所有文件的信息q ctlinfo 1----显示文件1的ctl文件信息set dfile 1---设置当前文件12、数据显示: set gxout printd var----在grads屏幕上显示变量var的值3、数据输出set gxout fwriteset fwri转载 2010-01-25 22:11:00 · 3590 阅读 · 0 评论 -
Grads转成Netcdf
zhanglj 的 Grads转成Netcdf(转动力论坛funny)用GrADS写NetCDF、 HDF-SDS或GRIB等格式文件,请使用lats4d.gs: 说明: 可对5种格式的数据进行转换:GrADS IEEE、GSFC Phoenix、GRIB、NetCDF、HDF-SDS。 输出格式有4种:GRIB、NetCDF、HDF-SDS、flat binary。 选项: -i filenam转载 2010-01-25 22:06:00 · 2975 阅读 · 0 评论 -
Grads中两种相关系数计算
1、显示空间相关系数:计算空间区域的相关系数(空间加权后计算相关系数) d scorr(var1,var2,x=0,x=360,y=-90,y=90)2、计算时间序列var1与空间场var2的相关系数(这里与正常意义的相关系数相似)tcorr(var1,var2,t=1,t=5)例:计算时间序列cc与空间场ps的相关系数空间场,即1维变量cc与空间lon*lat个格点变量ps转载 2010-01-25 22:10:00 · 5354 阅读 · 0 评论 -
province-basemap.gs文件
province-basemap.gs文件(见下面的帖子的附件,下载后拷贝到pcgrads/lib下)和32个省份的边界文件(见下面的帖子的附件)以及1个重新修改了的全国边界文件。实现这个需要的文件除了本省的数据文件外,还要有边界文件和地图文件(各省的地图可从aq9807制作的全国各省地图,见http://www.lasg.ac.cn/cgi-bin/fo ... um=3&topic=4转载 2010-01-25 22:14:00 · 3306 阅读 · 0 评论 -
grads处理多个ctl文件
用grads处理多个相同格式的数据时若单个单个处理非常麻烦,当文件非常多的时候是单个处理是不实际的。下面介绍一种方法; 第一步,在这种情况下可以重新写一个ctl描述文件,其文件变量都和已知的ctl相同,若原来的n文件只是时间不同,那么新描述文件的时间维数是所有原文件的时间的和。同样,若其他维数不同时也用同样的方法处理。 第二步,在第一行之后添加一行: options templa转载 2010-01-25 22:17:00 · 5291 阅读 · 0 评论 -
OpenGrADS气象绘图增强软件
What is OpenGrADS?OpenGrADS是一个致力于改进GrADS(格点分析和显示系统)的界面和扩展性的项目,它是由Brian Doty 及其合作者从20世纪80年代开始发展起来的。GrADS是一个交互式的桌面工具,它用来和容易的进行地球科学数据的接入、操作和可视化呈现。数据的格式可以是二进制、GRIB、 NetCDF、或者 HDF-SDS(Scientific Data S转载 2010-01-25 22:57:00 · 4505 阅读 · 0 评论 -
grads使用笔记
求导——— define dv = cdiff(v,x) define dx = cdiff(lon,x)*3.1416/180 define du = cdiff(u*cos(lat*3.1416/180),y) define dy = cdiff(lat,y)*3.1416/180 display (dv/dx-du/dy)/(6.37e6*cos(lat*3.转载 2010-01-25 23:04:00 · 9822 阅读 · 1 评论 -
画多种形状的点以及台风路径的gs
reinit open E:/eddy-statistic/ssh/test-ssh-2008.ctl t1=29 set t t1 q time utc = subwrd(result,3) enable print E:/eddy-statistic/PRC/2008/2008-ssh-argo-weekly-utcmark.gmf xs转载 2010-01-25 23:06:00 · 3417 阅读 · 0 评论 -
grads批处理ctl文件和用cmd生成idx文件的步骤
一:grads批处理涉及idx文件的grib资料时,编辑grads批处理所需要的ctl文件:(一般可以把已存在的第一个时次的ctl文件进行更改后另存为一个新的ctl文件即可,比如本次批处理的原来的ctl名字都是200206xxxx.ctl形式,则可以新建一个200206.ctl作为grads将来要调用的ctl文件))grib200206%d2%h2 时间循环 opti转载 2010-01-25 23:15:00 · 3850 阅读 · 2 评论 -
GrADS 画断面图例子
reinitopen tslv94.ctl enable print ttver3.gmfset lon 179set lat -79 90set t 3set font 5set vpage 1.5 9.5 4 7set z 1 18set yflip onset xlab offset grid offset grads offset yl转载 2010-01-25 23:44:00 · 3014 阅读 · 0 评论 -
grads中画站点图
1.把不规则的数据按照其经纬读写成stream型data(其fortran程序如下)parameter(n=160)character*8 stid(n),sti(n)real rlat(n),rlon(n),rval(n)!************************************************read rlat rlonopen(9,file=D:/stais转载 2010-01-25 23:50:00 · 3664 阅读 · 0 评论 -
生成GrADS绘图软件站点二进制文件的C语言函数库
由于一些Fortran编译器(ibm xlf90)不支持form=Binary方式,因此经常采用form=unformatted,access=sequential方式写二进制文件。但是该方式生产的二进制文件,GrADS绘图软件经常打不开。因此写了生产GrADS绘图软件站点二进制文件的C语言函数库,因为考虑了Fortran参数传递采用引用传递,所以可以直接在Fortran程序中调用。s转载 2010-01-25 23:53:00 · 1549 阅读 · 0 评论 -
GrADS
GRADS(Grid Analysis and Display System)是当今气象界广泛使用的一种数据处理和显示软件系统。该软件系统通过其集成环境,可以对气象数据进行读取、加工、图形显示和打印输出。它在进行数据处理时,所有数据在GRADS中均被视为纬度、经度、层次和时间的4维场,而数据可以是格点资料,也可以是站点资料;数据格式可以是二进制,也可以是GRIB码,还可以是NetCDF,从而具转载 2010-01-25 21:40:00 · 771 阅读 · 0 评论 -
Grads函数
GrADS的函数分两类, 一类是对格点/站点数据执行运算的,这一类我们姑且称之为分析函数; 另一类是脚本编程(gs)时使用的, 这后一类我们姑且称之为脚本函数. 看本帖的, 可参看下帖关于GrADS变量的内容:http://www.lasg.ac.cn/cgi-bin/forum/topic.cgi?forum=3&topic=6543 第一类分析函数又分为格点分析和站点分析两原创 2010-01-25 22:08:00 · 3652 阅读 · 0 评论 -
GRIB2 资料处理
说明:本人对GRIB数据格式还不是非常熟悉,这里只是初步了解后的一个整理内容,日后还会更新。 Windows环境: 1、保证在CMD窗口中可直接运行wgrib命令,否则在系统环境变量PATH中添加上GRADS的win32e目录的路径,重起电脑。Wgrib程序在GRADS的win32e目录中自带。 2、将grib2ctl.exe放在win32e目录下,以方便在CMD中运行。直转载 2010-01-25 23:09:00 · 8602 阅读 · 0 评论