r - Extracting values from nc files in a for loop -



r - Extracting values from nc files in a for loop -

i have 4 netcdf files in folder , want extract values these files. str of file is:

[1] "file c:1.dbl.nc has 2 dimensions:" [1] "lat size: 1" [1] "lon size: 1" [1] "------------------------" [1] "file c:\\users\\data.nc has 3 variables:" [1] "short so[lon,lat] [1] "short il[lon,lat] [1] short fg[lon,lat]

my loop is:

a<-list.files("c:\\users\\data", "*.nc", full.names = true) for(i in 1:length(a)){ f <- open.ncdf(a[i]) = get.var.ncdf(nc=f,varid="so",verbose=true) b <- get.var.ncdf(nc=f,varid="il") c <- get.var.ncdf(nc=f,varid="fg") write.table(t(rbind(a,b,c)),file="output-all.txt")}

there no errors there 1 line of results in output text file:

"a" "b" "c" "1" 500 200 300

which should many number of files in folder (4), instance:

"a" "b" "c" "1" 500 200 300 "2" 500 200 300 "3" 500 200 300 "4" 500 200 300

it seems loop not taking place here.

move write.table command out of loop. like:

a<-list.files("c:\\users\\data", "*.nc", full.names = true) d<-matrix(na,length(a),3) for(i in 1:length(a)){ f <- open.ncdf(a[i]) = get.var.ncdf(nc=f,varid="so",verbose=true) b <- get.var.ncdf(nc=f,varid="il") c <- get.var.ncdf(nc=f,varid="fg") d[i,]<-t(rbind(a,b,c))} write.table(d,file="output-all.txt")

r for-loop netcdf

Comments

Post a Comment

Popular posts from this blog

xslt - DocBook 5 to PDF transform failing with error: "fo:flow" is missing child elements. Required content model: marker* -

xslt - DocBook 5 to PDF transform failing with error: "fo:flow" is missing child elements. Required content model: marker* - i've inherited doc publishing process takes docbook , builds html , pdf output, using apache fop both. project had been started, never completed. lots of tweaking, i've able of doc sets build (20+), except handful, , fail severe code 1 of these 2 conditions: "fo:flow" missing kid elements. required content model: marker* "fo:block" not valid kid of "fo:root" the xml output comes out of docbook , docbook-fop checks out well-formed. the fop error in 1 book, "endecagloss", example, points end of string: <fo:flow flow-name="xsl-region-body" start-indent="4pc" end-indent="0pt"/> here’s log output in case: [java] severe: javax.xml.transform.transformerexception: file:/scratch/publishing/hudson/jobs/build-main-endeca-documentation/workspace/serverdoc/w...
Read more

mediawiki - How do I insert tables inside infoboxes on Wikia pages? -

mediawiki - How do I insert tables inside infoboxes on Wikia pages? - i want insert tables in infoboxes. example, have proteinbox i've developed using infobox template (http://health-and-medicine.wikia.com/wiki/template:proteinbox), , fasta field big fit comfortably in proteinbox. hence collapsible table within box. this tried: {{infobox | row 1 = {| class = "mw-collapsible mw-collapsed wikitable" |- | <!-- table content --> |} }} your problem is, putting nested table within template. pipe characters ( | ) in table syntax collides pipe usage in templates. the commonly used hack around create template called template:! , containing pipe character, , utilize when need set tables, parser functions or other stuff using pipe characters, within templates. table (with every | replaced {{!}} ): {{{!}} {{!}}- {{!}} a1 {{!}} b1 {{!}}- {{!}} a2 {{!}} b2 {{!}}} ...the equivalent of {| |- | a1 | b1 |- | a2 | b2 |} furthermore, hav...
Read more

Local Service User Logged into Windows -

Local Service User Logged into Windows - i created windows service tracks app server user connected through load balancer. it retrieves user logged windows on machine, cross references active directory (ad) business relationship homecoming first name, lastly name , e-mail (for notification through website created project) there 1 scheme on network users remote system, vista system, reason returning values "local service" logged machine, it's weird , security guys impossible "local service" business relationship logged windows. here how user logged windows managementscope ms = new managementscope("\\\\.\\root\\cimv2"); objectquery query = new objectquery("select * win32_computersystem"); managementobjectsearcher searcher = new managementobjectsearcher(ms, query); foreach (managementobject mo in searcher.get()) { _username = mo["username"].tostring(); } // remove domain part username string[] usernameparts = _...
Read more