fortran脚本在Unix下编译出现问题

MacTaylor

我之前有一个脚本，在win下做东西都没问题，在Unxi下编译出现错误提示：
PGF90-S-0201-Illegal I/O specifier - FILE (uq.f: 23)
  0 inform,   0 warnings,   1 severes, 0 fatal for MAIN
在我自己机子上的Linux下的gfortran编译出现以下提示：
open(unit=11,file=filename(M,:),form='unformatted',            
                           1
错误： FILE 标记在(1)处必须是标量

似乎是读入语句和输出语句有问题，是因为recl=4在不同环境下的设定的问题吗？我确认此程序在win下没有问题，请教大家是怎么回事？谢谢！

脚本如下：
        dimension uq(256,128,2400),uqt(256,128,2400,8)
         real lon(256),lat(128)
         integer M,L
         character*7 filename(8,7)
         
         filename(1,:)='uq1.bin'
         filename(2,:)='uq2.bin'
         filename(3,:)='uq3.bin'
         filename(4,:)='uq4.bin'
         filename(5,:)='uq5.bin'
         filename(6,:)='uq6.bin'
         filename(7,:)='uq7.bin'
         filename(8,:)='uq8.bin'

         do i=1,256
          lon(i)=0.0+(i-1)*1.40625
         enddo
         do j=1,128
          lat(j)=-90.0+(j-1)*1.40625
         enddo

         do M=1,8
         open(unit=11,file=filename(M,:),form='unformatted',
     &    access='direct',recl=256*128*2400*4)
         read(11,rec=1)uq
         close(11)
         
         
         do i=1,256
          do j=1,128
           do k=1,2400
            if(abs(uq(i,j,k)).ge.100000) then
             uqt(i,j,k,M)=0
            else
             uqt(i,j,k,M)=uq(i,j,k)
            endif
           enddo
          enddo
         enddo
        enddo

        do i=1,256
         do j=1,128
          do k=1,2400
           s=0
           do L=1,8
            s=s+uqt(i,j,k,L)
           enddo
           uq(i,j,k)=s
          enddo
         enddo
        enddo

        open(unit=12,file='uq.bin',form='unformatted',
     &    access='direct',recl=256*128*2400*4)
        write(12,rec=1)(((uq(i,j,k),i=1,256),j=1,128),k=1,2400)
        close(12)
                 
        end

残月

do M=1,8
         open(unit=11,file=filename(M,:)。。。） 在win下是对的？??
你这一个设备号连接了 n 个文件啊....
话说 filename为什么不定义成filename(8)  ，还要加个 ‘：’？？不是很懂你的意思。。。从程序来看定义为一维数组就可以了。

ogions

定义数组是不是应该定义成： character*7 filename(8)

MacTaylor

ogions 发表于 2012-6-14 11:19

                               
登录/注册后可看大图

定义数组是不是应该定义成： character*7 filename(8)

我也试过这样定义，但是错误信息是：
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_myprocnum':
initpar.c:(.text+0x2): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_lcpu' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_ncpus':
initpar.c:(.text+0x12): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_tcpus' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_getioproc':
initpar.c:(.text+0x22): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_ioproc' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_is_ioproc':
initpar.c:(.text+0x32): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_ioproc' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
initpar.c:(.text+0x38): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_lcpu' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_abort':
initpar.c:(.text+0x5f): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_lcpu' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_abortp':
initpar.c:(.text+0xeb): relocation truncated to fit: R_X86_64_PC32 against symbol `__hpf_lcpu' defined in COMMON section in /panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o)
/panfs/panasas/opt/pgi/linux86-64/9.0-3/lib/libpgf90.a(initpar.o): In function `__hpf_initarg':
initpar.c:(.text+0x127): relocation truncated to fit: R_X86_64_PC32 against `.bss'
initpar.c:(.text+0x151): relocation truncated to fit: R_X86_64_PC32 against `.bss'
initpar.c:(.text+0x17b): relocation truncated to fit: R_X86_64_PC32 against `.bss'
initpar.c:(.text+0x18b): additional relocation overflows omitted from the output
而且定义成二维在win下是调通了的，以前没在Unix下做过for，大家帮帮忙！谢谢！

MacTaylor

残月 发表于 2012-6-14 11:12

                               
登录/注册后可看大图

do M=1,8
         open(unit=11,file=filename(M,:)。。。） 在win下是对的？??
你这一个设备号连接了 ...

数组定义成二维在win是对的，之前一直能运行成功的，而且我也试过filename（M），但是也错误，错误信息请看我下面的回复。

godenflame135

确实由于系统架构不同，fortran使用不同编译器在读取二进制文件时有*4或不*4的问题。pgi就不用，ifort要用。由于此问题的存在，建议使用如下方式读取文件。
program read_uq
   implicit none
!-----------------------------------------------
dimension uq(256,128,2400),uqt(256,128,2400,8)
integer*4       ::   i ,j,k
!------------------------------------------------
...

     do M=1,8
         open(unit=11,file=filename(M,:),form='binary', access='sequential')
         read(11)  ( ( (  uq(i,j,k),i=1,256),j=1,128,),k=1,2400)
         close(11)
    end do
...
stop

end program read_uq

这样就不会因为架构不同，而有不同编辑方法的差别了。

MacTaylor

godenflame135 发表于 2012-6-15 13:02

                               
登录/注册后可看大图

确实由于系统架构不同，fortran使用不同编译器在读取二进制文件时有*4或不*4的问题。pgi就不用，ifort要用。 ...

十分感谢！学习了

zhuxueshuai

我也是，同求解。。。。

MacTaylor

zhuxueshuai 发表于 2013-4-12 19:59

                               
登录/注册后可看大图

我也是，同求解。。。。

试试改成:

character*7 filename(8)*7
         
filename(1)='uq1.bin'
......

zhuxueshuai

MacTaylor 发表于 2013-4-13 04:42

                               
登录/注册后可看大图

试试改成:

character*7 filename(8)*7

好的，谢谢啦，亲