FVCOM学习：SMS插入水深

写在前面：上一篇文章简单介绍了SMS打网格的过程。在我们打好网格之后，需要将水深数据导入SMS中，将水深数据插值到每个网格中，这样才可以进行下一步的计算。我们可以在保存目录中找到GRD格式文件，首行两个数字分别代表节点数和网格数，我的区域有36621个节点，18779个网格。其他行就是每个点的经纬坐标，因为没有插入水深，可以看到第三列数据都是0。FVCOM中水深为正，因此在SMS中插水深前，也要将水深改为正值。在插入水深之前，我们还需要关注一下开边界。开边界至关重要的，因为我们将从开边界点输入潮汐强迫等数据，是我们的模型不是一团死水，而是与实际海域海况相似。

Ironyz

2022/03/25

首先我们需要找出开边界点，SMS中操作如下。

1.创建一条节点串.点击create nodestrings，再开边界从第一个蓝点开始从左到右（从右向左也可以）选择开边界点（选择后会变红），按住shift可同时选多点。因为我们需要知道开边界点的数量，所以连线时最好数一下开边界到底有多少点。当最后一个蓝点被选择变红后，双击这个红点，就可以创建一条红色的nodestrings了。然后我们点击select nodestrings按钮，我们就可以看到在红色的线中间有一个小黑块，右键它，renumber，点击保存，SMS会在GRD文件中将开边界点排在前面。

2.插入水深

在SMS中，打开水深文件，我的水深文件是txt格式，在这里选择xyz格式，点击下一页，完成。有时会提醒在多少行有‘X’,插入错误，可以将后面的空格删除，即可导入成功。

导入后我们可以看见左边主栏出现scatter data。点击scatter data.界面上面功能栏出现scatter,点击它，在下拉菜单中选择interpolate to mesh.在interpolation中选择插值方式，笔者一般采用第二种inverse distance weighted，读者可以试着尝试其他两种方法比较插值效果，记住勾选Map Z

插值成功后，在左边主栏中Area Propetry中出现”elevation_interp”. ”elevation_interp”其实是我们在上一步操作中默认的New data set name,可以修改。

3.查看水深图

点击小电脑，点击2D mesh，勾选contours，再上面的点击contours，contour method中选择color fill.点击ok就可以看到水深云图了。最后我们保存，查看GRD文件发现第三列有水深数据了。