File: DATA33.FC of Tape: Various/ETH/fc2

(.FC,.FO Focal format converted to 'WRITE' listing) 

---

C-DATA  .TM OCT-74

01.01 C DATA33.FC
01.02 C 10 DEC 1974
01.04 C FUER FILES AB M30 BIS M33
01.05 C FUER DATEN BIS M29: DKB1:DATA29.FC
01.06 T !"CURRENT FILE:   "
01.90 F I=1,6;S Z=FOUT(MD(I))
01.95 Q

02.01 C AUT . PEAKS
02.05 E
02.10 D 3
02.11 S $=2000;S BE=0;S PS=1;F I=1,4;S SW(I)=0
02.13 S SI=0;S X=16;D 5.14;D 5.2;D 5.3;D 18
02.15 F I=1,4;A ?NU(I)?
02.16 G 2.2
02.17 D 2.13
02.20 S K=1;I (-FABS($-211))2.25,2.9
02.25 I (NU(K))2.26,2.17,2.26
02.26 I (SW(K)),2.5,2.3
02.30 S NU=NU(K);S PH=PH(K);S AF=-AF(K);S TA=TA(K)
02.40 D 13.8;D 21.3;D 21.32;S SW(K)=0
02.50 S NU=NU(K);D 10.05;D 9;D 20
02.60 D 13.8;D 21.3;D 21.32;S SW(K)=1
02.70 S NU(K)=NU;S PH(K)=PH;S AF(K)=AF;S TA(K)=TA
02.80 S K=K+1;G 2.25
02.90 F I=1,4;T %5.03,!,NU(I),PH(I),AF(I),TA(I)
02.95 S Z=FIN();D 3;A "'GROUP 8 OK?",Z;I (-FABS(Z-0YES))2.99;D 8
02.99 Q

03.01 C DATAFILE
03.02 T "FILE:";F I=1,6;S MD(I)=FIN()
03.04 O O STEUER;T !"3.1 O I DKA1:";D 1.9
03.05 T !"8.1 O O DKA1:";D 1.9
03.06 T !"G 3.1"!
03.07 O C
03.08 O I STEUER;Q
03.10 O I DKA1:M32D03
03.12 F I=1,6;S MD(I)=FIN()
03.20 A D,N;D 13.8;F I=0,N-1;A X;S Z=FCOM(I,X);D 13.5
03.25 F I=N,255;S Z=FCOM(I,0);D 13.5
03.30 A X;I (X-20)3.35,3.4,3.35
03.35 T "FEHLER"!;D 3.8;Q
03.40 A DK,HR,HI,HF,MI,PS,PT,RT,TM,TV,SH,MP,Z,Z,CH,GD,P0,P1,P2,P3
03.50 O I TTY:,E
03.60 A !"READ IN NUS",Z;I (Z-0YES)3.9,3.65,3.9
03.65 O R I
03.70 F I=1,4;A NU(I),PH(I),AF(I),TA(I)
03.80 O I TTY:,E
03.90 T "   ";D 1.9;T " IN"!;Q

05.01 C FOUR
05.05 A ?SI,SF,?
05.10 S X=(SF-SI)/16;I (X-FITR(X)-.5)5.14;S X=X+.5
05.14 S SS=FITR(X)/16;S S2=SI+SS*128;S SF=SI+SS*256
05.16 T %3,"'"SS*32"'"
05.20 D 13.9;D 12
05.30 S Z=FOUR(SI,SS);S ZM=FX(1,760)

06.01 C SR(FILTER)
06.10 S Y=G(21)*FCOM(K)
06.20 F J=1,20;S Y=Y+G(21-J)*[FCOM(K-J)+FCOM(K+J)]

07.01 C FILTER
07.05 I (-G(21))7.1;D 29
07.10 O I M24D19
07.15 A D,D,N
07.20 D 13.8;F I=236,510;A X;S Z=FCOM(I,X);S Z=FX(1,500)
07.30 F I=0,117;S K=256+2*I;D 6;S Z=FCOM(I,Y);D 13.5
07.40 F I=0,38;S Z=FCOM(197+I,FCOM(472+I))
07.50 D 7.2
07.60 F I=0,137;S K=217+2*I;D 6;S Z=FCOM(I+118,Y);D 13.5
07.70 F I=551,N;A X
07.72 A X;I (X-13)7.74,7.76,7.74
07.74 T "FEHLER: X=",%8.04,X;G 7.8
07.76 D 3.4
07.80 D 3.8;S N=256;S HI=1/(1/HI+21*DK*1E-7);S DK=2*DK
07.90 O I TTY:,E;Q

08.01 C STORE
08.05 D 1.9;A " OK",Z;I (Z-0YES)8.06,8.1,8.06
08.06 D 3.02;D 3.04;D 3.05;T !"G 8.1"!;O C
08.07 O I STEUER;Q
08.10 O O DKA1:M32D03
08.20 D 1.9;T %2,D,!,%3,N,!;F I=0,N-1;T %5.03,FCOM(I),!
08.30 T "20",!,DK,!,HR,!,HI,!,HF,!,MI,!,PS,!,PT,!
08.40 T RT,!,TM,!,TV,!,SH,!,MP,!;F I=1,8;T 0,!
08.50 O I TTY:,E;O O TTY:
08.52 A " OLD VALUES",Z;I (-FABS(Z-0YES))8.55,8.7
08.55 S I=1
08.60 T %1,"I=",I,"   ";A ?NU(I)?
08.62 I (NU(I))8.64,8.68,8.64
08.64 A ?PH(I) AF(I) TA(I)?!;S I=I+1;I (I-5)8.6,8.7,8.7
08.68 F J=I,4;S NU(J)=0;S PH(J)=0;S AF(J)=0;S TA(J)=0
08.70 T !;O R O
08.80 F I=1,4;T %5.03,NU(I),!,PH(I),!,AF(I),!,TA(I),!
08.85 F I=1,4;T 0,!,0,!,0,!,0,!
08.90 S Z=FOUT(140);O C
08.95 D 1.9;T " OUT"!;Q

09.01 C TAU + NUMAX
09.10 S X=256+(NU-SI)/SS
09.20 I (FCOM(X+1)-FCOM(X))9.25,9.25;S X=X+1;G 9.2
09.25 I (FCOM(X-1)-FCOM(X))9.3,9.3;S X=X-1;G 9.25
09.30 S XM=X;S AH=.5*FCOM(XM);S NU=(XM-256)*SS+SI
09.40 S X=X-1;I (FCOM(X)-AH)9.45,9.45,9.4
09.45 S X1=X;S A1=FCOM(X);S A2=FCOM(X+1);S X=XM+1
09.50 I (FCOM(X)-AH)9.55;S X=X+1;G 9.5
09.55 S X3=X-1;S A3=FCOM(X3);S A4=FCOM(X)
09.60 S X1=X1+(AH-A1)/(A2-A1);S X3=X3+(A3-AH)/(A3-A4)
09.70 S TA=(X3-X1)*SS/2
09.71 S TA=TA*N/256
09.72 I (TA-.95)9.74;S TA=TA^4-.75;S TA=FSQT(FSQT(TA));G 9.8
09.73 I (-TA)9.74,9.74;S TA=.1;G 9.8
09.74 I (TA-.885)9.9;S TA=FSQT((TA-.885)/.28)
09.80 S TA=TA*256/N;R
09.90 T "!TAU!";S TA=0.01

10.01 C PEAKEVAL
10.04 A ?NU?;I (-FABS(NU))10.05;A ?I?;S NU=NU(I)
10.05 S SI=NU-8;S X=1;D 5.14;D 5.2;D 5.3
10.10 D 9;D 20;T %5.03, "NU",NU,"  PH",PH,"  AF",AF,"  TA",TA,!
10.20 S NY=NU;F I=152,155;D I/10
10.30 A !"STORE AS NU(I)? ",?I?;I (-FABS(I))10.4;R
10.40 S NU(I)=NU;S PH(I)=PH;S AF(I)=AF;S TA(I)=TA

12.01 C FOURSCALE
12.10 S X=194;S A=SI/100;S Y=FITR(A)
12.12 I (-Y)12.14;S Z=FDIS(X+2,-96);G 12.16
12.14 S L=12.2;G 12.52
12.16 D 12.8;I (-Y)12.18;S Z=FDIS(X+1,-96);G 12.8
12.18 D 12.9;G 12.8
12.20 S X=201;S A=S2/100;S Y=FITR(A);S L=12.3;G 12.5
12.30 S X=212;S A=SF/100;S Y=FITR(A);S L=12.34;G 12.5
12.34 R
12.50 I (-Y)12.52,12.54
12.52 D 12.9;G 12.56
12.54 S Z=FDIS(X,-96);S X=X+1
12.56 D 12.8;D 12.9;D 12.8;D 12.9;G L
12.80 S A=(A-Y)*10+1E-4;S Y=FITR(A)
12.90 S Z=FDIS(X,Y-80);S X=X+1
12.92 G 12.2

13.01 C DIS SWITCH
13.40 S Z=FDIS(0,0)
13.50 S Z=FX(0,500)
13.60 S Z=FDIS(1,760)
13.70 S Z=FDIS(A,FABS(FDIS(A)));S Z=FDIS(B,-FABS(FDIS(B)))
13.80 S A=4;S B=5;D 13.7
13.90 S A=5;S B=4;D 13.7

15.01 C PERIODS  INCL ALIASES
15.05 A ?DK?,!
15.10 A ?NU?;S NY=NU
15.20 T "   P=",%4.2
15.30 T 512*DK/NY
15.40 S NY=512-NU;D 15.3
15.50 S NY=512+NU;D 15.3
15.90 T !;G 15.1

16.01 C PLOT
16.02 C WITH AXES AND TEXT: PA=0, ELSE PA=2
16.05 S PA=2;A ?PA?
16.10 S PF=4;I (FDIS(5))16.2;S PF=5
16.20 F I=2,19;S Z=FDIS(I,-FABS(FDIS(I)))
16.70 S X=PF+PA;D 16.9;S Z=FDIS(21,1)
16.80 D 16.9;S X=PF;D 16.9;S $=1024;Q
16.90 S Z=FDIS(X,-FDIS(X))

17.10 S Z=FDIS(0,0)

18.01 C SQ.ROOT OF POWERSPECTRUM(FCOM(1))
18.10 F I=256,511;S Z=FCOM(I,FSQT(FCOM(I)));S Z=FX(1,760)

19.10 A !,X;T FCOM(X);G 19.1

20.01 C FOUGNAU
20.10 S N1=NU*PI/256;S A=0;S B=0
20.20 F I=0,N-1;S A=A+FCOM(I)*FCOS(I*N1)
20.25 I (NU*(NU-256))20.3,20.4,20.3
20.30 F I=0,N-1;S B=B+FCOM(I)*FSIN(I*N1)
20.40 S AF=2*FSQT(A^2+B^2)/N
20.50 I (A)20.7,20.6;S PH=0;G 20.8
20.60 S PH=90*FSGN(B);R
20.70 S PH=180*FSGN(B)
20.80 S PH=PH+180*FATN(B/A)/PI
20.95 S AM=AF*TA*PI*N/(256*(1-FEXP(-TA*N*PI/256)))

21.01 C SUB PEAK
21.10 D 13.8
21.15 A ?NU?;I (-FABS(NU))21.2;A ?I?
21.16 S NU=NU(I);S PH=PH(I);S AF=AF(I);S TA=TA(I);G 21.3
21.20 A ?PH AM TA GD MP?
21.30 S NO=NU*PI/256;S PO=PH*PI/180;S TO=TA*PI/256;S A1=HI*DK*1E-7
21.31 S A=AF*TO*N/(1-FEXP(-N*TO));S H6=MP*NO*TM*.0048/A1
21.32 F I=0,N-1;D 21.35;D 21.4;D 13.5
21.33 R
21.35 S J=I*(1-GD*A1*I);S H1=1+J*A1
21.40 S Z=FCOM(I,FCOM(I)-A*FSQT(H1)*FCOS(J*NO-PO)*FEXP(-J*TO)/(1-FEXP(-H1*H6)))

23.01 C SUB LEGENDRE
23.10 A ?P0 P1 P2 P3?
23.15 D 13.8
23.20 F I=0,N-1;S J=I/N;D 23.3;D 23.4;D 23.5
23.25 Q
23.30 S Z3=20*J*J*J-30*J*J+12*J-1
23.40 S Z1=2*J-1;S Z2=6*J*J-6*J+1
23.50 S Z=FCOM(I,FCOM(I)-P0-P1*Z1-P2*Z2-P3*Z3);D 13.5

29.01 C FILTERKOEFF
29.10 O I FILTER
29.20 F I=1,21;A G(I)
29.30 O I TTY:,E

---