合成UV风我们用到的是大气探测学的方法,在我之前的文章里面已经介绍过了,这里不做多于赘述。本文主要是提供Fortran程序的实现和matlab实现的对比。
我们看下面这两组数据。
现有南海区域内(115-120˚E,5-10˚N)分辨率为2.5˚×2.5˚的850hPa水平风场U、V逐月资料,时间为2002年5月和6月。这两组数据提供的九个个点2个月的平均风向。
那么风速风向合成的实现就和之前文章里面很类似了。
风向的处理相对于风速会更加复杂,因为uv风向的合成需要判断uv的正负情况来加减arctan(u/v)。这样我们就会频繁用到IF语句,这里我们需要用八个IF判断来判断风的实际方向。原因是要判断4个象限与4个坐标轴。
matlab实现方法:
clear;
U=[-0.7 -1.1 -1.2;-2.9 -2.7 -2.1;-3.4 -3.0 -2.5;4.4 3.5 3.2;5.8 4.5 3.6;6.6 4.9 3.3];
V=[1.6 2.0 1.8;1.4 1.6 1.1;1.3