ROMS简单应用——绘制区域温度图

以jersey_lake的例子展示如何用roms输出的数据绘制一个给定区域的温度分布图。

 

1、配置输入文件

这一部分可以参考之前写的博客中关于变量设置的部分

Linux系统下ROMS安装流程与例程运行_SStegosaurus的博客-CSDN博客

如果之前已经完成了配置，这里就不需要修改my_build_paths.sh

2、运行ROMS

运行以下两个脚本之一完成安装
./build_roms.sh
./build_roms.sh -j 4 

修改VARNAME等参数，并运行mpirun指令，得到lake_jersey_avg、lake_jersey_his、lake_jersey_rst三个nc文件。

3、matlab绘图

作为刚接触ROMS的初学者，即使已经获得了输出文件，在提取数据和绘图方面也会遇到一些困难。用以下链接提供的绘图工具包，可以更方便地完成绘图。

GitHub - johnwilkin/roms_wilkin

ROMS Matlab tools TiddlyWiki - John Wilkin's notebook describing some Matlab tools for viewing and manipulating ROMS outputs and inputs

 代码：

filename='D:/tmp/lake_jersey_avg.nc';
grdfile='lake_jersey_grd_a.nc';
%  nc_inq(filename)
% nc_dump(filename)
info = ncinfo(filename);		%读取nc文件信息


nc_var = ncread(filename,'temp');
ocean_time=ncread(filename,'ocean_time');
grd = roms_get_grid(filename,filename);
s=.1;

[Data,Grid,Han]=roms_zview(filename,'temp',2,-12,grd,3,s,'k');