> gfortran -c makeplot.f90
makeplot.f90: In function 'MAIN__':
makeplot.f90:103: internal compiler error: in expand_expr_addr_expr_1, at
expr.c:6517
Please submit a full bug report,
with preprocessed source if appropriate.
______________
gfortran -v
Using built-in specs.
Target: i386-apple-darwin8
Configured with: ../gcc-4.2-20060610/configure --prefix=/sw
--prefix=/sw/lib/gcc4 --enable-languages=c,c++,fortran,objc,java --with-gmp=/sw
--with-included-gettext --host=i386-apple-darwin8
--with-as=/sw/lib/odcctools/bin/as --with-ld=/sw/lib/odcctools/bin/ld
--with-nm=/sw/lib/odcctools/bin/nm --mandir=/sw/lib/gcc4/share/man
--infodir=/sw/lib/gcc4/share/info --disable-multilib
Thread model: posix
gcc version 4.2.0 20060610 (experimental)
___________________________________
PROGRAM makeplot
!
! Take a lightcurve qdp file and try to configure a plot suitable
! for the WWW pages to be set up by Wayne
!
IMPLICIT NONE
!
INTEGER, PARAMETER :: max_val=4096, max_timescales = 10
CHARACTER*128 :: cmdline
CHARACTER*64 :: filroot,outfil, label
INTEGER :: len
REAL*8, DIMENSION(max_val,max_timescales) :: day_no, day_err, rate,
rate_err, signif, weights
INTEGER :: nn(max_timescales)
INTEGER :: nts
CHARACTER*(128) :: line
INTEGER :: ios
INTEGER :: i, j, k, n
LOGICAL :: ok1, ok2
INTEGER :: j_use, count
REAL :: max_y
REAL :: wmean
!
!---
!
day_err = 1.0
rate_err = 1.0
CALL GETCL( cmdline )
len = LEN_TRIM(cmdline)
OPEN (1,FILE=cmdline,STATUS='old')
DO i = 1,3
READ (1,'(A)') line
END DO
j = 1
ok2 = .TRUE.
DO WHILE ( (j <= max_timescales) .AND. (ok2) )
i = 1
ok1 = .TRUE.
DO WHILE ( ( i <= max_val ) .AND. (ok1) )
READ (1,'(A)', IOSTAT = ios ) line
IF (ios == 0 ) THEN
IF (INDEX(LINE, ' no ') > 0 ) THEN
ok1 = .FALSE.
! write (*,*) '1',trim(line)
ELSE
READ (line, *, IOSTAT = ios ) day_no(i,j), day_err(i,j),
rate(i,j), rate_err(i,j)
IF ( ios /= 0 ) THEN
ok1 = .FALSE.
! write (*,*) '2',trim(line)
ELSE
i = i + 1
END IF
END IF
ELSE
ok1 = .false.
IF ( ios < 0 ) ok2 = .false.
END IF
END DO
IF ( i > 1 ) THEN
nn(j) = i - 1
WRITE (*,*) j, ' Read ',nn(j), ' entries'
nts = j
j = j + 1
END IF
END DO
IF ( j > max_timescales ) THEN
WRITE (*,*) ' Buffer overflow - too many timescales ', j
STOP
END IF
signif = rate / rate_err
weights = rate_err**-2.0
! find shortest timescale on which we have at least 4 points at 2 sigma
j_use = 5
DO j = nts, 1, -1
count = 0
n = nn(j)
DO i = 1, n
IF (signif(i,j) > 2.0 ) count = count + 1
write (*,*) rate(i,j), rate_err(i,j), rate(i,j)/rate_err(i,j) ,
signif(i,j), count
END DO
write (*,*) j,count
IF (count >=4 ) j_use = j
END DO
WRITE (*,*) 'Use ',j_use
i = INDEX(cmdline,'.qdp')
filroot = cmdline(1:i-1)
DO j = 1, nts
n = nn(j)
wmean = ( SUM( rate(1:n,j) * weights(1:n,j) ) / SUM(weights(1:n,j) ) )
WRITE (label,'(I4,A)') 2**2.0
OPEN (17, FILE='hdr.pco', STATUS='unknown')
OPEN (18, FILE='qdp.mod', STATUS='unknown')
WRITE (18,*) 'CONS'
WRITE (18,*) wmean, ' 0.0 0.0 0 ! 1 CO'
WRITE (17,*) 'SKIP SING'
WRITE (17,*) 'TIME OFF'
WRITE (17,'(2A)') 'MODEL @qdp.mod'
DO k = 1, nts
IF ( k /= j) WRITE (17,*) 'COL OFF ',i
END DO
WRITE (17,*) 'R X 1440 1750'
WRITE (17,*) 'LOC 0 0 1 1'
WRITE (17,*) 'LAB 1 LIN 90 0.0199999996 CEN Bot'
WRITE (17,*) 'LAB 1 POS 1461 0.00200000009 "Jan-1"'
WRITE (17,*) 'LAB 2 LIN 90 0.0199999996 CEN Bot'
WRITE (17,*) 'LAB 2 POS 1551 0.00200000009 "Apr-1"'
WRITE (17,*) 'LAB 3 LIN 90 0.0199999996 CEN Bot'
WRITE (17,*) 'LAB 3 POS 1642 0.00200000009 "Jul-1"'
WRITE (17,*) 'LAB 4 LIN 90 0.0199999996 CEN Bot'
WRITE (17,*) 'LAB 4 POS 1734 0.00200000009 "Oct-1"'
WRITE (17,*) 'LAB 5 LIN 0 0.0500000007 JUS Lef'
WRITE (17,*) 'LAB 5 POS 1730 0.00100000005 " 1 mCrab"'
WRITE (17,*) 'LAB F ',label, ' day bins'
WRITE (17,*) 'LAB T ', filroot
WRITE (17,*) 'LAB X Swift Days'
WRITE (17,*) 'LAB Y counts/s'
WRITE (17,*) 'R Y -0.0004 0.00300000003'
CLOSE (17)
CLOSE(18)
CALL SYSTEM ('qdp '//TRIM(cmdline)//' < qdp.in')
CALL SYSTEM ('mv pgplot.ps '//TRIM(filroot)//TRIM(label)//'d.ps')
END DO
END PROGRAM makeplot
--
Summary: gfortran crashes and asks for bug report
Product: gcc
Version: 4.2.0
Status: UNCONFIRMED
Severity: normal
Priority: P3
Component: fortran
AssignedTo: unassigned at gcc dot gnu dot org
ReportedBy: skinner at milkyway dot gsfc dot nasa dot gov
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=34837
