MATLAB合成风场并生成dfs文件

个人看到文献说再分析风场和模型风场叠加会比较好，就上网查了下，源程序运行和输出有些不满足我需求，就改了下。我也是菜鸟一个，大神略过吧。送给有需要的人吧。

%************程序说明*****************
%由于原来网上下载的程序因为数据结构改变，运行不了，进行了修改，并增加了其他部分。

说明如下：
%（1）改写了风速的数组结构，文件存储格式变更导致，由1440*720*4标成了新格式4*1440*720，将新数组结构改为老结构
%（2）根据新命名格式更改了file文件名称
%（3）部分区域添加了绝对路径，保证运行不出错
%（4）增加了合成风速、风向、气压的输出。
%（5）增加了dfs文件生成所需表头文件，输出结果可直接用于dfs文件生成
%（6）输出结果将时间调整为北京时间，可直接用于数模计算，不用再进行其他转换

直接上代码吧：

 clc;clear all;
%=========通用基本参数======================
pa=1.2;%pa为空气密度
pn=1010*100;%外部气压，千米换算成米
R=6371.004*1000; %地球平均半径
%========================================
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%需修改的参数开始%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fid=fopen('F:\study\test\wind.txt','w');%                                                        +
                                      +
path='F:\study\test\data\';%                                                                     +
%-------------先计算模型风场----------------------------                                          +
%确定模型范围和分辨率                                                                              +
lons=105;%起始经度                                                                                +
lone=123;%结束经度                                                                                +
lats=4;%起始纬度                                                                                  +
late=24;%结束纬度                                                                                 +
dsx=0.125;%x向分辨率                                                                              +
dsy=0.125;%y向分辨率                                                                              +
                                                             +
%输入台风基本参数                                                                                  +
%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&台风路径所需excel格式为：（全部为数值格式）&&&&&&&&&&&&              +
%表头：日期  时间  时间序号 经度  纬度                                               &              +
%  20090806  6     6      123.5  24.8                                              &              +
%时间序号从0开始，后面按照时间加，如时间0,6,12,18,0,6序号为0,6,12,18,24,30            &              +
%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&              +
input=xlsread('F:\study\test\lujing.xlsx');%修改路径