search for: nlat

...;d appreciate your assistance. TIA, Tom Roche <Tom_Roche at pobox.com>---2 IDL routines follow to EOF--- from ftp://gfed3:dailyandhourly at zea.ess.uci.edu/GFEDv3.1/Readme.pdf ;++++++++++++++++++++ idl code to generate daily emissions +++++++++ nlon=720 ; number of grid points by longitude nlat=360 ; number of grid points by latitude gfedmly=fltarr(nlon,nlat) ; array containing monthly emissions gfeddly=fltarr(nlon,nlat) ; array containing daily emissions ; You must read monthly emissions to generate daily fluxes. ; For example, if you want daily emissions for January 21st, 2004, ; you n...

...t;,lat=73,long=144) { #d?finition de la fonction library(maps) a <- read.table(nomfichier) #lecture du fichier z <- array(0,dim=c(long,lat)) #cr?ation d'un tableau de 0 z <- t(a) #transposition de la matrice a en z z1<-rbind(z[65:long,],z[1:64,]) #centrage de la carte sur greenwich nlat<-seq(-90,90,2.5) #d?finition des latitudes / ordonn?es de 2,5? en 2,5? nlong <- seq(-180,177.5,2.5) #d?finition des longitudes /abcisses de 2,5? en 2,5? image(nlong,nlat,z1,col=rainbow (100),xlab="longitude",ylab="latitude") #tra?age map(add=TRUE) #superposition de la...

...fname <- paste(ncname, ".nc", sep = "") dname <- "prAdjust" # open a netCDF file ncin <- nc_open(ncfname) print(ncin)# get longitude and latitude lon <- ncvar_get(ncin,"lon") nlon <- dim(lon) head(lon) lat <- ncvar_get(ncin,"lat") nlat <- dim(lat) head(lat) print(c(nlon,nlat)) # get time time <- ncvar_get(ncin,"time") time tunits <- ncatt_get(ncin,"time","units") nt <- dim(time) nt tunits # get variable preci.array <- ncvar_get(ncin,dname) dlname <- ncatt_get(ncin,"prAdjust&...

...ed the code I use to build the netcdf file and 'ncdf' object below. Thanks in advance, Simon Knapp #--------------------------R code starts here-------------------------- # create a netcdf file holding a 'temporal stack' combined.ncdf <- with(new.env(), { require(ncdf) nLat <- 7 nLon <- 8 nTime <- 9 missing.val <- -999 filename <- file.path(tempdir(), 'combined_d6fg4s64s6g4l.nc') lat.dim <- dim.def.ncdf('lat', 'degrees', 1:nLat) lon.dim <- dim.def.ncdf('lon', 'degrees', 1:nLon) ti...

...quot;") > dname <- "prAdjust" > # open a netCDF file > ncin <- nc_open(ncfname) > print(ncin)# get longitude and latitude > lon <- ncvar_get(ncin,"lon") > nlon <- dim(lon) > head(lon) > > lat <- ncvar_get(ncin,"lat") > nlat <- dim(lat) > head(lat) > > print(c(nlon,nlat)) > # get time > time <- ncvar_get(ncin,"time") > time > > tunits <- ncatt_get(ncin,"time","units") > nt <- dim(time) > nt > tunits > # get variable > preci.array <- n...

.../* returned by malloc in ncread */ double ofs, scl; /* offset and scaling factors */ double *ptr; SEXP field, dim, dimnames, lon, lat, tim, missing; printf("ncfile= %s, vname=%s\n",fname,vname); printf("\nfield->%x\nlon -> %x\nlat -> %x\ntim -> %x\nmissing -> %x\n", field,lon,lat,tim,missing); /* Call the function dealing with the netCDF data */ /* The data is placed in memeory */ status=ncread(fname, vname, &nx, &ny, &nt, &ptrs, &scl, &ofs); /* Allocate space for poi...

...6288e+12 2005-05-17 -0.7524182200 230 1.116374e+12 2005-05-18 0.0271098670 Here is the code wich forms input matrix to make_stars function: transpose_stars = function(lat_b, long_b,lat_t,long_t,step){ matrix <- readSO2Data(lat_b, long_b) // read an arbitrary point for calculating matrix size nlat <- (lat_t-lat_b)/step+1 nlong <- (long_t-long_b)/step+1 ntime <- length(matrix$value) list_lat<-list() list_long<-list() list_time<-list() for( m in 1:nlat){ list_lat[m]<-toString(lat_b + step*(m-1)) } for( m in 1:nlong){ list_long[m]<-toString(long_b + step*(m-1)) } for(...

...directory mask = 0775 create mask = 0775 read list = virusbuster [vizsga] path = /vizsga browseable = Yes directory mask = 755 create mask = 755 write list = root _________________________________________________________________ 711 ?ll?saj?nlat k?z?tt biztosan tal?lsz olyat, ami Neked is megfelel! http://ad.adverticum.net/b/cl,1,6022,105302,170442/click.prm