LMM = LM - 1 



PI = 3.K+159 



READ (b»20) R»P0»PN»VF»S£.C1»C2.C3»WKC0R 



READ (5»2U) WXrRFF 

 ?0 FORMAT (10F7.2) 



READ (5»15) (SA(i\J) .N=l»r4M) 

 C IF IDATA = 1 .WIiMO STRESS» ASTRONOMICAL SETUP, PRESSURE SETUP AND WiNn 

 C STRESS COEFFICIENT DATA SUPPLIED BY AN EXTERNAL SOURCE* OTHERWISE 

 C DATA GENERATED INTERNALLY FROM WIND* RADII » AND THETA CURVES. 



IF (IDATA .EQ. 1) GO To 220 



READ (5flO) IR»IW»IT 



READ (5»30) (XR(I)f I=1»IR) 



READ (5»30) (XW(I) rl=lf IW) 



READ (5»30) (XT(I) »I=lr IT) 



READ (5»30) (VRd) »I = 1,IR) 



READ (5»30) (VW(I) f I='1»IW) 



READ (b»30) (VT(1) »I=lrIT) 

 30 FORMAT (10F8.2) 



XRE = XR(IR) + 1.0 



IND = XRE 



DO 35 1=1, IT 

 35 VT(I) = (VT(I) * Pl)/160.0 



IK = 1 

 C IN THE FOLLOWING THkEt LOOPS ALL ORDINATE VALUES ON THE CURVES ARE 

 C INTERPOLATED FOR EACH UNIT OF DISTANCE. 



DO 60 1=1, IR 



RS(IK) = VRll) 



II = XR(I) 



IKN = XR(I+1)+1 



III = IKN - 1 



IF (II .LQ. Ill) GO TO 50 



IK = IK + 1 



RATIO = (VRd+l) - VR(1))/(XR(I + 1) - XR ( I ) ) 

 UO RS(IK) = RS(IK-I) + RATIO 



IK = IK + 1 



IF (IK .Gt. IKN) GO TO 60 



GO TO 1*0 

 50 IK = Ik + 1 

 60 CONTINUE 



IK = 1 



UO 90 1=1, IW 



W(IK) = VW(I) 



II = XW(I) 



IKN = XW(I+1) + 1 



III = IKN - 1 



IF (II .EQ. Ill) 60 TO 80 



IK = IK + 1 



RATIO = (VW(I + 1) - VWd) )/(XW(I + l) - XWd)) 

 70 W(IK) = tJ(IK-l) + RATIO 



IK = IK + 1 



IF (IK .GL. IKN) GO TO 9U 



GO TO 70 

 80 IK = IK + 1 

 90 CONTINUE 



IK = 1 



KII = U 



DO 120 1=1, IT 



IF (KII .ta. 1) CO TO 93 



T(IK) = VT(I) 



IF (VT(1+1) .ST. 0.0) GO TO 9b 



KII = 1 



47 



