I read a ncdf file but I can't get the correct time info and I don't
know where is my error. Any help please. This is my code.
|# rm(list=ls())library(ncdf)nc=open.ncdf("tmp.1979-2014.nc")nc
[1]"file tmp.1979-2014.nc has 3 dimensions:"[1]"lon Size: 55"[1]"lat
Size: 41"[1]"time Size: 52596"[1]"------------------------"[1]"file
tmp.1979-2014.nc has 1 variables:"[1]"float tmp[lon,lat,time]
Longname: Temperature Missval:1e+30"lon=nc$dim$lon$vals
lat=nc$dim$lat$vals lev=nc$dim$level$vals time=nc$dim$time$vals
units=nc$dim$time$units start=as.POSIXct('1979-01-01
00:00',tz='UTC')end=as.POSIXct('2014-12-31
18:00',tz='UTC')times=as.POSIXct(ncin$dim$time$vals*60*60,origin='1957-01-0100:00:00',tz='UTC')t1=which(times
==start)t2=which(times ==end)t1 [1]369t2 [1]51504|
t1 and t2 must be respectively 1 and 52596 according to my file.
1. What is wrong in my code?
2. How can I calculate daily mean of my hourly file?
transform hourly data matrix to daily array
1 message · Sancta Vega