利用中央差分函数计算变量时出现的数组合并的问题

独一吖

现有air（time,level,lat,lon）,具体維数（120，12，144，73），需要求air对x的偏导，由于x随纬度变化，又不想用四重循环跑，因此想用NCL自带的中央差分函数center_finite_diff_n,通过纬度的循环变化，求出90N、87.5N……-90S下的偏导值的三维数组，并合并。

问题描述：在运行过程中，报错fatal:Assignment type mismatch,right hand side can't be coerced to type of left hand side.
分析：代码中运用center_finite_diff_n时，temp是不可以提前new的，但是不同纬度值下三維数组temp合并維到四维Q1_21时，temp是需要提前new或者明确具体維数的，因此无论是否预先创建temp，都会报错，也无法继续运行。

求教各位大佬，在沿用该思路的情况下，如何解决？

独一吖

或者说，对于由temp=center_finite_diff_n（……）得出的数组temp（不可以预先new），我如何声明其維数呢？

独一吖

好尴尬，分析的错的。我知道问题出在对Q1_21的创建上了，因为用的EC的数据，因该是double，因为之前跑NCEP是float，应该直接typeof(air)就好啦。不知道怎么删贴，大家见笑了。

圆梓君

哇，正准备求中央差分对x和y方向的函数