import xarray as xr import numpy as np import pandas as pd # === 文件路径 === swe_path = 'E:/1961-2020_swe_mm/ten_year/clip/1961-1970_clip.nc' rain_path = 'F:/gldas_10year_0.1_sep/clip/1961-1970_clip.nc' output_path = 'E:/RP_melt/rp_melt_month/1961-1970_rp_melt_month.nc' # === 加载 SWE 数据 === ds_swe = xr.open_dataset(swe_path) swe = ds_swe['sd_mm'] # 实际变量名 swe = swe.rename({'valid_time': 'time'}) swe['time'] = pd.to_datetime(swe.time.values) # === 加载降水数据 === ds_rain = xr.open_dataset(rain_path) ra