R/SaveRatios.R

In MC2toPath: Translates information from netcdf files with MC2 output into inter-PVT transitions.

SaveRatios = function(rawValues, outFile) {
   nRaw = length(rawValues)
   stopifnot(nRaw>=1)
   
   meanVal = mean(rawValues, na.rm=TRUE)
   stopifnot(meanVal!=0)
   
   ratios = rawValues/meanVal
   
   appendFlag = FALSE
   for (ndx in 1:nRaw) {
      if (is.na(ratios[ndx])) ratios[ndx] = 0
      cat(c(ndx - 1, ", ", ratios[ndx], "\n"), file=outFile, append=appendFlag)
      appendFlag = TRUE
   }
   
   return(ratios)
}