Re: How can generate Uniform Cartesian grid?
 

Rt, Have you sent me your simple-version code? I guess you miss ed my last message. My email address is s.b.kuang@gmail.com. Thanks for your help. Cheers, Harry

Re: How can generate Uniform Cartesian grid?
 

I am also interested in the changes done for compiling the code with gfortran. I am struggeling with that.

Re: How can generate Uniform Cartesian grid?
 

For these errors:

1)

IF(FLG_OCTREE) CALL MAKE_REMAINED_LEVELS

1 Error: ELSE IF clause at (1) requires a scalar LOGICAL expression In file CartGen.f90:139 you have to add ".EQ.1" to the if condition, that should be as IF(TRUE) ....

2) DO WHILE(FLG)

1 Error: Exit condition of DO WHILE loop at (1) must be a scalar LOGICAL expression In file CartGen.f90:1298

for this the same as before.

3) delete line 129 and 141 for TIMEF()

4) opening STL file change ACCESS and FORM:

!OPEN(UNIT=1,FILE=FNAME(IFF),ACCESS='DIRECT',FORM= 'BINARY',RECL=2) OPEN(UNIT=1,FILE=FNAME(IFF),ACCESS='STREAM',FORM=' UNFORMATTED',RECL=2)

5) the same for this line:

OPEN(UNIT=1,FILE=FNAME(IFF),ACCESS='STREAM',FORM=' UNFORMATTED',RECL=2) !OPEN(UNIT=1,FILE=FNAME(IFF),ACCESS='DIRECT',FORM= 'BINARY',RECL=2)

6) change line 1707 and 1855 :

!open(unit=ivtk,file=VTKFILE,form='binary',convert ='BIG_ENDIAN') OPEN(unit=ivtk,file=VTKFILE,ACCESS='STREAM',FORM=' UNFORMATTED')

but when I generate the binary VTK file I'm not be able to read it with paraview. I don't know if this is due to this last modification. Reading gfortran documentation CONVERT should be available, but only deleting it I compile the code.

My e-mail is: iota2@libero.it