& 
50 
si 
60 
COMMON /INZ BM(120) oB1 (120) »VELIN 
COMMON /U0UT /NPRINT »NPLOT »END 
FUNCT 
FUNCT 
COMMON /WAVE/ R(120) sPT (120) o> 2MAgZWMAgEMAS,ZZWMA, ZWEMA, Z2WMA,E2MAZF UNCT 
9ZWOOT (120) 
e 
COMMON /INTER/ II oKTT (10) oDIFF (10) 
COMMON /SEAWAVE/ START »RISE,RAMP 
COMMON /TEST/ VMA 
#2 2 #@ #@ » # @ INITIALIZE INTEGRAL SUMS 
0.0 
(1) ®SIN (XK (6)) oX (2) @CUS (X (6) ) 
Sx6 = SIN(x(6)) 
CX6 = COS(X(6)) 
WO = K#C 
# @ # @ # # SET UP THE FUNCTIUNS FOR THE INTEGRALS 
DO 90 I=1l»NUM 
IPART (1) =E (1) #E (1) #MA(T) 
QPART (I) =E (1) #MA(T) 
ZWOOT(I) = -RU#WO#SIN(PT (I) ) #RAMP 
U = X (1) #CX6—-X (2) *5X6e¢ZWDOT (1) #SX6 
VEL = VPART=X(3) #E (1) -ZwDOT (1) #CX6 
N 
N 
a 
= 
> 
1 0 tb 
Z1(1T) = MA(I)*#ZWOOT(T) 
Z2(1) = -Ma(1I)*#COS(PT(I)) #RAMP 
Z3(1) = E(1) #221) 
24(1) = E(T) #2101) 
Z5(1) = U#z2(1) 
Z6(1) = E(7)#25(1) 
Z7(1) = MA(I) #VEL#U 
IF (VELeLE.0.) GO TO 60 
IF (BI(1) LE. 0.0) GO TO 50 
DRDT = ZwOOT (1) #(X(L) Cox (3) ®(N( 1) @CX6=E (1) #SX6) ) /C 
DI(1) = — VEL®BL (1) #(X (2) x (3) # (CXO*E (1) #SX6#N(T)) 
GO TO 51 
DI(1) = 0. 
CONTINUE 
p2¢1) = €(1)*D1(1) 
Cl(1) = VEL#VEL#B(I) 
C2(I) = €(1)#C1 (1) 
GO TO 6 
OAS) en Oa 
D2(I) = 0. 
C1(1) = 06 
C2(1) = 0. 
60 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
(PAGE 4 OF NOFUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
=-DRDT) FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
FUNCT 
