子程序用open和read写入数组出现问题

rachel2002

大家好~有个小问题调了好几天，百思不得其解，希望大家帮帮忙，不胜感激！！
用usdfld子程序定义场变量，需要用到一个81*2维数组，用写入的方法，可是只要一用哪个open，read的就不算了，不知道问题出在哪里。如果不用这个写入的方法，直接用手动输入时可以的（orz）。
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
     5 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
     6 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
      INCLUDE 'ABA_PARAM.INC'
      CHARACTER*80 CMNAME,ORNAME
      CHARACTER*3  FLGRAY(15)
      DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),
     7 T(3,3),TIME(2)
      DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
以上是固定的title
real base(81,2)
real xdis,dis,ydis,A0,nox,oxtime,deltah,toler, deltah0
integer nsmall   

  open(unit=290,file='E:\USER\dhx\aba\20100818\222.TXT')
ydis = 1.0d10
toler = 1.0d-4
A0 = 0.4818
nox = 0.4007

这以下5行

do i=1,81
    read(290,*) (base(i,j),j=1,2)
     write(6,*) (base(i,j),j=1,2)
enddo   
close(290)

以下是判断场变量

if (abs(kstep-1.0) < toler) then
  oxtime = time(1)/3600.d0
endif   

if (abs(kstep-2.0) < toler) then
  oxtime = (time(2)-time(1))/3600.d0
endif   

deltah = (A0*oxtime**nox)/1000.d0
xdis = abs(coord(1) - 1.8465)

if ((xdis .le. deltah) .and. (coord(1).le.1.8465)) then
        field(1) = 1.0
else
        field(1) = 0.0
endif

      RETURN  
END
还请大家帮我看一下，谢谢。。我问了好些人了，都没发现错误，块逼疯了。。

rachel2002

1# rachel2002 下面是附件

cdstudio

你的程序在我机器上运行得完全正常, 不知道是不是你的222.txt的路径是错的.


下次请注意将运行所需要的文件, 包括子程序一起上传, 避免浪费其它人的宝贵时间

rachel2002

你的程序在我机器上运行得完全正常, 不知道是不是你的222.txt的路径是错的.


下次请注意将运行所需要的文件, 包括子程序一起上传, 避免浪费其它人的宝贵时间
cdstudio 发表于 2010-8-23 15:19

谢谢版主~首先非常不好意思，下次我会注意的。
在你的电脑运行正常？文件路径是写绝对路径对吧？你运行的时候也输出那个数组222了么？非常感谢！！

ycroy1986

楼主，我现在读入数据也出现问题，总是跳出内存不能为read，文件路径也不存在问题，不知道是怎么回事。你的问题是怎么解决的啊？

zhltang22

是不是你把文件名222.TXT弄成了222.TXT.txt了呢？你把文件的后缀名也完全显示出来看一下，windows的文件名还是要注意的

zcs220

      CALL GETJOBNAME(JOBNAME,LENJOBNAME)
      CALL GETOUTDIR(OUTDIR,LENOUTDIR)
!     OPEN(UNIT=290,FILE=TRIM(OUTDIR)//'\'//'222.txt')
      OPEN(UNIT=290,FILE=TRIM(OUTDIR)//'\'//TRIM(JOBNAME)//'.txt')

找到当前路径或文件。