NILU
Teknisk notat nr 22/72 Referanse: IO 000871 Dato: Februar 1972
PROGRA M MULREG Midlertidig program- og bruksbeskrivelse (CDC 3600)
av
Frederick Gram
NORSK INSTITUTT FOR LUFTFORSKNING Postboks 15, 2007 Kjeller
Norge
- i -
INNHOLDSFORTEGNELSE
Side 1
2
FORMÅL
3
...
OPPBYGNING .
2 .1 Input data .
EKSEMPLER
1 1 2 3
3.1 3. 2
Ana·1yse av temperatur- og vinddata fra
Nedre Telemark fra vintrene 67 - 69 . Analyse av forurensningsdata fra Oslo
vint·eren 70 - 71 .
3
8
- 1 - PROGRAM MULREG Midlertidig program- og bruksbeskrivelse (CDC 3600) 1 FORM ÅL
Det er ofte ¢nskelig å finne ut hvorvidt det ad statistisk vei kan utledes en sammenheng mellom en valgt variabel
(predikand) og en rekke andre variable (predikatorer). Dette programsystemet analyserer data for forel¢pig 20 variable
(kan lett utvides), og uttrykker predikanden som en lineær funksjon av henholdsvis en, to og tre predikatorer.
2 OPPBYGNING
Programmet består av hovedprogramm et MULREG, rutinene VAR, REG, REG2, REG3 og MATINV, samt input-rutinene INDAT og INLES.
Hovedprogrammet administrerer de ¢vrige rutinene. Navnene på de variable leses inn i VAR. Data leses inn i rutinen INLES, som må tilpasses hver enkelt kj¢ring. Det er prinsipielt
likegyldig for systemet hvorledes data foreligger, men ved hjelp av INLES bygger man opp et data-array Y (250,20).
Regressjonsanalysen forutsetter at man har fulle datasett for samtlige variable. Hvis det mangler data for en parameter i et observasjonssett, må hele observasjonssettet forkastes. En slik test kan passelig inkorporeres i INLES. Ved store datamengder kan array-grensene enkelt forandres, eller innlesningen kan foretas i flere trinn (konferer kort 1). Nedenfor er det vist eksempler på INLES-rutiner.
I INDAT overf¢res data fra Y til et tilsvarende array X (250,20) i den rekkef¢lge man ¢nsker. Det kan således være ¢nskelig å regne med flere forskjellige predikander eller å regne med forskjellig utvalg av datasett. Den videre analyse foretas på arrayet X, mens de innleste data stadig er uforandret i Y.
(Dimensjonene av X og Y kan lett forandres).
- 2 -
Rutinen REG foretar en korrelasjonsanalyse med alle variable
som er med i X, og skriver ut korrelasjonsmatrisen, middelverdier og standardavvik for de variable. Variabel nr 1, predikanden, uttrykkes ved de andre variable på formen
x 1 = A. X. l l + B. , l
og korrelasjonskoeffisienten skrives ut. F¢rste gang analysen foretas skrives det også ut de inntil 10 beste ligninger mellom de andre variable,
X-= A .. x. + B .. ,
l l] ] l]
sortert etter tallverdien av korrelasjonskoeffisienten.
Rutinen REG2 foretar en lineær regressjonsanalyse med to variable, og skriver ut de inntil beste relasjoner, vurdert etter den
multiple lineære regressjonskoeffisient Rl•ij:
X 1 = A. X . + A. X. + B ..
l l ] ] l]
Rutinen REG3 foretar en tilsvarende analyse med tre variable:
x 1 = A.X. + A.X. + AkXk + B. "k
l l J ] l]
2 .1 Input data
1) NPAR, IREG, NP, NTR (414)
NPAR - antall variable som skal være med i denne omgang.
IREG - antall trinn i denne regressjonsanalysen.
NP - det totale antall variable som skal leses inn til Y.
Hvis NP=O, regnes det med NP=NPAR
NTR - aata leses inn i NTR trinn. Hvis NTR er blank, settes NTR=l.
NPAR<O- ferdig, går til slutt.
NPAR=O- ferdig med dette datasettet, leser ett nytt kort
type 1 osv.
- 3 -
2) INST(I), I= l,NPAR (2014)
Angir rekkef¢lgen av de variable i X, i forhold til den rekkef¢lge de leses inn i Y.
INST(l) angir predikanden.
3) STAS(I,l), STAS(I,2), ST(I) (2A8,A4)
Variabel nr I i arrayet Y (250,20) heter STAS(I,l) - STAS(I,2) (16 bokstaver) og forkortes til 4 tegn i ST(I).
Variabelkortene legges altså i den rekkef¢lge de leses inn i Y fra INLES, i alt NP slike kort.
Deretter f¢lger innlesningen av data, avhengig av INLES. Det f¢lger så nye kort av type 1 og 2 for hver gang man ¢nsker en annen rekkef¢lge på de variable.
3 EKSEMPLER
Nedenfor vises forskjellige kj¢ringer med rutinen INLES, utskrift av datakort og resultatet.
3.1 Analyse av temperatur- og vinddata fra Nedre Telemark fra
vintrene 67 - 69, behandlet av KE Gr¢nskei i NILU rapport nr 26/71.
I rutinen INLES leses det inn data for 18 variable, ett kort pr dag.
INLES:
su iaout rus INI FS(t-lOAS,NP) C I J\I L F ~; T T L M UL ~ F G H Y Il R O,.. DA T A
CrJ,lMlll\! /8/ Y(;..>50,?0)
DI'lF,J<:;I()iJ XX(?O)
1Q1 FOf!M/\T CI7,181='5,1)
1 P 2 F OFl M /1 T C I 6 , H'l r 4 • 1 )
1 0 3 FO I~ M /\ T ( 2 J 8 )
1 0 4 F O tH 1 i\ T ( * Mf R nJN ? 5 0 n A GE p • q F 3 ~, t R ME!) D E 2 5 0 Fø Rs T f: * )
I:1
1 RF:An 10::>,IDAG,(XxtJ),J=l,NO)
P R r "1 T 1 n 1 , I ri A r • c Y. x c J > , J = 1 , \J P >
IF cIOAG,LE,0) Gn TO 3 IF c I,GT.~50) GO TO 1 DO~ .1:1,NP
2 Y(I, J>=XXCJ)
rr <I.EQ.,) I!"'ST:!DAG
I 1) S = I n Ar, ; I : I + 1- ; GO TO 1.
3 NO..:JS=I-1
PRINT to3,IJ)ST,In5
Ir C I • GT. ? 5 0 > p Rn· T 10 4
RE: T lJ II r-.1
Ef\' 1
- 4 -
INPUT DATA:
11' 3
1 7 13 2 o T ~ABR 07 TF07 T PIPE 07 TP07 T rORSBERr- 07 TT07 T LAKOLLE~ 07 TL07 U LAKOLLEN 07 UL07 V ~AKOLLE~ 07 VL07 T FABR 13 TF13 T PIPE 13 TP13 T rORSBERr, 13 TT13 T i AKOLLE,• 13 TL13
U (>KOLLEr,, 13 UL13 V LAKOLLEN 13 VL13 T FABR 19 TF19 T eIPE 19 TP19 T roRSBERr, 19 TT19 T ~AKOLLEN 19 TL19 U ~AKOLLE~ 19 UL19 V ~AKOLLEN 19 VL19
1 167 10 -8 -5 -20 2 167 -40 -59 -50 -40 3 167 5 -3 -10 -40 4 167 -60 -1a -45 -55
5 167 -25 -4n -35 -5o 6 J67 -90-li~-100-110 7 167-125-14~-135-130
a 167-115-13n-135-120 9 167 -75 -93 -65 -70
14
-20 -16 -70 -42 -50 -49 -28 -32 -45
3 9 15 4 10 16 5 11 17 6 12 18
-0 -10-::3-2----=-i:5 -2n ;.- 0- 20· =--4~- ,-_-6_2 40 -30 -- _- 1_ 0_2_ 3_
-o -15 -33 -5 -25 -60 -o -15 -25 -10 -40 -60 -o -o -5 -15 -10 -3o -60 -o -2~ -35 -20 -15 -60 -o 42 -55 -75 -5o -50 -35 35 -65 -80 -60 -60 -19 -46 -0 -1:> -35 -30 -45 -40 -o -45 -65 -40 -70 -46 19 49 -90-110 -95-110 -32 32-lOr-120-110-120 -25 25 11 -90-105-110-10~ -10 -O-lOn-123-115-115 -0 20 32 -90-100-100 -80 -35 35 -5n -45 -10 -3o -7o -o -0 -80 -98 -90-100 -25 25 -8r-100-100-100 -0 30
19 269 -50 -65 -70 -90 -42 -42 -35 -50 -50 -5(1 -35 -35 -4o -55 -55 -75 -55 -23 20 269 -55 -73 -75 -95 -45 -0 -55 -75 -70 -6<; -21 -21 -65 -80 -90-105 -42 -17 21 ?.69-110-130-125-145 -60 -0 -60 -78 _75 -8(1 -28 11 -55 -70 -70 -80 -40 -0 22 269 -50 -63 -60 -80 -23 10 _35 -48 -50 -55 -13 5 -25 -35 -40 -50 -32 -13 23 ?69 -20 •33 -30 -50 -11 -28 -10 -23 -25 -3o -21 -21 Q -15 -20 -40 -64 -64 14 269 0 -1a -20 -35 -49 -40 0 -lll -40 -20 -27 -65 0 -15 -40 -20 -4 9 -49 25 ?69 -30 -45 -80 -60 _49 -49 -20 -35 _35 -35 -55 -23 -20 -30 -60 -40 -55 -23 26 269 -60 -75-100 -85 -65 -27 -50 -68 -70 -65 -55 -23 -45 -60 -80 -65 -74 -31 27 ?69 _35 -5o -65 -55 -28 -11 -2:> -38 -50 -4n -23 _55 -15 -28 _55 -40 -46 -19 2.8 ;,69 -65 -78 _75 -60 _45 -0 -:> -5 _20 -2q -15 _37 -.5 -15 -40 -30 -21 -21
Hl 3
2 8 14 3 9 15 4 10 16 5 11 17 6 12 18 1 7 13
18 3
3 9 15 4 10 16 5 11 17 6 12 18 1 7 13 2 8 14
-999
..
' ..
~;. • • .;,· ,__;RESULTAT 5 (side 1 av 3)
•
'1 '
''
. J~;·'"i'!
·'} t .
I
v-... ... •, ... •
' 1 ·,~ ...
• ~ ~ ~ j
j t"
• ~ J
..
r .
<
lt •
·. 1 ',\.
j
... ,-..o,.
0 .-, V __, --, V
... "'q
'
I"") 0, U'\ 0,. tt\ -0 r-fr-f M .-I..._ :I) I-CV ~0('\
... > ••
f"') "7 C l"'>
' '
r-.o
...,...,.-,r-.
~:::, CO -tN,f'\
i-tr U"I .J.-.N ... > ••
-:r ..,. ., ...,
I I
,,.,,,.., r-,.. ""-0-
.-40, C 00 :C ::_ 0, ,-.. .J..-10 1- •• > .•
N -1 OP'}
I
..., ,... f\,, O, tf\ -:r"
- - .. .. . f'\ :'> J :"'. I"') l,l' _,__.V I•,-,,:._..•. NI":
' '
, .. 0 •:>
,,,,,,r,:-
0 ,0 N .-4!(\:.t\
(_ - _; r.r, , ...
,
::::..
"::' If' ~ "'\.' I
I'
V N
l,,. -0 tV ,-.. ,., -0 -• l."'I ·~ c,,-..,,
i... N .... .Jr'1C- '- N .,.., . ._, • ..,
..
::)..
I I
.,., --1 ti' U'tVV
...,r-. "')r-f,(\r-.
I... t- •
se'"'
•I- •• ..J N l V"1 •:, <7 V
I I
,-.. °"
7 P"? 0, -0 C" .... .:::: ... ,-.. I"".... ,0
cc_,,...~
.... . •I- • • • M ,:r --, ~
I I
. .
0
'°,....,...
00--00 or-.,o...J,
• - ·>
... 0 0
0 lf'\ 0, .-1 0,.
ON.-1...,.-t 0.-IOO_J
• - • ·::>
.-1000 I OC\,Q,0;,,..1-') 0 :"") lt\ U'\ N r-f
~ ;,,,,. O O O J
• • • • • =>
.-t ~ O O -:>
I o 0 r-,. I'? tn CO r-f ,-..
01N .-1 .0 "") co
::)-t.VOOO.J
.•• - - ·=>
.-l ":1 0 0 0 0 I I I O ""l -C -0 a'.. N .-1 0, O "C _,. .-, r-f['\. ...._ .-4 O "') M l"'\ ..-t r-f _, ...J
••• - •• •I-
.... o-::i = ::: o,=, .:=t V O N l"") .O r-.. --l M :>C\a,...,,.,C,.,,J',. ...
~0'- "':"':~~~";~
.-IOOOOCOC
gj~~6~~~~~6
.-, '1:: '() V N .-1 O:::, .-t ...J
• • • • • • • • • t- ... co o -:, -:io c o
I
'= ::n 0- ... co :\I -0 :::c a) 0,
°'
C "", -:,,'J\ a;,-..,<,-,. rl -cJ -i C -..0 ,;; :,,. ,.., (\/ ; ,J ':) ...+ .-! I-
• •••••••• •0-
.-t )0000000·:>
O J\ .,.. 0,. 0,. ...+ N O- -O C"> O, :,,') o-,__.lf'-::J,-...o:,-01"=L"l0'.-!
-:-:i C- r.- o, '.O l") ~ N -1 .-! .-! I- . . • . • • • • • • •I- -! (.!•: eie ~ t": -=, -.. --: ,, C ".ro.;":'\:"': "--: -:0 "": rl -:0.-1,-..
'=> T 'f' .-t O V 'C-,. Lo""") :'\J 0 0 X' r-. !.~ ,:: r-.. :,,') .-! rl O O -t 1-
•
....
-:cc-:-:o'::!o:>o--co
'
.::, '1 .._ - I 4 ~-, "T
~ ... ,. .,_,:,,.
~::-1,-v- T\-<.,...~,0-:,.,...-,,.-1 - C •r, ... ... "'l ":"- "1 "'J "IJ -:,, -:, -:, 1
.-, ,,,..._ ~ -t .f' C\: -t i"") "T) l,C'\ V tf" V "':
-:, -. -:i "" .... _, ::) - ,.. ,., 0 ... .... ") ....
cccca::~c-'t...,. ... -::i C"'" .- :!.
• •••••••••••• •0-
..-1r.or:i~ ~ coo-:,,cooc:
'
..., .,,..II' - ,... r- ... - r (',:.,... "7" ,. ..,. - ,...
_.. ... er-rit·-:: •. , c r",... ...,_ .-. .,.,l c-,.... r
.., ,:- ,.... ~ --:-- r,... - .... ... .., ... .-4 ... '::) 0 Q.
, 1 C" ~ C ':> C:- ,~ 0 C C' C" C r=i -:'J
' '
cr~c~~~ ... ~~~~•~ ... ~~
~"~~ ... ~ ...
~~aO-GMOV~"oa~~~~c,.-~nnNNn~,~
.-(0'=>00-:'.)0C"•Ot:>OC">OOOO
"oo~o~v--, ... ~o"n~~~~-~
0 0 \\I a..t ,u ..-i 1'1 ,.., Q >CJ :;, '(} ,..._ 0 /,,J ' " c., .-4 000-0-ro~~oc-~cn~.-10~.-1~
..-1occoco:conoooooo
'
oc,.m"'~OD~O"G-G~V~NG~
~~~~G~N~-,o~N"wm-0nn
~~~~~~~~~~~~~~~~~~~
"~c~~~-~c~o~ao=~~=
·~Mø--P~-~-M~~-ø-~M~
0..-1..-10-..-10--t=-.-10rlrlOrl.-!
~~~Q.Q.Q.1-0-1-JJJJJJ~~~
,_ _ ,_ ,_ .... ,_ ,_ l- _ _ ~1-~~~>>>
- ---
"
er
N a:, N c,.
II ti 0:
c,.
rl•
__,
rl__,
C)__,
;' ;> 1;- ,.., ,.,,
~ ~
-~o N
0 I
"'}00<00-~000l.00000,0 NC'OV-"IJNOl"'>N000..,0
I I I I I I I I
rl-C-O-vNor-..oc,oN,-...0- 1 I I I I I I I I I I I I I
N
...
7~ -;-'fr;'~~7~7~fr---,oc;-i
'' '
c=oooc.r.0001.t'looo .. '\
"-!NOV-trl.-.fOJ',::,vnN.-t I I I I I rl.-4 I rl I I I
I I I
rv r
t\O~C..,~C.J"'X);(;N-.-4 C!'J"'>CC-ONNVO..,- .-4!'\J-ONO I I I I I -'..-4 I--◄
I I II I I
O.t\0:1'1.00000':>'::>0:f\:t'tt'I '<T.-4"\.'0.;J'"::.::::,~•Cl"':N-,-..."\lc
I I I I I_. _. I I I I
' '
0 0 .-, '(\ 0 ~ ,::: t\ ·t\ "') 0 l""J r:, 0 U'\
"IJ-: ~..,,:,--, ":> _,, "\. ' '\! ~ f"\.' ':'I,-:, N
I I I I I I I I I
C:-C:,Cl")..,-~,..-11"')f\:"J':f'\.-l:C-::N I I I I I I I I I I I I I I
.,
'-' J 2 z
.,,
.u
__,
.D
<
""
<
l
0 w >('J !'\J "1 L" "T rl o CC O .,. o "1" "1' .-.f o I I I I I -f-t I .-.f I I
I o o
l!".r,C")OO,t'IOOOOn:J"IOOII'\
..,.,=, -tU"\~O- .-
◄
n::,, '.J'I "IJ .., .... :'\J .-4 I I I I! I .-.f_..,1 I I I I '\l.,.,.l'\J",,J\01.C'\''
.CCC C ..,,.::CJ'~~7 ~77~~:;rirr-: -;--c;-i,.,.,~
o I I
0 "\ "\ tt'\ !"'1 ':> ':t °=" ~ "" O r"\ ":l tf\ r"\
.-.f-fO :'\-f'=° 0- ,,.. ,:t,·""': N -:> "' J<:J" -' I I I I I I I I I I I I
0 I
0 I
"
er er "
.,,
c,." __, __,
> >
.
;:;
0- M"' "'
0 0 0
I ' I
. . . .
er o: er
II II 0:0 ... M
,.., _, C'
C, ~ C'
..,. """! .r,
V q
"'
0-
"
er
. . . . . .
. . .
~~~~6~~::~~~~~~~~~
~-~~~-~-~~__,~~n~--'---'
~l-1-~=>~~--~>---~~>
,---,,~~~---,-??
r~arc~ ... _ ... ~ ... ~..,
~zz~ ~~zz a•~-
o:ww~
~.u~w o:·uuw*.jjj ~jjj ~jjj
=~~=o,z~~,~?,cu·~?cc
:, !l... r c: :i:: --:: u n.. r: e. ~ -c: :, :.L r -:: -c:: ~ ::o,
c-no•<~-~••••-~••~
~i-~J~~~~J~~~Q.-JJ~~- --J >--1- - J>1- J>
. .
t I I I I I I I I t I I I ~ -~~,;:rU'\O~~~O .... "\l~V~Or-.n J ~
~
rl.--trlrlrlrlrl.-4rl
...
S>rl__,__,__,__,__,__,__,__,__,__,~__,__,__,__,__,__,__, w ww~w~wwww~w~~w~~ww ~,-...
ø=~=~~~I~~~I~~~z=~ ~
<C<<<<~~<~<~~¢<4~< W
--- ~
~«~x2a%xx««x~x«x~«
•••••c<ccc•cc••c•c
►
>>>>>>>>>>>>>>>>>>
z w __,
CD
<w
►
2'1"'>0-~~~r-..~o,r-.l")~~~o-r-..~~~oo~~~.-1..,~.-1~n~~-f -:, -I -f -:>-f -f-, _. .-4 J .-4 -1 ":> -f-' ":.I - -f I I I I I I I
~~~--~~ ... --~--'J--J~
,-.,-.,-.,-.~--~-,-.,-.,-.~~~>>>oG~No~mN~~oo~n~
N -' 1'- "i."••·: t' . .., -:- ..,. ..c-, I"\ ,t• ,:..
~ - ~ ~ - - - ~ - ,-... - ~ r,... ~ 0 ~ ~ - L I I 1 I I I I I I I 1 t I
~__.-r--~--~--t~-'-Orl-'
t'-:,OC''.t\'OO~r?'=-0':::l:>.-.~-:,
~~~~z~z~z~=~NP~~"'~M~~GM~~n-, tz:: '.!'" r' LI.' 'LI I..;.. W LU :..!.! U.: :J.: l!J I I I I I - ,.-! .-4 I I I I I
W~W_J_J_J_JJ_J_J_J_J I Io
mm~__,__,,__,__,__,__,__,__,.,,no~~o"'"'~n~~o~a
zer~~,,~~~~r~-~??~o- .
!J !J Il .L l'1.. :l. ,. C ?!'. !I:: 'C 'I! V: ·C: 'C'. C::: :C ',!.., 1"\-' ,r "')·~ ">•" I•(.' '1J ") •¥ 1•C• , ...,
"'! 1:-:t---~-::,-=:-:t<t<-..C<r<C<C"<:I I I 11.-4.-4.-11 I I 11 I
'"-~~~~~,-.-~ _ __,__,__,__,__,__,__,J
o IoU:t'C'l"'>r:)01"')1,(\0V}4,.,.,C'l"'),.,.J"\
1-1- ~-~---- -~-j JJ>>~ .. .•..• ••••....
O lf\ C' r-.. «r .-I «TI"') :-i- r-.. <"') '<T t"\ q- rr I t f I I I I I I t I I I I I I I I I I I I I - - - I I I I I 1
I I I
~Nl"')«rm~-~~o.-1~~«rn.o,~o~lf\~~olf\t\nno~.J"\n~
I
i" ~-" .
-''<TO-ONC,,,N-'~.l'\C'.•....,<"')~"\J JJJJJJ..JJJ..JJ~JJJJ~J I I I I .--t.--t I I I I I I
w~wwwwwwwww~~~~~w~ , ,
O~CD=~~~~~~~X~~==~~-,~---~,~ft~~,,~
···••<•<•cC<O~GGGGOC~~-~Gr~
- - - - - - - - - - - - - - - - - - - - - .-4 ..-4 - -t - ... - - .~ .... -
~~«~«~~~xx«x~«2xx~c~oooccoo~=o:>o~
••c••<•••cc••••••c-NM•"'e~ocøo-~M•"'
>>>>>>>>>>>>>>>>>> rl.-.frl.-4.-1-f
U"\U"\,-...Ol"')~U"ll"')OOONl"'),-...c,l"'>'<TO-l"'>.-40-rlO
~~~~~~~~~;~;~~~~~~~~~~~
I I I It I I I I I I I I I I I I I I I I
M~NO~~"'~OO~N~NONWC,.~VG ... 0
~~;;~~~~~~~~~~;~~;~~;~~
t I I I I I I I I I I I I I I I I I I I I I I 01.C'\OOU"llf\OOOOOOlf\U"I~~~~~~~~~
;;~~~~~~~~~~~~~I.C'\~NV,0..,-~0 I I It .-(,-4.-1.-4 I - I I I .-1 I I I I I I I I I
I I I I I I
t(\l,C'\01.C'\lf\OO~O~OU"\lf\OOCOOOOI.C'\OC
~~~~~~;;~;;;;;~;;;~;~;~
I I I I rl -t ..-i rl I I I I I I I I I I I I I I I I I I I
l"'>U"llf\~U"llf\lf\lf\lf\O">Olf\00~~:f'\O'=>rolf\O
~~~~~~;~~~~~~~~~~~~~~~~
I I I I I rl.--t.--t I .-4 I I I I I I I I I I I I I I I I I
0~"'~0000000~0~"'"'000~~~0
~~;~~~~~~;~~;~~~~~~;~~~
I It I I -t.-4rl I I I I I I I I I I I I I Noro~N~~~~ro .... r-..U"\M.-4~.-4U'\~~U"lr-..o
;~~~~~~;~~~~~~~;~~~~;~;
I I I I I I I I I I I I I I I I I
øtt\mtt\Ntt\~~ .... ro~tt'l~ .... ml"'>.-4,-...tt\~l"'>tt\=>
~~~~~~~~~~~~~~~~~~~~~~~
I I I I I I I I I I I I I I I I I I I I I 00.,,00~~0000~0"'0"'00"'~000
~;;;~~~~;~;~~~~~~~~~;N~
I I I I rl .-4 .-1 rl I .-.f I I I I I I I I I I I I I I I I I I
Olf\0~01,C'\U"\Olf\001.C'\OO~OU"\O~OQOO
~~~~~~~~~~~~~;~~~;~;~~~
I I t I rl--1.-1 I -t I I I I I I I I I I I I I I I I I I
~~l"'>I.C'\~OOCOCU"l-,01.t'l~CC~~~~~~O
;~~~;;;~~;~~~~~;~;~~~~;
I I I I I .-4.-4rl I .-4 I I I I I I I I I I I I I I I I I
ooooo~oooooo~"'o~oooc~~a
~~~~~~~~~~~~~~~~~~~~~~~
1 1, I 1 -t --t .-t I rl I I I I I I I I I I I I
I I I I
~~~~~~o-~arPNrO"K•~~-~c
~;~;~;;~~;~~;;;;~;;~~;;
1 1 t I I I I I I I I I I I I I I I
No~onN"'a:>~a-GN~OMrlO-~~a:,"'o
;;~~~;~~~~~;;;~~~;;~~;~
I I I I I I I I I I I I I I I I I I I I I
onaao"'"'"'o"'ooo"'"'oo~ormoo
~;;~;~~;~~~~;;;~~~~~~~~
I I I I ,--t.-.f.-4N.-4-1.-4 I I I rl I I I I I I I I
~U"II.C'\U"\6~66666~ott\~OOOOOU"IU"l0
... : ~ ,..._·,.,; s.: "' . .;
.-4··"\I· C:-0"1,-...· ,...,_·""· ~,.,;.,;,,;.;; .;,,....·,;1 t I It I rt rl '\l ,.... .-4 M I I• I .-4 I I I I .-4 I I I CCl"'>Ott\6~~~~~~11'1")61"'>~~U'\~OroO
~~;;~~;;;~;~~;~;~~;;~;~:
1 1 1 t ,....,...,... I .... .... I t I .--t I I I I I I I I C\I
I I I I I I C.
~~~~~~~~~~~~~~~~~~~~~~~:
OM-0~,0l"')~~c:,o~,;:rtt\l,C'\.-.ftt\NO~~~-oc 1 I I I M....t.-4..-irl I I I I,.... I I I I I I I
I I I I I I I"--
~ ~ ~C,. ~~~, ~ C,. C,.~ ~~C,.0- C,., C,.~ C,. 0 C~
~~'°'°'°~~~~,o~,o~,o,0<0-o-o-o-o-oc1.-4
NN~N~NNNNNNNNNNNNNNNNN C
0000000000000000000000
~~0-0....tN~Vtt\-01"--~0-0rlN~•I.O-O~~
rlrl.--t~,...--..-t.-1.-4NNNNNNNN~
,,
O C,.
0 ..0
~~ ...
0 0 ., 0 0 ..0
~o:t:
0 0, cc,"
0 0 (0 0
~°'. ,..._~
..0 0 0
0 ,0 -0 ,0 0, 01"")01"") ....
0 0 ..-4 O J
•.• ·>
.-t O O O 00-N1"') ...-!1"') O .... ,O N V .-1 OCOOOO...J
. . . . ·>
T""I O C O CI
' '
o-ør,..oo- o,.._
0 0- ,0 -:0 ,0 ,0 0
~";~":<=:~;;' .... 00';"7 ~
OLl"I 0- .... V :X: V,-, ':::)f\.!,-11")<1" ":)Y',-1 o .... ooorlN..J
...
•::)....io 0 0 0 ::- -::i
'
O N O ,o,-..II' :0 ...-,I"') 01">1.t'I :r.N .... ,-.. ..0 ....
0,-...00-::,.,....,.N..J
~,;~~~,;,;..;=>
''
O r-. ") U\ ~ .... ,01""),-.. r-,.
ON.-..01""':.0N..CO-O o,ov,::,e-~•..,...,-,..J .... oc~o":>eo-:>
' ''
~ _,, ,O O CO N .... .-f .-t 0 0- 0 <.0 .... .... .-f['\.,....J" 0 ~ .-, c:>•')1"")1"')-1 .... .-tC':J.C..J
O V O N ,.., ,O r,.. .-1 .0 0 .... ,"") ONr---0,.01""),I\0,0,0i."\...➔ C'">.,.l'"'lNN ... rl,...._'t>:C..J
. . . . . . . . -~
.-t O O C O O 'j O C'I C' ~
O o:f" tC l"'l ::> .t:; ,-.. .... a-GO- ,. .. 4:-...
O tt: N O -:r O'- C'\, '<!' N l' O V O O:OCX::VN .... O::.irl..CO-~..J
••••••••••• •1- .-IOOO':::)":><:,O::>':>OO
'
-, 0-,.. r, 0 '1 ::> 0 'O.,.. -., _. ") -,.
O M O- ..l\ '0 ,-.. £'1 0- "->"">_._.
0 "'O "'O -,. _,, :"J '.'J O -1 ,0 0 ::,. '-
OJ"lr-'>0-0-.-tNO-DO?-O',...:,-.r/"'7 O":J-f~C-...-:=>':lf\-Y..C'\1-... •
COO O CO ....
t'\;".,...,... - ,-..
0 C: ,._O""':N...._"-1,_,.-irc,~,-1"t" .... Tv ..._..._
C"I T ')' lr, •-. .C O:: ø:' .:- lt" f" C\J ~ ,-1 .... (:' 0 CD r-.. 0 :X-. ,... "'• .-1 ,-1 o C .... ?- 7' :r '-
.-IOQOOOC0:::>0000".:>l:)
'
oN,...0--1...-,vovo-,...o-~~~c~
ooo-v~ro ... 0-,00~~~,o,... ....
amvo-mm~•~Nao~--~~
..
,·-
.-. 0 0 0 0 0 0 ,::, 0 0 ,:, 0 0 0 0 0 O tt'l ,-.. v rl tt'I N rl ~ ~ ~ v co ~ ~ o ,-.. ,.., ooom.-.-u-,o,"u-m~u~o- 0 0- 0- ~ co o- o- co M N rl O O rl m o- o- ~
... ·-
rlOOO OOOO~ C~OOOOO O
'
o ~ ~ .-4 ~ ,-.. c,. ~ ,-.. N rl v ~ o- o, ~ ~ lf\ ,.._
Orl~,()lf\Vrl,00,-..,.._0"0,-..,-..NrlO o o- ~ o, co ,-.. o, ,-.. ,-.. ~ .-4 ~ M o o o- o- ~ ~
••••••••••••••••• •1-
~co~o~~~~oo~~~,::,~oo
' '
""'•""'°'"~•,-..,D.-.~o.-.n•
OrlrlOrl~OrlrlOrl~~rlrlOrl~
~~~l-1-1-..J...J.JJJ_J_J_JJ~~~
... ~jj>>>-1---
. .
.
:c
:r"
r.
I '::;"
:c
...
r:-
uN
"·
"
:,
.,
0
a ..J
..,
°'
,n
°' "' 6
0
'
7
'
0
'
0
" '
"'
0
'
"
"'
(side
"
a:
"
0
....
0.., '
"
0 ,-. 0... .... .... . .
,-.
,,,
.,, 0
N ,-.
.,,
,-."
0
....
0..2 av
. i
3)
... .,
N o
.,, 0
0 .,,
:,
'
0 0
' '
,-.
" "
0 O 0
~ ~ ~
~ co. .
0-.
,-.
., "'
• 0 c,.
., c,. .,
., .,
-
..... ...
...J.... .
" "
"' .
. .,
N ..0
0
....
0..:,
' .... .
0 0 N c,.
.. .
,.., '°
0, 0,O 0- ,.._ CO
Øl ,.._ ,.._ C0
°'
0,. 0,°'
II II
a:
a: :c
.,
,-.... 0
... ... .... .... .
,-. ,-.
... "' " ,,,
"" c,.
" c,.
:,
-
... ..J
.... .
u "
"' 0
"' "'
,-.
...
,-.
0 0
0.. 0..
...
rl :=i tJ"I J
r') N ,c s:
,() ,() ,0 C:
N N .-4 C' 0 0 0
' ' '
,-.
,,,
0 ...
... ... .
., .,,
"' 0 c,. N
"' "'
. .
. .
,-.
0 0..
....
0"
....
0.. 0"
0.
....
".,, 0 0 ,-. "' ..J .,, C,.
=>
n,;~
'
"
a:
0
'
0
'
a
'
0
'
,-. u
N 0
,. , ....
:;; .
0 0
0
...
:, 0
0
...
"
0
...
N
"'
"'
V
C,
'
0
'
"
,-.-• n 0..
:- ~
"'
0 -
.,
0
0
-
,.._ ,.., o, ,..,. ,.., c,. ,-.. ,.., c,. ,-.. n o- ,.._ ,.., o, ,-.. n o,.
~ rl ~ ~ ~ rl ~ rl rl ~ rl rl ~ ~ .-4 ~ .-4 rl oan1-1-1-..JJ...JJ...JJ...JJJ~~~
... =>=>=>>>>1-1-1-
c:;c,r.;,zzz;--zzzzz
•.-n::c:wwwwwwwwLL!
UJW UJ J ..J .J .J .J ..J ...J J .J lD lll :!l .J ...J ...J ..J ..J -1 ..J ..J ...J
ww~~~~ooooaoooo~tt~
~o..~•~•••c•••~c•••••
---~o~<4<44<<<<~<<
~ ~ ~ - 1- 1- ...J ...J ..J ..J ...J ...J ...J ..J ...J ~ ~ ~ ... - .... ~..J~>>>l-1-1-
1 ff I f f f f f I I I f I I I I I
" N ~ V lf\ -0 ,-.. CO O- O rl ~ ""'l V IJ"l -0 ,.._ CO rl rl rl T"4 .-4 .-4 rl .-4 rl
...J J ...J ..J ...J ...J ..J ...J ...J ..J ..J ...J J ..J ..J .J ..J ..J
wwwwwwwwww~w~wwwww m~x~~mææm~~=m~~~m~
<<<<~<
◄
<<<<<<<◄
~<~~~~~~~~~~x~~tt~~~~~
<<<<~<<~<<
◄
<◄
< <◄◄
<>>>>>>>>>>>>>>>>>>
,,
0"'
f"') 0- -00--0N rl ,.._ l"")rllf\lJ"I
~,.._ q.::<'-f'(J.
"") V N ...-
0 '
,.._ ,.._ 0,""') r-40, 0 .-4 V .-4 ,.._ M
~ ~ ~~"'- ~
tt'I '<:' N q-
0 '
0, ,;;r l(\,-..O,V .-4 C ,;;r e r,. CO t- N ,.,, L.. -0 CO I- • ·t- ••
¥) q- ,.., V
' '
I"'),:,,. IJ"IO-U'\-0 _.. _. "'l
--.r--
"O- C"" .... ...JC'{'\
I- • ·>.
I"? "T Cl";
' '
,.._ 0 .-t M rl r-- 0 OCOrlC'\;ll' ':,i: ... 1/'I J\ ..J .... C' t- • ·> .
> q- ..:- C ""
~ I I
...
r=,, 0, ,..., ,..,, ,.._ ,0 {l :t: .-to- 000-:t ,c( a_ 0,. r,.. ..J rl C C I- • ·> .
4:. "' ..,. ~I"
-,
'
':'Jr,.."/U'tt'II .-4 _....,-.-tlf\1 Q.. M 0, ..J rl '
~ . ·== .
") V '\I '' '
,.._ 0 ,() ,., 0 ., ,0 :'\/ ... './\
Q. ,..1-::, t- ::, ..
~ ,n N
' '
'"'RE•SULTAT: 7 (side .,
0:
....
0
"'
00 I 0:
N
"'
N C>
3 av
"
:r
- . . .
. .
<C
...
0
"'
u
'
3}
....
0
:
..J';
:·: t ::
. "'
. ...
. . . . .
.., .,.
...
..J'; .
...
li'\::."\
"' ... ....
.,_
0"' .
"
0X,
....
0
. . .
"
:r
:r "
er"
0
' ...
0 0 ..J ..J
'; ';
... ...
.... "' ....
...
N C>
C>
:r
"'
...
0... "'
N N
.... "'
0 ...
!>.. c..
~ ~ . .
.... "'
0 "'
"' ...
.,. 0
·~ ....
... 0
..J ,_
~ .
...
<>,..,
,,," 0 0 "
0
'
0
'
0
'
. . .
I'
OC>
O .-<
oa.
0"'..,
0 0 .-<0
.
C> ... ~
OU'\,.._ ,,0 0- 0 .... 0 ,.._ .-4
~°:°:O:~
.... oc. o 000-0U"'M OON<0-0.-4 0 0- 0- 0, (I)~
•• • • •I-
V
....
ro 0
:r
':J'" I(\ 'D -I'/'\"'
~.-< V :::,
..
N.-, I
0 0 0 I
....
0
..
0....
0. 0 0..... ....
0 0. 0 0....
0.....
0:
,.._~~,...~~,-..~~,...M~,...~~,-..M~
-~-~~~,--~--~--~-- ---LU~Q~~
-1--1-1-1-~=J>>>i-1-1-1-1-1-
,~-~~O,...-0,~c ... ~~,...~~
~--c--o--c~-~--c-~
~:, ,- :..= •••• z. :: :~ "? ~ ;;· :...:
'.C IC o"'. u, - U.. 11 U..' .lJ :..U i- W :.. • ..t· .lJ ..J -· .... _J ..J _J .J ...J ..J r...::: c.. ...J _, J .; ...J _, ...J _,J ...J
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ :r ~ ~ ~ W .I.<
- ~ ( 'C -t". '! C "C < ~ < 'C ".J U ~ Q. •.l.. t1.
~~?••···•·---
---.J..J....J...J..J..J..JJ_J~~~~~~.-.N~v~o,o-,.:::,.-.N~vn<>,.._m
rl.-4 .... -1 -I ..., .-1 ...
...J ..J ..J ..J ...J ..J J ..J J ..J ..J ...J ...J ..J J ..J _J ..J
w~w~~~~w~~w~~ww~w~
-~~~----T~~z-=~~~m
4 •••••••••••••••••
---
u,
;, u
-, v, .
>
w
0:VI u:
,.,
.,.
0=>
L
w
>
....
..J"'
C
MO CC:>'..,,.._
-' J\ U"l .-4 :J- -:) ....J CO ~-:,,. ,.._
~ • •I- ••
"'N
n.:v' '
,.._ ~ ,c...,
r-. N,., ,.._ -0 .-4 ..-1 -T ...J _,,,.:,,:;..MO- :, N "-'
. --
"') ~..
' '
:,- t\J <::" ,._ C t:;
.... :: ... C) c-.,~
.J " J "'J Cl.~ .--4
- • •I- ••
""'> '"' 0!:' o N
-~ .... ...., -1"
L"I- ,.... ~ 1 •• {'\, <q
,....!"--
C" ..., .-1 0, ,:'> ... q'" ... :-,. ...,~..,,. •~ - V :;
,_
.. - ..
o,~u-,,,-..ø-q, .--4 0 ,;:- O ..._ a) I- N ""'> - ,o .J)
,_ ...,..,.., . ·- ., .. . ..,.,
' '
,., V, 1,n V, Ul '4' ... ,., ... ,.... (0 ,- C> V -' C.. N I- • ·> ..
,., -:i- 0""}
' '
,.... 0 .... ,., .... ,....
~ =><0..-1NII\
1- 1n \f\ .J ""'4 N
- . ·> ..
<q V 0 ,.,
0 0, t0 II\ \I\ -0 r,..
~r:,:;o,-...-ONO oo- r,,,.o, cc,.... 1.o.
••••• •l- ---1 0 0 0 0 0 0 0 .c-0-:,,.-:: 0 0- ,-,,,: r:, ... ,.... ... .,, - 0 ~ .--4 0 0 .-4 0 _.J
...•. >
.... ':) 0 "':) 0 0 0 00-'\ll').-10,.:00,.1") 1':1 .-4 .C N ST ":>J'l...-4 ...-4 c cc o o o o o o .J
... ·>
--: 0 0 0 0 0 ,:, 0
' ''
o ,o r,.. o o, c ,o -:::r ...._ ,-...
co -oc:r.o,o ... -or'.>
or-,. ,0 .-4 0 0 .-. oo ...J
....•.•.. >
.--4 ": 0 0 0 ::::. "= :,o I I I I I O;(\ 0, ..-1 V :0 V~ J\ 0, 0, O C\. ..-1 ,-, ST O !""; O -C ,C .-.
':>.-I:::,. '::> ':> .--4 ('"~ .... .-4 ('",: ...J
••••••••• "='
..-1 0 0 0 0 0 0 0 0 0
'
O N O .C. r,,,.t{\a)...., ..-1 ,0 V I") r-, M !f\ \l\ r'\J _. ,-.. ,o r,.. ,..JQ, -1
o,-...ooo
...
... -1N-1NNJ ::::, ..-1 0 0 •.:, 0 0 0 0 0 0 -::i' '
=
,-...,-.;l!\r.o.-.,o,.,,-..
Nr":: o,-...C. "J.--4 -0 ""} 0 N -O ~ r--. o, .... -::, 0 .O ST 0 0 0 1"} ,-, l") l"') I"': ,q J
..••••••. - • ·=>
-1 0 0 0 0 0 0 0 0 0 0 0
'''
-:> "':' 0 ." ; "0 ",,1 ..-1 .-4 -f :> ,.... .... ~ ,,..
0 ":: .-1 .-4 .... 0\lr,.. l{\ 0 1'\'.) r.l V -t .-l
o "'> ""')..., .-4 .-4 .-4 0 :0 '2:1 ,-... 0:,. .J
...
.-•<:: C: C-·c.C· C c,c C, C c, c.
Os:f"ON,...,..Or,....-1,00.-40,.('J"')-, --::iv,-,..:,....,n-,..-=:io"'o-:,0-1 C O"l f\ f"\i-..- .-.,-..a: r:: ,-.. 0- CC_,
•• . • • • • • • . • • - ·l-
.... ':> ~ 0 0 0 0 0 0 0 0 0 0 0 0"." et: M 0 -0 ,-.... .-4 0 c:) 0- .... 0-:t'I_.,..._
C: .r f\. ,:: , '<:I' :) (\. ,,: N ,-.... C -v· .-1 .-4 <q O O 'T"l «"' ~ N - 0 0 .-4 '%) ~ <0 0, 0, (0 ...J
... - . ·-
.-I C> O U O O O O O O O O O O O
'
C :0 0, ..., T N -O <O C0 0, 0 ---1 ...., ,-.. .-4 ? 0,.
0 .-, ? ~ ~ ,-. G D .-< -O D .-, _ q ,-. q .-<
o m ~ ~ ~ N N o .-< .-< -O ~ ? ,-. m ? ~
011\M~O,~N0,,000,0,"')(0~v,-..~
o o ... ~ ru ,-. v o o ~ ~ • N ,.._ ~ ~ ~ - o~o-o,m~NN..-1.-1.-1,-...0,coco~co~
.... - ...
.-.oooocooooooococc 0"')~,-..~~""'><01"') .... CO .-t :O v ,-.. .-4 , N ,-..
:>t00-~ ... -0~00V\l\~NN.-1.-1-0 on
~~~~~~~~~~~~~~~~~~=
... ooo~c~~o~oooooo~o,.._~~~-~~~~,.._~o,-.~• .... ~~
~~~~~~~~~~~~~~~~~~
"
0:
....
0
- 8 -
3.2 Analyse av forurensningsdata fra Oslo vinteren 70 - 71. I rutinen INLES leses det inn d¢gnverdier for S, SO2, sot, relativ f.uk t i.ghe t og middeltemperatur på Blindern, og nye parametre SO~, SO2 • sot og SO2 • fuktighet genereres. Ett kort pr dag.
INLES:
SL~ROLTlNE lNLFS(NQ~S,NP>
C 11\ I. f S T I L t·
1LL R 1: G S _ u AT A OS L 0
CCMM □~ /8; Y(2~o,2n>
I:O ; F=0.1
Pn If• T 101
101 Fcru-1td (1HO)
1 REA[ 100,IP,IM,IY,~,SO?,SOT,HH,TT
100 FCR~AT (3I?,F4,1,3F4,n,F4,1) IF ciu,L~.G) Gn TO ~99
IF (S,L~.o) Gn TO 1
l = f + 1
SC?=~u?*F , SOT=SOT*t
Y<I,1)=S ; Y<I,?>=~02 ; Ye l,3>=SOT i Y(I,4>=RH Y(I,5l=TT : Y(f ,b):So2~s02 ; Y(I,7l=S02*SOT Y(I,o>=SO;,:>*RH
PF;If\T 3J0, ID, !"1, IY, <Y< I ,J) ,J=1,8)
Jon FCH~Af" <2x,3I3,10FQ,2)
Ge ro 1.
999 NCHs=I
HETt:1-d-.
E: f\ [)
9
INPUT DATA:
·8-;.·::J-:-1 ' -,,~
·1 ... 2 - -"~-~ .... 4 5
S B,IsKEcHY,so?. ..
ser,
·;~._. .ELJNCERN"'·
Jl"El"1tJDERN
·sc2·•so2
SC2•SoT
·s,sc?.
· - ··:S1,:1u 101270
• 141270 151270 lb1<70 171270 221<70 -~--23-12-'io-
241270 2
1.71 '3 171 4 171 5 171 6 171 16 1 71 17 171 29 171 2 271 3 27t
4 2715 271 S 271
7 n18 271
·i
:>71 1 ') <'71 17
l7.i.1, <7l
1 ;; c. 7 l2J '71
:?
1 l7 l~? t.7 1
? ~ C 7
j
,? .. 271
J71
96 ~71 1tt4
rl?
,4487 131 133 90 17
58 148 1401,
-1?7J 120 118 O?
-~l83 033 66 95 -t?
'l5 '1"10
r
io
01 1J60 104 55 ,1·-42 -27-1,n·-15--77. _.,s
51 16A 55 d4 -R4 51 •91 33 94 -05 96 17/l -!33 96 -,2 43 86 55 70 -73 78 170 118 67-ltj 77 <'37 85 a9 -78 71 115 r,q ~4 -,2
51 65 44 97 -11 55 70 77 75 -72 34 95 66 56 -22
sn 120 125 76 ?2
58 i51 140 A6
?I70 i72 14J 7~ ?7
~5 149 125 6~ 2n
43 115 66 ~Q ,~17 32 22 7J nl
o?
92
~Q ~?-?S
9~ 65 :55 ~.; -13
·i:
J9 ~~ 1,1 nn34 3·1 22 ~G nl
11:.,3 :n
.;ii -1
o48 jl 22 ~~ -,u
54 ~ 3
0·1 ,,,::-;s4 71
,,
5-:; 7~ -,d<1 46
:;;: ~£'1
s sn2
·SOT run TEMP S022 sso1 5/52
(.':'.
:"Vi :; 1 - "!:)
·,r-,_
l •r;5 i71 J; 1i2 ~; ~~
-7h!J ~lt .:~ 59 '3 -t.C:
J71 ~~ l~~ 77 7' )
-~ J7.i. ~,: o J
5!
,t;0
~c~:2.~ro;vE1 ,.:-.~ □ 2
.,. ..
soj: .,., .:.- . ;.;;. "' soT
Re• EL I NDERN FUKT'
T C<ISL'l NDERN TEMP -
se 2, s02 ~- S022
s02,soT ssoT
SISC2 s1s_ 2 __
. 24 "27i- 22 .. ~2r 26 50 -no
·ro J71 38 53 7 55-lOO 5 371 16 47 21 53 -70 .l71 33 . 35 12 a:1 -4<'
371 32 76 26 70
o Ia
1 3 4 5 6
10 1,r7.-ro3 69 94 -5<
17 171
8635 133 97 -1 l
;,8 171 42 123 101 66 -4 4
2·;171 20 86
t)4 7'}-7<
2 271 42 156 155 56 -:>2 3 271 il4 170 31c 76 ?<
;c
4 ·211 •82 193'294 66 ;, 7 . 5 ··271 147 270 406 7u 21 .. -·6 271 70 177 190 66 20 7 271 54 186 164 69 l,
8
271 ·27 78 64 73 n<
9
271 61 128 107 02 -;,3 10 271 152 98 114 98
-1 J17 271 ?.9 61 45
96 nd1a 271 45 75 48 90 0l 19 271 85 110 89 90 -16 20 211·~ 84 73 48 96 -1u 21 271 - 65 113 74 95 n2 22 271
-·43 106 _126
7? .-..QO23 271 39 75 ·79 8? -no
;>4