Converting a variable's units from gigatons to teratons
ncvar_get() should return a data frame, so just get2.teratons <- get2 / 1000 Even if it's a higher-dimensional array, that will work. You could then write it back with ncvar_put(), subject to the data type of that variable in the NetCDF file. Sarah On Thu, Mar 28, 2019 at 11:55 AM rain1290--- via R-sig-Geo
<r-sig-geo at r-project.org> wrote:
Hi there,
I am currently working with a variable whose units are expressed in gigatons. However, I would like to convert all of these values in this variable to teratons. Effectively, I would have to somehow have R divide all values (90 values) within the variable by 1000.
So far, I have read in the variable (called CanESM2) as follows:
ncfname1 <- "cumulative_emissions_RCP45.nc"
Model3 <- nc_open(ncfname1)
get2 <- ncvar_get(Model3, "CanESM2")
Is there a way to accomplish this conversion?
Thanks, and any assistance would be greatly appreciated!
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Sarah Goslee (she/her) http://www.numberwright.com