ncl中要计算正滞后和负滞后，必须进行两次计算。
mxlag = 9
x_Lead_y = esccr(x,y,mxlag)
y_Lead_x = esccr(y,x,mxlag) ; switch the order of the series
ccr = new ( 2*mxlag+1, float)
ccr(0:mxlag-1) = y_Lead_x(1:mxlag:-1) ; "negative lag", -1 reverses order
ccr(mxlag:) = x_Lead_y(0:mxlag) ; "positive lag"
❓q1 谁超前谁？
esccr(x,y,mxlag)
计算序列x最右边维及序列y最右边维之间的交叉相关系数。mxlag表示x超前y的时长，建议 0 <= mxlag <= N/4，N为数组x或数组y最右边维的序列长度。
eg. 给定数组x(M,N)和数组y(K,N)，计算x与y的交叉相关系数。
cor=esccr(x,y,10)；返回值为cor(M,N,0:10)，
其最左边11个数表示x超前y从0至10个时次的交叉相关系数。