编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

牧师

标题

身体

有人知道如何计算垂直集成水分通量在WRF使用NCL。最好是如果你能分享一个NCL脚本。这是我的代码,我用这段代码来计算平均水分通量在一个月内。但我困惑的单位最后的输出。另一个问题是,我不确定我需要输出除以9.8吗?;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;dir1 = " /用户/ gerry /桌面/ WRF_EXTPCP / wrftest / FILES1 = systemfunc (“ls”+ dir1 + wrfout_d01_ *) a1 = addfile (FILES1 +“.nc”、“r”); - - - - - -计算总层水汽通量- - - - - - - - - - - - - - - - - - - - - - - - - - - - psfc = wrf_user_getvar (a1,“psfc”, 1) * 0.01;表面的压力,转向hPa (;单位在WRF PSFC Pa)总统= wrf_user_getvar (a1, "压力",1);气压hPa(时间、级别、纬度、经度)temp = wrf_user_getvar (a1, tk ", 1); K rh = wrf_user_getvar(a1,"rh",-1) ; % ua = wrf_user_getvar(a1,"ua",-1) ; m s-1 va = wrf_user_getvar(a1,"va",-1) ; m s-1 ; printVarSummary(psfc) ; (time,lat,lon) ; printVarSummary(pres) ; (time,level,lat,lon) ; printVarSummary(temp) ; (time,level,lat,lon) ; printVarSummary(rh) ; (time,level,lat,lon) ; printVarSummary(ua) ; (time,level,lat,lon) ; printVarSummary(va) ; (time,level,lat,lon) psfc1 = dim_avg_n_Wrap(psfc,0) pres1 = dim_avg_n_Wrap(pres,0) temp1 = dim_avg_n_Wrap(temp,0) rh1 = dim_avg_n_Wrap(rh,0) ua1 = dim_avg_n_Wrap(ua,0) va1 = dim_avg_n_Wrap(va,0) ;https://www.ncl.ucar.edu/Document/Functions/Built-in/mixhum_ptrh.shtml shum = mixhum_ptrh(pres1,temp1,rh1,2) ;units: (hPA,K,%,kg/kg);calculate the specific humidity(2),(1)means mixing ratio printVarSummary(shum) uq = ua1 * shum ;vapor flux vq = va1 * shum copy_VarCoords(psfc(0,:,:), psfc1) copy_VarCoords(ua(0,:,:,:), uq) copy_VarCoords(va(0,:,:,:), vq) copy_VarCoords(pres(0,:,:,:), pres1) ;;--vertical integrated moisture flux ptop = 100 ; top layer pressure pbot = 1100 uqflux = vibeta(pres1(lat|:,lon|:,lev|:),uq(lat|:,lon|:,lev|:),2,psfc1,pbot,ptop)/9.8 ;should divide 9.8?? vqflux = vibeta(pres1(lat|:,lon|:,lev|:),vq(lat|:,lon|:,lev|:),2,psfc1,pbot,ptop)/9.8 ; 1 means linear interpolation, 2 means log interpolation qflux = sqrt(uqflux^2+vqflux^2)

标签

编辑总结

取消

1

\ begingroup美元 esrl.noaa.gov / psd /数据/复合材料/天/ calculation.html你可以尝试使用,页面上的信息。告诉我们你在哪里卡住了,然后我们可以试着插入。 \ endgroup美元
- - - - - -gansub
2018年12月22日在3分
\ begingroup美元你将如何定义一层?表面一层不是压力。通常是一些两个表面的平均数值。该链接提供给了我一个想法如何计算层 \ endgroup美元
- - - - - -gansub
2018年12月22日,在16
\ begingroup美元是的,你是对的。实际上,我已经搜查了邮件列表,没有发现任何有用的信息给我。 \ endgroup美元
- - - - - -Love_qq_xq
2018年12月22日在4:17
\ begingroup美元我使用了NCL vibeta函数得到整个层的水分通量从1100年到100年。 \ endgroup美元
- - - - - -Love_qq_xq
2018年12月22日在4:19
\ begingroup美元不,不是一个定义的一个数值大气的层。一层可以被定义在许多方面。通常在考虑表面然后上方和下方。取平均的值这三个表面。然后整合跨层或做一个辛普森的求和跨层集成或其他数值积分方案 \ endgroup美元
- - - - - -gansub
2018年12月22日,在43

| 显示3更多的评论

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复