爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5458|回复: 2

求计算台风GPI的MATLAB程序

[复制链接]

新浪微博达人勋

发表于 2018-4-4 14:10:56 | 显示全部楼层 |阅读模式
MATLAB
问题截图:
问题概况: 或者请大神能帮忙解释下这个,根据这个ncl的脚本我需要输入哪些数据,引入的外部函数也必须是ncl文件吗

function cal_vpot(T[*],R[*],PSL[1],SST[1],Plev[*]) ;; temperature(C),mixing ratio(g/kg),Psfc(hPa),SST

begin

; SUBROUTINE PCMIN(SST,PSL, P, T, R,NA,N,PMIN,VMAX,IFL)

; see readme.txt r r rd rd rd i i r r i

; emf::PCMIN(SST,PSL, P, T, R,NA,N,PMIN,VMAX,IFL)



;inputs

if (any(T.ge.100))then

Tc = T-273.15

Tc@units = "degree C"

else

Tc = T

end if

Tk = Tc+273.15

if (.not.(ismissing(SST).or.ismissing(R(0))))then

if (any(SST.ge.100))then

SSTc = SST-273.15

T@units = "degree C"

else

SSTc = SST

end if

else

resu = new(1,"float")

resu@capee = resu

resu@capem = resu

resu@cape_ast = resu

return resu

end if

SSTk = SSTc+273.15



if (max(Plev).ge.10000)then ; if Pa then to hPa

P = Plev*.01

else

P = Plev*1.

end if

if (PSL.ge.5000)then ; if Pa then to hPa

PSL = PSL*.01

end if



if (isMonotonic(P).eq.1)then ;; from bottom to top

Tc = Tc(::-1)

Tk = Tc+273.15

R = R(::-1)

P = P(::-1)

end if



NA = dimsizes(P)

N = NA
我看过提问的智慧: 看过
自己思考时长(天): 30
系统平台:

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2018-4-4 17:02:09 | 显示全部楼层
大哥,这是ncl啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-4 17:14:39 | 显示全部楼层

我知道啊,我只是在求GPI的MATLAB,同时这个ncl是我在网站上找到的希望能够知道其中的意思
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表