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

lee309

NCL脚本中需要对数据进行分段计算，但是有些数据比如温度是三维数据（经度，维度，高度），直接把分段函数加进去会报错，维度不一致？？？
      if(R.le.-3.7) then
      gama=0.05
      end if
      if(-3.75.lt.R.le.-0.12) then
      gama=0.065+0.004*R
      end if
      if(-0.12.lt.R.le.0.14)  then
      gama=0.109+0.367*R
      end if
      if(R.ge.0.14)           then
      gama=0.155+0.021*R
      end if
上面脚本中R包含经度和维度的，因此报错
fatal:Conditional statements (if and do while) require SCALAR logical values, see all and any functions
这个问题怎么解决？

longlivehj

用where函数

lee309

longlivehj 发表于 2014-12-11 16:09
用where函数

能详细说下吗，对ncl的函数不了解

lee309

怎么去除经纬度信息，直接得到变量值？

longlivehj

看下官网where函数帮助吧，这个是最基本最常用的。

longlivehj

用(/.../)来去掉变量元数据。楼主最好抽空看看ncl 基础教程。

lee309

longlivehj 发表于 2014-12-11 18:36
看下官网where函数帮助吧，这个是最基本最常用的。

你说的where函数相当于C 里面的三目运算符，不是我需要的。
我的问题是判断语句中的量要求是标量，但是我不知道怎么把数据中带有经纬度的变量转化为标量。
比如温度T ：     （35,46）      50       ;    (lat,lon)    温度值
因为if（T>0）中T要求是标量，所以我只需要提取后面的温度值，有没有这样的函数？

longlivehj

lee309 发表于 2014-12-12 21:42
你说的where函数相当于C 里面的三目运算符，不是我需要的。
我的问题是判断语句中的量要求是标量，但是 ...

提取二维数组T里面的值还要函数么？直接用下标不就可以了！
非要用if语句么？你给出的代码，4个where函数就可以搞定！

lee309

longlivehj 发表于 2014-12-12 23:17
提取二维数组T里面的值还要函数么？直接用下标不就可以了！
非要用if语句么？你给出的代码，4个where函 ...

嗯，我尝试用where写了下，发现where只能有两个结果，比如where（a,b,c）结果是b/c
而我有四个if语句，有四段呢，能不能帮忙写下，谢谢

lee309

longlivehj 发表于 2014-12-12 23:17
提取二维数组T里面的值还要函数么？直接用下标不就可以了！
非要用if语句么？你给出的代码，4个where函 ...

请版主大人帮忙写下怎么用where写我的四段函数吧，这个问题怎么都解决不了。谢谢！