10
\ begingroup美元

我正在学习Wavewatch III波模型(http://polar.ncep.noaa.gov/waves/index2.shtml)。我困惑的稳定校正因子,名叫ASF的代码。它的目的是什么?与物理有关,还是与数值方法和CFL条件?我不能立即算出来的代码。

\ endgroup美元
2
  • 1
    \ begingroup美元 你能点中的特定部分网站或代码,你是指? \ endgroup美元
    - - - - - -等密度线振荡
    2015年1月3日在5:56
  • 1
    \ begingroup美元 @IsopycnalOscillation没有在网站上被讨论,但我将试着拼凑代码段使用ASF,然后回来报告。 \ endgroup美元
    - - - - - -boxofchalk1
    2015年1月3日,6

1回答1

12
+ 100
\ begingroup美元

校正系数的稳定ASF有关大气稳定性的影响(浮力和剪切功能)波的增长,并实现在Wavewatch3杜尔曼Chalikov 1996输入源项。修正的代码可以在发生w3updtmd.ftn:

!5。稳定校正(/ STAB2) !初始设置:!!SHSTAB = 1.4 !OFSTAB = -0.01 !CCNG = -0.1 !ccp = 0.1 !FFNG = -150。 ! FFPS = 150. ! !/STAB2 STAB0 = ZWIND * GRAV / 273. ! !/STAB2 DO ISEA=1, NSEA !/STAB2 STAB = STAB0 * AS(ISEA) / MAX(5.,U10(ISEA))**2 !/STAB2 STAB = MAX ( -1. , MIN ( 1. , STAB ) ) ! !/STAB2 THARG1 = MAX ( 0. , FFNG*(STAB-OFSTAB)) !/STAB2 THARG2 = MAX ( 0. , FFPS*(STAB-OFSTAB)) !/STAB2 COR1 = CCNG * TANH(THARG1) !/STAB2 COR2 = CCPS * TANH(THARG2) ! !/STAB2 ASF(ISEA) = SQRT ( (1.+COR1+COR2)/SHSTAB ) !/STAB2 U10(ISEA) = U10(ISEA) / ASF(ISEA) !/STAB2 END DO

修正基本上占大气稳定基于海气温差、应用和计算风速在十米级高度“有效”,从输入风速。注意,这段代码只在很活跃/ STAB2启用开关(可选),只有工作/ ST2开关(杜尔曼和Chalikov 1996源项包)。它的实现和测试结果中描述杜尔曼2002。

这也是记录在手册的版本4.18中,见方程(2.86)-(2.89)在部分2.3.7。

引用:

  • 杜尔曼,h·l·d·诉Chalikov 1996:第三代风浪模型中的源项。期刊。Oceanogr。,26, 2,497–2,518.

  • 杜尔曼,h·L。,2002: Validation of WAVEWATCH III version 1.15 for a global domain. Tech. Note 213, NOAA/NWS/NCEP/OMB, 33 pp.

\ endgroup美元
3
  • 1
    \ begingroup美元 感谢@milancurcic详细的回复和引用!是的,我检查了,我确实使用STAB2ST2开关——他们都包含在默认开关文件。从手动,似乎ASF是蒙混因素提高波增长过程的准确性。为了确保(我没有发现本手册中提到的),我仍然应该输入10米以上的风速表面,而不是在其他一些高度,对吧?我认为U10就是WW3接受,因为它提到“请注意,这有效风速是派生的风在10米高。”手册。 \ endgroup美元
    - - - - - -boxofchalk1
    2015年1月3日,在二十29
  • 3
    \ begingroup美元 很高兴有专家在特定的模型这里谁能回答这些问题:-) \ endgroup美元
    - - - - - -半日西蒙
    2015年1月4日,在43
  • 1
    \ begingroup美元 @boxofchalk1是正确的,风输入字段保持在10米高。 \ endgroup美元
    - - - - - -milancurcic
    2015年1月4日22:04

你的答案

通过点击“发布你的答案”,你同意我们服务条款,隐私政策饼干的政策

不是你要找的答案?浏览其他问题标记问你自己的问题