Help with readBin

-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf
Of kapo coulibaly
Sent: Thursday, May 03, 2012 10:57 AM
To: r-help at r-project.org
Subject: Re: [R] Help with readBin

 I believe here is the structure of the file I'm trying to read:
record marker (4 bytes), 2 integers (4 bytes each), 2 doubles (8 bytes
each), one string (16 bytes or 16 characters), 3 integers (4 bytes each), 1
record marker (4 bytes) and a big array of doubles (8 bytes each).
Everything in the file is read correctly except for the doubles.
If any indication, I've read similar file before with readBin the only
difference is this one was created with a code compiled with gfortran in
linux 64 bit. I was able to read the same output binary file when the
fortran source code was compiled in windows xp 32 bit. The values I'm
expecting should be between 0 and about 32.




The code I used is:



# Loading Required libraries
library(tcltk)

# Tk inputbox function
 inputBox<-function() {
  tt<-tktoplevel()
  Zmin<-tclVar("0")
  Zmax<-tclVar("0")
  dZ<-tclVar("0")
  entry.Zmin<-tkentry(tt,width="20",textvariable=Zmin)
  entry.Zmax<-tkentry(tt,width="20",textvariable=Zmax)
  entry.dZ<-tkentry(tt,width="20",textvariable=dZ)
  lbl.Zmin<-tklabel(tt,text="Number of layers")
  lbl.Zmax<-tklabel(tt,text="Number of Stress Periods")
  lbl.dZ<-tklabel(tt,text="dZ")
  tkgrid(lbl.Zmin,entry.Zmin)
  tkgrid(entry.Zmin)
  tkgrid(lbl.Zmax,entry.Zmax)
  tkgrid(entry.Zmax)
  #tkgrid(lbl.dZ,entry.dZ)
  #tkgrid(entry.dZ)

  OnOK <- function()
  {
      # NameVal <- c(tclvalue(Zmin),tclvalue(Zmax),tclvalue(dZ))
        tkdestroy(tt)
    }
  OK.but <-tkbutton(tt,text="   OK   ",command=OnOK)
  # tkbind(entry.Name, "<Return>",OnOK)
  tkgrid(OK.but,columnspan=3)
  tkfocus(tt)
  tkwait.window(tt)
  res<-as.numeric(c(tclvalue(Zmin),tclvalue(Zmax)))#,tclvalue(dZ)))
  return(res)
}

########################################################################
########
# Main program
########################################################################
########

# Model Parameters input (number of layers and stress periods)
param<-inputBox()

# Select and open Modflow Binary file for reading
fich<-tclvalue(tkgetOpenFile(title="Modflow Binary File",filetypes="{{hds
binary Files} {.hds}} {{All files} *}"))
zz <- file(fich, "rb")

# Cycling thru time steps and layers
for (k in 1:param[2]) {
  for (i in 1:param[1]) {
    readBin(zz,what="numeric",n=1,size=4) # record marker typical of
fortran access="sequential" in gfortran
    readBin(zz,what="integer",n=2,size=4)->N1
    readBin(zz,what="double",n=2,size=8)->N2
    readChar(zz,16)->txt1
    print(txt1)
    readBin(zz,what="integer",n=3,size=4)->N3
    tnber<-N3[1]*N3[2]
    readBin(zz,what="integer",n=1,size=4) # record marker typical of
fortran access="sequential" in gfortran
    readBin(zz,what=real(),n=tnber,size=4)->N4
    readBin(zz,what="integer",n=2,size=4) # record marker typical of
fortran access="sequential" in gfortran
    print(N4[1:10])


    }

}

close(zz)

On Thu, May 3, 2012 at 1:26 PM, Duncan Murdoch
<murdoch.duncan at gmail.com>wrote:

On 03/05/2012 12:41 PM, kapo coulibaly wrote:

I'm trying to read a binary file created by a fortran code using readBin
and readChar. Everything reads fine (integers and strings) except for
double precision numbers, they are read as huge or very small number
(1E-250,...). I tried various endianness, swap, But nothing has worked so
far.
I also tried on R 64 bit for linux and windows (R 2.14) and R 2.11 on
windows XP 32 bit.
Any help would be appreciated.

As I wrote to someone else with a similar problem a couple of weeks ago:

You need to see what's in the file.  The hexView package can dump it in
various formats; see example(viewFormat) for a couple.

Duncan Murdoch

	[[alternative HTML version deleted]]