On Fri, May 6, 2011 at 10:11 AM, francoise orain
<francoise.orain@meteo.fr>wrote:

> Hello
> I am a new user of R .
> and I ve problem with R and netcdf .
> I succed installation , I could use all examples .
> But when I take my netcf it is different .
>
> I want to do statistic on this kind of file .
>
> 1)
> first calculate mean .
> my data is like that
> through ncdump -h test.nc
>
> netcdf test {
> dimensions:
>       lat = 301 ;
>       lon = 401 ;
>       time = UNLIMITED ; // (80 currently)
> variables:
>       float lat(lat) ;
>               lat:long_name = "latitude" ;
>               lat:standard_name = "latitude" ;
>               lat:units = "degrees_north" ;
>               lat:valid_range = -60., 60. ;
>       float lon(lon) ;
>               lon:long_name = "longitude" ;
>               lon:standard_name = "longitude" ;
>               lon:units = "degrees_east" ;
>               lon:valid_range = -100., 45. ;
>       double sst(time, lat, lon) ;
>               sst:long_name = "sea surface temperature" ;
>               sst:standard_name = "sea_surface_temperature" ;
>               sst:units = "K" ;
>               sst:_FillValue = -32768. ;
>       double time(time) ;
>               time:long_name = "time" ;
>               time:standard_name = "time" ;
>               time:units = "seconds since 1981-01-01 00:00:00" ;
>               time:ancillary_variables = "sst_data_processed_flag" ;
>
> I succeed to read it with R, there is only one variable sst   80 days and
> lat = 301     lon = 401 (this file is already a concatenation with nco
> fonctions .
>
> Here is what I do :
>
> nc<-open.ncdf("test.nc")
> summary (nc)
> data<- get.var.ncdf(nc)
> print (data)
>
> answer extract :
>
> [331,] -32768.00 -32768.00 -32768.00 -32768.00 -32768.00 -32768.00
> -32768.00
> [332,] -32768.00 -32768.00 -32768.00 -32768.00 -32768.00 -32768.00
> -32768.00
>         [,295]    [,296]    [,297]    [,298]    [,299]    [,300]    [,301]
>  [1,]      0.10      0.10      0.12      0.15      0.15      0.14      0.15
>  [2,]      0.12      0.12      0.13      0.15      0.15      0.14      0.15
>  [3,]      0.13      0.13      0.13      0.16      0.14      0.14      0.14
>  [4,]      0.15      0.13      0.11      0.16      0.15      0.15      0.15
>  [5,]      0.17      0.16      0.15      0.15      0.16      0.17      0.15
>
> #_fill_value is at  -32768.00
>
>
Hello,

I'm not sure what the problem you are experiencing is. As a practical
matter, I suggest you try something like this:

nc<-open.ncdf("test.nc")
summary (nc)
data<- get.var.ncdf(nc)
data[ data < -32000 ] <- NA

Regards,

--Dave

-- 
David W. Pierce
Division of Climate, Atmospheric Science, and Physical Oceanography
Scripps Institution of Oceanography
(858) 534-8276 (voice)  /  (858) 534-8561 (fax)    dpierce@ucsd.edu

	[[alternative HTML version deleted]]

