NCL脚本中怎么加入分段函数公式

longlivehj

lee309 发表于 2014-12-13 10:04
请版主大人帮忙写下怎么用where写我的四段函数吧，这个问题怎么都解决不了。谢谢！

gama = R
gama = where(R .le. -3.75, 0.05, gama)
gama = where(R .le. -0.12 .and. R .gt. -3.75, 0.065 + 0.004 * R, gama)
gama = where(R .le. 0.14 .and. R .gt. -0.12, 0.109 + 0.367 * R, gama)
gama = where(R .gt. 0.14, 0.155 + 0.021 * R, gama)

lee309

脚本中的R为格点数据，不是一个标量
(41,31)        -0.087815
(41,32)        -0.08835597
(41,33)        -0.08903541
(41,34)        -0.09113897
(41,35)        -0.09109438
(41,36)        -0.093316
(41,37)        -0.09435109
(41,38)        -0.09539092
(41,39)        -0.0966844
(41,40)        -0.09728384
(41,41)        -0.09863391
难道if语句不允许出现这样的数据么，哪位可以指点下？

longlivehj

lee309 发表于 2014-12-13 21:58
脚本中的R为格点数据，不是一个标量
(41,31)        -0.087815
(41,32)        -0.08835597

http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml#If

if(scalar_logical_expression) then
        statement list
end if

lee309

longlivehj 发表于 2014-12-13 21:54
gama = R
gama = where(R .le. -3.75, 0.05, gama)
gama = where(R .le. -0.12 .and. R .gt. -3.75, 0. ...

谢谢版主，已经成功了，以后多多指点