R/config_read.R

In nauta008/VIC_setup:

CONSTANTS <- list(routing=list())
#CONSTANTS$grid_mapping <- list()
#CONSTANTS$grid_mapping$lambert_azimuthal_equal_area <- list(name="lambert_azimuthal_equal_area",alias="laea")
#CONSTANTS$grid_mapping$latitude_longitude <- list(name="latitude_longitude",proj4_params="+proj=longlat +datum=WGS84 +no_defs", EPSG_code="EPSG:4326")
CONSTANTS$grid_mapping$lookup <- list(
  #aea = "albers_conical_equal_area",
  #aeqd = "azimuthal_equidistant",
  laea = "lambert_azimuthal_equal_area",
  #lcc = "lambert_conformal_conic",
  #cea = "lambert_cylindrical_equal_area",
  longlat = "latitude_longitude"
  #merc = "mercator",
  #omerc = "oblique_mercator",
  #ortho = "orthographic",
  #stere = "stereographic",
  #tmerc = "transverse_mercator"
)

CONSTANTS$grid_mapping$default <- list(
  grid_mapping_name = "latitude_longitude",
  semi_major_axis = 6378137,
  inverse_flattening = 298.257223563,
  longitude_of_prime_meridian = 0
)

routing.origin=list()
routing.origin[['LISFLOOD']] <- list(outlet_val=5)
routing.origin[['RVIC']] <- list(outlet_val=9)
routing.origin[['ARCMAP']] <- list(outlet_val=9)

CONSTANTS$routing[['origin']] <- routing.origin
rm(routing.origin)

missing.vals <- list(float=1e20, integer=-.Machine$integer.max, double=1e20, short=-32,767,byte=-127)
CONSTANTS$missing_vals <- missing.vals
rm(missing.vals)


VICSetup <- new.env()

config.read <- function(file){
  config.restore()
  VICSetup$config <- read_yaml(file)
  defaults.set()
  log_info(sprintf("Using the following configuration: \n %s", as.yaml(VICSetup$config)))
  domain.init()
  # create setup output
  if(!is.null(VICSetup$config$output)){
    if(!is.null(VICSetup$config$output$domain)){
      domain.create()
    }
    if(!is.null(VICSetup$config$output$routing)){
      defaults.routing.set()
      rout.params.create()
    }
    if(!is.null(VICSetup$config$output$params)){
      vic.params.create()
    }
  }
}

config.restore <- function(){
  rm(list = ls(envir=VICSetup), envir = VICSetup)
}