批量读入数据后，选取其中一段时间的数据，求得其最大值？

tiantianyixin

批量读入nc数据后，选取某一时间段，求这时间段内数据的最大值，已完成的脚本如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin                              
all_files = systemfunc ("ls /home/chenyg/ncl/air.*.nc")
fall      = addfiles (all_files, "r")   ;读入的数据个数为2
ListSetType (fall, "join")
T    = fall[:]->air  
tt =T(:,0:89,{20:70},{30:130})
？？？
end
;选取每个数据前90天，这样得到的tt维数表现为（2,90,26,54），我想求出各个格点上这180天（2*90）中的最大值？
我的想法是把前两个维数合并，但是实现不了，还请大神指点！

sun_shine_Xia

你新定义一个数组第一维是180，分别把前后各90天的数组放进去就好了~

tiantianyixin

sun_shine_Xia 发表于 2014-11-23 15:08
你新定义一个数组第一维是180，分别把前后各90天的数组放进去就好了~

谢谢你哈！只有这一种办法吗？有没有什么函数可以实现呢

sun_shine_Xia

有的吧，你试试dim_max_n，我还没用过这个函数，不知道能不能实现~

sun_shine_Xia

上面那个函数是可以直接找最大值，reshape这个函数可以调整数组维数~

芦苇草

{:5_213:}

longlivehj

tt（2,90,26,54）前两维最大值
dim_max_n(tt, (/0, 1/))

tiantianyixin

longlivehj 发表于 2014-11-24 23:22
tt（2,90,26,54）前两维最大值
dim_max_n(tt, (/0, 1/))

多谢多谢大神！

tiantianyixin

sun_shine_Xia 发表于 2014-11-23 15:20
上面那个函数是可以直接找最大值，reshape这个函数可以调整数组维数~

好的，我试试，多谢啦

姨妈血红

也可以分别读两个数组，再把两个数组拼在一起，然后用dim_max_n_Wrap