1. 下載數(shù)據(jù)和代碼

下載TPXO9數(shù)據(jù)：https://www.tpxo.net/tpxo-products-and-registration

下載潮汐預報代碼 TMD（基于matlab），或者用OTPS做預報（基于fortran & Linux）：https://www.tpxo.net/otps

如果用TMD做預報，下載TPXO9數(shù)據(jù)時要選擇binary格式

下載?tpxo_atlas2local.m 代碼，用于將TPXO9生成的各分潮文件合并，合并的文件能用于TMD預測

2. 用 tpxo_atlas2local.m 處理TPXO9數(shù)據(jù)

該函數(shù)需要調(diào)用TMD里的程序，文件路徑為：你的路徑/TMD/FUNCTIONS/。所以，可以把?tpxo_atlas2local.m 放在該路徑下，或者在?tpxo_atlas2local.m 里添加TMD函數(shù)的路徑：

????????addpath('D:\Code\TMD\FUNCTIONS')

tpxo_atlas2local.m 的調(diào)用方式：

????????function []=tpxo_atlas2local(atlas_modfile,out_modfile,lat_lims,lon_lims);

atlas_modfile指向TPXO9數(shù)據(jù)的路徑，例如，我的數(shù)據(jù)放在 ‘D:\Code\TMD\DATA\TPXO9_atlas_v5’ 文件夾下，為方便起見，我在‘D:\Code\TMD\DATA’ 文件夾下保存atlas_modfile文件（放哪都行）。創(chuàng)建atlas_modfile文件：新建一個txt文檔，輸入三行：

????????D:\Code\TMD\DATA\TPXO9_atlas_v5\h_*_tpxo9_atlas_30_v5*

????????D:\Code\TMD\DATA\TPXO9_atlas_v5\u_*_tpxo9_atlas_30_v5*

????????D:\Code\TMD\DATA\TPXO9_atlas_v5\grid_tpxo9_atlas_30_v5

h是潮位文件，u是流速文件，grid是模型網(wǎng)格，*可以歷遍所有相應文件。TPXO9有不同版本，最后的部分需要根據(jù)你的文件名修改，比如我的文件名帶有v5。保存txt文件，然后去掉文件后綴，我保存的文件名是下圖的 “Model_tpxo9_atlas"

接下來編輯out_modfile文件，指明tpxo_atlas2local.m生成的文件需要放在哪個路徑下，這里我簡單采用了函數(shù)提供的默認寫法。首先，新建txt文件，寫入以下三行：

????????D:\Code\TMD\DATA\MyArea\h_MyArea

????????D:\Code\TMD\DATA\MyArea\uv_MyArea

????????D:\Code\TMD\DATA\MyArea\grid_MyArea

它表示tpxo_atlas2local.m生成的文件將儲存在 ‘D:\Code\TMD\DATA\MyArea’文件夾下，生成的三個文件分別是：h_MyArea，uv_MyArea，grid_MyArea。我把txt文件命名為Jiangsu_coast，并去掉后綴名。注意：你需要先建好MyArea文件夾

注意：上述提到的文件路徑，盡量不要有中文，也避免空格、橫線。我曾使用以下的文件路徑，導致運行失?。?/p>

????????D:My - code\Code\TMD\DATA\MyArea

完成上述準備工作后，可以在命令行窗口運行函數(shù)：

????????tpxo_atlas2local('DATA\Model_tpxo9_atlas','DATA\Jiangsu_coast',[32.4,35],[120,124]);

最后兩個參數(shù)是你需要的經(jīng)緯度

3. 運行TMD進行預報

直接運行TMD.m文件，會提示你選擇輸入模型，這里我直接點擊 ”Jiangsu_coast"（見上圖），函數(shù)會自動定位到‘D:\Code\TMD\DATA\MyArea\’文件夾并讀取上一步生成的h_MyArea，uv_MyArea，grid_MyArea文件。運行函數(shù)后得到GUI界面：

這里我把所有分潮都點上，左下角可以選擇predict tide（做預報）或者extract tidal constants（提取分潮常數(shù)）。你可以用鼠標在界面上點擊選擇想要研究的站位（不推薦），我會提前生成一個文件，從文件里讀取經(jīng)緯度坐標和時間（見下一小節(jié)）：

注意：全球潮模型用的是格林威治時間，如果你要預報北京時間（GMT+8），需要算好時差。比如你想預報北京時間 2018/10/15 14:00:00某個站位的潮位，起始時間應選為2018/10/15 06:00:00

上圖中的文件夾“D:\Code\TMD\LAT_LON”里有TMD自帶的坐標文件，你可以隨便選一個玩玩

4. 創(chuàng)建坐標文件，用于TMD預報

下載OTPS或者OTPSnc：https://www.tpxo.net/otps

解壓后在文件路徑 ‘D:\OTPSnc\matlab’ 下可以找到mk_lltime.m。函數(shù)內(nèi)容如下：

lat和lon是匹配的經(jīng)緯度，上圖顯示我有5個站點，其中一個的坐標是（34.4136，120.745）。dt是時間間隔，單位是小時。d1\d2分別是起始時間和結(jié)束時間，注意用的是格林威治時間。以上就是你需要根據(jù)自身情況修改的內(nèi)容，直接運行即可得到五個文件，每個文件對應一個站點：

文件里長這樣（左邊的行序號是軟件顯示的，文件里沒有）：

生成的這些文件放在“D:\Code\TMD\LAT_LON”里，可在TMD的GUI頁面里直接調(diào)用。

5. 用OTPS或OTPSnc做預報

這兩個和TMD本質(zhì)上是一樣的，只不過這些是基于Fortran和Linux，OTPSnc會生成NetCDF格式的文件，需要 netcdf 和 netcdff。如果你會使用LInux系統(tǒng)，但不懂MATLAB，可以考慮用這些軟件做預報。由于本人沒成功運行OTPSnc，就不過多介紹了，主要是一直報錯說我找不到 lnetcdf 或 lnetcdff。原因是這樣的，我用的是學校的HPC系統(tǒng)，它已經(jīng)安裝了Netcdf和netcdf-fortran，lnetcdf.so在netcdf文件夾下，lnetcdff.so在nefcdf-fortran文件夾下，但我不懂怎么在makefile里同時添加兩個文件夾：

如果你不需要用到NetCDF文件，那么可以用OTPS做預報！