ncl调用自定义function

王磊

最近想把以前的程序改写下，所以想到了ncl自带的function功能，做了一个简单的测试程序。开始放置于主程序的路径下，但怕找不到这个函数，就在$NCARG_ROOT/lib/ncarg/nclscripts/csm路径下又放置了一份子程序的代码。但从运行结果看，似乎还是没能找到这个函数，我参考了一下调用contributed的函数调用方式，调用方式应该是没有问题的。所以请教该如何修改程序。

子程序myfunc.ncl的内容为
undef ("foo")
function my_func(x:numeric,opt:logical)
local dimx,rankx,xx
begin

    dimx = dimsizes(x)
    rankx = dimsizes(dimx)
    if(typeof(x).eq."float")then
        if(opt.and.isatt(opt,"scale"))then
            xx = x*opt@scale
        else
            xx = x
        end if
    else
        return(x)
    end if
end


主程序main.ncl的内容为：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/myfunc.ncl"
begin


opt = True
opt@scale = 10
opt@add = 1000
opt@wgts = (/0.25,0.50,0.25/)
;opt@a3d = array_3D


data=(/1.0,2.0,3.0/)
var=data

print(var)
;printVarSummary(opt)
scale_var= my_func(var,opt)    ;调用子程序,但这一行出问题了，为第20行
print(scale_var)
end


错误提示为：

fatal:Illegal right-hand side type for assignment

fatal:["Execute.c":8578]:Execute: Error occurred at or near line 20 in file main.ncl

zongqing123

楼主的问题解决了没有，我也想调用一下自己的函数，但是也出现问题

zongqing123

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"         
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"            
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "./ftest.ncl"

王磊

修改后的代码如下：
（1）建立一个主程序，名为try.main.ncl,内容如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "./try.func.ncl"
begin


opt = True
opt@scale = 10

data=(/1.0,2.0,3.0/)

;printVarSummary(opt)
unpacked_values= my_func(data,opt)
;myfunction is used to multiply  the data by 10(opt@scale)
print(unpacked_values)
end

(2)建立一个子程序，名为try.func.ncl,内容如下：


load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;you don't need to do tihs（load libraries）if your don't need to use the  Built in function of  NCL


undef ("my_func")  ;undefine one function name in one time
undef ("my_func2") ;you don't need to do tihs if your function name is not a Built in function.


;====================================================================================================
;this function canbe used to get the true unpacked value
;(unpacked values) = scale_factor*(packed value)


function my_func( x[*]:numeric , opt :logical)  

;define the parameters you need,in which dimension must be declare
;the declare should include all the varibles(include demension and type) you need to pass
;numeric contains double,flost,long,integer,short,and byte et. al.

local dimx,xx     ;declare the local varibles you want to use(not a must,but recomended)
begin

    if(typeof(x).eq."float")then
        if(opt.and.isatt(opt,"scale"))then
            xx = x*opt@scale
            return(xx)        ;note return is a must,because function procedure need return values
        else
            xx = x
            return(xx)        
        end if
    end if
end

王磊

关键点有3：
其一：实参和虚参的维数必须一致，否则出现：“Illegal right-hand side type for assignment”的问题提示
其二：必须要保证function具备返还值，否则无法调用，因为你的格式是：reture_value=function(a,b)
其三：你所书写的function只需要指定正确的路径即可，无需放在特定的文件夹下。

王磊

此外，关于return values：
(1)you can reture a varible which ony contains value, like return(a),
(2)you can retrun a array, like retrun(a(:)) or retuan(a), where a=a(1:10)   
(3)you can return a varible or arry, which contains some meteinformation ,such as units

此外，你可以在一个文件里面同时定义主程序和子程序，并执行之。

807727227

学到了！！！

luwen

环境变量中加载export NCL_DEF_SCRIPTS_DIR=/your/scripts/path