APSIM 模型运行参数设置与结果保存

本教程展示如何使用 R 语言设置 APSIM 模型运行参数，运行模型并保存结果。

1. 设置模型运行参数

sim_parameters <- list(
  simulation_name = 'my_simulation',
  override_met = TRUE, # 是否使用默认气象数据
  datapath = '.', # 数据文件夹路径
  soilpath = '.', # 土壤文件夹路径
  version = '7413.3', # APSIM版本号
  start_date = as.Date('1990-01-01'),
  end_date = as.Date('2000-12-31'),
  run_type = 'auto', # 运行类型（auto或batch）
  num_workers = 1 # 并行运行的进程数
)

2. 读取土壤参数

soil_df <- read.csv('E:/weidan/SOIL.csv')

3. 循环运行模型

for (i in 1:nrow(soil_df)) {
  site <- soil_df$site[i]
  lat <- soil_df$lat[i]
  lon <- soil_df$lon[i]
  
  for (depth in c(5, 15, 30, 60, 100, 200)) {
    # 设置土壤参数
    soil_parameters <- list(
      Site = site,
      Latitude = lat,
      Longitude = lon,
      Depth = depth,
      LL15 = soil_df$LL15[i],
      DUL = soil_df$DUL[i],
      SAT = soil_df$SAT[i],
      BD = soil_df$BD[i]
    )
    
    # 运行APSIM模型
    file_apsim <- 'E:/weidan/region1wheat/SHIYONG.apsim'
    writeLines(apsim, file_apsim)
    system(paste(apsim_exe, file_apsim), wait = TRUE, ignore.stdout = TRUE, show.output.on.console = FALSE)
    
    # 读取输出结果，并保存为CSV文件
    output_file <- paste0(result_folder, '/', site, '_depth', depth, '.csv')
    output_df <- read_APSIM_output(output_file)
    write.csv(output_df, output_file, row.names = FALSE)
  }
}

注意： 代码中的 read_APSIM_output 函数需要用户自行定义。该函数用于读取 APSIM 模型的输出结果。

4. 定义结果保存路径

result_folder <- 'E:/weidan/region1wheat/SHIYONG'

# 创建文件夹
dir.create(result_folder, showWarnings = FALSE)

5. 错误处理

代码中的错误提示 Error in is.data.frame(x) : object 'output' not found 是因为在之前的代码中没有定义 output 对象。需要先定义该对象才能运行后续代码。

6. 总结

本教程介绍了如何使用 R 语言设置 APSIM 模型运行参数，运行模型并保存结果。希望本教程对您有所帮助。