NCL 为什么数据维度都是一样的，我加起来取平均就显示数据维度左右不匹配呢

ststst

f111        =addfile("/mnt/d/test/pr_climate_BCC-CSM2-MR_historical_r1i1p1f1_gr.nc","r")
f222        =addfile("/mnt/d/test/pr_climate_CanESM5_historical_r1i1p1f1_gr.nc","r")
f333        =addfile("/mnt/d/test/pr_climate_CESM2_historical_r1i1p1f1_gr.nc","r")
f444        =addfile("/mnt/d/test/pr_climate_CNRM-CM6-1_historical_r1i1p1f1_gr.nc","r")
f555        =addfile("/mnt/d/test/pr_climate_CNRM-ESM2-1_historical_r1i1p1f1_gr.nc","r")
f666        =addfile("/mnt/d/test/pr_climate_FGOALS-g3_historical_r3i1p1f1_gr.nc","r")
f777        =addfile("/mnt/d/test/pr_climate_HadGEM3-GC31-LL_historical_r1i1p1f3_gr.nc","r")
f888        =addfile("/mnt/d/test/pr_climate_MPI-ESM1-2-LR_historical_r1i1p1f3_gr.nc","r")
f999        =addfile("/mnt/d/test/pr_climate_UKESM1-0-LL_historical_r1i1p1f3_gr.nc","r")

pr31        =f111->pr(:,{17:55},{72:136})
pr32        =f222->pr(:,{17:55},{72:136})
pr33        =f333->pr(:,{17:55},{72:136})
pr34        =f444->pr(:,{17:55},{72:136})
pr35        =f555->pr(:,{17:55},{72:136})
pr36        =f666->pr(:,{17:55},{72:136})
pr37        =f777->pr(:,{17:55},{72:136})
pr38        =f888->pr(:,{17:55},{72:136})
pr39        =f999->pr(:,{17:55},{72:136})
;printVarSummary(pr39)
;将缺测值设为0
pr31=where(ismissing(pr31),0,pr31)
pr32=where(ismissing(pr32),0,pr32)
pr33=where(ismissing(pr33),0,pr33)
pr34=where(ismissing(pr34),0,pr34)
pr35=where(ismissing(pr35),0,pr35)
pr36=where(ismissing(pr36),0,pr36)
pr37=where(ismissing(pr37),0,pr37)
pr38=where(ismissing(pr38),0,pr38)
pr39=where(ismissing(pr39),0,pr39)
;printVarSummary(pr39)

pr3=(pr31+pr32+pr33+pr34+pr35+pr36+pr37+pr38+pr39)/9
printVarSummary(pr3)


Variable: pr39Type: floatTotal Size: 933888 bytes            233472 valuesNumber of Dimensions: 3Dimensions and sizes:   [month | 6] x [lat | 152] x [lon | 256]Coordinates:            month: [   1..  12]            lat: [17.125..54.875]            lon: [72.125..135.875]Number Of Attributes: 6  long_name :   Monthly Precipiation Climatology  standard_name :       precipiation_flux  units :       mm  coordinates : lon lat  missing_value :       ( -, 9, 9, 9, 9, f )  _FillValue :  -9999


我认认真真对过，每一个数据都是一样的维度，但就是显示维度不匹配。。。。
报错：fatal:Dimension sizes of left hand side and right hand side of assignment do not match

HANSEN

哪一行的维度不匹配？把对应行的左右两边的数组print出来看看就知道了，NCL报错一般不可能出错的，一定是你的维度有问题

ststst

HANSEN 发表于 2024-3-28 08:40
哪一行的维度不匹配？把对应行的左右两边的数组print出来看看就知道了，NCL报错一般不可能出错的，一定是你 ...

就是pr3这里，我打印了pr31,pr32,pr33,pr34,pr35,pr36,pr37,pr38,pr39，都是6*152*256的，我再把这些数据加起来除以9，就会显示维度不匹配

KasuganoSora

如果确定数组没问题的话，pr3这个变量在之前是否定义过

t_imba

大概率你的PR3之前有定义过了，你整个程序查一查看看

HANSEN

ststst 发表于 2024-3-28 09:42
就是pr3这里，我打印了pr31,pr32,pr33,pr34,pr35,pr36,pr37,pr38,pr39，都是6*152*256的，我再把这些数据 ...

那你pr3的维度呢？右边其他都是6*152*256

ststst

KasuganoSora 发表于 2024-3-28 11:51
如果确定数组没问题的话，pr3这个变量在之前是否定义过

的确被定义过，我换了个变量名，就不报错了

ststst

HANSEN 发表于 2024-3-28 15:01
那你pr3的维度呢？右边其他都是6*152*256

已经解决啦，我要计算的就是pr3，打印不出来。后来我换了个变量名，就不报错了

ststst

HANSEN 发表于 2024-3-28 15:01
那你pr3的维度呢？右边其他都是6*152*256

已经解决啦，我要计算的就是pr3，打印不出来。后来我换了个变量名，就不报错了

ststst

t_imba 发表于 2024-3-28 14:35
大概率你的PR3之前有定义过了，你整个程序查一查看看

是之前定义过了，不过为什么之前定义过的就会报错呢，我前后的pr3都是一个用途