查看ww3_case_01运行结果
写在前面
话接上文(WAVEWATCHⅢ6.7.0安装记录),按照上文操作,最终会生成以网格形式给出的计算结果,计算结果文件为GrADS格式的,可以用GrADS软件打开,所以这一篇就记录一下用GrADS软件查看计算结果的过程。
1. ww3_case_01运行结果
运行结果在以下文件夹中(其实还有很多中间结果被删除了,后面会写一篇记录一下案例运行过程的分析)
/home/zxs/WW3/cases/work_case_01
# 切换到该目录查看目录中的文件
cd /home/zxs/WW3/cases/work_case_01
ls
# 其中表示网格输出结果的文件为以下这些
grd1.ctl ww3.grd1
grd2.ctl ww3.grd2
grd3.ctl ww3.grd3
grd4.ctl ww3.grd4
grd5.ctl ww3.grd5
# 左边一列为输出结果的描述文件,右边一列为数据
# 以上文件直接打开会乱码,应该下载GrADS软件然后在软件中打开
2. 安装GrADS
对于该软件的介绍就免了,如果是由于我们的系统的ubuntu所以安装十分简单,只需要打开应用商店,搜索GrADS,然后等待安装完成即可:
安装完成后在终端输入以下命令查看GrADS是否正确安装:
grads
# 如果出现以下结果,表明安装成功,输入yes则直接打开软件
Grid Analysis and Display System (GrADS) Version 2.2.1
Copyright (C) 1988-2018 by George Mason University
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information
Config: v2.2.1 little-endian readline grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile
Issue 'q config' and 'q gxconfig' commands for more detailed configuration information
Landscape mode? ('n' for portrait): yes
# 随后会出现绘图窗口,不过此时窗口内一片漆黑,在终端输入以下命令退出软件
quit
3. 利用GrADS查看计算结果
GrADS的相关指令及详细教程可以查看官方用户手册:GrADS官网,这里仅仅使用最基本的查看和绘图指令。
# 切换到结算结果所在文件,并打开GrADS
cd /home/zxs/WW3/cases/work_case_01
grads
yes
# 输入以下指令打开grid1网格的计算结果
open grd1.ctl
# 得到以下输出
Scanning description file: grd1.ctl
Data file ww3.grd1 is open as file 1
LON set to -85 0
LAT set to 0 70
LEV set to 1000 1000
Time values set: 2006:1:1:0 2006:1:1:0
E set to 1 1
# 输入以下指令查看输出数据包含的信息
q file
# 得到以下输出
File 1 : WAVEWATCH III gridded data
Descriptor: grd1.ctl
Binary: ww3.grd1
Type = Gridded
Xsize = 86 Ysize = 71 Zsize = 1 Tsize = 25 Esize = 1
Number of Variables = 8
map 0 99 grid use map
dpt 0 99 Water depth
uwnd 0 99 Wind U (m/s)
vwnd 0 99 Wind V (m/s)
ice 0 99 Ice Conc. (-)
hs 0 99 Wave height (m)
fp 0 99 Peak Freq. (Hz)
dp 0 99 Peak Dir. (rad)
# X刻度86个,Y刻度71个,数据维度均为1,时间步25个,E维度为1
# 一共包含8个变量,我们比较关心的是风速(uwnd、vwnd)和波高(hs)
# 输入以下命令设置图形为填色图而非等高线图
set gxout shaded
# 输入以下命令设置查看第25个时间步的数据
set t 25
# 输入以下命令设置查看的变量为波高
d hs
# 输入以下命令将图形输出为svg矢量图
gxprint grid1hs.svg svg
# 由于优快云 Markdown不好插入矢量图,所以以下为截图
4. 生成其他数据图片的命令
- grid5波高
open grd5.ctl
q file
set gxout shaded
set t 25
d hs
gxprint grid5hs.svg svg
5. GrADS常用指令及解释
# 打开file文件,如果file文件不在当前目录
# 需要切换到当前目录,或者使用全部路径
open file
# 查看file文件中包含的信息
q file
# 设置绘图模式为填充图
set gxout shaded
# 设置绘图数据的时间维度为25
set t 25
# 绘制hs变量的图形
d hs
# 将当前图形保存为名为grid5hs.svg的svg矢量图
gxprint grid5hs.svg svg
# 关闭之前的绘图以进行新的绘图
reinit
6. 在Windows下安装GrADS
- 从以下网址下载安装包
- 点击安装,一定要选择默认位置,否则运行过程可能会出现如下报错
Open Error on Map Data Set: /usr/local/lib/grads/lowres
Error opening stroke character data set
Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
-
安装过程中所有可以打钩的地方全部勾上,免得自己设置环境变量
-
找到桌面上的三个图标,打开第三个所在的位置
-
记录下opengrads.exe可执行文件的全部路径
C:\OpenGrADS-2.2\Contents\Cygwin\Versions\2.2.1.oga.1\i686\opengrads.exe
- 在终端切换至需要打开GrADS文件的文件夹下,执行第五步的可执行文件从而进入GrADS,然后正常使用GrADS