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向分辨率                                                                              +