Table G-7. Lis,ting of the computer program NIADSEN. --Continued 



END 

 1 bllHpH'ITINF INTfR(^,TiL.HU»A|NUiriTI,HIiWt»lFUAG) 



L0rMuN/'^4nS/F-STt<>ill)iKUT(9,intRT(l7f n)iT)((9,10)|RX(9ilO) 

 DIH^i'-siOM TS( 10) •"StlO) 

 5 HfAl, flKL.l l^U.K^«LAMtlDA,N 



SS=C.M/0.0S)*»2 



10 L A>^Bi)A=l , 



no. 



KC=, 7n, 



IC = n 



IS K = I1.M 



U=A»&iJRT(<J,80/Hn)/(l,*l.AMBUA) 



f<r;ai)»n/MJ 



F=N'CKO*L)»fSQHT(l .♦(l,»RC/Rl/)«iJ6.*BETA«**L/(5,«3,10l59»M0) ))•!,) 



LAMnLAcKO»L*f-/{?,*N) 

 20 ir( TC.GT.IO) GO TO 5 



IF C f AHStfN-fp /F) .Gr .0,02) GO TO 2 

 ■5 T 1=1 ,/{ 1 .♦LAhHOA) 



KI = l <iManA/( 1 .♦LAMBDA) 



KS=F/SS 

 25 C "KlTtfbiJ^T) F.FSiU.HD 



597 (DWMATt?i.ix, (F,FS»H.R0s(,uEi5,5) 



IftKJKL.GT.n,?) IFLAG31 



If (M'^L.GT.O.'?) RtTuHN 



If (H<l.,LT,[i.n RETURN 

 30 If'CFb.GT.JS.) FS=J5, 



J = ^ k, L ♦ 1 , 



1*FS 

 C iMFf'PUiATf MApSEM CURVES 2 i'^D i 



UC 1 -^stflO 

 35 hSC'O=RX(J,M)*(RX(j*l.M)-Rx(j.>i))*(MM-n,l*j)/0,l 



I T3(''i) = TX(J,H) + (TX(J*l,M),Tx(J.,''))«(NiKUi>n,l»J)/0.l 



lF(FS.LF,t,n) TIiTSfn*AUOGlO(FS)*(TS(ln)"TS(n) 



IMrS.LF.l.n) RT = W5Cl)*ALOG10CFS)*(RS(lO)-HS(l)) 



IMrS.GF.lO.) TIsTS(l(1)»(3S.»FS)/2sI 

 UO If (FS.GF ,1 0.) RTsRS( 10)*(l.»KS(105)»(f'S„l0t)/?5. 



lf(rb,Lf.l,0,nR^FS,Gt,10,0) RETURN 



KI=R5(I)+fRS(Ul)-RS(I))»(AL0G(FS)«AL0G(I*l.))/t*L0G(I+l,)»*L0C(l« 

 ♦1.)) 



TIsTS(n*f TS(l*l)»TS(I))*(ALOG(FS).ALOG(I»l,))/(*l,OG(I*l,)«*tOG(l« 

 «5 ♦! ,)1 



"F Tul'N 



t ^D 

 1 SUbwUUTTNE ltNr.T( DLOiDL) 



•<EAl LDil.nNjEw.LOO 



L03, ,0/nLn 



LOD-LO/DLO 

 S N=l 



Rl=5,iy tS9 



1 ARG-c".0*PT/LO 

 I.ONf« = lO0*TANHCaRG) 



10 DIFF=AyS(l.pN£w,uD) 



lF(fj.?0O) j.U.a 

 5 IF(nlFF-0,0O0S) 2»2t5 



S Ln=f inNFw*Ln)/2.o 



GO 7 1 



*RlTt(6.10P) riLOtfL 

 loo F0RM4T(uaH SUBROUTINE LENGTH DID N07 CONVERGE!) ri/UO ■ jFtCSs 



1 8nI'/L » tFlO.S) 



2 ULBl,n/LON[w 

 20 t<fTuRN 



EMD 



187 



