Dear David,
Thanks,
I could read and open .nc file in R, but now how to plot a simple filled
color.
Data and variables of .nc after reading in R are as below. Data have time
steps of 99 (monthly mean data for 8 years).
If I want to fill / plot co2 variable at time=1 then how to do it. Also,
how can select lat lon
over a particular region. For example, fill co2 [lat=0:40, lon=50:90,
time=1] and then overplot a continental map.
how to write this fill command in R,
Thanks, Yogesh
> library(ncdf)
Warning message:
package ncdf was built under R version 2.14.2
> file.nc=open.ncdf("AIRS_2002_10_added-all-years.nc")
> print(file.nc)
[1] "file AIRS_2002_10_added-all-years.nc has 3 dimensions:"
[1] "time Size: 99"
[1] "lat Size: 91"
[1] "lon Size: 144"
[1] "------------------------"
[1] "file AIRS_2002_10_added-all-years.nc has 1 variables:"
[1] "float co2[lon,lat,time] Longname:co2 Missval:1e+30"
On Mon, Apr 2, 2012 at 9:20 PM, David William Pierce <dpie...@ucsd.edu>wrote:
> On Mon, Apr 2, 2012 at 1:54 AM, Yogesh Tiwari <yogesh....@googlemail.com>
> wrote:
>
> > How to read netcdf files in R ?
> > Which packeges do we need to install for this,
> > and what commands are used for reading netcdf files.
>
> Typical code would look something like this:
>
> library('ncdf'')
>
> filename <- 'data.nc'
> varname <- 'Temperature'
>
> ncid <- open.ncdf( filename )
> data <- get.var.ncdf( ncid, varname )
>
> If you need to see what variables the file contains to get the
> variable names, just do a "print(ncid)" after you open it.
>
> If it is a very large file, too large for the data to fit in memory
> all at once, you can step through and read each timestep in the file
> separately using the "start=" and "count=" arguments to the
> get.var.ncdf() call.
>
> Regards,
>
> --Dave
>
> --
> David W. Pierce
> Division of Climate, Atmospheric Science, and Physical Oceanography
> Scripps Institution of Oceanography, La Jolla, California, USA
> (858) 534-8276 (voice) / (858) 534-8561 (fax) dpie...@ucsd.edu
>
--
Yogesh K. Tiwari (Dr.rer.nat),
Scientist,
Centre for Climate Change Research,
Indian Institute of Tropical Meteorology,
Homi Bhabha Road,
Pashan,
Pune-411008
INDIA
Phone: 0091-99 2273 9513 (Cell)
: 0091-20-25904452 (O)
Fax : 0091-20-258 93 825
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
