fluxUnit suggestion
Created by: emorway-usgs
For your customers with native data units in ng/L, it would be helpful to have a fluxUnit option of kg/yr. After pulling the data from NWIS, I do the recommended conversion of:
Sample$ConcLow <- Sample$ConcLow / 1000000 Sample$ConcHigh <- Sample$ConcHigh / 1000000 Sample$ConcAve <- Sample$ConcAve / 1000000
to get the data into the expected units of mg/L. After progression through the WRTDS workflow and running tableChange(), I ignore the warning message at the end of the screen output regarding units of ng/L in INFO, but am hopeful that a units value for these smaller annual fluxes could be provided? If you would prefer users do the conversion themselves, no problem, only a suggestion.
tableChange(eList, fluxUnit=8, yearPoints=c(1998,2005,2013))
... Flux Trends time span change slope change slope 10^3 kg/yr 10^3 kg/yr /yr % %/yr 1998 to 2005 0.00076 0.00011 10 1.5 1998 to 2013 0.0052 0.00035 70 4.7 2005 to 2013 0.0044 0.00055 54 6.8
or changing fluxUnit=3 yields
tableChange(eList, fluxUnit=8, yearPoints=c(1998,2005,2013)) ... Flux Trends time span change slope change slope kg/day kg/day /yr % %/yr 1998 to 2005 0.0021 3e-04 10 1.5 1998 to 2013 0.014 0.00095 70 4.7 2005 to 2013 0.012 0.0015 54 6.8