source code unknown compiler
bm304_at_cam.ac.uk
Date: 01/31/05
- Next message: Jean-Baptiste FAURE: "How to speedup program compiled by G95 ?"
- Previous message: David Ham: "Re: handing over a parameter to a program"
- Next in thread: Klaus Wacker: "Re: source code unknown compiler"
- Reply: Klaus Wacker: "Re: source code unknown compiler"
- Maybe reply: robin: "Re: source code unknown compiler"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
Date: Mon, 31 Jan 2005 12:50:48 GMT
I am doing research and found a source code file and the author dosent say what compiler
was used.
i tried to compile the code with g77 and it did not work.
the source code is below.
any help or advice is appreciated.
Bryan
-----------------
C* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* * * * * * * * * * * * * * *
C PROGRAM YVAR. FOR *
C Calculates the Aperture Admittance of an Open-Ended
C Coaxial Line radiating in free space under *
C Variational formulation *
C *
C * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
COMPLEX WK,CEPCI,AREA1,AREA2,AREA3,ADM,ADMN
COMPLEX Ell,F22,F,COEF
OPEN (90, FILE='OUTYn',STATUS='NEW' )
OPEN (91, FILE='OUTYa' , STATUS='NEW' )
WRITE (*, *) 'HOW MUCH IS THE INNER RADIUS IN cm ? '
READ(*, *)AA
WRITE ( * , * ) 'HOW MUCH IS THE OUTER RADIUS IN cm ? '
READ(*,*)BB
WRITE(*,*)'WHAT IS THE DIELECTRIC CONSTANT er '
WRITE (*, *) ' INSIDE THE CABLE = ? 1'
READ (*, *) ER
WRITE(*,*)'HOW MUCH IS THE FREQUENCY in GHz = ? '
READ (*, *) FREQ
PI = 4. * ATAN(1.)
A = AA / 100.
B = BB / 100.
OMEGA = 2.*PI*FREQ
WKL = OMEGA*SQRT(PI*4.*ER*0.88542)
CEPCI = CMPLX(l.,0.)
WK = OMEGA*CSQRT(PI*4.*CEPCI*0.88542)
COEF = CMPLX(O.,1.)*CEPCI*WKL/(PI*ER*ALOG(B/A))
YCHARA=2.*PI/(100.*SQRT(PI*4.0/0.88542/ER))/ALOG(B/A)
FRNORM = 2. * PI * FREQ * B * sqrt(ER) / 0.299792
CALL SIMSN3 (WK, A, B, AREA3)
ADMN = COEF*AREA3
ADM = COEF*AREA3*YCHARA
WRITE (*, *) ' freq. in GHz=' , FREQ, ' norm. freq=' , FRNORM
WRITE (*, *) 'YL admittance normalized =' ,ADMN
WRITE(*,*) 'YL admittance =',ADM
WRITE (90, * ) ' freq. in GHz=' , FREQ, ' norm. freq=' , FRNORM
WRITE(90, *) 'YL admittance normalized =',ADMN
WRITE (91, * ) ' freq. in GHz=' , FREQ, ' norm. freq=' , FRNORM
WRITE(91,*) 'YL admittance =',ADM
STOP
END
SUBROUTINE SIMSN3(WK,A,B,AREA3)
COMPLEX AREA3,ENDS,EVEN,0DD,AREA,F22,WK
TEST = 0.0009
LIM = 20
NOI = 0
INT = 1
V = 1.
ODD = CMPLX(O.,O.)
EVEN = CMPLX(O., 0. )
AREA3 = CMPLX(O.,O.)
ENDS = F22 (WK,A,A,B) + F22(WK,B,A,B)
219 H = (B-A) / V
ODD = EVEN + ODD
X = A + (H/2.)
EVEN = CMPLX (0., 0.)
DO 319 I=l,INT
EVEN = EVEN + F22(WK,X,A,B)
X = X + H
319 CONTINUE
AREA = ( ENDS + (4.*EVEN) + (2.*ODD) ) * H / 6.
NOI = NOI + 1
R = CABS ((AREA3-AREA)/AREA)
IF (NOI-LIM) 3119,3219,3219
3119 IF (R-TEST) 32l9,32l9,4l9
3219 RETURN
419 AREA3 = AREA
INT = 2*1NT
V = 2.*V
GO TO 219
END
COMPLEX FUNCTION F22(WK,RHO,A,B)
COMPLEX WK,AREA2
CALL SIMSN2(WK,A,B,RHO,AREA2)
F22 = AREA2
RETURN
END
SUBROUTINE SIMSN2(WK,A,B,RHO,AREA2 )
COMPLEX AREA2,ENDS,EVEN,ODD,AREA,F11,WK
TEST = 0.0009
LIM = 20
NOI = 0
INT = 1
V = 1.
ODD = CMPLX(O.,O.)
EVEN = CMPLX (0 ., 0. )
AREA2 = CMPLX(O.,O.)
ENDS = Fll(WK,RHO,A) + Fll(WK,RHO,B)
21 H = (B-A) / V
ODD = EVEN + ODD
X = A + (H/2.)
EVEN = CMPLX(0.,0.)
DO 31 I=1, INT
EVEN = EVEN + Fll(WK,RHO,X)
X = X + H
31 CONTINUE
AREA = ( ENDS + (4.*EVEN) + (2.*ODD) ) * H / 6
NOI = NOI + 1
R = CABS((AREA2-AREA)/AREA)
IF (NOI-LIM) 311,321, 321
31l IF (R-TEST) 32l,32l,4l
321 RETURN
41 AREA2 = AREA
INT = 2*INT
V = 2.*V
GO TO 21
END
COMPLEX FUNCTION F11(WK,RHO,RHOP)
COMPLEX WK, AREA1
PI = 4. * ATAN(1.)
CALL SIMSN1 (WK,O.,PI,0.001,20,AREA1,RHO,RH0P)
F11 = AREA1
RETURN
END
SUBROUTINE SIMSN1(WK,X1,XEND,TEST,LIM,AREA,RHO,RHOP)
COMPLEX AREA1,ENDS,EVEN,ODD,AREA,F,WK
NOI = 0
ODD = CMPLX(O.,O.)
INT = 1
V = 1.
EVEN = CMPLX(O.,O.)
AREA1 = CMPLX (0.,0.)
ENDS = F(WK,X1,RHO,RHOP) + F(WK,XEND,RHO,RHOP)
2 H = (XEND-X1) / V
ODD = EVEN + ODD
X = Xl + (H/2.)
EVEN = CMPLX(0. , 0. )
DO 3 I=l,INT
EVEN = EVEN + F(WK,X,RHO,RHOP)
X = X + H
3 CONTINUE
AREA = ( ENDS + (4.*EVEN) + (2.*ODD) 1 * H / 6.
NOI = NOI + 1
R = CABS((AREA1-AREA) /AREA)
IF (NOI-LIM) 31,32,32
31 IF (R-TEST) 32,32,4
32 RHOMl = RHO + RHOP
RHOM = 4.*RHO*RHOP/(RHOMl**2)
CALL ELPINT(PILP,RHOM)
AREA = (2.*PILP/(RHO+RHOP)) + AREA
RETURN
4 AREA1 = AREA
INT = 2*INT
V = 2.*V
GO TO 2
END
COMPLEX FUNCTION F(WK,X,RHO,RHOP)
COMPLEX ARGU, WK
R = SQRT(RH0**2+RHOP**2-2.*RHO*RHOP*COS(X))
ARGU = WK*R*CMPLX(O.,-1.)
IF ((RHO .EQ. RHOP) .AND. (X.EQ. 0. )) F=WK,CMPLX (O.,-1.)
IF ((RHO. EQ. RHOP) .AND. (X .EQ. 0.)) RETURN
F = (COS(X)*CEXP(ARGU)-l.)/R
RETURN
END
SUBROUTINE ELPINT(ELP0,ELPC)
ELPI = 1. - ELPC
A0 = 1.3862944
A1 = 0.096663443
A2 = 0.035900924
A3 = 0.037425637
A4 = 0.014511962
B0=0.05
B1=0.12498594
B2=0.068802486
B3=0.0033283553
B4=0.0004178701
AL=((((((((B4*ELPI)+B3)*ELPI)+B2)*ELPI)+B1)*ELPI)+B0)
ELP=((((((((A4*ELPI)+A3)*ELPI)+A2)*ELPI)+A1)*ELPI)+A0)
IF (ELPC .GE. 0.999999) ELPC = 0.999999
ELPI = 1. - ELPC
X1 = 1. \ ELPI
ELPO = ALOG (XI) * AL + ELP
RETURN
END
- Next message: Jean-Baptiste FAURE: "How to speedup program compiled by G95 ?"
- Previous message: David Ham: "Re: handing over a parameter to a program"
- Next in thread: Klaus Wacker: "Re: source code unknown compiler"
- Reply: Klaus Wacker: "Re: source code unknown compiler"
- Maybe reply: robin: "Re: source code unknown compiler"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
Relevant Pages
|
