Table n-7. Listing of the computer program MADSEN. --Continued 



20 »HTS PSOGWAH IS(,F6,2»( SEC FOR THIS CONDlTlONf) 



LSUs0.7'?9 

 lOS I=(l SL*10.+ 1 ,) 

 C INTeKPOi.ATE INPUT TABLE FOR THIS LSL VALUE 

 11*1 SL*?0,+ 1 • 



25 00 5 j=i.n 



FSSfJ)3FST(I,J)+(PST(l+liJ).FST(I,J))*(LSL-{I-n»0,|)/0.l 



ftUSf J3=HUT(IiJ)*CRUT(I*l»J).RlJT(l,J))*(LSL''(I»'n*0.n/0.l 

 1 «S( j)=RTcil,JufHTciI*l.J)-RT(II,J))*(USt-(II»n»0,05)/'0,05 



C GllfeSS PhI AMD ITERATE 

 30 P^-ls5,0 



6 J«PhI 



MCs(ALnGfPHl»i.)»ALUC(J*.l,))/(4L0Gf J*i»,)'ALOGtJ*l,)) 



Fs=pss{j*i)* FAC»(^ss(J♦^)•fss(J♦^) 



35 Hi)=ouSf J*i )♦ FAC*(KUS(J*25»RUS(J*1)) 



KUsWS(J+1 )*(RSf J*25"«S(J*n)*FAC 



ARGsO,?''*f 0/Hn)«*p,?»{RU«?.*A/(HO*TANB)1**0,5*FS 

 PhlMin,S*ATAMfA»G)*57,29b78 



nO DtLsAHSrPHIN-PHT) 



IF tM.GT.en) GO TO 9 



Pwi-Phini 



If (PMl.LT.n.oi) PHUO.ol 



IF IP-HI.GT.9.V<») PH1«9,9<) 

 4S I^lntl .GT.o.OS) GO TO 6 



9 KIl3rtII*CF 



WETllH^ 



1 SlIHRQUTINif REAOT 



C^YMU^l/MA^)s/^STf'»•^)lF<UTt9,ll)tRTC17•ll)lTx(9t^O),RXt9«lO) 



177 FOWMATC^Xi 17FU.?) 



PU 1 Ms 1,11 



S 1 KfAnCS,)77) CFST(N,M) ,H=1»9) 



DO S Msl,l 1 



2 HfAntS.177) (RUT(N,M),N»1,9) 



DO 5 MslHl 

 ■i WtAn(S,177) (RTfN,M),Nal ,17) 



10 DO u Ms 1,10 



« REAn(S,177) (TxrN'M) »Nsl,9) 



00 S Ms 1,10 

 S ReAn(5,177) (Rx(M»M)|Nil,9J 



RETiiHN 

 15 EMO 



1 SlIbRUllTINE EORwf0H7.,0HT,LE,H0,HS,TANB.N9,0H,TOPW) 



coMMON/MAPisi/MM,NL,o(in,Ncin» Ldjf U)»TM(in 



OIMpNSION Hf TA(1 I) ♦0M( 11) 

 f^t AU M.L«1 E .f^R 

 5 '"«=«. aj5 



BETaHs2,7»(1 ,-NR)/(MR»»3*0R) 

 on ■>! tsliNM 

 ?1 BMA(n=2.7*(l,-NCn)/(N(I)*»3*0(I)) 

 TMlrO. 

 10 11^230, 



on a J=i .''^L 

 Twi=THi*THtJ) 



OH{J)sTHtJ)/MO 

 15 IF(TM1 .GT.HO) OH(J)8(HO"TH2)/HO 



IP(tHi ,GT,hO) GO to S 

 U TH^5Th2*TW( J) 



S SUm2=0. 



00 16 J=1,NVL 

 20 5UM,rt), 



t'O i7 l3l,NM 

 17 SUM, =sum1*BETA(T)/8ETAR*L(I,J) 

 16 SUM?=suM2+dm(J)/(SQRT(SUM1)) 

 LF=i./(SUM2*«2)«DHE/DHT 



186 



