StackOverflow2013

Note that there are some explanatory texts on larger screens.

plurals

POwriting a loop for upscaling precipitation for USA
text
Body
copied!<p>I am writing a code to calculate the mean amount of precipitation for different regions of conterminous USA. My total data has 300 times 120 (lon*lat) grids in Netcdf format. I want to write a loop in R to take the average of each 10 by 10 number of grids and assign that value (average) to all of the grids inside the region and repeat this for the next region. At the end instead of a 120 by 300 grids I will have 12 by 30 grids. So this is kind a upscaling method I want to apply to my data. I can use a for-loop for each region separately but It makes my code very huge and I don’t want to do that. Any idea would be appreciated. Thanks. P.S: Here is the function I have written for one region (10by10) lat*lon.</p> <pre><code>upscaling <- function(file, variable, start.time=1, count.time=1) { library(ncdf) # load ncdf library to manipulate ncdf data ncdata <- open.ncdf(file); # open ncdf file lon <- get.var.ncdf(ncdata, "lon"); lat <- get.var.ncdf(ncdata, "lat"); time <- get.var.ncdf(ncdata, "time"); start.lon <- 1 end.lon <- length(lon) start.lat <- 1 end.lat <- length(lat) count.lon <- end.lon - start.lon + 1; # count number of longitude count.lat <- end.lat - start.lat + 1; # count number of latitude dat <- get.var.ncdf(ncdata, variable, start=c(start.lon, start.lat, 1), count=c(count.lon, count.lat, 1)) temp.data<- array(0,dim=c(10,10)) for (i in 1:10) { for (j in 1:10) { temp.data <- mean(dat[i,j,]) } } } </code></pre>

Querying!

Guidance

An individual column

Larger individual text columns get their own page to allow for proper reading.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload