07c-dataPrep_fS.R

## NOTE: 07a-dataPrep_2001.R and 07b-dataPrep_2011.R need to be run before this script
do.call(setPaths, dataPrepPaths)

gid_fSsimDataPrep <- gdriveSims[studyArea == studyAreaName & simObject == "fSsimDataPrep", gid]
upload_fSsimDataPrep <- reupload | length(gid_fSsimDataPrep) == 0

fSdataPrepParams <- list(
  fireSense_dataPrepFit = list(
    .studyAreaName = studyAreaName,
    .useCache = ".inputObjects",
    fireYears = 2001:2020,
    igAggFactor = 10000 / unique(raster::res(simOutPreamble[["rasterToMatch"]])),
    sppEquivCol = simOutPreamble[["sppEquivCol"]],
    useCentroids = TRUE,
    usePCA = FALSE,
    whichModulesToPrepare = c("fireSense_IgnitionFit", "fireSense_EscapeFit", "fireSense_SpreadFit")
  )
)

simOutPreamble$rasterToMatch <- raster::mask(simOutPreamble$rasterToMatch, simOutPreamble$studyArea)
fSdataPrepObjects <- list(
  .runName = runName,
  cohortData2001 = biomassMaps2001[["cohortData"]],
  cohortData2011 = biomassMaps2011[["cohortData"]],
  historicalClimateRasters = simOutPreamble[["historicalClimateRasters"]],
  pixelGroupMap2001 = biomassMaps2001[["pixelGroupMap"]],
  pixelGroupMap2011 = biomassMaps2011[["pixelGroupMap"]],
  rasterToMatch = simOutPreamble[["rasterToMatch"]],
  rstLCC = biomassMaps2011[["rstLCC"]],
  sppEquiv = simOutPreamble[["sppEquiv"]],
  standAgeMap2001 = biomassMaps2001[["standAgeMap"]],
  standAgeMap2011 = biomassMaps2011[["standAgeMap"]],
  studyArea = simOutPreamble[["studyArea"]]
)

invisible(replicate(10, gc()))

ffSsimDataPrep <- simFile(paste0("fSsimDataPrep_", studyAreaName), Paths$outputPath, ext = simFileFormat)
if (isTRUE(usePrerun)) {
  if (!file.exists(ffSsimDataPrep)) {
    googledrive::drive_download(file = as_id(gid_fSsimDataPrep), path = ffSsimDataPrep)
  }
  fSsimDataPrep <- loadSimList(ffSsimDataPrep)
} else {
  fSsimDataPrep <- Cache(
    simInitAndSpades,
    times =  list(start = 2011, end = 2011),
    params = fSdataPrepParams,
    objects = fSdataPrepObjects,
    paths = dataPrepPaths,
    modules = "fireSense_dataPrepFit",
    .plots = NA,
    #useCloud = useCloudCache,
    #cloudFolderID = cloudCacheFolderID,
    userTags = c("fireSense_dataPrepFit", studyAreaName)
  )
  saveSimList(fSsimDataPrep, ffSsimDataPrep, fileBackend = 2)
}

if (isTRUE(upload_fSsimDataPrep)) {
  fdf <- googledrive::drive_put(media = ffSsimDataPrep, path = as_id(gdriveURL), name = basename(ffSsimDataPrep))
  gid_fSsimDataPrep <- as.character(fdf$id)
  rm(fdf)
  gdriveSims <- update_googleids(
    data.table(studyArea = studyAreaName, simObject = "fSsimDataPrep",  runID = NA,
               gcm = NA, ssp = NA, gid = gid_fSsimDataPrep),
    gdriveSims
  )

  source("R/upload_fSDatPrepFit_vegCoeffs.R")
}