计算PMcoarse

;============================= ; 计算 PMcoarse = PM10 - PM2.5 ;============================= begin version = "MEICv2.0" in_year = 2024 out_year = 2024 mech = "SAPRC07" sector = "5" in_dir = "/....../" out_dir = "/....../" months = (/ "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12" /) departments = (/ "agriculture", "industry", "power", "residential", "transportation" /) ; === 根据你的 ASC 文件头行数修改：6 === header_lines = 6 ; 输入输出路径 in_file_dir = in_dir + version + "_" + in_year + "_" + mech + "/" out_file_dir = out_dir + version + "_" + in_year + "_" + mech + "/" ;============================= ; 循环每个月、每个部门 ;============================= do q = 0, dimsizes(months)-1 do w = 0, dimsizes(departments)-1 month = months(q) department = departments(w) in_asc_pm10 = in_file_dir + in_year + "_" + month + "_" + department + "_PM10.asc" in_asc_pm25 = in_file