使用iris的merge命令时,有时会出现由于两个或多个cube的坐标完全一样,Python无法区别,因此会认为是同一个cube重复append而无法合并的问题。这时需要在append之后给每个cube加一个index再合并。
import iris
......
all_cubes.append(cube)
cube_list=iris.cube.CubeList(all_cubes)
for index,cube in enumerate(cube_list):
cube.add_aux_coord(iris.coords.AuxCoord(index,long_name='new_time',units='no_unit'))
# long_name可以是任何自己想取的坐标名字,如果改成standard_name,则坐标名字只能是time,latitude,longitude等cube 固定用的坐标名。
new_cube=cube_list.merge_cube()