合并多年的逐小时nc文件案例，求指教改进

主啊

# 将所有的nc文件进行合并，生成一个nc文件
# 多年的nc数据

import xarray as xr
import netCDF4 as nc
import os
import numpy as np
import pandas as pd

# 外置函数
# 列出文件下所有文件名，返回输入文件路径+文件名
def getALLFiles(targetDir):
    listFiles = os.listdir(targetDir)
    for i in range(len(listFiles)):
        listFiles = targetDir + "/" +listFiles
        #print(listFiles)
    return listFiles


# 获取文件夹下的年份
targetDir = "E:/xinzheng2021-2022/xinzheng"
years = os.listdir(targetDir)

# 年份循环
# 逐小时数据
file_new = []
for everyyear in years:

    # 当年nc路径构造
    ncfilepath = targetDir + "/" + everyyear
    files = getALLFiles(ncfilepath)
    #当年逐日数据读取拼接
    for file in files:
        timefiles = getALLFiles(file)
        for timefile in timefiles:
            # print(timefile[-16:])
            data = xr.open_dataset(timefile)
            data["Temperature_height_above_ground"] = data["Temperature_height_above_ground"] - 273.15
            file_new.append(data)


da = xr.concat(file_new, dim='time')
print(da)
print(da.time.values)
da.to_netcdf('E:/xinzheng2021-2022/mulyear.nc')

万里路始于足下

楼主这个代码为啥在帖子上这么漂亮指点指点{:eb302:}

主啊

万里路始于足下 发表于 2023-12-1 23:55
楼主这个代码为啥在帖子上这么漂亮指点指点

直接粘贴的