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岁,2497 - 2518。

  • 杜尔曼,h·L。第三,2002:验证WAVEWATCH 1.15版本为全球域。213年科技。注意,NOAA / NWS / NCEP / OMB, 33页。

\ 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

你的答案

通过点击“发布你的答案”,你同意我们服务条款并承认您已阅读并理解我们的隐私政策的行为准则

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