Frequentiedemodulatie, analoge numerieke integratie, hybriede banddoorlaatfiltering, meer-kanaals correlatie en energiedichtheidsspectra

(1)

ANALOGE DATA PROCESSING

Frequentiedernodulatie, analoge numeriekeintegratie, _hybriede banddoorlaatfiltering, rneer-kanaa is correlatie en

energiedichtheids spectra.

Ke'

to)

Rekencentrum T.H. Deift

augustus 1973

Afdelin Hybried Rekenen _{Ir. F.J.Pasveer}

(2)

INHOUD

DEEL I

Samenvatting

i

Inleiding

2

Frequentiemodulatie en - deinodulati

hij

analoge instrumentie recorders

₅

Frequentiedemodulaties chake ling

6

Analoge nunierieke integratie

9

Hybriedo banddoorlaatfilter

10

Covariantiefuncties en energiedichtheidsspectra

12

Appendix

14

DEEL II

Bandnelheid bij opname en weergave

Specifieke tijdsintervallen hij de dat.averwerking

18

Dernodulatieschakeling

Timing - en Control circuit

19

Frequcntievrrnenigvuldiging in hot golfversnellingskanaal

24

Organisatie software

Organisatie COrIMON

25

AD4 interrupthcwaking en - afhandeling.

30

Referentiebewaking

(

stthroutine REFBE)

31

Bij lagen:

listings

strooms cherna' s

overdrachtsfunctie bancìdoorlaatfilter

voorbeeld presentatie van de uitvoer

(3)

Samenvatting

In samenwerking met het laborat:oriuri voor Scheepsbouwkunde is op de Centrale Hybriede Rekeninstallatie een procedure ontwikkeld voor de parallel verwerking van recistraties van scheepsbewegingen _{op analoge}

tape.

De dataverwerking ornvat de deeltaken: frequentiedemodulatie, _analoge numerieke integratie, hybriede banddoorlaatfiltering, berekeningen _van covariantiefunctics en energiedichtheidsspectra.

Door een passende taakverdeling over analoga en digitale gedeelte werd een aanzienlijke reductie verkregen ìn de totale verwerkingstijd tussen aanbod van informatie vanaf analoge tape en presentatie van de

eindresultaten.

(4)

-1--Inleiding

Tabel i geeft een overzicht. van de verlangde resultaten _{en specificaties} van de geregistreerde metingen op analoge tape.

De gehele procedure van dataverwerking hestaat uit 3 afzonderlijke

hoofciDrogramlna

Is.

initialisatie

parallel processing

spectrale berekeningen, presentatie van resultaten

Tijdens de initialisatie vindt de invoer van gegevens plaats. De tweede fase omvat de daadwerkelijke data processing wolke o,a. bestaat uit

berekening van auto- en kruiscovariantiefuncties. In de derde _{fase worden} powerspectra uit de covariantiefuncties berekend en worden overdrachts-functies bepaald. De eindresultaten worden in de _{vorm van grafieken met} hijbehorende getaiwaarden op de regeidrukker gepresenteerd.

De parallel processing omvat naast digitale covariantieberekeningen een aantal voornamelijk analoge voorbewerkingen, te weten:

frequentiedemodulatie

analoge numerieke integratie

ci. hybciede banddoorlaatfilteren

ada:

Vanwege de ineonstantheid van de harxdsnelhetd gedurende de _{opnarnen kan} bij de verwerking geen gebruik worden gemaakt van de interne frequentie-dernodulatoren. Daartoe zijn aparte frequentiedemodulatieschakelingen _op de analoge computer geprogrammeerd die gestuurd worden door een afzon-derlijk _{eferentiesignaal van constante frequentie dat tegeiijkertijd} _met de meetregistratics op tape is meegeschreven.

ad b:

Ter verkrijeing van de vereiste verplaatsingssignalen dienen de versnellings-registraties tweemaal geintegreerd te worden. Dit _{integreren vindt eveneens} plaats op de analoge computer volgens _{eon op de rechthoekregel voor}

digi-tale nunierieke integratie ciebaseerd orincipe. De integratieschake lingen staan in serie met de frequentiedemodulatoren van de versnellingsregis-traties.

zie pag. 3

(5)

-2--registratie

Tabel 1: Overzicht veriangde resultaten en relevante specificaties

-) houdt integratie in

)

afgeleid uit 50 punts covariantiefunctie

°) qebaseerci op 4x versnelde weergave

U (t) U) U) (I) 4_J t

-4-i t-1 U

rlN

tie

--ti

oX

U)>

rj4J

ct-i

ow

ti (ti (D t--i o

roihoek

X X

0,4

13,3

5400

-r-4

0-1,2

roerhoek

X X

0,4

13,3

5400

01 , 2

slirigerhoek

X X

0,4

13,3

5400

0-1,2

s tarnphoek

X X X

0,4

13,3

5400

0-1,2

domoversnell.

X X X

0,4

13,3

5400

(riaar ampi.)

-'/c)

-,2

golfversnell.

(naar ampi.)

k) X X X

0,4

9,2

500/

1000

-1 2

(6)

13,3 kc

i,3kc

13, k k

20 kd

k 2 k3

k4

ED.

F.D.

FD.

-4-ad C:

Ter onderdrukking van off set.s en andere ongewenste _{frequentiecomponenten}

tijdans de integratte worden banddoorlaatfilters in _{serie met de} integra--tieschakelingen toeqepast. D _{f ilterberekeningen vidèn plaats in de} digitale computer, daar bet onwaarschijnhijk is dat een dergelijk smal-bandig filter (doorlaatband van 0,2 tot 1,2 Hz) _{eenvoudig met analoge} componenten gerealiseerd kan worden

Een overzicht van de hybriede dataverwerking is in _{f iguur i gegeven.}

'1'

L5

;1J

1333 k

ttiii

& CGflkoL

Fij.1:

bLokd.iaq'ern

£.

vn

_{tLLi!. dk.ve'k3.}

J

L;

mpt.Lrg

i'tuptte ut

hirt

ops1

uLtti

lo +

coviartie

btekcni.

2°

FLteft op

Lntc'vtup.-'t C t

btuLng ckodvLtûn.

I

6aootLaatf1Le.

PD.

_{/tcFcntk DodLaL.}

f

integ

ta

'toce

f'tq.

-dc(e

1Su.

A D C o r V e ¿t e 'i,

ki

1,3 k

(7)

Frequentiemodulatie en -demodulatie bij analoge instrumentatie recorders

Bij frequentiemodulatie wordt _{een frequentieverandering ten opzichte van} een constante frequentie veroorzaakt door de momentane waarde van het te registreren signaal. Dc conotante frecjuentie _{noemt men daaggolffrequentie} en de maximale frequentieverandering de frequentiezwaai. _{Omdat uitsluitend} het aantal nuldoorgangen per tijdseenheid maatgevend is en niet de vorm, wordt eon fmsignaal als blokvorm _{op tape geregistreerd.}

Net fm sigraal kan dan ecnvoudig beschreven _{worden met (f +f (t)) , waarin}

f _{de draaggolffrequentie is}

en tf(t) de, met bet uitsturend. signaal even-redige frequentieverandering. Zie f iguur 2 voor een principeschets.

p1a Sqn.

_-)-

F'tyenttz

-5-pe

0Fra t j

Ikop

F.Z: p'tcp

pramc

'-

ue't.avc.

hçj cer anLje

Voor betrouwbare opname en weergave is _{een constante bandsnelheid vereist.} Tijdcns de opnalnen wordt t.g.v. wisselende _{bandsnelheid een andere} frequen-tie op tape geregistreerd dan die _{door de modulator wordt afgegeven.}

Hierdoor ontstaan bij weergave frequentieafwijkingen die _{d.c demodulator} ten gevolge van fluctuaties in de handsneìheid niet kan _onderscheidn van het qeregistreerde signaai

De inv].oed van wisselende bandsnelheid

kan uitgeschakeld worden door een apart signaal van constante frequontie _{direct op tape mee te schrijven}

niet

de frequentiegemoduleerde _{rieetregi.straties}

Nu is in principe wel een betrouwbare _{frequentiedemodulatie uit te voeren,} omdat de onderlinge tijdsrelatie _{tussen de fm reqistratie en het} referen-tiesignaal niet verandert daar beide op de zelfde wijze door wisselende bandsnelheden worden beïnvloed.

M9

ri

de oduLafo

(8)

5knaaL.

(fc+L)

Frequentiedemodulatieschakeling

De frequentiedemodulatieschakelingen voorzien de digitale computer en de analoge numerieke integratieschakelingen van hemonsteringen van de analoge meetregistraties. In het ene goval worden rechtstreeks covariantiehere-keningen uitgevoerd. In het andere gc-val vindt eerst dubbele integratie plaats van de versnellingsregistraties alvorens covariantieberekeningen worden uitgevocrd.

De bemonstcrfrequentie wordt evenals de logische besturing van de demodu-latoren geheel uit bet referentiesignaal afqeleicl. Figvur 3 geeft het principe van de demodulatieschakeling.

H

3rdp

be s tu nq

od4LCe.

Lr

-6-+ t

.3:

_cpc

f't

te- .c ckmo

cLii Lcka

Uit de referentiefrequentie wordt een subdraaggolffrequentie

5b

van 13,3

kHz afgeleiä, overeenkomend met de ingestelda draaggolffrouentie bij op-name. De pulsreeksseparator leicit uit bet meetsignaal en f twee

puis-reeksen af. In geva:l van positief uitsturencl meetsignaal _{op tape worden} "positief gaande modulatiepulsen" via electronische schakelaar S aan de integrator aangeboden. In het andere geval "negatief _gaande modulatie-pulsen" via S

iaa'i.

diU1e,

(9)

nure-

-/-Het aantal pulsen wordt in dc integ.cator "geteld' _{en aan het einde van} een bemonsterinterval is de uitgaxiqsspnning een maat voor de waarde van het geregistreercie signaal.

Indien f _{golijk is aan f} _{is het versehil in frequentie tussen}

meet-sub c

registratie en f5 gelijk aan de modulrende frequentie Lf.

Indien f verschilt vari f _{ziet de demodulatieschakoling dit als een}

sub c

modulerende frequentie. Di.t ri'anifest.ecrt zich aan de integratoruitgang als een constante offsctspanriing. Deze offsetspannirig wordt m.b.v.

cofficint K

gecompenseerd, zodat amplifier A steeds het juiste

demo-dulatieproduct aanhicdt.

Met coëfficiènt K wordt dc hoogte van de analoge telpulsen _ingesteld die door S en S _{uit de logische pulsen van de pulsreeksseparator} _zijn verkregen.

Voor het uitgangssignaal van de demodulator geldt

V=(f +Lf(t)-f ).x.T.k .k. ₍₁₎

c sub c i

waarin

ff(t)

frequentie gemoduleerd signaal, f draagaolffrequentie,

Af(t) frequentiezwaai overeenkomend met het tijds-afhankelijke vurschijnsel

subdraaggolffrcquentie, afgeleid uit de referentie X _{"teltij&, tijdsinterval tussen 2}

bemonstertijd-stippen

T _{breedte van de logische pulsen uit de} pulsreeks-separator

k _{analoge pulshoogte}

k. _{integratiesneiheid van de demodulatieinteqrator}

De invloed op de u.itgang van de dernodulatieintegrator bij _ongelijke waarden van f en f voigt uit (1)

c sub

Indien f'

5b

gcldt bij afwezigheid van uitsturend signaal:

_/

V L L

).X.L.

off c sub c i

Coèfficiënt K dient dus op de waarde _V0ff inqesteld _{te worden.} Met behuip van cofficiént K _{wordt de maximaal toelaatbare spanning} V _{bij maximale modulatiediepte ¡f} _afgeregeld.

(10)

Hiervoor geldt: V max

K=

c if .X.T.k. max

Tijdens de initialisatiefaso worden uit de per kanaal opgevraagde f en

f de waarden van K eli K _{berekend en aan het begin van de}

dataver-max C o

werking op de analoge computer ingesteld.

Figuur 4 geeft in principe weer op welke wijze de bemonsteringen door de frequentiedemodulatieschakeling uit een analoge meetregistratie _worden achterhaald.

Xse

I

Th

-8-mpL coirnnd

fl?'

die. epuEe

u ti dernodLELc. O voLt. ouEpu

dcodIa.e.

F

g .4:

moL.iin

ri

meL.n

_L

(11)

Analoge nurnerieke integratie

Beide versnellingsregistraties dienen _{t.weemaal geintegreerd te worden orn} over verpiaatsingssignalen te beschikken waarop covariantieberekeningen uitgevoerd kunnon worden. Qwdat de frequentiodemodulatieschakelingen reeds bernonsteringen afgeven kan er geen continue integratie plaatsvin-den, zodat een discrete integratieprocedure gevolgd moet worden.

Uitgangspunt voci: de integratieprocedure is de rechthoekregel

I(n)h.(f +f +...+f +f

1 2 n-1 n

waarin

Deze procedure wordt m.b.v. twee teruggekoppelde _{integratoren en}

bijoassende logische besturing volgens f iquur 5 op de analoge computer geprogrammeerd.

F.5: anaLoge riumetieka

LriU.

stapgrootte

samples van de te integreren tijdsfunctie integraal na n samples

(12)

Hybriede banddoorlaatfilter

In serie met het riumeriek integreren van de versnellingssignalen worden banddoorlaatfilters opgencmen orn offsets

en

(zeer) laagfrequente

drift-verschijnselen te onderdrukkcn. Zio fiuur 6.

-lo-Fi..6:

6anddoL

Ltts

in

se't

rnc

nLce

urn-tke.

De spocificaties van het banddoorlaatfilter zijn.

w=5 rad/sec (f-O,79 Hz)

centrale frequentie

Wbr 3,8 rad/sec

_{br°'6 Hz)}

halve bandbreedte

Het is niet waarschijn]ijk dat een dergelijk smalbandig _filter _{( van 0,19} tot 1,4 Hz) eenvoudiq met analoge componenten gerealìseerd _{kan worden.} Met cen digiaa1 filter is dit echter _{zeer wel mogelijk.}

In principe komt digitaal f ilteron neer o _{een benadering van de} convo-lutieintegraal, dic de relatie beschrijft _tusseri

in- en

_uitgang

_van

_een lincair systeem i.c. filter.

Indien x(t) de ingang, y(t) de uitgang en h(t) de puisresponsie van een filter is, volgt voor y(t)

(causaal systeern; h(t)=O voor t<0) Discretisatie leidt tot

m y.=/t.'h.

.x

J _i i j-].

Do gewichtsfactoren h. worden vooraf herekend uit de pulsresponsie hbd(t=iAt) van het banddoorlaatfilter.

x) _{Het is eon ervaringsfeit dat het spectrum van oceaangolven zich}

uit-sbrekt tot ca. 2 rad/sec. Bij verwerking _{wordt 4 maal versneld afgespeeld,} waardoor, met inachtname van eon gewenst spergebied

_rond

_{0 Hz bovenvermelde} specificaties voor

hot

_{banddoorlaatfilter volgen.}

t-x-,

r dcle o Laa

-fL.

(13)

Tijdens het filteren wordt op jeder hcmonstertijdstip _X. _{ingelezen en} worden de voorafgaande samples (daaroe in het geheugen opgeslagen) na verrnenigvuldiging met h. gesoiameerd en via _{een DAC naar de analoge com} puter teruggevoerd.

De functie hbd(t) wordt gevormd uit het verschil van de puisresponsies van twee causale laagdoorlaatfilter met ongelijke afsnijfrequerìties. Het f ilterprograruma wordt geactiveerd door interrupt requests vanuit de analoge computer.

Een uitgebreide toelichting op de berekening _{van hbd(t) is in de} appendix opgenomen.

(14)

Covariantiefuncties en eriergiedichtheidsspectra

De digitale benaderingsformule van do acvf (autocovariantiefunctie)

voor ceri eindig aantal bemonsteringn worc3t gegeven door

N

c(j)= ---Xx(i)x(j-i)

N aarital samples x(i)

i

en kan herleid worden tot

c(j,t)=c(j,t1)+x(t).x(t-j)

10r1,..,m

c(j,O)=O

n=l,2,... ,N

x(t -j)=O indien t

_n _n

j

m gewenste anntal lags

Uit deze formule blijkt dat na iedere ingelezen sample alle bijdragen

tot de acvf bij de voorgaande bijdragen worden gesommeerd. Na af loop

van een analoge tape run vindt deling door N plaats.

Op overeenkomstige wijze wordt de kruiscovariantiefunctie bepaald voigens

en

k_xy

(j,t )=k

(j,t

) +x(t).(t-i)

n

xy

n k

(j,t )=k

(j,t

)1y(t ).y(t -j)

yx

n

yx

fli

n n k

(j,O)=k

(j,O)=O

xy

yx

x(t_j)=Y(t_j)=O indien t

j

j=O,1,. .

.

n=I,2,..

.

De kcvf dient zowel van x naar y als van y naar x bepaald te worden orn ceri

volledig kruisspectrurn te kunnen berëkenen.

Het autospectruin is de fourier getransformeerde van de acvf en wordt

benaderd met:

f

m-i

C(j)=2At( c(0)+2'c(i).cos(irji/m)+c(m).cos(rj)

i

Evenzo worden co- en kwadratuurspectrum uit de kcvf bepaald volgens

CO(j)= At) (k

(rn)+k

(m)cos(j)+2k

(0)±21(k

(i)+k

(i)

1

cos(rji/rn)i

\ X

X

_/

_\ _xy

yx

_/

J

m-1

KW(j)=2t(k

(i)-k

(i)

sin(ji/m)

Lj\XY

X _/

(15)

-1

3-Na bepaling van daze "raw" spectra vindt "smootling" plaats m.b.v. het Hanning window volgens:

SP (1)=F(1)±F(2)

SP(j)=¼F(j-1)+½F(j)+F(j+1) 2j&rn-1 SP (m)=½F (m-1 )+½F (ni)

waarin F(j) het 'raw" spectrum en SP(j) het "smootled" spectrum voorstelt.

Opmerking.

De accurnulaties tot de covariantiefuncties worden, orn reden van sneiheids-winst, in 32 bits integers uitgevoerd. Hierdoor is de maximale waarde van in te lezen samples beperkt, daar anders de woordiengte overschreden kan worden.

De dernodulatoren worden dan ook op 10 Volt maximale uitsturing ingeregeld, zodat de ingelezen samples tussen -1000 en +1000 blijven.

(16)

Appendix

De puisresponsie is de inverse fourier qetransformeerde van de complexe overdrachtsfunctie H(w) van het filter.

H(w)=A (w)exp (-jO (w)) A(w) amplitudegedrag

O (w) fasegedrag

-

11-Terugtransformatie naar het tijddomeim van H(w) vindt plaats volgens

h(t)=(2)'fH(w)exp(lwt)dw

Uitwerking voor eon ideaal iaagdoor-laatfilter A(w)=1 w O (w)wt levert op h(t)=0 voor tKO te worden voldaan.

Aen doze voorwaarde kan worden voldaan door vermenigvuldiging van h(t) met een eenmalige blokfunctie die symmetrisch ugt t.o.v. to.

Bij voldoende faseverschuiving inì het ideale filter - waardoor do top van h(t) voldoende ver verwijderd ugt van de oorsprong - zal afsnij-ding hot gedrag van h(t) slechts gering hemnvloeden. Dus zal bet

ampli-tudespectrum van het causale filter "niet veel" afwijkn van het ideale

filter. w _(û ii (t-te) -w o C

De functie h(t) bereikt de maximale waarde w/Tr voor t=t0. Nuldoorgangen van h(t) liggen op afstanden rr/w vanaf tc

Bij practisch realiseerbare filters dient aan de causaliteitsvoorwaarde

-1 -1 sin (w (t-t0)

h(t)=(2) (exp(jw(t_t0))dw(1r) _{j cos(w(t-t0)d} C

(17)

Vermenigvuldiging van h(t) met de eennìalige blokfunctie p(t) _{geeft de} causale tijdsfunctie h (t).

C

Uit de relatie

h(t)=h(t) .p(t)

in het tijddomein voigt. do rejatie

H (w)H(w)

P()

C

in hat frequentiedonein voigens de convolutiestelling voor fourier trans-formaties. Voor de in bijgaande figuur getekende functie p(t) en zijn getransformeerde P() geldt

P(w)p(t)

vol gens P(u)=2sin(wt0)exp(-ujt0)/w

-15-K->

;-If1C .

Grafisch is het resultaat van H (w) oak wel _{"aan te voelen' indien} C

men bedenkt dat convolutie de opperviaktebepaling _{inhoudt van het} product van twee functies, waarhij de ene functie steeds op zijn pleats blijft en de andere langs de horizontale _{as wordt verschoven.} De uitkomst van het opperviak is de waarde van de convolutie behorende bij die verplaatsing lanqs de horizontale as.

(18)

-1G-De overdrachtsfunctie van een banddoorla.ttfilter is te beschouwen als bet vorschil van twee laagdoorlaatfilters met ongelijke afsnijfrequenties.

Hbd (w) -H1 (w)

op grond van de lineariteitseigenschap veor fouriertransformaties voigt hieruit

hbd(t)=h2(t)_hl (t)

zodat

hbd(t)=(sin(w2(t_tO)_sin(wl(t_tO))/(rr(t_tO)) (in beide oiderdoor1aten identiek fasegedrag veri

\Jke met de substituties ondersteid)

w =(w +u ) en

w =(w -w

) w centrale radiaalfre-c

21

br

21

c -quentle; 0br halve bano-breedte is te herleiden tot h (t)=2sin(w (t-t ) .cos(w (t-t )/rr(t-t bd br O c O O

Bovenvermelde methode van causaal maken van een ideaal filter ieidt in het amplitudespectrum tot. sterk aangroeiende _{oscillaties in de omgeving} van de afsnijfrequentie. Dit is het gevoig van de abrupte overgangen in de afsnijdende blokfunctie. Bij gebruikmaking _{van cosinusvormige} afsnij-ding echter, zullen de oscillaties grotendeels _verdwijnen.

(19)

-17-Irnmers, in de fouriergetransformeerde van een eenmalige cos-functie zijn

de zijiobben in beduidend mindere mete vertegenwoordigd dan in die van

de blokfunctie.

Voor de berekening van de filtergewichten vari de banddoorlaatfilters worclt

dan ock uitgegaan van

hbd(t)=

_ir(t-t0T

SÍfl

\h(tt0

Ow(t_t0)» (

1+cos(r(t_t0)/t0))

waarin de laatste factoren de cosinusvormig afsnijdende functie

vertegen-woordigt.

Oprnerking.

De berekening van de juiste f iltergewichten heeft aanvankelijk tot eniq

onderzoek geleid. flet ging er voornamelijk orn dat de overdrachtsfunctie

van bet digitale banddoorlaatfilter maximale verzwakking in de omgeving

van 0 Hz opleverde, zodat offsets en zeer laagfrequente drift eon

mini-male invloed zou uitoefenen op de numerieke integraties achter de filters.

Er is clan ock m.b.v. eon apart programma ruimschoots

geéxperimenteerd

met het aantal gewichten en de tijdvertraging t0, waarin juist de invloed

van het causoal maken op de overdrachtsfunctie rneespeelt.

Het filtergewichten bepalende programma I3FGEW vraagt orn onderstaande

parameters over de typewriter cp:

NH

even aantal f iltergewichten

DELTT

tijdsinterval

t (sec)

TN

tijdvertraging t0 (sec)

WC

centrale frequentie (rad/sec)

(20)

(21)

-18--Bancisneiheid bij opname on weergave

Specifieke tijdsintervallon hi..j de dataverwerking

Van bet gebruikte type recorder zijn do specificaties bij de gekozen opnanle- en weergavesneiheid van 3 resp. 1.5 ips:

speed f f +Lf f -tf Af

c C max C max max

(ips) (kHz) (kHz) (kHz) (kHz)

3 3,375 4,725 2,025 1,350

15 13,500 18,900 8,100 5,400

Op real time basis is een bemonstertijd van 1,6 sec gekozen. Dit leidt tot een toelaatbare spectrale breedte van O tot 0,31 Hz. Bij een

cavan-antiefunctie van 50 lags laat deze bemonstertijd een oplossend vermogen toe van 0,006 Hz in het energiespectrurn.

De teltijd van de demodulatieintegratoren is vastgcsteld op 1,4 sec.

De dataverwerking wordt 4 maal versneld uitgevoerd, zodat de bemonster-tijd 0,4 sec en de telbemonster-tijd 0,35 sec wordt.

Op deze tijdsintervallen is de software geschreven en do analoge patching orìtwikkeld. Herziening van die afspeelsnclheid waarhii f en Af _{niet in de}

c max

omgeving van 13,3 kHz en 5,4 kHz liggen, leidt tot vnij ingnijpende modi-ficaties in de analoge schakeling en _{op diverse plaatsen in de software.}

(22)

Demodulatieschake ling Timing- en Control circuit

Hot ingangsgedeelte van de demodulator bevat een aparte flip flop waarmec de blokvorm vanaf tape in de AD4 timing wordt overgenomen.

Met behuip van differentiatoren worden op de opgaande flanken in het uitsturend signoal eri in de sbreferenti enkelvoudige pulsen gemaakt en aan de puisreeksseparator aangeboden, waarin de scheiding plaats vindt tussen positief en negatief uitsturend signaal. Deze enkelvoudige pulsen worden m.b.v. eon monostahiele multivibrator, bestaande uit flip flop en 1- decade teller, omgezet in logische pulsen van constante breedte (20 psec). En vervolgens via electronische schakelaars als analoge telpulsen aan de demodulatieintegrator aangeboden.

In plaats van iedere f lankovorgang in referentie- en fm signaal met een analoge telpuls te honoreren worden uitsluitend die flankovergangen in het referentie- en fm signaal doorgelaten, die rechtstreeks bijdragen tot hot uitsturend signaal. Op deze wijze wordt eon aanzienuijke reductie ver-kregen in bet aantal schakelovergangon van de electronische schakelaars. Het principe van de pulsreeksseparator is in figuur i toegeiicht.

(C o c.o.

I

19-FLJtJTLJ1 rLrLJ .fÇt

fl..

I

i

_F

_fsub r

I'

I r I i

.J.

I I I I t

L..

O3íÇcLf(4)

i

J

O15UI

j r r

It

r Ir

II

I

_i

L_J

L

J

L_J

O Fj

__

I

th5).

Y

fL

Ji_____ ahaL.!uLs

2014.LSec

/

teLpu1sn

_¿

.

/

_ULtnfJ

p 'incipe put

eksetfoi.

_{ddLrn p.ekssepai.ato't}

F.1:

u

kssepo.

_MC

mccuLa

Bij f =13,33 kHz, f _{=5400 Hz on een teltijd van 0,35 sec ugt het aantal}

C max

schakelovergangen tusseri 0,35. (13333+5400) en 0,35. (13333-5400) in geval van honorering van iedere f lankovergang. Bij toepassing van de pulsreeks-separatoe bedraagt het schakelovergangen maximaal 0,35.5400.

(23)

-20

In het timing- en control circuit wordt uit do referentiefrequentie van 20 kHz afgeleid:

teitijd van 0,35 soc voor besturing demodulatoren samolingtijd van 0,4 sec

subreferentiefreauertie van 13,3 kHz

De frequentic van 13,3 kHz is eigenhijk de frequentie waarmee in de demo-dulatoren uit (f +f(t)) de uitsturende zwaai Af(t) wordt uitgeteld. Er is voor 13,3 kHz gekozen omdat de hijbehorende oeriodeduur van 75 psec eenvoudig uit de halfperiodeduur van 25 psec van de referentie van 20 kHz m.b.v een schuifregister samengesteld kaki worden.

Bovendien ugt de waarde van 13,3 kHz dicht in de omqeving van de draaggolf-frequentie hij de weergavesnelheid, zodat slechts een geringe off setspanning gecompenseerd behoeft te worden.

Na deling van de halve periode van de roferentie ( 25 psec) door 1000 wordt

een periodeduur van 25 msec verkregen. Met behulp van vier tweedelers en enkele gates worden de tijden van 0,4 sec en 0,35 sec hieruit afgeleid. Figuren 2 en 3 geven de schakeling en het tijdiagram van het timing- en

(24)

RE F 2okt1 2 dec;c]er 13.c.D. dac.

1.C.D.

i3; kH

dernoduLtic:

cyEc- enabLircj

r

der'oduL4.e-con.

5)

ere

nt

í4IOO

TSH-bstutLn.

L_i_J i

L

I

.1

J

L

F

L______

J

L

tJJ

ce'ta(e ffcd

I

1OD

TSH meE ck J

r

i

F.3: {jdd

an

_'1iiri

_{0,35 s}

_o

_c

_{o k1z}

_'t4.

F2

FI

9i

I 2 I

jo

r'0

-Fi

I

ctoL'cuib.

I I

LLL

_I

L LL

I

L

i

F2 _J

I

L_t_J

I I I 92 F3

II

I I I I I F

l_

Ji

j

ioo

ec.

icr.ri ¡oo

I

sarLc

voot

_fi

(25)

-22-Dc wijze waarop de subreferentie van 13,3kHz uit de van tape afkornstige referentie van 20 kHz worcit afgeleid is ii het tijddiagram van figuur 4 aangegeven.

2j

c at.Ç ie 'tLod L.n o k F 'ec (. V3

n

______r

fcct

H

L

_j--Fi.: afLcLd.n

ubie

_{'i.ntL. vzin 13 ktj..}

Volledigheidehalve is in figuur 5 eon tijddiagram gegeven waarin de sarnenhang tussen de diverse stappen: frequentiedemodulatie, filteren en analoge nuinerieke integrat:ie wordt aangegeven.

(26)

cLt.pvkcr i800

OAçec sarnpLLr sbuLircJ

fderow

-Oup. deod.

TS-4

-.---fit

hurnt.eke

j

itaLe)

-Ç-: Lntzto't.

;o J:

Ln 4 r-&c wo'cLt de

vc.1.

I

1ck.

5Ò't'

[1

1:1.?ck I

çi [

FLg.5:

dc

deut

e

ek

i £

' 2

-c.

0V I

i

oub dc.

ti

+ oV

(27)

Frequentievermenigvuldiging in liet golfversneliingskanaal

Tijdens de opnamen was hot orn techniscie redenen niet mogelijk de afge-sproken draaggolffrequentie en maximale frequentiezwaai van 13,3 kHz resp. 5400 Hz te handhaven.

BIj de verwerking diende dan oak rekening gehouden te worden met een afwij-kende f van 9,2 kHz en een Lf van 503 Hz.

C _max

De op 13,3 kHz ingestelde demodulator ziet eon f van 9,2 kHz als een constante frequentiezwaai van 4,1 kHz. Dit heeft bij eon Af van 5400

max

Hz een offset van ca 80% tot gevoig.

Door toepassing van freguentievermenigvuldiging met eon factor 4/3 wordt de afwijkende draaggolffrequentie getransformeerd tot 12,2 kHz, waardoor de offset geerduceerd wordt tot ca 10%.

Figuur 6 geeft het principe van de frequentievermenigvuldiging.

íJ

L

c2 L t t

-

L t I

FI J

L_J

J

_J1

L L I I

4-- -.

c 3 j_

I J I I I I I _______L__.

/f sec

uit

i I I I L

F7:

24-pncLpe

teuervtricìuLcLLjn

Pc=3

De uit Gi afkomstige blQkvorm wordt op de op- en neergaande flanken gedifferentieerd, waardoor de frequentie wordt verdubbeld. Met deze pulsen wordt flip flop Fi geset en door de teller gereset. Fi geeft dan een blokvorm van de dubbele frequentie. Via twee differentiatoren komt dan uit G3 cen pulstrein van vier maal de ingaande frequentie. Tenslotte vindt rn.b.v. een op 3 ingestelde C precount) teller doling door drie

plaats. Zie f iguur 7 voor het tijcldiagram.

(28)

-25-Organisatie software Organisatie COMMON

De software voor dataverwerking bestaat uit drie hoofdprogramma's

INIFD iriitialisatie

MUCFD eigenlijke dataverwerking C verdeeld over 2 runs)

RESFD presentatie van resultaten

In de hoofdprograxnma's wordt qebruik gemaakt van subroutines voor diverse deelfuncties.

INIFI on line filterprogramma op interrupt requests van de AD4 INTEN interrupt enabling

t.b.v. banddoorlaatfiiteren INTDI interrupt disabling

J

REF]3E referentiebewaking

NEGSP negeren enkele sense pulse (synchroniseren digitale programma) WRTDI schrijven van COMMON naar file op disk

REFDI lezen van COMMON van uit file op disk

MUACC p- en e accumulatie, bepaling lengte meetregistratie AUKRC auto- en kruiscovariantie bepaling

HUARZ alle huip arrays op nul zetter

Bij de presentatie van resultaten wordt nog qebruik gemaakt van de subroutines

SKALE scaling real array in integer array C t.b.v scpoe display)

DISP1 display integer array op scope

DAUPW digitaal autopower uit autocovariantiefunctie DKRSP _{idem CO- en KWaäratuurspectrum uit kruiscov.} HANWI _{smoothing raw spectrum door Hanning window} PRPLT print plot van spectrum op regeidrukker

(29)

-26-Opdrachten in INIFD

- op nul zetten van arrays en ISCAN (bepaling lengte meetregistratie)

- inlezen coèfficintadressen van sonskaart

- inlezen van filtergewichten vari ponskaart

- over typewriter opvraqen van

- aantai te verwerken analoqe databloks

- f

en xf

per kanaal, oer analoog datablok

o max -.

COMMON wegschrijven naar f ilr op disk

Opdrachten in MUCFD

- over typewriter opvragen van

- nummer analoog datablok

_{( voor eigen administratie)}

-

vragen op MUACC reeds is gedaan

zo neen

_{: dan eerste run voor p- en}

G

accumulatie en

telling van het aantal scans

zo ja

: dan tweede run voor covariantieberekeningen

na verwerking ana bog clatablok COMMON naar file

_{op disk}

Opdrachten in RESFD

- opvragen nummer analoog datablok

- print out van

- coêfficiëntinsteliing

- draaggolffrequentie

- maximale frequontiezwaai

- aantal scans

- opvragen opnamegevoeiighoden van

- vaarsneiheid (in/sec)

- invaishoek

(0)

- berekening en print out van

- auto- en kruiscovariantiefuncties

- hijbehorende spectra

- golfhooqte-,dornphoogte- en stamphoekspectriun

_{op w}

_schaal

- ovardraclitsfuncties

- stamphoch/golfhoogte

- domphoogte/golthoogte

(30)

-27-Specifieke gegevens en of tussenresuitaten van een meetregistratie worden in COMMON opgcslagen. File handling van en naar COMMON vindt plaats met twee aparte routines WPTDI en REFr)I.

Organisatie in COMMON met. toelichting op de grootheden:

FC(6) ,FzW(6) ,FR,CON IZWCF(6) rIZWWA(6) MUCOF(6) ,IGEM(6) NtJICF(12) ,NUIWA(12) IGAF 6) INFEX (40) I SCAN SMSAM(6) ,SMKWA(6)

CORi (50) ,COR2 (50) ,COR3 (50)

COR4(50) ,COR5(50) ,COR6(50)

CRCO(50) ,CRKW(50) 11(50)

12(50) ,13(50) ,14 (50) ,15 (50)

I6(0) ,ICO(50) ,IKW(50) SXKWA,SYKWr IFI1 (74) ,1F12 (74) ,1F13 (74) 1F14(74) ,1F15(74) ,1F16(74) IHGEW(74) ISE1 t/m ISE6O inelezen f en ¿1f c max f _{en hulpeonstante} re f

normeer coëf, en waarde voor elke ¿1f max offset cof en waarde

numerieke integratie coëf en waarde gem. waarde analoog datablok

( te verdisconteren in IGEM)

extra gereserveerde geheugenruimte aantal scans (lengte meetregistratie) scm samples, som kwac3raat samples accumulaite arrays en huip arrays voor covariantieberekeningen (6x auto-en ix kruiscov)

en Ey2 extra grootheden hij kruiscovarian tiefucntie

huip errays voor 6 banddoorlaatfilters

array met filtergewichten ( ¿1t

ingecelcu-leerd t.b.v. snelle convolutieberekening) Verpiichte hulpgrootheden in COMMON

(31)

-28-De routines waarrnee de file handling wordt verzorgd zijn WRTDÏ (IFIL) en REFDI (IFIL) voor resp. sc1ìr:jven naar en lezen vanaf disk. De parameter IFIL is het file nummer dat overeenkomt rnt hct nummer van de meet-registratie.

Intern wordt de inhoud van COMMON over 4 files verdeeld. Dus uit het

enkele nummer IFIL worden 4 files geadresseerd. De rela en integer getallen worden afzonderlijk ter file weggeschreven orn te voorkornn dat de inhoud van een file gernengd zou zijn. Voor de gebruiker is dit echter van geen enkel belang. De gebruiker dient er van uit te gaan dat de software voor-ziet in 15 files op disk voor het opslaan van resultaten.Dit houdi dus in dat er 15 analoge databoks ( uiteraard verdeeld over een aantal analoge tapes) in ken keer verwerkt kunnen worden.

Onderstaand overzicht geeft wear welke grootheden in COMMON wel en niet aan de file handling deelnernen.

(32)

niet noodzakelijk naar disk, wel steeds clearing

in MUCFD

hulpgrooth. voor file handling,

J

mogen niet naar disk

totaal integers voor

file

handlirig,

kan in één sector )

Daartoe de DEFINE FILE

DEF FILE i (1,169,u,IsEl)

2(1 ,300,UISE2)

3(1 300 ,TJ,15E3)

4(1 ,256,U,ISE4)

etc.

In WRTDI(IFIL) en

REFDI(IFIL) worden

naast dc file

IFIL

ook nog dc

volgen-de 3 files gebruikt voor datatrans3rt van en naar disk.

INTEGERS PEALS MUCOF(6) 6

FC(6)

x

IGEM(6) 6 FZL(6)

x

NUICF(12) 12 FR x NUIWA(12) 12 CON x IZWCF(6) 6 SMSAN(6)

x

IZWWA (6) 6 S1KWA(6)

x

IGAFvJ(6) 6 SXKWA X INFEX(40) 40 SYKWA X I SCAN

11(50)

x i CRCO(50)

x

CRKW(50)

x

12 (50)

13(50)

x x CORi (50) x COR2(50)

14(50)

X COR3(50)

x

15(50)

x COR4(50) x

16(50)

ICO(50)

x COR5(50) x COR6(50)

z

IIKW(50)

IFI1(74)

128 reals

1F12(74)

1F13(74)

150 reals

1F14(74)

1F15(74)

i50

reals

1F16(74)

IHGEW(74) 74 +

Alle

reals verdelcn over 3 sectoren

169 in resp. lengten:

ISE1 x

2.128=256

2.150=300

ISE2 X 2. 150=3flfl

Voor significante inhoud van COMMON

1SE59 I

X naar disk 4 files à i sector

reserveren:

ISE6O x in édn file alle integers

(33)

-30-AD4 interruptbewaking en -aihandeling

Interrupt requests van de AD4 (t.b,v. digitaal f ilteren) dienen volledig

gescheiden gehouclen te worden van die van IBM devices, omdat het huidige

Operating System niet zondor inner voorziet in correcte afhandeling

van

AD4 interrupts.

0m deze reden wordt er veer gezorqd dat f ilterinterrupts

_uitsluitend

optreden in lint "analoge data verwerkende gedeelte" en gemaskeerd zijn

tijdens informatietransport naar en van disk.

De interruptbewaking vindt plaats in én routine met de twee entry pòints

INTEN

_{enabling AD4 interrupt}

INTDI

_{disabling AD4 interrupt}

Daze routines worden in MUACC (li-a accuinulatie) en AUKRC (auto- en kruiscov.)

viak voor en direct na hybriede verwcrking van cnn analoog datablok

aangeroepen.

De convolutieberekeningen vinden plaats in het via interrupt branches

bereikte gedeelte van het programma INIFI. Ilet

FORTRAN- aanroepharc

gedeelte in INIFI verzorgt de parameterinlezing.

CALL INIFI(IHGEW,74,IFI1)

IHGE

_{eerste adres van array met filtergewichten}

74

_{aantal gewichten}

IFI1

eerste adre

_{van de cerste van zes direct achter alkaar}

voigende arrays IFI1 tim 1F16 welke als back up arrays

voor de 6 kanalen worden gebruikt voor ooslag van delayed

samples

Opmerking

(34)

Referentiebewaking (subroutine REFBE)

Daar zowel het opkornen als het afvalinn _Cil _{afwezig zijn van} referentie-frequentie indicaties bevat voor de voortgang in _{de software is voor de} bewaking van het referentiesignaal _{cen afzonderlijke routine gemaakt.}

REFBE (IPAR)

IPAR=0 _{verblijf in REFBE tot referentie "op'} _is

IPAR=O _{entry in REFEE tijdens referentie 'op" (anders}

foutmelding) wacht tot referentie "af" is

exit indien referentie weer "on" is

Figuur 8 geeft de bij REFEE behorende schakeling.

R R

c k.

Fi

-31-Fi. 8

_j

F2 krijgt na elke SQ 3lsec een setpuls. _{Intern is op verschillende plaatsen} een wachttijd van 2 msec aanwezig. Bierdoor wordt pas na 40 setpulsen

intern steeds clearing:) geconcludeerd _{dat de referentie inderdaad "op'} is, zodat eventuele valse starts _{tot een minimum beperkt blijven.}

SLR (jr bu)

(35)

O

C)

;k cr. '...a _{_a ...e*... ....&a,,} _... _.

t)

O

Bij lagen.

II JO3

/1 FOR INIFD

::rl.IST SUUkCE PÜGRÁM

PUG.AM.

::fNt

;nRD ItTEGERS

::rIDCS C TYPEWRITER) :1 OCS ((EYBOARÜ ) :r _{( c} _{RD )} :: _{( D I S K )}

;:Irc5N443 PRIMIER)

CÜti1DN IZWCF() , !ZW!A(6)

, 1GAFWC} ,MUCOF(f) ,!GEMC6)

V { 53 ) ,

I 3 50 )

I 4 ( 50 ) I 5 ( 50 ) , I 6 1 50 ) , I CU I O ) I }K ( 50 ) I F _I

( 7'

} , I F I 2

V , I F I 3 ( -7'.s ) I F I ' I 74 ) , I F I ( / ) , I F I C 74 ) , i HU t i 4 ) , F C ( ) , FL W C o) , FR,:

Ví(; ( ) , C 6, ,

XI'A , SYKit , CRC[J( 5U)

CRK (f0 )

_{CD1 ( 50 )}

_{, C Lt2 C}

D ) ,

V3 (50)

_{LiPid bU) ,COr550) CUô( 50) ,JSilSE2, 15L3, JSt4,}

CiiMIIIJN

VE9IE40, IS F4iiSF'+2, .SE3, I SE4', J SF45, 15E66, i

_{SF47, 15F48, I}

VE50,I5EI,ISE52jISE53piSEI5E55ISE55J5E57,ISE5D;I5EY9,I5FbO

([JMMEIN 1N1-Ex(40hNUICF(l2},NUIWp(l2),C0N

DEFINE FILE 1(1,IS9,U,IE-i)

DEFiNE FILE 2(1,300U,iSE2)

DEFINE FILE 3(i?300,U,15E3)

DEFINE FILE ti(1,256,U,15E4)

DEFiNE FILE 5(ii69,h1SE5)

DL FINE FILE 5(1 3O0,UI5E6)

DEFiNE FILE 7I,3O.,U,ISE7)

DEFINE 1-ILE 8(1,256,U,ISE3)

DEFINE FILE 9(1,ì69,UI5E9)

DEFINE I:TLE 10(i300UISEI0)

DEFINE FILE 1l(1,300,U,1SEIÏ)

DEFINE FILE 12(1256,UISEL?)

DEFINE ElLE l3(1I69,U,ISE1)

DFFIE FILE 1t(l,300,UISEi4)

DEFIN-E FILE 15(1,300,U,ISE15)

DEFINE FILE 1S(1,256UISEI5}

DEFiNE FILE 17U,169,U,1SE17)

DEF-IFE FILE lall,300,U,ISEJ8)

DEFINE FILE 19( i9300U, 15E19)

DEFINE FIL.E 20(i,25b,L.I,15E20)

DEFINE FILE 2i(I,iô0,tJISE21)

DEFINE FILE 22(1,300,U,15E22)

DEFINE FILE 23(1,300,U,i5E23)

DEFINE FILE 24(IT25ÔPUtISE2'+)

DEFINE F1L

_{25(1,169,LI, 15E25)}

DEFINE FILE 26(1,300,u,1sE26)

DEFIÑE FILE 27(1303,U,ISE27)

DEFI'E FILE 28(l,256,UiSE2)

DEFINE FILE 2(1,Iô9,U,ISE29)

DEFINE FILE 30(1,300,U,J5E30)

DEFINE FILE 31(1,3()ü,U,15E31)

DEFiNE FILE 32(i,256,U,1SE32.)

DEFINE FILE 33l,i69,U,1SE33)

DEFINE FILE 3l(ì,30G,U,I5E3t)

DEFINE FILE 35(1,30091J,IE35)

[)EF I NE FI LE 3(( i , 256U, i E3(

DEFINE FI

- '7 (0 H

; )

'I

C)

(36)

H

Ç)

Q

o

Th

r

;j112

z:?

j

DEFINE FILE 38(1,3ÜO,UtSE33)

DEFINE FJLE 39(1,3Of)U,ISE39)

DEFINE FILE

_{0(1r?.ó,U,J.SE4Ü)}

DEFINE FILE 6111,1&9U,!SE41)

DEFINE FILE 42(13ODU,)SE42)

DEFINE FiLE

43(1,300,J,1SE43)

DEFINE FILE 44(1 ,25o,U, JS4)

DEFINE FILE 45(A,iô92UE45

DEFINE FILE 4o1,30OU,ISE46)

DEFINE FILE 67(l,300,U,ISE47)

DEFINE FILE 48Lt,256UEJSE48)

DEFiNE FILE 4$(i,I69,UI5E49)

DEFINE FILE 5O(i,30OUISEO)

DEFINE FILE 5H1,3OU,UJSES1)

DEFINE FILE

_{2(1,25(7UiSE52}}

DEFINE FILE 53(1169,U,15E53)

DEFINE FILE 54(1,300,UiSE54)

DEFINE FILE 55(l,3oQ,U1SE55)

DEFiNE FILE

_{6(1,25,U,ISF.56}}

DEFiNE FiLE 57(1,1C9UI5E57)

DEFINE 1-ILE 58(1 3OOU, 15E58)

DEFINE FILE 59(1,300,U71SE59)

DEFINE FILE £O(i,25,U,!SE6Q)

I SEl =1

15E2=l

15E3=1

15E4=l

ISE=l

T f

r i

J U - i

I SE7=i

iSEd=l

15E9=1

I S E 10 = i 15E 11 1

I 5E12l

J 5E 1 3 = i I S F 1 ' i

I 5E151

1 5 El 6 1

I SE 17=1

I 5E 18=1 15E 19..: 1

I SE2O=1

1 5E21=1

I 5E22=1 I SF2 3=1 I SE2'.=1 I 5E2 5l E. r: L .J L. L 'J -15E27=1 I 5E28=1

ISE9=1

ISE3O=1 I SE31 =1

1SE32=i

I 5E33=I

ISE3tf=1

15E35=I

I SE361

I SE37= I

(37)

Lì

í)

n

pm

I S E 39 = i I S E 40 = 1 ISE41=1 I SE42=1. I SE43i I SE'.4=1 I S E 45 1

i SE461

I SE 47 1 ISE4B=1 I 5E49=1

i SE01

1S51=1

I5E521

15E53=1 I 5E54- 1 I SES5=1 I 5ES ¿S i I SE571 I 3E58=1 I SE 59=1 I SF60 1 F R 13333 33 (T [i N O O

DU i J15O

COR1( J) =0 0 CUR2(J ) =C)JJ

CUR3 (J) =0.0

C OR 4 C J) O . O COR5 (J Y=0.O CUR6( J) =0.0 CRCO(J ) =0,0

..1 CRK(J)=OO

S X K W A O O S Y K a A O O

RLAD{2,3) IZWCF

3 FUMAT(6(I41X))

RED(2,3) MUCUr-.

REAUC2,4) UU1CF

4 FORMAT (12(14e IX))

DU 7 Jz16

IL1!CF(J)=IDEOC(IZWCF(J))

7 M1C0F(J)DEDC(UCD(J)

[)0 9 J'112

9 NU1CF(J)=IDEOCU.UICF(J)

1W IlI\ (11=4000 NUIWAC 21=10000 NUI A (3) =10000

HUIA(4 1=4000

NUIWA( 5)-100û0 1W i ( ¿S):: 13000 'UI1A( 7)=4000 (8) =10000 NUI OA (91=10300 F!UIWA(1U)=4300 raJrwA(li =i0000

NUJW( l)=10000

REAl) (2 II) IHGEW

li Fl)1JÏ(13JS)

DU 12 J=140

(38)

L)

C?

(

k

LI

ç)

\j1

o

40 R1TE(1

'1 FORMAI ( 'GEEF At TAL

AN4LD&E

DATAOLOKS OP IN 12',!, 'CALI8RATIE5

VT METFLLEN ')

CALL TYPFO

REAù(6,62) LADAB

42 FORMAT(I2)

Jf IADAB-1) 45,454:

43 WRITE(1,44)

44 FORMATI '':

ERROR

1, 'MAX AANT NIET GRUTER DAN 15')

OU TO 40

45 WRITE (1 js6)

46 FORMAT (

GEEF AANT KAN OP TAPE IN 11')

CALL TYPED

RLAD(6,67) IANT

47 FORMAT il)

IF(IANT-6)

050,4$

48 WRITE(1,49)

49 FLJRMA1 I 'ERROR MAX IS 6')

GD TL1 65

50 INFEX( 1)IANT

IDBTLl

54 WkITF(153) 1DBTL

53 FORMATI 'GEEF V DATAOL ',12,' ACÌITEREENV FC IN

_IS')

IZWT=1

57 URITE(1,55) IZWT

55 FORMATI 'VAN KANAAL ',12)

CALL TVPFÜ

READ(6,56) 10F

56 FURiAT(I5)

FC(IZWT)::JDF I zol =1 ziíi

IF(IZWT-6) 57,5758

58 WRITE(1,63) IDBTL

63 FORMATI 'GEEF V DATAOL ',12,' ACHTEREENV FZWAAI IN

_14')

I Z W T = i

67 WRITF( 1,65) IZWT

65 FORMATI

'VAN KANAAL ',12)

CALL TYPÍÍJ

READ(6,66) ¡DF

66 F[RAT ¡'t)

FZI( 1ZWT)1D1

I iWT=I ZUT+1

¡FI ¡ZUT-6) 67,6/,68

63 DU 69 J=1,6

CCAkR=0.1E0I(7.0:FLU(fl)

GGEM( FR-FC (J)) CCARR:O00Ü7

1ZWÎA(J1=CCARR

69 IGFVfJ)=GGFM

CALL

_{WRTDI(IDE3TL)}

IDBrL=IDOTL+1

IF(IDBTLIADAB) 54,54,77

77 CALL EXIl

END

FEATURE.S SUPPORTED

NONPROC F SS UiE dflRD INTEGERS

Ocs

-

44

,-'-

k -. -

#U.

o

Ç)

(39)

i J i I

L)

n

L)

L

C)

i

n

Li

Ç)

n

C?

Ç)

1)

II FO

DU? FUNCTION COMPLETED

STURECI

¡NIFD ¡I:

;FILESC 1,FIL1,O) F ¡LES (2, F I L2, 0) F I L E .5 C 3, F I L 3 , 0) 'F I LES C 4, FI L4 0 F I L E 5 (5 , F I L 5, 0) :::F J LES C 6, F I L 6,0) F-ILES(/,FiL7?D) F I LES 3 ,F ILO 0) lLES (9 ,FIL9 0) F I LES (10, F J LiD 0 F J L E 5 (11, F I L1 1 ₀₎ F I LES (12, F I Li 2 D F I L ES 13 , F I L ). 3 , 0)

FILES(14,FIL14,0)

1- 1 LES I i S F I L 1 5 0)

FJLE5(16,FILi6,0)

;rFILES(17,FIL17,Q) F I LESI 18 ,FI L1R 0) FJLESC1C),FIL19,0)

FILE5(20,FIL2Q,0)

'FILES(21,FIL21,D)

:-lLES(?2,FIÌ27,O)

FiLESI23,FIL23,O)

FILE-SC24,FIL2t0)

F1LES(25,FIL25,D)

F I LES (26, F IL 26,0)

FILES(27,FiL27,0)

' FI LES (28, F I L28 0) :F I LESI 29, F ILD ?O)

.F-ILES(30,FIL30,Ò)

FIL(31,FIL3i 0)

FILES(32,FIL32,0)

F I LES (33, F1 L33 0) F I LES (34, FI L34 0) ::F I LES C 35, Fi L35 , 0) F I LESI 36 F J L 36,0)

FILES(37,FIL37,D)

F I LES (38, F I L38 0 F I LES (39, r-i L39 0)

;F I LESI 40

,FI L40 0) FILES (41., FIL41 0)

FILES(42FIL42,0)

F I LES (.43 EI Lc+3 0) F I LES (4'+, F I L44 0 F J LES (45, F I L45 0)

F1LES(46,FIL46,0)

F I LES (47, FI L47,O F ILLS ( 4 , F! L48 0)

lLES (49,

F1L49 0) -F J LES (50 ,FI LSD 0) F I LES (51 F1 L51 Û)

FILES(32,FIL52,O)

1- I LESI 53,Fi L53 0) 'F I LESI 54,FI L540 ).

-I

-tZ r

----. -.-

-.

--.-

-,--.-.--i

_{END UF CLIMPILATIUN}

(40)

- . r I L £ (

9, i

î L J i

'i_J

FILES(6,FIL6O,D)

3:t C C E N D

CLB, BUILD INIFO

CL[, INIFO Lù XC

DUP FUNCT I CN CONPLETED

r)

_{k j} A)

n

_u,

I

-J-.)

t.

I

o

C)

_u

ç)

C)

LI

C..

Li

Ç)

C)

(41)

o

C)

/1 Jû13 li F!TR

UCFD

L1ST SURCE PRUGRMl

SS PR[JGAM

LiNL

CRD INTtGERS

:: ₎ _{S { T Y P E W R I T E R)} :'I flcs (KEY3OARÜ Y i: C? C. S ( L) I S K)

DIMENSIUN IFDÁC(6)

CLi1[N

I ZWCF ( ( ) , I Zt'!WA C } IGAFW( ) ,NUCOF ( E) , I GEM C 6) , I SCAt, 11 _SC')

VC 50) , 13( 50) , Jt( 50) , I5( 53) , 16( 50) ICU (53 ) , (50 ) , IFI i (74 , I FI2{

V , I F I 3 ( ) , I r: I ' C 74 ) , I FI 5 C 74) , I FI ô ( 7) , IHGE4( 74) , FC t ) , FZ ) , FR,

V3(5:j),CLR(Do3rLfl5(0)rCOR6(5o),ISE1,ISE,ISE3,ISE4,ISE5,JEb

CuMMrt

CUMMDt IFEX(4D)UICF(1?)i1tJIW112),CON

DEFINLT FILE i(1,i69U,ISE1)

DEFINE FILE 2(1,3OOU,IE2)

DEFiNE

FILE 3(13ÜD,L,iE3)

DEFINE FILE 6(i,?5úU,IE4)

DEFINE FILE 5(ï,16,U,ISE5)

DEFINE FIL.E 6(1,3ODU,ISE6)

DEFINE ErLE 7(1,300,U,15E7)

LiEFJNt F1LL- E(1,25f,U,IE8)

DFINE FILE )(1,i69,U,I:9)

DEFINE FILE

1O(1,3ÛO,U,IE1O)

DEFINE FILE ii (1 ,300,U1 SEII

DEFINE FILE 12(l25óUISEi2)

DEFJNc FILE 13(1,169U,1SE13)

DEFINE FILE 14(1 ,3OO,U ISEI4)

DEFINE FILE i(I,3OO,UISE1S}

DEFINE FILE l6(1,256,UlSE16)

DEFINE FILE i7(l,1,U?ISE17)

DEFINE FILE 13(1,303,'U,15E18)

DEFINE FILE 19(1,3O,u,ISE19)

DEFINE

FILE 2O(1,255,UIS2Q)

DEFINE FILE 2i(1,1$9,U,

15E21)

DEFINE FILE 22(i,300,U,ISE22)

DEFINE FILE 23(l,300,U,iSE23)

DEFiNE FILE 24(125b,U,1SE24)

DE-FINE FILE 2(1169,U,ISE2E)

DEFINE FILE

26(1,3OO,U,ISE2)

DEFINE FILE 27(1«OC,U,ISE27)

DErINE

FILE ?8(1,25,U,1SE28)

DEFINE FILE 29(1,169,U,15E29)

DEFINE FILE 3O(13O3,U,ISE33)

DEFINE FI Lt 31 (1,30u,U, iS&31 J

DEFINE FILE 32(1,256U,I5E32)

DEFINE FILE 33(119,U,iSE33)

DEFINE FILE 36(1930Ü,U7ISE34)

DEFINE FILE 35(1,3E09u,15E35)

DEFINE FILE

3o(1 ,25,U, 15E36)

(42)

--C)

()

o

'O

(j

o

(Th

o

-

.* *

t

C_a-4 DEF¡NE O E F I F. E DEFI F D EF i N E O E F I N E D E F S N E D E F I N E D E F i NE D E F I N E D E F I N E D E F I N E DEP I NE D E FI N E DEN tiE D E F i N E D EF I NE D E F I N E D E F I N E D E F I NE D E F I E D EF I N E D E F I N E

ISEl

=1 r s E i t Cr_{.J L}

-I 5E4=1 i SE5=l I 5E 6- 1 I SF7 =1 I SE&= i

15E91

15F IC::1 ISEi 1=1 I

SF12::

i ¡SEl 31 r s E 1 4 i I SE151 I S E 16:: 1

¡SEl 7i

J 5E J. F i I SEI 9 i I SF20:: i I SE21 =1 I SF22:: i I SF23:: I

ISE241

I SE 2.5=1 I SE6=1 I 5E27- i I 5E

2E::

1 I SE29=i 1SF

30

1 I SF31 =1 I SF32:: i I

3E33=i

ISE34i

I 5E35--1 I 5E36=l I SF37= I

15E 381

F I LE FI LE F I LE I: i L E

FILE

FILE F I L E FI LE F I LE

FILE

F ILE F I LE FILE f:J L E FI LE

FILE

FILE FILE FI LE

FILE

FILE Fi LE

39(1,300,

U9 I SF39) 40(19 2u 0 J SF40

41(1 169,0,1 SE1)

142(l,3OL;U 15E42)

(f3(

_{130O,U,i5E43)}

44 (1, 2h ,

U, ISE4')

45(1, 169,0, 15E45)

66

1 :300, U IS

F66)

4

7(1

300, U, I SE-47

43(1,±5b,U,JSE43)

49(l,69,U,1SE49)

50(1

,300,U, 15E5Û SiCi 300,0, 15E51

5211,

25o , U, 15E-52)

53(1, 169,Up 15E-53)

54(1,300,IJ,IS[:56)

55(i,300,U,1SE5S)

56(1 ,256U, 15E56) 57(1,169, U, ¡SF57) 581 1,300,0, ¡SF50) 59(1 ,300,U,1SE59)

60(1,25óU, 15E60)

(Th

G

U. o

r'

(43)

*

I SE'iO= i JSE41=1 I S E ' 2:: 1 ISEA3=1 ISE/*4::1 I SE 4 5:: 1 I 5F46 i 15E47::1 15 [48:: 1 ISEL+9::1 I SE5O1 I SF51:: i I 5E52=1 J5[53::1 I SE5'+1 15E55=1 1SE56::1 I 5F57 ISEDB=1 J SESO = i I SEbO i

CALL IIYTSTC3,GJ

CALL I Nl TACI E, O)

CALL REF[)I(1)

CALL STGLK( NUIC F, NUl WA 12)

IFOAC(1)::JDEOC(23001

IFDAC(2)lDEOC(3O2)

\

I

Y

IFOAC(3)=IDLOC(2304)

1FDAC(4)::IDEDCL336)

J

/

I FLAC (5) JDEOC (2310) FLAC (6) ::jDC (¿312)

I AANT5O

IKA;N=6

i ERITEC i 2)

2

FURi'AT('GEEN AUTCINAT lNREG,/,Z.ELF STEEDS REC UP JL'ISTE POSITIE

VETTEN 1 'LT DATASW O 4F EN START 3800 DOOR IN PAUSE 7777

I, 'L

VSW U IS NUODSTDP IN (1UACC')

PAUSE 7777

9 WRITE(1,3)

3 FüRAT C 'GEEF F(L1

ANAL DATAL IN 12 ,DATA5

O AF ')

CALL TYPFU

READ(6,4) 1ACAB

1 FOTfl)ATC 12)

CALL REFDI( IADAB)

CALL INIFI C IUGEW,7 1 Eli)

WRIT E (1 , 5

5

FGRI4AI C 'NUACC AL GEUAAN JA=1 NEEO'

CALL TYÍ'FU

REI\D(6,o) Iô.SL

6 FD;14T (Il

Ir-( IBSL)' 23r10,20

1G WRITE(1,11)

ii FURAT C 'SÌART EEC [I'4

1800 iN P7777

CAl_L STisLK (I ZWC F, I ?.WWA, 6)

CALL STLLK(MUCOF,TGEM,6)

CALL HUARZ

CALL STbLR(JFDAC,IFI1,b)

PAUSE 7777

CALL IN1UX( IE,i)

(44)

3. C)

D

p

o

Ö

o

(Th

t. .

CALL TY?FU

CALL MLtACC WR 1 T E (1

14)

14 F[JRNAI ( EXIT FRUM

UACC

GO TO 30

20 4RITE( 1,11)

CALL

STLK( IZWCFrIZWA,5)

CALL SThLK(vUCDF, IG'S

CALL HUARZ

CALL STPLK( JFDAC,IFl1,6)

PÍ\U.'E 7777

CALL IN1UX(IE,1)

WRI TE (1,21)

21 FÛRMAT( 'ENTPY IN AUKRC'

CALL TYPED

CALL AUKRCiIKANN,IAANT,

ISCAN,C0RlI1,CflR2I2,CDR3,I3,CDR4,l4,Cfl

VI,EÛR6,16,SXKWA,SYKWA,CRCO,ICO,CRKtIK)

WRITE (1,22)

22 FUSAT'EXIÏ FRON AUKRC)

30 CALL WRTDI(IADAB)

WRITE (1 31)

31 FURMAT('DUDRGAAN JA1 NE[=0')

CALL TYPED

REA(66) IB.SL

IF(IDSL) 9,34,9

34 CALL INITA(IE,D)

CALL EXIT

END

EtA ruREs SuPPORTED

NO N P. DC F S S

UNE WRD IN1EGL-R5

hics.

CORE REU1REMENT5 FUR NUCFD

COMMUN

1930

INSKEL CD4MUN

O

VARIABLES

378

PR0GRAi

672

END UF CUMPILATION

'i tic F D

DUP FUNCTION COMPLETED

51UECI

MUCFD MUCFD

:F I LES (1,F IL 1,0)

F I LE (2 :1 L2 ri)

FILES(3,FIL3,O

J LE

(4,

F I L4 O)

FILES(5,FIL5,O)

F ILE S (6 Ei L6 0) F I LE S C , F I L 7, 0 ) F I LES (3, F I L , U) F1LE..5(,FIL9,O)

FJLES(10,FILi0,D)

F I LES (11, r I Lii 0) F I L5 (12, FIL 12,0)

::rF I LES (13, ri

Li 3,0)

:F I LES (14 FI L14 3)

I LE (15,1 LiS 3

(45)

R

-Th

o

ç

C)

o

T)

F I L E S C 19 , F I L i C) , O F I LES (20, F I L20 F I LES (21, FI L21 0

FILES(22,F1L22,0)

FILES(23F1L23,0)

F ILE S (2, FI L24 , 0)

:FILES(2,F1L25,0)

::F I LES (26, FI L26 0)

FILFS(27FIL27,3)

F I L E S ( 23 , F I L ¿8, 0) - I LES (29, F I L2C)

0)

lLES ( 3D, FIL3O 0) ::FjLE5(31,FJL31,9) I LES (32, F1 L32 O) F1LLS(33,fiL.i3,D)

:FILE5(34,FjL34,D)

1LE(3,FIL35,D)

I LES (36, F I L36 0)

'FILE5(37,FIL37,0)

-b ILE 5 33, F1 L8 , 0) F ILE 5 (39, F I L 39,0) Fi LES ( O F I L+O 0

FILES(41*F1L'i3O)

F i LES C f I L't2 0) F I LES (t3, FI L43 0 F i LES (44,1 L44 , O ::F I LES C 45,1 L45 0)

FILE5(66,FiL46,D)

::F I LES 47 FI L47 ,û vE I LES t 48 F I L48 0) FILES(49,FIL'+9,C))

FILFS(5O,FIL50,3)

F I LE.S (51 f I L51 .0) F I LES (52, F 11.52 0)

FILES(53,FIL53,0)

F1LES(b4,FIL54,O)

FILE5(55,FIL55,0)

F I LES (56,FI L56, O F I LES C 57,1 L57 0)

ILES(5,FIL58,0)

I LES (59 ,FI L.59 0) F I LES (O , F1600 CC E N D

CLB, BUILD MUCFD

(LB, MUCFD LO X

(46)

)

:2

T)

tC.c-S.L. a.tQtMt. »... V tS .4S.S... ....%4.» -+ - ,- - - - +: +.. - .V...t. - _. Ws. /1 JiiB . 1/ Ffl RE.SFD

-VLIST SUUCE PCSRAM

:NìPrUiCES5 °RDGRAM

;:3J[

_{':D INTEGERS}

::

ircs

_{t TIPLWR I TER)}

:: iC5 C KEYcCl/ki) )

i::cs (CARL))

VI)CÇ(DJK)

;:VIflCS ( :

43 PF-INTE)

COMMPN IZCF1ó) ,IZWA(6) , IGFW(6) ,MUCUFC6) ,

GEM( 6) ,

SCÁN Il (50) ,

V ( 50 ) , 13 ( 50 ) , I 4 ( 3D ) 15 (50 3 , I ( 50 )

. 1CJ (

) ,

IK; ( 50 )

i F I I C 7i I , I F I 2( 7

V ,

I FI 3( 7) , J FI £( ¡() , I

I5 ( 74) , J F-I6( ?

)

IFT&E( 7) , FC C 6 ) 1FLI(b ) ,FR,Sr

CuMfl--CJ:'iEi

DEFINE NLE !H,162,UI5E1)

DEFINE FILE 2(i,300U,ISC2)

-DEFINE FILE 3(l,3L)0,UISE3)

- +

DEFINE FILE 4(l,26U,ISE+)

DEFINE FILE- 5(1,169,lJ,I.SE)

DEFUE FILE

(1,300,U,tSEó)

-

-DEFINE FILE 7(1, 00,U, ISE?)

-DEFINE FILE 8(1,2S6,U,13E8)

DEFiiE FILE 9(I,169,U,ISE9)

- DEFINE FILE 10(1 !300,U ISE1U)

DEFINE FILE U(i,30ù,U,ISEI1)

DEFINE FILE 12(1,25b,U,LSE12)

DEFINE FILE 13(1169,UI5E13)

DEFINE FILE 1-(1,3DÜ,U,iSE14)

DEFINE FILE 15(i3OU,U,IS15)

DEFINE FILE 16(1,256,U,15E16)

DEFINE FILE 17(1-,169,U,15E17)

+

DEFINE FILE I(I300r1J,ISE18)

-DLFINE FILE 19(i3OJ,U,I5E19

-DEFINE FILE 20(1,25h,U,ISE2O)

DEFINE FILE 21(1,169,J,15E21)

DEFINE FILE 22(1,3CV0,U,1SE22)

DEFINE FILE 23(1,300,J,I.SE23)

-DEFINE FILE 24(125LV,U,I5E24)

DEFINE FILE 2(11ó9,U,1SE25)

DEFiNE FILE 25(i3OU,U,ISE26)

DEFINE FILE 27(i3Oü,ULS21)

DEFINE FILE 28(1,2.ó,U,ISE28)

-DEFINE FILE 29(ìl69,U,ISE?9)

--DEFINE FILE 30(1,flL),UISE33)

--DEFINE FILE 31(1rÛU,U,IE31)

DEFINE FILE 32(1,256,U,15E32)

O

DEFINE FILE 33(1,16U,1S33)

DEFINE FILE 34(1,30u,U,ISE3()

-DEFINE FILE 35(1,Ûì3U,ISE35)

DEFINE FILE 36(1,256U,I'SE6)

(47)

C)

L)

n

DEFINE FILE 33(.1,OJ,U,1SE3)

DEFI' ÍIL

3J(1,3D3,ULSE39}

DEFISE

_{FILE LD(1,2'6,U,iSE4.))}

EFINL :TLE (1 (1, 19,J

ISL-+1)

DEFIE FJLE 62(1,3O32UI5E42)

DEFiNE FILE

't3(1,3CDU,.SE.'.3)

DEF I NE EI LE '+4(1 ,25 U I SE4'+

DEFINE FILE '+5(1,ì$U,1SE45)

DEFINE FILE 4ó(173ÇU,I5E()

DEFINE FILE '+7(1

,3flO,U

15E47)

DEFINE FILE 4h(1,25h,UDISE)

DcF I NE EI LE 4')( 1 i63 ,U , I .SL49)

DEFINE EILE 53(1,30u,U,ISE5Q)

DEFINE FILE 51(1 ,30u,U, 15ES!

DEFINE FILE 52(ì9256,U,iSE52)

DEFINE FILE 53(1 ,169,U, 15E53)

DEFiNE FILE

(i?3OO,UISE54)

DEFINE EILE 55(i,3OO,U,iSE5)

DEFINE

FiLE

56(i,256,UjSE5)

DEFINE FILE 57(1, 16),U 1SE57)

DEFINE FILE 5b(l,30U,U,15E58)

DEFINE FILE 59(13ÛO,L',I5E59)

DEFINE FILE

U(1725(,U,15E6O)

I SE! = 1 I 5E2=1 I E3 =1 I S E 4 1 I 5E5=1 ISE6=1 I 5E7= i I SEcS 1 I SE i I SEl i =1 I SE12= i I SE 13 1 I SEl 4=1 (SF1 5=1 15F! 6=1 15E17=1 ISE18=i I SE19=1. I SE20=1 ISE21=1 I 5E22=1 I S[23=I I SF24 i I SF25 1 I S[26=1 15E27=l I sE28= i I SE29=1 I S E. D i ISE31=1 I SE-32 =1

I SF33 i

ISE4=i

I 5E35=1 I S E 3&=i I sE371

- _a..- ..a.. _{...t.ei .*.. .}

i

L

L)

)

(48)

t)

C)

i)

D

p

C)

..

...,_ ,.at.flUaJ ts,.

..,.

*_ L

-1

I .

31 r SF4] =1

I SE'+21

IE+31

ISEY+1

I sE45=1

ISL'61

I S E 7 1

I SEC1

ISE+91

I SF5G

I

ISEi11

15E52=1

ISE531

I SE541

i

F55=1

I S F 56 1

I SE371

IS581

15F59=1

I S F6 O

i

CALL riYTST(3,O)

CiLL IN1TA( JEDO)

t: R I T F C

i , I)

FURi'IAT C

tSCJPE DI SPL CU.

_{EN PLThER VIA DAC 314 EN DAC 316 tJEL CL}

VEIl 30?Y)

10 CALL R51

CALL RES2

CALL RES3

R I i E (1, 2 )

2

_{FLU'tAT( 'DUflRGAA}

JA=i NEED)

CALi. TYPFU

REAL)(6,3)

fl3

3 FÚRMtT (il)

IFUD-1) 11,1O11

11 CALL EXIT

E N D

FEATU<ES SUPPL1TED

N CN P DC F5 S UNE

DR9 INTEGERS

I C5

C[jE REJI REENTS FUR RESFO

COMMUN

193

INSKEL CDMMD O VAR lADLES

362

PRO(RA

3?.

END OF cOPILTION

RESFD

DUP FUCTIFN CDM'LETED

/1 F(h<

:LIST SLiURCE I'JSRtM

SURÜUT1E

ES1

DIMESILlN ALFA(b)

C

HlflN

I Z

C FI 6) ,

I Z.

A 16) 1

AF; (6) , MUCUF (6)

, i Ç.E'

t 6), 1 5CA

, Il C SU)

V t O ) ,I C 50) , 14 C _{0 )} , I 5 C 5 i) , 1 6Ç 0

, I

C BC 50 , I K C 0 ,

I F13 ( 7

3, iF I

Ci

(49)

Li

ç. -C1J'iiU

E17,I5E1g,IsF1Q,ISFo,r21,i5Ez2,1SE23,ISr2',,IE25,Io,ISE27,

C[iMMiTi' 2 WR I i E i , 3)

3 FDRNAT( 'CEEF Nfl ANAL DATA BL IN I2)

CALL TYI'FU

EAD(6,) IADAB

6 FURIAT (12)

IF(1ADA-15) 66,5

5 RITE( 1,7)

7 FLRAT( 'ERRUR MAX IS 15')

C.D ID 2

6 CALL REFD1(1DAU)

IANT It!FEX( 1)

RITE(3,9) IADAB

W R I T E ( 3 , 6 0

40 FUAT

C 'OINGESTELÛE CUEFF WAARDEN' I

WRITE(3,41) IL:A

41 FURMAT(a017CF UP

',o16)

tFI1E(3,L+3) IGE

43 FDRAT( '0IUCDF UP

',616)

WRI FN 3 ,t4) IGAEN

66 FURIAT C OI N "OCDE DE GEW! D WAARDEN 61 ,

REEDS VERREKENU

WRITF(3,45)NUIWA

15 Fr1R1Ar( 'flNJICF UP

'

₁₂₁₆₎

VRII[-(3,9) 1A)AB

9 FUR'AT( '1

RESULTATfl4 ANAL DATî

bLUK ',IZ

' GFEN

UTuNATISCH 11PE

VLEN',í,

,I)

DU

0 Jl,IANT

ICAFzFC(J)

ICAKDFLW(J)

1(WR1TC(3,li) J,TCARF

CARD

11

FDR.'AT(' KA'AAL ',Ii,' DRAAGÇUI.F 'I5,' HERZ ''

AX ZWAAJ

',l'.,

VHER1' )

WAITE(3,i2) ISCAN

12 FùR:1ttT(

DAANIAL SCANS

ID)

DU 13 J=1,5

13 ALFA(J)=0.O

lk1 TE(3,17)

17 FURMA1HOUPGEGE¼iE

UPNAE GEVUELIGHEDEN',/,Q')

t

IT Fil 1 .)

16 FuAT 'GE

GPi GEVH PER KAN VULDENS 1, '+XXX XXX')

DU 2

J1,1ANT

WR.tTF(1,15) J

15 FURiIAT C 'KANAAL ,Il)

CALL TVFD

REAU(6,16) GK

16 FLR;1T(F.3)

ALFA(J)GH

:krrE(3,u) J,GH

1 Fi]tAT(' VALI KAi'IAAL ',I1,X7Fd.3)

20 CUNIINHE

DUiL=1000.D

RILJ1U003O0.O

(50)

C)

i)

DUP FUCT I 1N Cfl(PLETE3

-

-.

S J GXSAf J)

S1CKXSA( J)

A=ALFA(J)

B A A S I G

S I GXA/L)UI7

S I

S I CKX3/R

_{I L.j} CL S I G / S C A N

('IïE(?',2l)

J 21 FU7SMA f OKANAAL '

i i ,/

O

.RITE(3,22) sic,x

22 FCR'IAT(' SIGMA X'r17X?L-15.5)

WRI ft (3,23) SICKX

¿3 FUk1AT(' SJLA X Ki1+X,155)

ERIIF(2,24) SIG

24 FtJR1AT(' SIG(tA

X

GECORRIGEERU ',E15)

FRITE(3,25) SICK

25 FUMAT(

SIG'A X KW GECURRIGEERD ',E15.)

27

_{.kITE(3,26) GEM}

26 FÛAT(' GEJDD WAARDE UECORR',4X,E155)

WRITE (1,30)

FuFNT( 'CLEF

V

E'i Mu IN GRADEN UP VULGENS ,/ , '+,\XX,xXX XXX

CALL TYFD

PEAD(6,31 I

V,MU

31 FURMAT(F83,1X,I3i

WITE(3,32) V,MU

32 FURNAT( 'O[JrGEGE VEN SNELdE ID

_{F33, '}

_M/SFC

_{INVALSHOEK ',13, '}

_GR

VER')

FCT1 ALFA (i) ALFAC 1) / (RMTLJSCAfl

FCT,ALFA(2);::ALFA(2)/(RILJSCAN)

FC) 3 ALF'A C 3 iL F A C 3) / C N J LJS LAN

FC T4A LEA C ) 'ALF A (4) / C M I U SCAN

FCr5ALFA (5)ALFA( 5)256 .U/ C RiiLJSCAN)

F

35 J3 ,U

CC1 (J)

CL1

1(J)

_FCT1

C UR2 (J) CU:2 (J) ::FCT2

CUR3I J)

CuR

(J)

FCT3 C Ufl4 C ) =CU

(J)

F C 14

CDR5 (J) =CÜR5 (J)

FCT

35 CUa( J)

CL].G( J) FCT6

F C T A LF A C 5 ) ALF A C o ) 2 5 O / CA 1 IL J S C A N i

DU 36 J-1,53

r'r-

.. UI

-.rf

.1\LL 3

(:kK.(J)=CkJFCî

RLTJN

END FEATURES SUPPORTED

uNE -OR

_1NTEG5

CL5F E REUI REMtNT5 FUR RESi

CDMUN 193's INSKEL COMMUN O _VARIADLES ₆₂ _PROGRA) ₉₃₃

EFC UF CUMPILATIJN

C)

Ç)

()

D

(.t

C)

T)

)

(51)

)

T)

C)

n

C)

::

r)

-

.__.__..__.S_. :...,a .. ,....a...,. ..4.a Í..t,b, ..W._.#.

SJSrUTE kES2

DIEtSIUN

1(iO) ,?2(0)P(5) ,?4(U ,P5C)

,Pó(503 ,PCfl(50) ,PKWt5

V)

c: u[N I Z:C F C b ) , i LW. & ) i¶.JF ( C) ) MUCLF ( 6 ) I GEM t ( , I 5CAi, I i ( 50 ) , I