使用conform_dims不足维数

chongzika

函数使用时，总是显示错误，看了半天官网还是不懂怎么使用的，具体维数怎么定义的搞不懂。files=systemfunc("ls ./1995-*.nc")f   =addfiles(files,"r")
E0=f[:]->DSTSFDRY
lat=f[0]->lat
lon=f[0]->lon
nlat = dimsizes(lat)
nlon = dimsizes(lon)

E=new((/96,144/),float)
E1=conform_dims((/12,96,144),E0,(/1,2/))
E=dim_avg_n_Wrap(E1,0)
E!0="lat"
E!1="lon"
printVarSummary(E)

longlivehj

chongzika 发表于 2014-4-28 19:31
可以解决，这个函数读取数据会直接增加一维，不过那个conform_dims函数还是没不明白，总之，感谢你解决了 ...

conform_dims，还有conform，根据官网的解释和例子应该都可以弄懂。继续你贴出的代码，E0为(/12, 96, 144/)的数组，E为(/96, 144/)的数组，现在要求每月与年均值的差值。显然，E0和E直接相减肯定不行，两者大小不一致；可以用循环，循环12次，作12次减法就行；比较符合ncl程序惯例的做法是，将E复制12次，变为和E0结构相同的数组，然后直接相减，这个复制的工作就由conform或者conform_dims完成。

E1 = conform_dims((/12, 96, 144/), E, (/1, 2/))
anom = E0 - E1

五味子

不懂，还是给楼主顶顶

longlivehj

估计是缺少这条语句：ListSetType (f, "join") 。
一些关键变量的信息最好贴出来，比如1995-*.nc是否每个代表一个月？DSTSFDRY的维数信息如何？

chongzika

longlivehj 发表于 2014-4-27 21:18
估计是缺少这条语句：ListSetType (f, "join") 。
一些关键变量的信息最好贴出来，比如1995-*.nc是否每个 ...

数据是这样子的每个都代表一个月，

还有一些量的相关信息，DSTSFDRY,DSTSFWET,DSTSFMBL是我需要的3个量，我想把他们求一下95年的平均

longlivehj

chongzika 发表于 2014-4-28 08:34
数据是这样子的每个都代表一个月，

还有一些量的相关信息，DSTSFDRY,DSTSFWET,DSTSFMBL是我需要的3个 ...

files=systemfunc("ls ./1995-*.nc")
f = addfiles(files,"r")
ListSetType (f, "join")
E0=f[:]->DSTSFDRY
lat=f[0]->lat
lon=f[0]->lon
nlat = dimsizes(lat)
nlon = dimsizes(lon)

E=dim_avg_n_Wrap(E0,0)
E!0="lat"
E!1="lon"
printVarSummary(E)

chongzika

longlivehj 发表于 2014-4-28 10:27
files=systemfunc("ls ./1995-*.nc")
f = addfiles(files,"r")
ListSetType (f, "join")

可以解决，这个函数读取数据会直接增加一维，不过那个conform_dims函数还是没不明白，总之，感谢你解决了我的燃眉之急。