- 积分
- 1096
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-3-18
- 最后登录
- 1970-1-1
|
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
|
|
窥视卡
雷达卡
发表于 2018-4-4 14:10:56
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
成长值: 32430
发表于 2018-4-4 17:02:09
楼主