《海洋数值模型的理论及应用》.pptVIP

水平对流项的计算 ■ 虽然模型采用有限差分计算格式，但是对于每个网格 对流项的计算都是按照有限容积的方法进行处理，即 温度对流过程可表示为 ■ 速度的对流过程与温度相似，可写为 其中， 是由z坐标转化为σ坐标后 产生的弯曲项 (28) (29) 垂直扩散项的计算 ■ 垂直扩散项的计算公式（第k层，1kkb-1） ■ 可以改写为 (30) (31) ■ 展开后合并同类项得到 其中 (32) (33) ■ 现在假定温度解的形式为 ■ 由(34)得到的 代入到(32)就可以得到 其中 的值由(33)求得， 的值由上 一层的 值确定 (34) (35) ■ 当k=1时，即表层海水，海水温度主要取决于海表面 温度通量。公式（34）的解可以近似写为 ■ 这样，表层海水的温度只能根据式（31）求解，即 ■ 上式可以进一步表示为 (36) (37) ■ 式（37）还可以进一步写为 ■ 再与温度的通解（34）比较 ■ 可以得到 (34) (38) (39) ■ 当短波辐射通量 的计算公式 ■ r，ad1，ad2 均为光辐射常数，根据不同水质取值如下 (40) ntp 1 2 3 4 5 Jerlov type I Ia Ib II III r 0.58 0.62 0.67 0.70 0.78 ad1 (m) 0.35 0.60 1.0 1.5 1.4 ad2 (m) 23.0 20.0 17.0 14.0 7.9 来源：Jerlov，1976；Paulson and Simpson，1977 ■ 当k=kb-1时，即底层海水，假定海底的热通量为0， 根据前面的推导方式同样可以得到底层海水温度的计 算公式 ■ 对于盐度方程垂向扩散项的计算与温度方程相同，只是不考虑太阳短波辐射这一项 (41) 时间步长的CFL限制条件 ● 由于模型水平方向采用显格式，因此时间步长的选取 必须要满足CFL稳定条件。 ● 对于外模态，时间步长的限制条件为 其中 ， 可能预见到的最大流速 ● 对于内模态，时间步长的限制条件较外模态的情形宽松很多，主要是速度较快的外重力波已经在外模态中考虑了。一般 取值30~50即可 (42) 时间步长的其它限制条件 ● 对于动量或标量还有其它一些时间限制条件 其中 或 ● 以及 其中 分别为地球角速率和地理纬度 (43) (44) 侧开边界条件（bcond） 陆地及岸线是由dum、dvm、fsm控制的，在陆地上这些变量的值设为0，有水的地方设为1 子程序bcond（idx），idx=1对应的水位边界条件；idx=2对应垂向平均流速边界条件；idx=3对应三维水平流速边界条件；idx=4对应温盐边界条件；idx=5对应垂向流速边界条件；idx=6对应湍流动能和湍流混合长度边界条件 模型对开边界条件的要求很高，而开边界条件本身又具有很大的不确定性，因此有必要对模型的外模态和内模态的开边界分别进行处理 数值积分流程图 Formula Boundary Code Inflow condition EAST uaf(im,j) = 2*bc(j)/(h(im,j)+elf(im,j) + h(imm1,j +elf(imm1,j)) elf(im,j) = elf(imm1,j) vaf(im,j) =set WEST uaf(2,j) = 2*bc(j)/(h(1,j)+elf(1,j) + h(2,j)+elf(2,j)) elf(1,j) = elf(2,j) vaf(1,j) = set NORTH vaf(i,jm) = 2*bc(i)/(h(i,jm)+elf(i,jm) + h(i,jmm1) + elf(i,jmm1)) elf(i,jm) = elf(i,jmm1) uaf(i,jm) = set SOUTH vaf(i,2) = 2*bc(i)/(h(i,1)+elf(i,1) +