http://www.nnseek.com/e/comp.lang.fortran/ Posts for comp.lang.fortran Wed, 23 Jul 2008 03:40:49 PDT http://www.nnseek.com/ http://www.nnseek.com/img/64.png 64 64 NNSeek http://www.nnseek.com/e/comp.lang.fortran/read_a_header_and_then_a_binary_file_109393410t.html http://www.nnseek.com/e/comp.lang.fortran/read_a_header_and_then_a_binary_file_109393410t.html
Im new to fortran and Im trying and (failing) to solve my little
problem. Any help would be much appreciated?!

Im trying to read a text header that I have appended to the top of a
binary file and then read the binary file

ie.

1
360
720
0
Blah blah blah

where blah is the start of the binary file.

I can't seem to work out how to read the binary part - i.e. skip over
the header bit. So far I have

PROGRAM read_binary
IMPLICIT NONE

INTEGER :: num_rows , num_cols , num_frames , sizeofbyte = 1
INTEGER, ALLOCATABLE, DIMENSION (:) :: pft
INTEGER :: junk(10)
INTEGER :: i, stdin = 10, hdr_length = 4
INTEGER :: ioerr = 0

OPEN ( stdin, FILE = '/dev/stdin', STATUS = 'old', ACCESS = 'stream',
&
FORM = 'FORMATTED', IOSTAT = ioerr )
IF ( ioerr /= 0 ) THEN
PRINT *, "Error opening file"
STOP
END IF

DO i = 1, hdr_length
READ ( stdin, *, iostat = ioerr ) hdr_stuff(i)
IF ( ioerr /= 0 ) THEN
PRINT *, "Error reading file header"
STOP
END IF
END DO

num_frames = hdr_stuff(1)
num_rows = hdr_stuff(2)
num_cols = hdr_stuff(3)

ALLOCATE ( pft(num_rows * num_cols * num_frames) )

PRINT *, num_frames
PRINT *, num_rows
PRINT *, num_cols

CLOSE ( stdin, IOSTAT = ioerr )
IF ( ioerr /= 0 ) THEN
PRINT *, "Error closing file"
STOP
END IF

OPEN ( stdin, file = '/dev/stdin', form = 'unformatted', access =
'direct', &
recl = sizeofbyte, iostat = ioerr )
IF ( ioerr /= 0 ) THEN
PRINT *, "Error opening file for bin"
STOP
END IF



END PROGRAM read_binary


As I said I think that the way I have approached it thus far is ok -
open file read header, close file, reopen file and try and read
binary. Only I have no idea how to read the binary part - I can't seem
to skip over the binary bit!!

Thanks.

Posted In: comp.lang.fortran

no comments

Reply

]]> Wed, 23 Jul 2008 03:40:49 PDT http://www.nnseek.com/e/comp.lang.fortran/china_sell_max_87_89_80_95_ltd_puma_ugg_boots_lacoste_109392386t.html http://www.nnseek.com/e/comp.lang.fortran/china_sell_max_87_89_80_95_ltd_puma_ugg_boots_lacoste_109392386t.html www.jumpmanlocker.com.cn
air max 87 89 90 95 ltd timberland jeans ugg boots lacoste sandals
hoodies,
t-shirts, mauri shoes, dsquared , hogan shoes, dunks , red monkey,
polo t-shirts,
evisu jeans, bbc jeans , dior, lv, dg, versace, coach puma shoes, nfl
jerseys
shox r2 r3 r4 r5 r6 tn tl1 tl3, sandals, nhl jerseys, mlb jerseys,
nba
jerseys
probowl jerseys,prada shoes,kobe james, hockey jerseys, nfl jerseys,
football
jerseys, baseball jerseys, jordan shoes, jordan fusion air force ones
25 years
basketball jerseys Men's women's shocks OZ NZ TL shoes Discount
Coach
Sandals,
Dior Sandals, Prada Sandals, Chanel Sandals, Versace Sandals, Crocs
Sandals,
Women's Sandals Men's Slippers From China jordan shoes, jordan fusion
air force ones www.jumpmanlocker.com.cn
Affliction T-shirts lacoste T-shirts Polo T-shirts Brand ShirtsGGG T-
www.jumpmanlocker.com.cn
shirts Designer T-Shirts Helen Coat burberry coat Jacket Juicy
www.jumpmanlocker.com.cn
Couture bbc hoodies bape hoodies Designer Hoodies NFL NHL NBA
www.jumpmanlocker.com.cn
MLB Jersey Lacoste Trainers, Prada Sneakers,

Posted In: comp.lang.fortran

no comments

Reply

]]> Wed, 23 Jul 2008 02:27:04 PDT http://www.nnseek.com/e/comp.lang.fortran/converting_fixed_form_f77_into_free_form_f90_109389058t.html http://www.nnseek.com/e/comp.lang.fortran/converting_fixed_form_f77_into_free_form_f90_109389058t.html
Question sounds maybe familiar, but can anyone help me to find a proper working download for
converting fixed form (F77) files into free form (F90) files.

Somewhere I found a file convert.f90 that does the job for you, but I am looking for a ready-to-use
download (binary) under Windows or Linux.

Some help would be nice.

Thanks a lot !

Regards,
P.

Posted In: comp.lang.fortran

2 Comments

Reply

]]> Wed, 23 Jul 2008 00:42:36 PDT http://www.nnseek.com/e/comp.lang.fortran/calling_c_from_fortran_109386498t.html http://www.nnseek.com/e/comp.lang.fortran/calling_c_from_fortran_109386498t.html

main()
{
ofstream myfile;
/********* PART 1: Checking machine bit *********/

FILE *stream;
char sys[6],*sysptr;
int bit;
sysptr = &sys[0];
sys[6] = 0x0;
stream=popen("/bin/uname -m","r");
fread(sysptr,1,6,stream);
pclose(stream);
//printf("%%s\n",sysptr);

if ((strncmp(sysptr, "i686",4)==0)||
(strncmp(sysptr, "i386",4)==0))
{
bit=32;
}
else {
bit=64;
}
myfile.open ("machine");
myfile << bit<<"\n";
myfile.close();
printf("This is a %%d bit machine \n",bit);
}

wose output is written in "machine" file. a f90 program is supposed to
read it as



open(2,file="machine",status="old")
read(2,*) bit
close(2)

I want to use the C code as a function of the f90 program.
can anybody tell me how can i do that?

Posted In: comp.lang.fortran

2 Comments

Reply

]]> Tue, 22 Jul 2008 22:59:49 PDT http://www.nnseek.com/e/comp.lang.fortran/segmentation_fault_in_calling_fftw_109385986t.html http://www.nnseek.com/e/comp.lang.fortran/segmentation_fault_in_calling_fftw_109385986t.html I have written a simple program and called FFTW in it.
but I envisaged Segmentation Fault after "CALL dfftw_execute(plan)"
Is there anyone who can guide me where have I made a mistake?
------------------------
program test
implicit none
include "fftw3.f"
integer::u,jatom,nat,h,functn,pp,ind1,ind2,ind3,nx,ny,nz,counter
real::i,j,k,rx,ry,rz,originx,originy,originz,te,tee
real,allocatable::struct(:,:),tau(:,:),r0(:,:)
complex,dimension(100,100,100)::dat,der
INTEGER*8 plan
read(88,*)nat
allocate(struct(nat,6),tau(4,3),r0(nat,3))
do pp=1,4
read(55,*)(tau(pp,h),h=1,3)
end do
originx=tau(3,1)
originy=tau(3,2)
originz=tau(3,3)
do jatom=1,nat
r0(jatom,1)=tau(jatom,1)-originx
r0(jatom,2)=tau(jatom,2)-originy
r0(jatom,3)=tau(jatom,3)-originz
read(666,*)(struct(jatom,h),h=1,6)
end do
nx=100;ny=100;nz=100
counter=1
i=0;j=0;k=0
do ind1=1,nx
i=i+0.1065
do ind2=1,ny
j=j+0.1065
do ind3=1,nz
k=k+0.1065
rx=i-originx
ry=j-originy
rz=k-originz
u=0
do jatom=1,nat
tee=struct(jatom,4) !*0.529177249
te=sqrt((rx-r0(jatom,1))**2+(ry-r0(jatom,2))**2+(rz-
r0(jatom,3))**2)
if (te.GT.tee) u=u+1
end do
if (u==3)then
dat(ind1,ind2,ind3)=(1.0,0.0)
else
dat(ind1,ind2,ind3)=(0.0,0.0)
end if
write(36,*)dat(ind1,ind2,ind3)
counter=counter+1
end do
k=0
end do
j=0
end do
CALL
dfftw_plan_dft_3d(plan,nx,ny,nz,dat,der,FFTW_BACKWARD,FFTW_ESTIMATE)
print*,"111"
CALL dfftw_execute(plan)
print*,"222"
CALL dfftw_destroy_plan(plan)
print*,"END"
end program test
------------------------------------
the program can print "111" , but it isn't able to print "222".

I'm anxiously looking forward your reply and guidelines.
with best regards,
Fatemeh

Posted In: comp.lang.fortran

1 Comment

Reply

]]> Tue, 22 Jul 2008 22:51:02 PDT http://www.nnseek.com/e/comp.lang.fortran/c_fortran_writing_big_binary_files_issue_109375746t.html http://www.nnseek.com/e/comp.lang.fortran/c_fortran_writing_big_binary_files_issue_109375746t.html

> What's going on ?
>
> What am I doing wrong in the fortran version ?
>
>
> TIA.
>
>
> Cheers,

Did you make the data files in either of the above syntaxes? If so, if you
posted that, others could replicate.

--
The cynics are right nine times out of ten.
H. L. Mencken

Posted In: comp.lang.fortran

no comments

Reply

]]> Tue, 22 Jul 2008 18:24:29 PDT http://www.nnseek.com/e/comp.lang.fortran/c_fortran_writing_big_binary_files_issue_109357826t.html http://www.nnseek.com/e/comp.lang.fortran/c_fortran_writing_big_binary_files_issue_109357826t.html wrote:


> !!! intel compilo
> open(unit=20, file='input.sep', form='binary')
> open(unit=21, file='a.sep', form='binary')
>
> !!! GNU compilo
>!!$ open(unit=20, file='input.sep', form='unformatted', access='stream')
>!!$ open(unit=21, file='a.sep', form='unformatted', access='stream')

Intel Fortran supports access='stream' - unless you are using a version
several years old.
--
Steve Lionel
Developer Products Division
Intel Corporation
Nashua, NH

For email address, replace "invalid" with "com"

User communities for Intel Software Development Products
http://softwareforums.intel.com/
Intel Fortran Support
http://support.intel.com/support/performancetools/fortran
My Fortran blog
http://www.intel.com/software/drfortran

Posted In: comp.lang.fortran

1 Comment

Reply

]]> Tue, 22 Jul 2008 10:50:25 PDT http://www.nnseek.com/e/comp.lang.fortran/c_fortran_writing_big_binary_files_issue_109355010t.html http://www.nnseek.com/e/comp.lang.fortran/c_fortran_writing_big_binary_files_issue_109355010t.html
wrote in <87tzeiuedp.fsf@free.fr>:

> Using this code as an example:

> int main(int argc, char *argv[])
> {
> int i, j, k;
> int nx=600, ny=800, nz=720;
> float x;
> FILE *file_in, *file_out;

> file_in = fopen("input.sep", "rb");
> file_out = fopen("a.sep", "wb");

> for (i=0; i> for (j=0; j> for (k=0; k> {
> fread(&x, 1, sizeof(float), file_in);
> fwrite(&x, 1, sizeof(float), file_out);
> }

> fclose(file_out);
> fclose(file_in);

> return(0);
> }

> takes 1 mn to dump the file.

I'm surprised, it looks horribly inefficient.

> Using this code to do the same thing:

> integer :: i, j, k
> integer :: nx, ny, nz
> real :: x

> nx = 600; ny = 800; nz = 720

> !!! intel compilo
> open(unit=20, file='input.sep', form='binary')
> open(unit=21, file='a.sep', form='binary')

> do i = 1, nx
> do j = 1, ny
> do k = 1, nz
> read(20) x
> write(21) x
> end do
> end do
> end do

> close(20)
> close(21)

> is... endless :-(

> What's going on ?

I'd imagine there's something different in the buffering.

> What am I doing wrong in the fortran version ?

Have you tried something like:

integer :: i, j, k
integer :: nx, ny, nz
real :: x(720)

nx = 600; ny = 800; nz = 720

open(unit=20, file='input.sep', form='binary')
open(unit=21, file='a.sep', form='binary')

do i = 1, nx
do j = 1, ny
read(20) x
write(21) x
end do
end do

?

--
Ivan Reid, School of Engineering & Design, _____________ CMS Collaboration,
Brunel University. Ivan.Reid@[brunel.ac.uk|cern.ch] Room 40-1-B12, CERN
KotPT -- "for stupidity above and beyond the call of duty".

Posted In: comp.lang.fortran

1 Comment

Reply

]]> Tue, 22 Jul 2008 09:52:00 PDT http://www.nnseek.com/e/comp.lang.fortran/f2003_and_f2008_in_windows_system_are_only_dreams_far_109316610t.html http://www.nnseek.com/e/comp.lang.fortran/f2003_and_f2008_in_windows_system_are_only_dreams_far_109316610t.html F2008, can be used after 5 years?

Posted In: comp.lang.fortran

2 Comments

Reply

]]> Tue, 22 Jul 2008 07:21:24 PDT http://www.nnseek.com/e/comp.lang.fortran/netcdf_fortran_and_cygwin_109291778t.html http://www.nnseek.com/e/comp.lang.fortran/netcdf_fortran_and_cygwin_109291778t.html binaries to work with gfortran. Feedback led me to install Cygwin and
try to get NetCDF to work there. There are still problems. I have some
tentative ideas about why, and wonder if anyone has any comments.
1) There are no native Windows NetCDF binaries for Fortran 9x, like
gfortran;
2) There are no Cygwin NetCDF binaries for 32 bit systems, like my 5
year old Pentium 4.
My options, therefore, are:
a) Try to get a (free) f77 compiler (f2c?);
b) Build a 32 bit NetCDF under Cygwin myself (gulp);
c) Buy a new 64 bit computer (double gulp).
I have contacted Unidata (NetCDF), but have yet to hear back. I would
appreciate feedback,especially about the complexity of option b). I am a
Meteorologist, not a computer specialist, and have never built anything.
Lou

Posted In: comp.lang.fortran

4 Comments

Reply

]]> Mon, 21 Jul 2008 19:27:17 PDT