R/Project.R

In SHUD-System/SHUDtoolbox: Toolbox For SHUD Model System

#' Prepare Input file of SHUD model
#' \code{shud.filein} 
#' @param projname Character, project name, default= PRJNAME which is a global variable.
#' @param inpath SHUD model input directory, default = inpath which is global variable
#' @param outpath SHUD model output directory, default = outpath which is global variable
#' @return Character of full path of input files for SHUD model
#' @export
shud.filein <- function(projname=get('PRJNAME', envir=.shud),
                        inpath = get('inpath', envir=.shud), 
                        outpath = get('outpath', envir=.shud)
                        ){
  dir.create(inpath, showWarnings = FALSE, recursive = TRUE)
  fn.mesh = file.path(inpath, paste0(projname, '.sp.mesh' ) );
  fn.att = file.path(inpath, paste0(projname, '.sp.att' ) );
  fn.pt = file.path(inpath, paste0(projname, '.sp.points' ) )
  fn.edge = file.path(inpath, paste0(projname, '.sp.edges' ) )
  fn.mdriv = file.path(inpath, paste0(projname, '.sp.riv' ) )
  fn.mdseg = file.path(inpath, paste0(projname, '.sp.rivseg' ) )
  
  fn.ic = file.path(inpath, paste0(projname, '.cfg.ic' ) )
  fn.para = file.path(inpath, paste0(projname, '.cfg.para' ) )
  fn.calib = file.path(inpath, paste0(projname, '.cfg.calib' ) )
  
  fn.soil = file.path(inpath, paste0(projname, '.para.soil' ) )
  fn.geol = file.path(inpath, paste0(projname, '.para.geol' ) )
  fn.lc = file.path(inpath, paste0(projname, '.para.lc' ) )
  
  fn.forc = file.path(inpath, paste0(projname, '.tsd.forc' ) )
  fn.bc = file.path(inpath, paste0(projname, '.tsd.bc' ) )
  fn.lai = file.path(inpath, paste0(projname, '.tsd.lai' ) )
  fn.rl = file.path(inpath, paste0(projname, '.tsd.rl' ) )
  fn.mf = file.path(inpath, paste0(projname, '.tsd.mf' ) )
  gisdir = file.path(inpath, 'gis')
  fns = c(
    inpath, outpath, gisdir, 
    fn.mesh, fn.pt, fn.edge, fn.mdseg, fn.mdriv, fn.att,
    fn.ic,fn.para, fn.calib,
    fn.soil, fn.geol, fn.lc,
    fn.forc, fn.bc, fn.lai, fn.rl, fn.mf)
  names(fns) = c(
    'inpath', 'outpath', 'gispath',
    "md.mesh","md.pt", "md.edge", "md.rivseg", "md.riv", "md.att",
                 'md.ic', 'md.para', 'md.calib',
                 'md.soil', 'md.geol', 'md.lc',
                 'md.forc', 'md.bc', 'md.lai', 'md.rl', 'md.mf' )
  fns
}