BLOCK DATA
C
C
C   BLOCK DATA AND
C   MAIN ROUTINE FOR STEADY STATE EDDY CURRENT PROGRAM
C
C   CONTAINS DIMENSIONS OF ALL COMMON BLOCKS
C
      DOUBLE PRECISION RS,T,WC,W,R
      INTEGER IRN,ITEST,NELD,MATE,IRNU
      COMMON/GENER/MGEN,ISOLN,NV,IVAL(50),AVAL(50)
      COMMON/INPUT$/IER,ICONT,NCHMAX,ICP,IPRINT,KLP,IGOTO
      COMMON/FEXT/HXEXT,HYEXT,FREQ
C   COMMON/INT/ CONTAINS SPACE FOR MATRIX AND LU DECOMPOSITION
      COMMON/INT/RS(6000,6,2)
C   COMMON/INT1/ CONTAINS SPACE FOR ROW NUMBER TABLES OF
C                MATRIX AND LU DECOMPOSITION
      COMMON/INT1/IRN(6000,6)
C   COMMON/INT2/ CONTAINS SPACE FOR RHS VECTOR/ SOLUTION VECTOR
      COMMON/INT2/T(6000,2)
C   COMMON/INT3/ CONTAINS SPACE FOR DIAGONAL POINTERS OF
C                MATRIX AND LU DECOMPOSITION
      COMMON/INT3/IP(6000)
C   COMMON/INT4/ CONTAINS REAL WORKSPACE FOR SOLUTION
      COMMON/INT4/R(6000,4,2)
C   COMMON/INT5/ CONTAINS REAL WORKSPACE FOR SOLUTION
      COMMON/INT5/W(6000,4,2)
C   COMMON/INT6/ CONTAINS COMPLEX WORKSPACE FOR SOLUTION
      COMMON/INT7/WC(6000,2)
      COMMON/MESHX/X(6000)
      COMMON/MESHY/Y(6000)
      COMMON/NODES/NNDS,NELD(6,6000)
      COMMON/TEST/ITEST(6000)
      COMMON/PCONST/IR,NEQ,NTAB(6000)
      COMMON/TEMPQ1/TEMPS(12000)
      COMMON/REGEL/NREGEL(57)
      COMMON/ELEMXK/XM(6000)
      COMMON/ELEMYK/SIGMA(6000)
      COMMON/ELEMP/PHASE(6000)
      COMMON/ELEMS/EK(6000)
      COMMON/ELEM/NEL,MATE(6000)
      COMMON/COILS/NC,NCMAX,PCOIL(60,20)
      COMMON/REGION/NR,NRMAX,PREG(60,75)
      COMMON/SOLUT/U(6000),UINV(6000)
      COMMON/MAXSIZ/NELMAX,NODMAX,NERMAX,NORMAX,NEQMAX,NIRMAX
      COMMON/CNVT/NVT
      COMMON/BOUND/UVAL(400),IUPT(400),NDU
C
C->   DOUBLE PRECISION RS,T,WC,W,R
      DIMENSION WORKIN(12000)
      EQUIVALENCE (RS(1),WORKIN(1))
C
C->   INTEGER IRN,ITEST,NELD,MATE,IRNU
C
      DATA NODMAX/6000/,NCMAX/20/,NRMAX/75/,ISOLN/
0/
      DATA KLP/6/,IPRINT/1/
C
      END
      PROGRAM PE2D AC
C
C   MAIN ROUTINE
C
      COMMON/GENER/MGEN,ISOLN,NV,IVAL(50),AVAL(50)
      COMMON/COILS/NC, NCMAX, PCOIL(60,20)
      COMMON/REGION/ NR, NRMAX, PREG(60,75)
      COMMON/FEXT/HXEXT,HYEXT,FREQ
      COMMON/RHSE/MAXW,ICNT,WORD(1)
      COMMON/GAUSS/WGG(80),NG
      COMMON/INT/RSC(1)
      COMMON/INT2/T(1)
      COMMON/INT3/IP(1)
      COMMON/PCONST/IR,NEQ,NTAB(1)
      COMMON/MAXSIZ/NELMAX,NODMAX,NERMAX,NORMAX,NEQMAX,NIRMAX
      COMMON/CNVT/NVT
      DOUBLE PRECISION RSC,T
      DOUBLE PRECISION T1,T2,CPULFT
      INTEGER*4 IFD3,IFD4,IFD1,IFD2
      DATA IFD3,IFD4/'#FRE','QUEN'/
C
      NG=7
      NELMAX=NODMAX
      MAXW=0
      NIRMAX=NODMAX*6
      IR=NIRMAX
      NEQMAX=NODMAX-1
      NERMAX=NCMAX
      NORMAX=NRMAX
C
      CALL LOAD(1)
C
      FREQ1=FREQ
      IO=1
      IF(NV.EQ.6)IO=2
C
    1 READ(5,500)IFD1,IFD2,NF
  500 FORMAT(2A4,I2)
      IF(IFD1.NE.IFD3.OR.IFD2.NE.IFD4)GO TO 1
C
      DO 7 IF=1,NF
      IF(NF.LE.1)GO TO 2
C
      READ(5,501)FREQ1
  501 FORMAT(G10.3)
C
    2 IF(IF.GT.1)GO TO 3
      FREQ=FREQ1
      CALL SETUP
C
      NEQ2=NEQ*2
      IR2=IR*2
C
      WRITE(11)(T(I),I=1,NEQ2)
      REWIND 11
      GO TO 5
C
    3 FF=FREQ1/FREQ
      DO 4 I=2,IR2,2
      RSC(I)=RSC(I)*FF
    4 CONTINUE
      FREQ=FREQ1
C
    5 IFLAG=1
      IF(IF.EQ.1.OR.FF.EQ.0.0)IFLAG=0
      CALL SOLVE(IFLAG)
C
    6 ISOLN=1
      IVAL(4)=1
      CALL UNITS(IF)
      CALL STORE(IF)
C
      READ(11)(T(I),I=1,NEQ2)
      REWIND 11
C
    7 CONTINUE
C
      STOP
      END