Issues with porting Fortran code

WD_star

I have a Fortran (F77) code which I used to use at my previous institution that I am having trouble with at my new job. I used to compile it with ifort on computers running Fedora (14 if I recall correctly). At my new job they don't have ifort so I have had to revert to the use of g77 (f77) on systems running Red Hat.

It took awhile, and a little help, to figure out that one of my codes would compile but crash once executed because g77 didn't automatically initialize all variables to zero. Luckily switches do exist to make g77 do this and I thought I was home free. This week, I tried running what is essentially the same code, but I get an ''illegal unit number error" when opening a file. I have Googled extensively to both understand the problem and find a solution but have ultimately been unsuccessful. I know illegal unit numbers can come from choosing things out of the range of 0 to 99, but in this case the unit number in question is 87!! Furthermore, the other code I had trouble with previously includes the exact same subroutine where that file is opened and it doesn't cause a problem!!!

Any and all suggestions are useful. Please let me know if I can supply any additional information that might help pinpoint the problem.

Thanks.

WD_star

It's a simple open statement:

Code:

      open(87,file='data/fort_HALPHA.22')

And yes I have tried adding a status='old' statement but I get the same error which is as follows:

open: illegal unit number
apparent state: unit 87 (unnamed)
last format: list io
lately reading sequential formatted external IO
Abort (core dumped)

gsal

Maybe the problem is not when opening the file...looks like that already worked...maybe the problem is at some reading statement somewhrw else where somehow the unit number is not correct. Is the code long? May want to post it

WD_star

The whole code is quite long, here is the portion where the file is read(it coincides with the end of the subroutine):

Code:

      nxbal=0
      read(87,*) nxbal,stnnec,stnchc,vneuc,vchac
      do i=1,nxbal
         read(87,*) xlbal(i),(plbal(i,j),j=1,NNMAX)
      end do
      close(87)
      stnnec=10.0**stnnec
      stnchc=10.0**stnchc
      iwarnc=0
      return
      end

There is also a related issue that I didn't mention in my first post because I found a way to solve it even though it shouldn't be a problem at all. There is an input file that the code cannot find (''file not found'' type of error). The name of the input file is read as a string in another input file. However, if I hardcode the filename in the open command, then there is no problem. The relevant code for that follows:

Code:

      
      read(iunit1,100) file
c      open(iunit4,file=file,status='old')
      open(iunit4,file='data/grid_DA_FINAL',status='old')

WD_star

Bump

Anybody have any suggestions??

gsal

You did not post much about the "open" statement for unit 87...maybe that statement is not working out too good...did you verify that it worked? You can use one of the arguments to the OPEN statement and inspect that upon return from the OPEN statment to see if it worked or not.

As far as that other file name being read from a file vs hardcoded....did you verify what you are reading from the file? You need to read the filename and write it back out and see what you are reading, THEN if that does match the actual filename.

Also, if the file you are trying to read is not in the exact same directory where you are running from, I wonder if you can work with relative path or you need full path, etc.

WD_star

gsal, thanks for taking the time to try and help me.

I have actually solved my problem and it wasn't related to Fortran (at least not exactly). I read your last post and started running some tests. In the process, I noticed that at the bottom of the emacs window for the input file where the other input filename was read, it was written 'DOS'. After a quick Google search, I understood that I had somehow inadvertently changed the buffer encoding system of that file from Unix to DOS!! (http://edivad.wordpress.com/2007/04/...nd-vice-versa/) I switched it back to Unix and presto, the code works! Clearly Fortran could not read the file with the DOS buffer encoding.

In any case, thanks again gsal!

Case closed!

Similar discussions for: Issues with porting Fortran code
Thread	Forum	Replies
fortran code	Engineering, Comp Sci, & Technology Homework	1
Having some issues with a piece of fortran95 code	Programming & Comp Sci	9
FORTRAN 90 iteration issues	Programming & Comp Sci	1
Help me to write fortran code	Programming & Comp Sci	1
Convert latex code to fortran code?	Math & Science Software	1

WD_star	It's a simple open statement: Code: open(87,file='data/fort_HALPHA.22') And yes I have tried adding a status='old' statement but I get the same error which is as follows: open: illegal unit number apparent state: unit 87 (unnamed) last format: list io lately reading sequential formatted external IO Abort (core dumped)

gsal	Maybe the problem is not when opening the file...looks like that already worked...maybe the problem is at some reading statement somewhrw else where somehow the unit number is not correct. Is the code long? May want to post it