FORTRAN SOURCE PROGRAM 



C STRESS IN CABLE LOWERING LOADS TO DEEP OCEAN. 

 C CAY=CONSTANT,K, UM=PARAMETER ♦ SFR=STEP IN FREQUENCY. 

 C FMAX=MAXIMUM FREQUENCY, STU=STEP IN U-ZERO, 

 C UMAX=MAXIMUM U-ZERO, FZER=INITIAL FREQUENCY. 

 1 READ 100,CAY,UM,SFR,FMAX,STU,UMAX»FZER 



U0=1.0 

 8 BETA=CAY*UO 



PUNCH 101 



PUNCH 102, CAY, UM, BETA, UO 



FREO=FZER 



PUNCH 111 , CAY, UM,SFR,FMAX,STU, UMAX jFZER 



PUNCH 109 



5 PHI=ATAN(FREQ/UM) 

 ALPH=2,*PHI 

 DELT=FREQ+PHI 



UONS=(COS(DELT) )**2/( 2.*BETA»»2*S I N ( PHI ) ♦♦2*S IN < FREQ ) ♦*2 ) 

 TERM={ l.+( (BETA**2*SIN(FREQ)**2*SIN(ALPH)**2)/C0S(DELT)»»4) ) 

 TERM=TERM»»0.5 

 TERM=TERM-1. 

 U0N=U0NS*TERM 

 U=UON«»0.5 



TPSI = ( (0.5*BETA**2.*UON)*SIN(PHI )/ COS (PHI ) ) - ( COS ( ALPH ) /S IN ( ALPH ) ) 

 PSI=ATAN(TPSI ) 



STRS=(1.+SIN{PHI ) /COS (PHI )*(SIN(PSI ) /COS (PS I ) + l./COS(PSI ) ) ) 

 STRS=STRS»U0N»(FREQ**2. ) 

 STR=STRS»*0.5 

 PUNCH 110, STR, FREQ 

 IF{FREQ-FMAX)3,4,4 



3 FREQ=FREQ+SFR 

 GO TO 5 



4 IF(U0-UMAX)6,1,1 



6 UO=UO+STU 

 GO TO 8 



100 FORMAT(2F10.2,2F10.8,2F5,2,F10.8) 



101 FORMAT(8X,9H CONSTANT ♦4X , 3H UM,10X,5H BETA,8X»7H U-ZERO) 

 10 2 F0RMAT(3X,F10.2,3X,F10.2,3X,F10.2,3X,F10.2) 



109 F0RMAT(7H STRESS » lOX , lOH FREQUENCY) 



110 F0RMAT(E15.8,3X,F8.4) 



111 F0RMAT{2F10.2,2F10.5,2F5.2,F10.5) 

 END 



68 



