[输入输出] 复数数据读取

PROGRAM FMMM
 IMPLICIT NONE
 INTEGER I,J,K,ERROR,BBB,B,AA,BB
 REAL P0(360000),P1(360000),P2(360000),P3(360000),P4(360000),P5(360000),P6(360000),LEG0(5000,72),LEG1(5000,72),LEG2(5000,72),LEG3(5000,72),LEG4(5000,72),LEG5(5000,72),LEG6(5000,72)
 COMPLEX HANKEL0(5000),HANKEL1(5000),H(7,5000),HH(7,5000)
     OPEN(UNIT=30,FILE="HHKR_0.TXT")
     READ(30,*) HANKEL0
     OPEN(UNIT=66,FILE="P0.TXT")
     READ(66,*) P0
     OPEN(UNIT=68,FILE="P1.TXT")
     READ(68,*) P1
     OPEN(UNIT=70,FILE="P2.TXT")
     READ(70,*) P2
     OPEN(UNIT=72,FILE="P3.TXT")
     READ(72,*) P3
     OPEN(UNIT=76,FILE="P4.TXT")
     READ(76,*) P4
     OPEN(UNIT=78,FILE="P5.TXT")
     READ(78,*) P5
     OPEN(UNIT=82,FILE="P6.TXT")
     READ(82,*) P6
     DO AA=1,5000
       DO BB=1,72
           LEG0(AA,BB)=P0(BB+(AA-1)*72)
           LEG1(AA,BB)=P1(BB+(AA-1)*72)
           LEG2(AA,BB)=P2(BB+(AA-1)*72)
           LEG3(AA,BB)=P3(BB+(AA-1)*72)
           LEG4(AA,BB)=P4(BB+(AA-1)*72)
           LEG5(AA,BB)=P5(BB+(AA-1)*72)
           LEG6(AA,BB)=P6(BB+(AA-1)*72)
       ENDDO
     ENDDO  
     OPEN(UNIT=90,FILE="L6.TXT")
     WRITE(90,*) LEG6
 PAUSE  
 STOP   
 END

real :: ra , rb
character(len=1) :: sign
read( 30 , "(f17.15,1x,a1,1x,f17.15)" ) ra , sign , rb
if(sign=='-') rb = 0 - rb
hankel = cmplx(ra,rb)