Dear R-users
I am working with ncdf data using the variables time (1-365), lon (longitude),
lat (latitude) and the Temperature variable daily). After setting the
parameters for the model, I am able to calculate the output for each lon-lat
grid point. The model works well including one ncdf file (TabsD29) (1). Now I
want to include more than one ncdf file (2). Alls the ncdf files have the same
variables with exception of the temperature variable. However, I get the output
"wrong arrays". I think may suggestion to calculate the mean is wrong?
Many thanks
Sibylle
nlon <- length(lon)
nlat <- length(lat)
nday <- length(TabsD29[1,1,])
Tmin <- 10.
GDDmax <- 145
DOYstart <- 1
(1)
Teffs <- pmax(TabsD29 - Tmin, array(0., dim=c(nlon, nlat, nday)))
(2)
Teff <-
pmax(mean(as.numeric(TabsD80+TabsD81+TabsD82+TabsD83+TabsD84+TabsD85+TabsD86+TabsD87+TabsD88+TabsD89+TabsD90+TabsD91+TabsD92+TabsD93+TabsD94+TabsD95+TabsD96+TabsD97+TabsD98+TabsD99+TabsD20+TabsD21+TabsD22+TabsD23+TabsD24+TabsD25+TabsD26+TabsD27+TabsD28+TabsD29))
- Tmin, array(0., dim=c(nlon, nlat, nday)), na.rm=TRUE)
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
