Question

我正在尝试使用7月的所有WRF观测值创建CSV文件（总共744观测值或744行）。但是我的代码一次只执行一个文件，而最终的output.csv只是一行。我在终端上运行我的R脚本。我记得在文件名的第一部分后使用星号，但是这次不起作用。例如，所有文件都以相同的名称“ wrfout_d04_2013-05-01”开头，而我只需要在其旁边放置（*），这样它就可以在所有以该名称开头的文件上运行。

 #!/usr/bin/env Rscript

 ## This script will read ncdf4 encoded WRF output files and extract
 ## variables T2, RAINC and RAINNC for pairs of Lat, Lon given in a 
 csv
 ## file.
 ##
 ## Usage:
 ##
 ##     $ ./extract_for_validation.R weather_stations.csv 
   wrfout_d04_2013-05-01_00\:00\:00 output.csv
 ##
 ## Where
 ## - "weather_stations.csv" contains pairs of Lat, Lon in rows
 ## - wrfout_d04_2013-05-01_00\:00\:00 is a ncdf4 encoded file 
  output by WRF
 ## - output.csv is where triads of T2, RAINC, RAINNC will be saved, 
    in
 ##   a single row, one for each pair of Lat, Lon given at input
 ##

 library(ncdf4)
 library(ncdf4.helpers)

  args <- commandArgs(TRUE)

  estaciones <- read.csv(args[1])
  climate_output <- nc_open(args[2])

  lon <- ncvar_get(climate_output, varid = "XLONG")
  lat <- ncvar_get(climate_output, varid = "XLAT")
  temp <- ncvar_get(climate_output, varid = "T2")
  pp0 <- ncvar_get(climate_output, varid = "RAINC")
  pp1 <- ncvar_get(climate_output, varid = "RAINNC")

  date <- tail(strsplit(args[2], "_")[[1]], n=2)[1]
  time <- tail(strsplit(args[2], "_")[[1]], n=2)[2]

   out <- c(c(date, time))

 for (x in 1:length(estaciones[,1])){
    column <- match(round(estaciones[x ,1], digits=2), 
    round(lat[1,], digits=2))
    row <- match(round(estaciones[x, 2], digits=2), round(lon[,1], 
    digits=2))
    out <- append(out, (c(temp[row, column], pp0[row, column], 
    pp1[row, column])))
    }

    write.table(t(unlist(out)), file=args[3], row.names=FALSE, 
    col.names=FALSE)

如何扩展它以同时在所有744个netcdf文件上运行，并创建包含所有744个观测值的单个output.csv？

如何将多个netcdf文件（每个文件是一个观察结果）提取到单个output.csv文件中？

0 个答案: