Die ausführliche Berechnungsmethode der sekundärenZerstreung in den Planetenatmospären durch dieRechnungsmachine

Sci. Bull. Fac. Educ., Nagasaki Univ., No.27, pp.133-148 (1976) 133

Die ausführliche Berechnungsmethode der sekundären Zerstreung in den Planetenatmospären durch die

Rechnungsmachine

Takao Sato

Fakultät von Erziehung, Nagasaki Universität

(eingenonnen Oktober 31 1975)

Ubersicht

In Bezug auf zwei Planeten d. h. Erde und Mars, hat der Verfasser nahe mit keiner Vernachlässigung die ausführliche Auswertung durch die Ruchnungsmachine veransch-

aulicht.

Einleitung

Der Verfasser hat im Falle der Erdatmosphäre über die Veränderung der Luftdichte mit Höhe, die von Humphreys angeführte Wert vorgebracht, wobei ihre Feuchtigkeit d.h.

Wasserdampf in Betracht gezogen worden war.

Dagegen im Falle des Marses hat er selbst die Formel dafür folgendermassen abge- leitet.

Bezeichnet man nun den Druck, Dichte bezüglich der Luft und Schwerkraft bezogen auf Boden mit po, Po und g, dann wird die Meereshöhe der homogenen Atmosphäre in

der From

H — P°

Pog ausdrücken läßt.

Mit Hilfe der Untersuchung mittels des Rocketzeuges kann man den Wert von H aus die Formel berechnen, zur Folge H = 8056 in Meter. Daraus kann man weiter die vertika- le Verteilung der Dichte mit der Höhe folgendermassen geben

P=Poe'H

wobei x die Meereshöhe in Meter ist

2. Problemstellung, Zweierlei Fortran Programme.

Die Verschiedenheit zwischen beiden Programmen zeigt sich ausschließlich in der aus

irgendeinem erster Zerstreung bestrahlenden Punkt, namentlich P in Literatur 1, gesehen-

en Richtung, in welcher liegt der zweite von der direkten Sonnenstrahlung bestrahlenden

Punkt, Q in dieselbe Lit.,

Es sei nun 0 der Mittelpunkt des Planeten, und sei 01 der W'nkel LPQ und sogar 0 der spezielle Wert von 0, wenn PQ gerade der Planetenoberfläche tangentiel berührt.

In Programm A hat der Verfasser als 01 die mit 5 teilbaren Werten zwischen 0°

und 180° gewählt

In Program B dagegen hat er mehr verwickelterweise die mit 5° teilbaren Werten innerhalb zwei Bereich 0°-80° und 90°-180° und noch hinzu 1(80+0)0,1(90+0) 2"2 innerhalb 80°-90° gewählt.

Diese zweierlei Programme sind folgendermassen gegeben.

(Unvollendet Fortsetzung folgt)

Literatur 1 : On the Polarization Angle in the Scattering Problem. Journ. Met. Soc.

Japan, Vol. 34, No. 2 (1956)

Die ausführliche Berechnungsmethode der sekundären Zerstreung in den

Planetenatmospären durch die Rechnungsmachine ₁₃₅

ABTEILUNG VON LF

AuSwEhT0,,c, VON NACH DLR rERECHNUNG VON DICHTE U.S. 44

!-,E1 ,=1,KI=1, Fut-0 Gt7SAMTE KPP

nAs HElzZT DIE vFRKL,FRLuNG VON CO0u44 IN KOKUDMI UNIV.

INIEGEK ZE(7)0.(3),,1T(42),AS,15(91),AF,ISSSt3,9,7.7,12)

REAL MmT(100) ,L1,L2,1.3.M101.2,M3,N1,N2,N3,PPOW(12) 1 DDLi (129).HF(129),F1(129),DLt:(65),DF(65),ALE(65),

1 LJDLS(124),-)DHS(129),!'2(129),DLS(05),DHS(65), PSKLJ( 7,9112.12), 21)(65),D1(u5,65),,(65,12),SK(12).SK7(91,12).PSKZ(91,12),SL(7,33,12

3),S(7,33,12),SSS(j,9,(.7,12)

4, mF(65),mb(8)) ,DHF(65) ,ALS(65),LE,LS ,CV(65)

5, DDL 1(129) .,..)HT 1129) ‚MT (65) 0-3(1.29) (6_5)

6, IR(A1J),HH(12),PK(12),VM(410) HO(42)

S(X)=SIN(A) C(X)=COS(X)

PAI=3.1415927

HAp,=( ii.o*P41)**1*6000.0 PAkID=9.0/PAP

RLAD(5.20u) (PP(L),L=I.12) 200

100

4567 300

66b 669 6666

111

1

FOKHAT(12ho.5)

READ(5,10D) (PK(L)•L=1,12) FOKMAT(6F_I1.5)

DO 45b7 L=1,12 PHwü(L)=PK(L)**2*'-,pr, CuNT1UE

eAD(5,300; (IR(.N=1.81U) FOkMAT(10I-L.(1)

Ir<(H01)=4u(

Vm(1)=o.0

*RITL(Attir,S) (Pk(1.),L=1.12) FOkmAT(1H +12F1u.5)

,RITr(6,6(-)9) (PK(L),L=1,1?) FuhMAI‘1H .6Ele.5)

wi-dTE(6,h6b6) (1ht(I),.1=1,810) FUkMAI(10F9.u)

A=6370,0 ii=6410.0 CC=1G37.4 PA'0,0

1)0 111 K=1,79q,z PA=PA+1R(K)+4.0*If,(1,4.04IK(K+2) KM.(K+0/2+1

Vm(KM)=PA*1,.05/(3Lu0.u*Cc) C00/1.uF

vM(402)=1.0 vM(403)=1•() VM(404)=1.0 Vm(4u5)=1.0

P=3.14154/1bn.'.

PImP IR^801)=40(

1)1) 1 I=M(12,R11 IR(I)=0 CONIINUt

00 9("1 m=1.3

"o0)=303Fm w=p*htnAl(H(m))

SH=3(w)

chi=c(,.)

601

60c

2

5

50 60

wRIFE(6.6n1) FoWMAT(1H1) b0 99 m,1=1,9

wg1TF(46>.-P0)

1=5*K1-4- hT(1)=1-1

H0(1)=A+FLUAT(1-1) DU 4 Kt•P=1\.7 AF =30*(KIJP-1) Al=P*FLOAT(Ah)

PSA=S(;(1) f,cA=C(xl)

J=1,7 dFC,J;=1,..•*( 1-1) Li=r*FLu,T(LF(J;)

Th,T0*, KSTrI=Siih) KIhi=v(Th)

1i=c)(1.E.(9.1.0.J.LE-;) (-,ü To 4(, 1F(KI.U9.9.ANi.,.J.,,F.',) CO TO 4C.

IF(J.LE.J) CO IO d IF(J•GT.i) GO Tu 3

LF=-SW-g(A**2-(o()(i)*S(2.1))**2)+HO(1)*C(Z1) GO iu 5

LF=+SQkl(**2-^.mn(1)eS(ii))*?)+HO(I)*C(71) 01=h-

f-A,. 6 K=1,33 upb-(()=IF*FL00(,-,.)/32.0

U0HF(K)=SLRT(H0(1)ilild+PnLF(K)**2-2.n*H0(1)*Dt,LF(K)*C(Z1)) u=2u."(DnHE(K)-A)4.1.0

n=r,) f=-FLOAl(\J)

,F(ki•FW•i) 10 5t 1F(N-80u) 0,0+.50.1)u F1(0.(1.(:-F)*IP('.)+Feih(k*1)

GO TO b

Fl(e)--=(1.0-F)*IP(r,)+r,*1(N41)+F*(E-1....,)*(1h(N+2)-1(N+1)-1F<(N)+IR 1(-1))/4.0

6

7

9999

CONT1NUL frF(1)=().0

Pm=0.0

00 7 K=3,33,2 Ph=Pm+F1(K-2)+4.,*Flv-1)+E1(K) KM=(K+1)/2

MF(KM)=PM*00LF(2/l300U.O*CO COh1lNut

IIII=0

IF(Km.t .W.3) 1111=1 TmACE(1111)*(9999.18) CONT1NUE

nü 8 1(=.1433, Km=(K+1)/2

DLF(KM)=DnLF(K) 011F(KM)=DnHF(K) DF(Km)=Fi(K)

G=DLF(Km)*S(Z1)/DHI-((m)

Die ausführliche Berechnungsmethode der sekundären Zerstreung in den Planetenatmospären durch die Rechnungsmachine

10 11

ALI-.(KM)=ATAN(GiSQi-0(1.0-("i**2)) Y1,ALF(011

hzALA:J(Y1) MCALaC(Y1)

(2=AlüHF(KM3 GZ=.AHSIN(0) T=AK;(bHF(KM)-) TT=ABS(DHF(KM)-h)

IF( T.LT.0.5)00f(j2000

IF(TT.LT.0.5) (20 i0 2001 K1=1

KKK=37 Go 10 2002 2000 K1=14

KKK=37 ()() TO 2002

2001 K1=1

KKK=19

2002 PI=P

DU 91 KP=1+12 AS=30*KP-30 X2=P*FLOAT(A)

SA1=S(X2) (A1=C(X2)

OU 42 jJ=K1000:

LSCJJ)=5*JJ-5 L2=P•FLOAT(ZS(JJ))

PS1h1=S(Z2) 1,CThi=C(I2)

999 93

- ._=-(kS1r-i*SP+K:Th*CH*r(A) Ri1=-PUTHieksA

fi\J=-(1,STM*C4-fTt,*Stli:f,cA)

RKS=-KCAL*hCA*hSTrii*CA1+;,,Sp*RSTN1*SA1+SAL*RCA*RCThl Elm=-KLAL*hSA*PSTml*CAl-RC*RS11-1*SA1*kSAL*PSA*i.eCTh1 zEj=RSkq*kSTF41*CA14NChL*PCTH1

(lAri=hKS*Ch+7E1*Sh*

KNAP=-RKS*SH+LET*Ui X4=GAM*PKAP Ye=F1A*RKAP i4=-1.°KAP**2 XD=GAL*Flk„

Y5=-1.0+Efr**2

i).1.0-PKAp,"2-(Rl*›4*0-"*Y4+RN*24)**2 Ub=1.ü-ETA**2-(PL*X5+km*Y7+RN*Y4)**2

POL=12+DO

If(72-v7) 10,10+11

Ls-m-swf(T(A**2-inHh(k'm)*S(L2))**2)+DHFam)*C(Z2) GO 10 1/4/1

(Km.k(v.17.Aivt u0 lo 949

LS= Sudh,T(B**2-(pHi (Km)*S(72))**2)+0HF(Km)*C(72) G( TO 93

LS=0.0 pl=r,

Du 12 L=1,33 DDL(L)=LS*FLOA1(L-1)/32.00

uDHS(L)=S0WI(UHMi**2+DDLS(L)**2-2.0*DilF(KM)*NAS(L)*C(Z2)) 0=20.0*(nbhS(1)-A)+1.0

Fr-,(9-f-lOAT (N)

137

13 14 12

15

3800 CvtLM)=C()*(C(Y1)*C(A1)*C(X2)*(Y?)-:701) 11)*C

1(Y2))+S(w)*(C(Y1)-KL(Yz)-0O2)*S(Y1)*S(Y2)) mT(Lm)'mm“LM)/CV(LM), Dr'-(KM,LhA)=mF(Km)+N(Lt,)+MT(Lm) GU TO 16

180U

2800 16

IF(N.FW.1) GG TO 13 IF(N-800) 14013,13 f2(L)*(1.u-F)*IH(N)+f.*IR(N+1) 60 TO 12

hd(L)=(1.0-r)*IP(N)+r*I+1)+F*(F-1.0)*(1R(N+2)-1R(N+1)-1R(N)+IR 1(N-1))/4.0

cONTINuE PM=0.0 MS(1)=0•0

DO 15 L=1,3312

PN=Pm,F2(L-2)+4.0*F2(L-1)+F2(L) LM=(L+1)/2

MS(IM)=h-M*DULS(2)7(j(00.0*CC) CANIIN01-

(10 16 L.11J3,2 IM=(L+1)/2

DLS(LM).DNLS(L) DHS(Lm) =Dnms(L) ns(LNI)=h2(L)

0=10.0*(011S(Lf-)-p)+1.0 F=Q-ELOAT(N)

MMI(Lm)=VV(N)*(1.0-F)*h*Vm(h+1) 0-1 u

G=1)L5(L)*S(72)/nHbtL^)

ALS(UM)=ATAN(G/SRI(1.0-G**2)) Y2=ALS(LM)

1.1=RCAI*QC4 L2=KCAL*RA

L3=-RSAL M1.-kSA M2=RCA N1=h§AL*PCA N2=hSAL*HSA

N3=PCAL

Xlw=DHS( (_r)2)*CA1 Y1(v=DHS‘ Lm),:,5(y2)*SA1 Z1u=DHS

yG=L2*X1w+m2*Ylu+N2*il'' Z(J=L3*X1(9+m3*Y1,>+,1.3*,

Xw1= Lh*Y(\)*SH41L

“'1=Yw 7w1=-Sh*Y4-(H*,

IF(XWi.LT.0.(,) Gu To 1800

'ti .M)=C()*(C(Y1)*C(A1)*C(X2)*(Y?)- .':(X114,5(),z)*;,(v2)+c(x1)*5(y

u1=SwkT(Ywl**2+Zul**2)

IF((Q1.1..T.A) 0 Tu 280ü GO TO 3800

Om(KM,Lm).0.0

DS(LM)=°.0

CON1INUL

00 17 W=1,17

Die ausführliche Berechnungsmebhode der sekundären Zerstreung in den

Planetenatmospären durch die Rechnungsmachine ¹³⁹

18 17

20 19

21 92

2 5 240 91

169

17U 801

2 8

2 7 8 171

172 9

602

30

1)0 18 L=1+12

h(N.L)=Db(N)*HP(L)**U'i(rr+M) COkTINUE

CONIINUI DO 19 L=1.12 pm=0,0

DO 20 N=1.15.2

PM=PM+R(N•L)+4.G*RiN+1,L)+W(N+2,1) CONTINUF

SK(L),,OLS(2)*PM/3000.0 CONTINUh

DO 21 1=1,12 SKZ(JJ.L)=SK(L)*b(L2) PSKL(JJ.L)=SK/(JJ.L)*POL

CONTINte CONT1Nuf DG 240 L=1.12

f,f4=0,0

KKT=KKK...'2 DG 25 JJ=KT,KKI%2

h'r1=PM+PSKL(JJ,L)+4.0*PSK7(JJ+let.)+PSKL(JJ+2,L) CON'TINUF

fSKLJ(KPp.KmeKP,L)=Pm*5.0*P/3.0 CONIINU

C°NTINUE ehITE(64600)

wRITE(6•10) (L.L=1.12)

FOF'MA1(1H0.4H H.4H HT .4H Ah•4H 2F .4H KM.41, KP•8X, 1)2k 6H PSKLJ.12),/)

DO b01 KP=1,12

,,'ITE(6.170) H(M) •HT (1).AF.ZF(J)004,1,4,.(PSKZJ((P1~.KM.KP.L),L=1.12 1)

FONMAT(1H 1614..12F9.2) CONTINGE

D() 27 L=1.12 PM=0.0

00 28 KP=1+12 -,,...-Pm+PSK/J(KPPosm,0,,L)

CONTINUE

sL(KPP.Km,L)=PM*30.u*P

sS(KPP.KmeL)=Ur(eM)*U(KPP,KM,L) CONIINUF

CONTINUE.

WRITF(6,171) (LIL=1'12)

FoRMAI(10,411 F•4H HT,4H AF•4Ei ZF,4H Km.3)c,12(7H SL•

112).1)

DG 29 Km=1,17

WhITE(6,172) ii(M)•HT(1)4AF,7F(J)'KM«SL(KPP+10-,L),L=1,12) FO4IAT(1H .514.3X,12E9.2)

CONTINUF wRITE(64602) F

OmMAT(1HÖ,//) Uu 3u Km=1.1,

^N,--ITELA.172) i(m;•UT(1),A.F.ZFJ),KN.(SS(KPP,Km,L),L'1,12) CukT1NuE

Uu 31 L=1,19 PH=0.0

t;0 3£ ^m=1,15,2(

Pm=Ell+SS:KPP,Km,L)+,,.n*SS(KFP,Km+7,L)+5S(KIP1•',Km+r,L)

32 CuNT1NuE

5SSCH,KI,KPH,J,L)=Ph,*PPtlas(L)*PLF(2)

ISSS(t.i0(1,KPP,J,L)=—I(JOO.C*ALUG1°(SSS(MeKI.KPP,J,L))

31 CONTINGE

WRI1F(tg1''1) (L,L=1,12)

181 FOHmA1(1H0,4H H,4H hT,4H 21-,4H KM,7x422(71 SSS,

112)./)

.hlit(,•164) h(M),hT(I),AF,ZF(J).(SSSmol,KPI,.J.L),L=1,12) 184 FOkmA1.1H •414,7X,12F9.3)

WwITE(7.1b5) (ISSS(m,KI,KPP,J,L),L=1112),H(m).hT(1)4AF.ZF(J) 185 FOkmA-I (1715,2),.,41L+, 2H t

40 CONTINuF

4 CONTINuE

99 CONTINuE

90 CUr.TINuf

sTuP

ENU

Die ausführliche Berechnungsmebhode der sekundären Zerstreung in den

Planetenatmospären durch die Rechnungsmachine ^14]

NEUE METHODE ZUR BERECHNUNG OrR SEK ZER. MIT BERUECKSICH

TIGUNG AUF DER TAr'4GENTIELLER RICHITUNG NACH DEM AUSWERTUNG

DER NOETIGEN DATA

REAL PPQU(12)

INTEGER ZF(7),H(3),HT(42),AF,A5

REAL LI,L2,L3,M1,M2,M3,N1,N2,N3,LF,LS,IR(810),PP(12),PK(12),VM(410 1),DZLS(33),DZHS(33),F2(33),MZS(17),ALZS(17),H0(42),DZS(17),MZT(17) 2,MS(17,40,17),DDLF(33),DM(17),DDHF(33),DHF(17),GZ(17),F1(33),01-(17 3),ZS(1/,40),Z55(40),Z2(17,40),R5TH1(17,40),RCTHI(17,40),T(17),TT(1 47),DDLS(33),DDHS(33),DLS(17,40),DHS(17,40,17),DS(17,40,17),MMI(17, 540,17),ALS(17,40,17),GG(33),FG1(33),ALF(17),POL(17,40),CV(17),MT(1 6/),R(17,12),SK(40,12),PSKZ(40,12),PSKZJ(12,12),SL(17,12),SS(17,12) 79S5S(1,1,5,3912) 9M1-(17)

S(X)=SIN(X) 7

C(X)=COS(X) PAI=3.1415927

PAP=( 8.0*PAI)**2*6000.0 PAPP=2„5*9.0/PAP

READ(5,200) (PP(L),L=1,12) 200 FORMAT(I2F6,5)

READ(5,100) (PK(L),L=1,12) 100 FoRmAn6E12.5)

DO 4567 L=1,12 PPOO(L)=PK(L)**2*PAPP 4567

300

668 669

111

CONTINUE

READ(5,300) (IR(N)0=1,810) FORMAT(10F6,0)

VM(1)=0,0

WRITE(6,668) (PP(L).L=1,12) FORMAT(1H ,12F10.5)

WRITE(69669) (PK(L),L=1,12)

FORMAT(1H ,6E12.5)

WRITE(696666) (IR(I),I=1,810) FORMAT(I0F9.0)

PA=0.0

DO 111 K=1,799,2 PA=PA.IR(K)+4.0*IR(K+1)+IR(K+2) KM=(K+1)/2+1

A=6370.0 B=6410.0 CC:1037.4

VM(KM)=PA*0.05/(3000•0*CC) CONTINUE

VM(402)=1,0 VM(403)=1,0 VM(404)=1,0 VM(405)=190 P=3.14159/180.0 PI=P

601

M=I H(M)=30*M

W=P*FLOAT(H(M)) CH=C(W)

SH:5(W) WRITE(69601) FORMAT(1H1)

KI=1

WRITE(6,600)

I=5*KI-4

HT(1)=1-1

H0(1)=A+FLOAT(1-1) KPP=1

AF =30*(KPP-1) X1=P*FLOAT(AF)

R5A=5(X1) RCA=C(X1) WRITE(6,600)

WRITE(6,600)

600 FORMAT(1H0)

L5=50RT(B**2-A**2) 00 6600 L=1,33 PZL5(L)=L5*FLOAT(L-1)/32,0 DZHS(L)=5ORT(A**2+PZLS(L)**2) 0=20.0*(DZH5(0-A)+1,0

N=0

F=C4FLOAT(N)

IF(N.E0.1) GO TO 6601 1F(N-800) 6602,6601,6601 6601 F2(L)=(1.0-F)*IR(N)+F*IR(N+1)

GO TO 6600

6602 F2(L):(1.0•.,F)*IR(N)+F*IR(N4.1)+F*(F.-1.0)*(1R(N+2)-1R(N+1)-1R(N)+IR 1(N.•.1))/4.0

6600 CONTINUE MZS(1)m0.0 PM=0.0

DO 6603 L=1,31,2

PM=PM+F2CL)+4,0*F2CL+1)+F2(L+2) LM=(.+1)/2 +1

teS(LM)=PM*DZL5(2)/(3000,0*(C) ALZ5(LM)=ATAN(DZL5(L)/A) DZS(LM)=F2CL)

0=10,0*(DV-15(LM)..A)+1.0 N=Q

F=0-FLOAT(N)

IF(N.EQ.1) GO TO 6605 1F(N..400) 6606, 6605,6605 6605 MZT(LM)=VM(N)*(100..F)+F*VM(N+1)

GO TO 6603

6606 MZT(0)=VM(N)*(1.0..F)+F*WA(N+1)+F*(F-.1.0)*(VM(N+2)...VM(N+1).-VM(N)

6603

6604

50

1 +VM(N...1))/41,0 CONTINUE DO 6604 K=1,17 MS(K,20,1)=0.0

DL5(K,20)=DZL5(2) DO 6604 L=1,17 M5(K,20,L)=MZS(L) 05(K,20,L)=DZS(L) MMT(K120,0=1.0,-MZT(L) CONTINUE

PO 6001 K=1,33 DPHF(K)=A+425*FLOAT(K-1) 0=20.0*(DDHF(K).A)+1,0 N=0

F:0-FLOAT(N)

IF(N.E0.1) GO TO 50 IF(N.-800) 60,50,50 FI(K)=.(1.0-F)*IR(N)+F*IR(N+1)

GO TO 6001

Die ausführliche Berechnungsmethode der sekundären Zerstreung in den

Planetenatmospären durch die Rechnungsmachine 143

60 Fl(K )=(1.0...F)*IR(N)+F*IR(N+1)+F*(-..1.0)*(IR(N+2)-IF(N4.1)-.1R(N)+IR 1(N•1))*0.25

6001 CONTINUE 00 6 K=1,33,2 KM=(K+1)/2

DHF(KM)=DDHF(K) G:A/PHF(KM) GZ(KM)=ARSIN(G) DF(KM)=F1(K )

Z5(KM,18)=(80,04,GZ(KM)/P)/2.0 i5(KM,19)=GZ(KM)/P

ZS(KM;20)=GZ(KM)/12

ZS(KM,21)=(90.0+GZ(KM)/P)/2.0

6 CONTINUE DO 901 JJ=1,40

1F(JJ.LE.17) GO TO 899 IF(J.J.GE.22) GO TO 898 GO TO 901

899 Z5S(JJ)=5*J.,45 GO TO 901 898 Z55(JJ)=5*(JJ^4) 901 CONTINUE

DO 897 KM:1,17 DO 896 JJ=1,17

Z5(KM,JJ)=Z55(JJ) 896 CONTINUE

DO 895 JJ=22,40 Z5(KM,JJ)=Z55(4J) 895 CONTINUE 897 CONTINUE

DO 894 KM=1,17

894

8000

8001 8002

8003 10

DO 894 J,=1,40 Z2(KM,JJ)=P*ZS(KM,JJ) R5TH1(KM,JJ)=5(Z2(KM,JJ)) RCTH1(KM,JJ)=C(Z2(KM,JJ)) CONTINUE

DO 800 K=1,3392 KM=(K+1)/2

IF(KM,E0.1) GO TO 8000 IF(KM.EU.17) GO TO 8001 KTml

KKK=40 GO TU 8002 KT=22

KKK=40 GO TU 8002 KT=1 KKK=22 P1:P

DO 902 JJ=KT,KKK

IF(JJ.E0.19) GO TO 8003 IF(JJ.E0.20) GO TO 902

IF(KM.EO.1.AND,JJ,EU.22) GO TO 97 IF(Z2(KM,JJ).GZ(KM)) 10,10,11

L5z5ORT(OHF(KM)**2-A**2) GO TO 94

KJ=-1

0=A

GO TO 93

11

93

999 97

94

95 96

IF(KM.E0,17.AND,JJ,E0.22) GO TO 999 KJ=1

D=B

L5=FLOAT(KJ)*SORT(D**2.(DHF(KM)*5(22(KM,JJ)))**2)+DHF(KM)*C(Z2(KM, 14J))

GO TO 94 LS=0.0 GO TO 94 L5=SORT(B**2-A**2) P1=1)

DO 12 L=1,33

DDL5(L)=L5*FLOAT(L1)/32.0 IF(JJ.E0.19) GO TO 95

OOH5(L)=5ORT(DHF(KM)**24DDL5(L)**2-2.0*DHF(KM)*DDL5(L)*C(Z2(KM IiJJ)))

GO TO 96

DDHS(L)=SURT(A**24(LS*FLOAT(33, .)/32.0)**2) 0=20 .0*(DDHS(L).A)+1,0

N=0

F;.0-,FLOAT(N)

IF(N.EQ.1) GO TO 13 IF(N-800) 14,13,13 13 F2(L)=(1.0.-F)*IR(N)+F*IR(N)

GO TO 12

14 1(N F2(1.),;(1.0....F)*IR(N)+F*IR(N41)4F*(F^1.0)*(1R(N.62)..IR(N41) -IR(N)41R .1))/4.0

12 CONTINUE

MS(KM,JJ,I)=0,0 PH5(KM4J11)=DDHS(1)

DS(KM,JJ,I)=F2(1) MMT(KM,JJ,1)=0,0 ALS(KM,JJ,1)=0.0 PM=0.0

DO 15 L=1,31,2

PM:PM.F2(L)+4,0*F2(441)4F2(1 .42)

LM=(1.+1)/2 41

NM=(L+1)/2 41

M5(KM,JJAM)=PM*DDLS(2)/(3000.0*CC) DL5(KM,JJ)=DDLS(2)

DH5(KM,JJ,LM)=DDH5(1.) D5(KM,JJ,LM)=F2(L) 0:10,0*(DHS(KM,JJ,LM)..A)4.1

,0 N =U

F;Q-FLOAT(N)

IF(N.E0.1) GO TO 7700 IF(M-400) 7701,7700,7700

7700 MMT(KM,JJ,01)=1,0.VM(N)*(1.0.F)..F*VM(N41) GO TO 15

7701 MMT(KM,JJ,LM)=1.0•-VM(N)*(1.0-F)-,F*VM(N+1) 4-*(F -1.0)*(VM(N+2)-VM(N 141)-VM(N)+VM(N.•1))/4.0

G=PDL5(L)*S(Z2(KM,JJ))/DH5(KM,JJ

,LM) ALS(KM

,JJ,LM)=ATAN(G/SORT(1.0-G**2)) 15 CONTINUE

902 CONTINUE 800 CONTINUE

DO 401 KM=2,17 DO 401 LM=1,17

IF(KM..1) 4001,4001,4002

Die ausführliche Berechnungsmethode der sekundären Zerstreung in den

Planetenatmospären durch die Rechnungsmachine 145

4001

4002 401

ALS(KM,20,0)=0.0 ALS(KM,19,LM)=0,0 GO TO '401

AL5(KM,20,LM)=AL5(KM,19,LM)+ALZ5(LM) CONTINUE

00 40 J=5,7 ZF(J)=30*(J-1) ZImP*FLOAT(ZF(J)) TH=-904ZF(J)

THzP*TH RSTH=5(TH) RCTII=C(TH) MF(1):°•°

PM=0.0

700

2000

2001 2002

DO 700 K=1,33,2 GG(K)=A*S(L1)/DDHF(K)

FGICK )=F1(K )/SURT(1.0-GG(K)**2) CONTINUE

DO 7 Kz1,33,2 KM=(K+1)/2

Pm=PM+FG1(K)+4,0*FG1(K+1)+FG1(K+2) mF(KM)=PM*1.25./(3000,0*co

ALF (KM) =3. 141592+4 -ARbiN(Gu(K) Y1=ALF(KM)

R5AL=S(Y1) RCAL=C(Y1) CA=RCA*RCAL CB=RCA*R5AL ETA11=-RCAL*RSA ETA32=R5AL*R5A RKS11=-RCAL*RCA RK532=R5AL*RCA LI._RKS11 L2=-ETA11 L3=-R5AL m1=-RSA M2=RCA M3=0,0 N1=RKS32

N2=ETA32 N3=RCAL

RL=.(R5TH*5H+RCTH*CH*RCA) RM=-RCTN*R5A

RN=-(R5TH*CH-RCTH*511*KCA) IF(Km„EQ.17) GO TU. 2001

IF(KM,E0,1 ) GO TU 2000 KT=1

KKK=40 GO TO 2002 0=22

KKK=40 GO TO 2002 KT=1

KKK=22 PI=P

00 91 KP=1,1-

A5=30*KP•30

X2=P*FLOAT(A5)

SA1=S(X2) CA1=C(X2)

XXQI=DHS(KM,JJ,LM)*CA1 YNI=DHS(KM,JJ,LM)*5A1 RK521=RKS11*CA1 RK522=RSA*SA1

3800

1800

2800 16

18

R031=RKS21+RKS22 ETA21=ETAII*CA1 ETA22=-RCA*sA1 ETA31=ETA21+ETA22 ZET31=RSAL*CA1 CE=RSA*SA1 DO 92 JJ=KT,KKK

RKS:RKS31*RSTH1(KM,JJ)+RKS32*RCTH1(KM,JJ) ETA=ETA31*RSTH1(KM,JJ)+ETA32*ROH1(KM,JJ) ZET=ZET31*RSTHI(KM,JJ)+RCAL*RCTHI(KM,JJ) GAM=RKS*CH+ZET*SH

RKAP=-RKS*SH+ZET*CH X4=GAM*RKAP

Y4=ETA*RKAP Z4=-1+0+RKAP**2 X5=GAM*ETA Y5=-1+0+ETA**2

-D=1 ,0-RKAP**2-(RL*X4+RM*Y4+RN*L4)**2 DD=1.0-ETA**2-(RL*X5+RM*Y5+RN*Y4)**2 POL(KM,JJ)=D+DD

DO 16 LM=1,17

XX01=DH5(KM,JJ,LM)*CAI r(01=DHS(KM,JJ,LM)*SAL Y2=ALS(KM,JJ,LM) X10=XXUI*S(Y2) Y10=YY01115((2)

Z1O=DH5(KMIJJ,LM)*C(Y2) XO=LI*X10+MI*Y10+NI*ZIO

YO=L2*X10+M2*Y10+N2*110 Z0=1,3*X1O+M3*Y10+N3*Z1U

XO1= CH*X0+SH*ZO YOI=Y0

ZU1=-SH*XU+CH*Z0 IF(X01,LT,0„0) GO TU 1800 CD=CA1*S(Y2)

CV(LM)=C(W)*(CA*CD-CE*S(Y2)+Ct3*C(Y2))+S(W)*(RCAL*C(Y2)-CD*RSAL MT(LM)=MMT(KM,JJ,LM)/CV(LM)

DM(LM)=MF(KM)+MS(KM,JJ,LM)+MT(LM) IF(PM(LM).LT+0.0) DM(LM)=0.0 GO TU 16

01=SORT(YO1**2+Zü1**2)

IF(01,LT,A) GO TU 2800 GO TO 3800

DM(LM)=0.0 DS(KM,JJ,LM)=0.0 CONTINUE

00 17 Nm1,17 DO 18 L=1,12 .

R(N,L)zDS(KM,jj,N)*PP(L)**DM(N) CONTINUE

17 CONTINUE

Die ausführliche Berechnungsmethode der sekundären Zerstreung in den

Planetenatmospären durch die Rechnungsmachine 147

DO 19 L=1,12 PM=0,0

DO 20 N=1,15,2

PM=PM+R(N,L)+4,0*R(N+1,L)+R(N+2,L) 20 CONTINUE

5K(JJ,L)=DL5(KM,JJ)*PH/3000.0 19 CONTINUE

PO 21 L=1,12

P5KZ(JJ,L)=5K(!J,L)*5(22(KM,JJ))*POL(KM,JJ) 21 CONTINUE

92 CONTINUE A/TE(6,600) DO 178 JJ=KT,KKK

WRITE(6,177, AF,LF(J),KM,A5,Z5(KM,JJ),(

177 FORMAT(1H ,4I4,F7.2,1X,12E9.2)

178 CONTINUE

WRITE(6,600) WRITL(6,173) 173 FORMAT(1H0,4HPSKLe/)

WRITE(6,175) (L,L=1,12)

75 FORMAT(1H0,4H AF,4H j,4H KM,4H AS,5H

SK(JJ,L),L=1,12)

JJ,3X,12(3X,4H R,

1I2))

PO 174 JJ=KT,KKK

WRITE(6,177) AF,ZF(J),KM,AS,Z5(KM,JJ),(P5KZ(JJ,L),L=1,12) 174 CONTINUE

DO 240 L=1,12

IF(KM-0 242,242,243 242 PM1=0.0

00 241 JJ=22,38,2

PM1=PMI+PSKZ(JJ,L)+4.0*P5KZ(JJ+1,L)+P5KZ(JJ+2,L) 241 CONTINUE

PM =PM1*5.0*P /3,0 GO TO 246

243 IF(KM..17) 244,245,245 244 PO 8 LK=1,12

P5KZ(20,LK)=P5KZ(20,UK)+P5KZ(19,LK) 8 CONTINUE

PM1=0.0

DO 2410 JJ=22,38,2

PM1=PM1+PSKZ(JJ,L)+4.0*P5KZ(JJ+1,L)+PSKZ(JJ+2,L) 24I0 CONTINUE

PM1:PM1*5,0*P PM2=0,0 09 25 JJ=1,15,2

PM2=PM2+P5KZ(JJ,L)+4,0*P5KL(JJ+1,L)+P5KZ(JJ+2,L)

25 CONTINUE

PM2=PM2*5-.0*P

PM3m(P5K2(17,0+4.0*P5K2(18,04P5KZ(19,L))*(GL(KM)..80,0*P) PM4=(P5KZ(20,L)+4,0*P5K2(21,L)+PSKZ(22,0)*(90.0*P..GZ(KM))

245

250

888

PM=(PMI+PM2+PM3+PM4)/3.0 co TO 246

PM2=04 DU 250 JJ=1,15,2

PM2=PM2+P5KZ(JJ,L)+4,0*P5KZ(JJ+1,01.P5KZ(JJ+2,L) CONTINUE

DO 888 LE=1,12

P5K2(20,LE)=P5K2(20,LE)+P5KZ(19,LE)

CONTINUE

246 240 91

169

170 26

28

2 7 7 171

172 29

602 181

30

PM22PM2*5.0*P

PealPSKZ(17,L)+400*P5KZ(18,L)+PSKZ(19,L))*(GZ(KM)..8040*P) PM4=(P5Ki(20,0+4,0*PSKZ(21,0+PSXZ(22,0)*(904,0*PGZ(KM)) PM2( PM2+PM3+PM4)/3.0

PSKZJ(KP,L)=PM ) CONTINUE

CONTINUE WfITE(6,e>")

WRITE(6,169) (L,L=1,12)

FORMAT(1H0,4H H,4H HT,4H AF,4H ZF,4H KM,4H KP,8X, 112( 6H PSKZJ,12),/)

DO 26 KP21,712

WRITE(6,170) H(M),HT(I),AF,ZF(J),KM,KPOPSKZJ(KP,O,L21,12) FORMAT(1H ,614,8X,12E8.2)

CONTINUE WRITE(6,600)

DO 27 L=1,12 PM=0.0 DO 28 KP=1,12

PM=PM+PSKU(KP,L) CONTINUE 5L(Km,L)=PM*30.0*P

55(KM,L)=5L(KM,L)*F01(4)

CONTINUE CONTINUE

WRITE(6,171) (L,L=1,12)

FORMAT(1H0,4H H,4H HT,4H AF,4H ZF,4H KM,3X,12(7H 1I2),/)

I12),/)

DO 29 KM:1,17

WRITE(6,172) H(M),HT(I),AF,ZF(j),KM, (sL(Km,L),L.1,12) FORMAT(1H ,5I4,3X,12E9.2)

CONTINUE WRITE(6,602)

hORMAT(IHO,//)

WRITE(6,181) (1.,1-21,12)

FORMAT(1H0,4H H,4H HT,4H AF,4H ZF,4H KM,3X,12(7H DO 30 KM=1,17

WRITE(6,172) H(M),HT(I),AF,ZF(J),KM, (S5(KM,L),L=1,12) CONTINUE

DO 31 L=1,12

SL,

55,

32 31 183 184 40

PM.40.0

DO 32 KM=1,15,2

PM=PM+55(KM,L)+4,0*55(KM+19L)+55(KM+2,L) CONTINUE

S5S(M,K1,KPP,J,L)=PMAPPUU(L) -CONT1NUE

WRITE(6,183) (L,L=1,12)

FORMAT(1H0,4H H,4H HT,4H AF,4H ZF,7X,12(7H 555,12 ),/) WRITE(6,184) H(M),HT(1),AF92F(J),(555(M,K/OCPP,J,L),L=1,12)

FORMAT(1H ,414,7X,12E9.3) CONT1NUE

STOP

END