R/readSWAT.R

In SWATmodel: A multi-OS implementation of the TAMU SWAT model.

readSWAT <-
function(outfile_type){
if(missing(outfile_type)){print(" 'outfile_type' is missing, should be rch, sub, or.. ask drf28 for other types")}
if (outfile_type=="sub"){
   varformat="x6,a4,1x,a8,1x,a4,a10,30a10"
   dataformat="x6,i4,1x,i8,1x,i4,f10,30f10"
} else if (outfile_type=="rch"){
   varformat="x6,a4,1x,a8,1x,a5,30a12"
   dataformat="x6,i4,1x,i8,1x,i5,30f12"
} else { print ("You need to add your file type to this function if it is not output.sub or output.rch")}
  vfrformat = unlist(strsplit(as.character(varformat), ","))
  dfrformat = unlist(strsplit(as.character(dataformat), ","))
  outvars=read.fortran(paste("output.",outfile_type,sep=""),vfrformat,skip=8,nrows=1)
  outdata=read.fortran(paste("output.",outfile_type,sep=""),dfrformat,skip=9,col.names=outvars)
  return(outdata)
}