fvcom 科氏力文件cor制作==基础文件制作
fvcom 科氏力文件cor制作基础文件制作
fvcom 科氏力文件cor制作基础文件制作20250705=
学习阶段慢慢来！一个一个弄懂了之后，可以代码一键制作。都写成函数调用即可。
目前已经制作的文件有：
fvcom 垂直坐标sigma制作
fvcom 水深文件dep制作
fvcom 网格文件grd制作
我们的目标：
图片
在 FVCOM（Finite Volume Community Ocean Model）中，cor 文件用于提供科里奥利参数（Coriolis parameter），即地转偏向力的计算基础。这个参数与纬度有关，是动力方程中不可或缺的一部分。
制作最重要的是理解它。
科氏力计算公式：
图片
而对于fvcom模式来说，我们只需要提供纬度。
那就简单了：
文件格式就变成了
lon lat lat
图片
制作方法：
还是三个代码：
第一个直接调用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_coriolis(‘yellowsea_cor.dat’)
if name == ‘main’:
multiprocessing.freeze_support() # Windows上防止多进程问题
main()
即可得到！
经过验证没毛病和示例一样的：
-------=那就不需要自己在做了有直接使用即可！=
想自己做也OK的，上次更新的函数，里面一改就可以了。因为上次最后一列是0 ； 这次只需要直接给真实数据水深即可。给写个吧。
图片
图片
图片
图片
学习阶段慢慢来！一个一个弄懂了之后，可以代码一键制作。都写成函数调用即可。
clear;clc;close all;
% use function ===
% 海洋与大气科学： 20250705
convert_2dm_to_cor(‘yellowsea.2dm’, ‘yellowsea_cor.dat’)
本次案例网格和水深展示
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边界潮汐高程的输入文件
引用：
引用：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制作
持续更新=20250705------=====