请问有大神会用python对nc数据进行切割么？

leeying422

我现在手里有一些nc数据，每一个nc数据都是很多天的，我现在想把它切割成每一天的数据单独成为一个文件，不知道有没有可用的脚本可以参考呢？或者是我下载的是ERA5的数据，请问有下载ERA5数据每天单独成一个文件的脚本么？多谢大佬们指教

wjy_ecnu

zhe难道不是多此一举吗

qzf688018

可能是大范围网格，层数也多的原因吧。试试循环吧，每个时间循环打开一个nc，然后把该时刻的数据写到这个nc里边。

bahzhong

# -*- coding: GBK -*-

import cdsapi
import calendar

c = cdsapi.Client()
for year in range(2000,2020):
        for mon in range(1,13):
                monthRange = calendar.monthrange(year,mon)
                for day in range(1, monthRange[1]+1):
                        print('------------  Download %s%s%s'% ("%04d" % year,"%02d" % mon,"%02d" %day))
                        c.retrieve(
                                'reanalysis-era5-pressure-levels',
                                {
                                        'product_type': 'reanalysis',
                                        'variable': [ 'geopotential','u_component_of_wind',
                    'v_component_of_wind', 'vorticity',],
                                        'pressure_level': ['700','850',],
                                        'year':str(year),
                                        'month':("%02d" % mon),
                                        'day':("%02d" % day),
                                        'time':['00:00','06:00', '12:00', '18:00',],
                                    'area':'40/95/25/111',
                                        'grid':'0.25/0.25',
                                        'format': 'netcdf',
                                },
                                'd:\\era5\\CDS%s%s%s'% ("%04d" % year,"%02d" % mon,"%02d" %day)+'.nc')

Lancelot

cdo