grads如何运行gs文件

运行gs文件，所安装的程序没有直接关系，用以下软件也可以直接打开。

可打开GS文件的软件: Verity Software House GemStone, Geosoft Oasis montaj.

gs格式：剧本写的谷歌Apps脚本，基于JavaScript和用于自动化整个谷歌产品的任务，在网上一门编程语言，用于创建在线Web应用程序，如Web界面，发送自定义的电子邮件程序或一个程序，自动填充和格式谷歌文档文件。

高盛的文件进行编辑与谷歌Apps脚本，运行在Web浏览器中的代码编辑器。

<H1 >其他GS格式： </H1 >与GemStone的创建文档，流式细胞仪基于概率状态建模系统上的数据分析程序，包含由流动细胞仪捕获的数据模型和分析，可以在宝石中以图形方式浏览应用程序。

有关GemStone的更多信息，请访问宝石的网站拆闹灶。

<H1 >其他GS格式： </H1 >用绿洲MONTAJ ，地理空间数据处旅扮理和制图应用程序的脚本文件，包含了一组类似的宏观指导用于记录一组重复的任务，并在软件自动化他们。

可以运行从GX →脚本→运弯银行脚本的脚本..菜单选项。在出现的对话框允许用户提供可选参数的脚本。

给你一个小程序，是用perl写的，然后funny转成了exe文件，你可以用它生成ctl，但生成的ctl文件还需要自己去掉这个程序强制添加上去的一些信息，然后，你用gribmap.exe生成index文件，就可以显示了。

这个zip包里有原始的perl程序，转好的exe文件和 gribmap.exe三个文件。

C:\drawing\ncep-monthly\ex>grib2ctl -i prs.grib.mean.y1980>y1980.ctl Using NCEP reanalysis table, see -ncep_opn, -ncep_rean options Using NCEP reanalysis table, see -ncep_opn, -ncep_rean options

C:\drawing\ncep-monthly\ex>gribmap -i y1980.ctl Open Error: Unknown keyword in description file -->The invalid description file record is:

-->this exe file was created with the evaluation version of perl2exe.

The data file was not opened. File name is: y1980.ctl

"Using NCEP reanalysis table, see -ncep_opn, -ncep_rean options" 就是说你应该用 -ncep_opn 或 -ncep_rean 的选项，具体看帮助

下面就不用说了，ctl都没有形成，自然不行的了

Hi,funny Thanks！

错误与“ -ncep_opn, -ncep_rean options”无关，是ctl中endvars后“-->this exe file was created with the evaluation version of perl2exe.”的这句话作怪，删掉后可正常得到idx。

＝＝＝＝＝＝＝＝＝

D:\PCGrADS\win32>grib2ctl E:\。。。\air2m.mon.mean.nc >air2m.ctl 显示的信息如下： Big problem:

#either e:\ncep\nc\air2m.mon.mean.nc is missing or not a grib file or wgrib is not on your path

可这个人家告诉我确实是格点文件啊银汪 ＝＝＝山禅＝＝＝＝＝＝＝＝＝＝＝＝＝

这个文件是格点文件，但不是grib码文件，直接用sdfopen 文件名就可以了。

#我在DOS里面用命令C:\Program Files\PCGrADS\win32>grib2ctl.exe grib_file 为什么不能生成grib_file.ctl?请各位前辈指点。 [这个贴子最后由tibet在 2004/06/30 04:13pm 编辑]

C:\Program Files\PCGrADS\win32>grib2ctl grib_file >grib_file.ctl or

gs>!grib2ctl grib_file >grib_file.ctl

#请教个问题，如何形成几个时次在一个文件中的ctl 和idex

请参考grads中文手册。 ctl文件中可以设置一个参数 options template。生成一个时次的ctl之后自己编辑修改ctl。不然，把所有时次都写一个grib文件中。

'reinit'

'sdfopen d:\zhj\slp.1948.nc' 'set lon 120 240' 'set lat 20 60' 'set lev 0' 'set t 289 366' 'set gxout fwrite'

'set fwrite d:\zhj\slp.grd' 'd slp'

'disable fwrite'

锋唯仔写成.grd文件再写一个.ctl文件就可以了 不过是不是有点麻烦

#各位前辈

我最近才接触到grads ,现在用的是.nc格式的ncep reanalysis资料 看了前面的贴子,还是有很多不太明白的,现在我用sdfopen可以打开资料 也可以画图,但是怎么能生成.ctl文件呢.希望能说的清楚一些 非常感谢

你这个是格点资料 直接用sdfopen读就可以了，nc文件自带描述文件的，不需要你编ctl。要是你把nc

资料转化成grd格式的再画图就需要编写ctl了

￥

D:\PCGrADS\win32e>gribmap -i hhca50.ctl Open Error: Data file type invalid

-->The invalid description file record is: -->options yrev

The data file was not opened. File name is: hhca50.ctl

我用girbmap命令生成idx时候，会提示这样的错误，是怎么回事啊？

1　GRADS有关的文件类型和维数环境

可以认为和GRADS有关的文件类型有：十进制的原始数据文件（＊.TXT）、二进制的数据文件（＊.DAT）、数据描述文件（＊.CTL）、批处理文件（＊.GS）和图形文件李戚（＊.GMF）。

气象业务中使用的站点资料或格点资料都是以十进制形式存放，而GRADS只能识别二进制的数据格式，所以使用GRADS之前的第一步就是转换数据文件；数据描述文件则是对数据文件进行说明，以便后续的 *** 作有对象；批处理文件是把进入GRADS绘图环境后所要输入的命令写成批处理格式，以便可以自动执行输入的各项 *** 作命令（批处理文件可单独建立，也可以略过不写，而在进入GRADS环境后一步一步地输入各项 *** 作命令）；图形文件是GRADS已经绘制好的图形，只能用GV打开浏览。

数据文件一般都是用Power Station或Visual Fortran来转换（也可以用C语言或其它工具）。数据描述文件、批处理文件可以在“写字板"中写好，只是在存档的时候，要把文件的后缀改为.CTL和.GS。

GRADS将每一个物理量场视为1个四维数据集，它包括空间三维和时间一维。维数环境的定义由SET LAT（纬度）/哪御陵LON（经度）/LEV（高度）/TIME（时间）来设置，也可由SET X/Y/Z/T来设置。

2　GRADS的使用方法

文中采用1982年1月～1994年12月共156个月气象拆扒场中的1 000 hPa的海平面气压、温度以及850 hPa、500 hPa、200 hPa的u、v风场的月平均数据和45°×15°的格点资料，来举例说明GRADS的具体使用方法（经度为40°E～150°E，纬度为10°S～25°N）。

2.1　数据文件的转换

格点资料用Power Station转换，程序如下：

PARAMETER(NX=45,NY=15)

REAL GRID(45,15)[格点资料数组大小]

CHARACTER＊60 AA(156)　[定义数组前的说明为字符串]

CHARACTER＊20 DNAME(8)　[资料文件的数量]

DATADNAME/'D:.DAT','D:.DAT',

＄'D:.DAT','D:.DAT',

＄'D:.DAT','D:.DAT',

＄'D:.DAT','D:.DAT'/

[须描述资料文件的列举]

DO 10 N=1,8</FONT>

OPEN(10＋N,FILE=DNAME(N))

10 CONTINUE

[依次打开各个资料文件]

OPEN(20,FILE='D:.GRD',FORM='BINARY',

＄ACCESS='DIRECT',RECL=NX＊NY＊4)

[把所有须描述的资料文件数据同时放入1个文件ALLDAT中，以便只须1次转换后便可以随意调取]

INUMBER=0

DO 100K=1,156</FONT>

DO 20 N=1,8</FONT>

READ(10＋N，50)AA(K)

READ(10＋N,40)((GRID(I,J),I=1,NX),J=1,NY)

INUMBER=INUMBER＋1

WRITE(20,REC=INUMBER)((GRID(I,J),I=1,NX),J=1,NY)

20 CONTINUE

100 CONTINUE

[本程序按经纬度顺序先读写第一时刻的SLP、SST、U850、U500、U200、V850、V500、V200，再读写第二时刻的资料，依次类推...]

40 FORMAT(15F4.0)

50 FORMAT(A60)

END

2.2　数据描述文件的编写

数据描述文件DAT.CTL的编写格式为：

DSETD:.GRD[DSET是给出所描述文件的文件名]

TITLEWEATHERDATA[该数据描述文件的标题]

UNDEF－9.99E33[缺省记录的标记]

XDEF45 LINEAR40 2.5[X方向共45个格点，起始纬度为40°，步长为2.5°]

YDEF15 LINEAR－10.0 2.5[Y方向共15个格点，起始经度为－10°，步长为2.5°]

ZDEF4 LEVELS1000 850 500 200[Z方向分为4层，即1 000、850、500、200 hPa]

TDEF 156 LINEAR JAN1982 1MO[时间共156个月，起始时间为1982年1月，步长为1个月]

VARS4 [共SLP、T、U、V 4个变量]

SLP 0 0 [SEALEVE LPRESSURE]

T 0 0 [TEMPS]

U 3 0 [UWINDS]

V 3 0 [VWINDS]

ENDVARS [结束变量说明标志]

2.3　批处理文件的编写

如果我们欲编制批处理文件时，则应注意在用写字板编写GS文件时，必须在各项命令两边加单引号。如绘制经度为40°E～150°E，纬度为－10°S～25°N，时间为1982年10月（第10时刻）的850 hPa等风速矢量图时，其批处理文件DAT.GS可以写成以下语句：

‘OPEN D：．CTL’　

‘SET LON 40 150’

‘SET LAT －10 25’　

‘SET TIME T 10’　

‘SET LEV 850’　

‘ENABLE PRINT D；．GMF’　

‘SET GXOUT VECTOR’　

‘SET CCOLOR 5’　

‘DISPLAY U；V’　

‘PRINT’

‘DRAW TITLE 10/1982 WINDS’　

‘QUIT’　

以上各项命令都可以根据自己的需要进行更改。比如：想要绘制经度为40°E～150°E，纬度为－10°S～25°N，时间为1984年6月的1 000 hPa气温场和温度场的合成流线图，就可以把以上的批处理文件改为：‘OPEN D：.CTL’，‘SET LON40 150’，‘SET LAT－10 25’，‘SET T 30’，‘SET LEV 1000’，‘ENABLE PRINT D；.GMF’，‘SET GXOUT CONTOUR’，‘SET CCOLOR 5’，‘SET CSTYLE 3’，‘DISPLAY T’，‘SET CSTYLE 1’，‘DISPLAY SLP’，‘PRINT’，‘DRAW TITLE 6/1984 WINDS’，‘QUIT’。

3　结语

当然，GRADS可以绘制的图形远不止这些，若固定经度、时间，让纬度、高度变化，则可以画出气象场的剖面图；若固定经、纬度和高度，让时间变化，则可以画出某气象要素场随时间的变化。另外，GRADS还可以通过调用函数名来直接绘制气象图形，如当设置显示命令“DISPLAY HCURL（U，V）”和“DISPLAY HDIVG（U，V）”时，就可以直接得出垂直涡度和水平散度图等。

---