fvcom 水深文件dep制作
fvcom 水深文件dep制作20250630

本次案例网格和水深展示

vv

image
Figure 1 Model domain

本次制作其它驱动文件的输入文件为yellowsea.2dm

格式2dm; 文件内容格式详细介绍参考：

https://www.xmswiki.com/wiki/SMS:2D_Mesh_Files_*.2dm

制作方法可以参考往期教学：

matlab读取shp文件做SMS的输入文件cst【FVCOM,SCHISM】

【教程制作】从google map提取作为SMS可使用的cst格式的岸线

【matlab教程】matlab画fvcom网格sms网格==仅是展示作用；

模式打网格不只SMS

模式网格制作：SMS===这里最后一步选择保存为2dm格式. 而不是grd格式

【再次感谢烟台海岸带所的学者来稿】制作FVCOM边界潮汐高程的输入文件

上节我们跑通了河流案例=fvcom 实操—河流算例

我们可以打开看看dep文件（tst_dep.dat）的格式：

第一行就是多少个格点数目；然后每个格点经纬度坐标水深

lon lat depth

图片
知道了格式那就好整了：
还是三个代码：
第一个直接调用fvcom 函数包：

PyFVCOM
pip install PyFVCOM

引用：Cazenave, P. W. et al. (2018). PyFVCOM (version x.x.x) [software]. Plymouth, Devon, United Kingdom: Plymouth Marine Laboratory. https://doi.org/10.5281/zenodo.1422462

运行如下命令：

from datetime import datetime
import PyFVCOM as pf
import multiprocessing

def main():

Define a start, end and sampling interval for the tidal data

start = datetime.strptime(‘2016-01-01’, ‘%Y-%m-%d’)
end = datetime.strptime(‘2016-12-31’, ‘%Y-%m-%d’)
interval = 1 / 24 # 1 hourly in units of days
model = pf.preproc.Model(start, end, ‘yellowsea.2dm’, sampling=interval,
native_coordinates=‘spherical’, zone=50)

Write out the files for FVCOM.

model.write_grid(‘yellowsea_grd.dat’, depth_file=‘yellowsea_dep.dat’)
if name == ‘main’:
multiprocessing.freeze_support() # Windows上防止多进程问题
main()

图片
即可得到！

经过验证没毛病和示例一样的：

图片
-------=那就不需要自己在做了有直接使用即可！=

想自己做也OK的，上次更新的函数，里面一改就可以了。因为上次最后一列是0 ； 这次只需要直接给真实数据水深即可。给写个吧。

图片
图片
验证成功：
图片

引用：

引用：Cazenave, P. W. et al. (2018). PyFVCOM (version x.x.x) [software]. Plymouth, Devon, United Kingdom: Plymouth Marine Laboratory. https://doi.org/10.5281/zenodo.1422462

之前处理复杂算例的matlab 代码：
matlab读取fvcom画站点温度剖面随时间变化（Hovmöller 图）
matlab读取fvcom结果并计算密度画出密度剖面分布图
matlab读取fvcom输出结果并画温度垂直剖面图
matlab读取fvcom结果并计算密度画出密度平面分布图
matlab读取fvcom输出结果并画温度水平分布图
matlab读取河流算例nc并画图展示

fvcom 网格文件grd制作

持续更新=20250630------=====