Cer eee22-222 2 TEST IF INTERVAL LENGTH CAN BE DOUBLED 
KUTME1O] 
100 IF (ERRUR#64.-EPSE(I)) 11051205101 KUTME102 
lol INC = 1 KUTME103 
110 CONTINUE KUTME 104 
C-- 2-2 - == - = = UPDATE T AND SOLUTION KUTME105 
ll T = TeHC KUTME106 
DO 112 I=1,ND KUTME107 
le YO(I) = v2(1) KUTME 108 
Ce - ee ee = = - = ~ GET SOLUTION IN NEXT INTERVAL KUTME109 
LOC = LUC+1 KUTME110 
IF (LOC-IPLOC) 12092105210 KUTME LLL 
120 IF (INC)210,130,210 KUTME112 
130 IF (LOC-(LOC/2) #2) 21051405210 KUTME113 
146 IF (IPLOC-1) 210,210,200 KUTME114 
C--+22e-+2+- - = = DOUBLE INTERVAL LENGTH KUTME115 
200 HC = 2,%HC KUTME116 
LOC = LUC 72 KUTME117 
IPLOC = IPLOC/2 KUTME118 
210 IF(IPLOC-LOC) 30,329,30 KUTME119 
329 BwACL = F0(2)-xACCL#F0(3) KUTME120 
CGACL = FO0(2) KUTME121 
RETURN KUTME 122 
END KUTME123 
END KUTME124 
SUBROUTINE DAUX (TIME »X»RHS) DAUX 2 
Cc DAUX 3 
C TIME TIME AT WHICH SYSTEM IS TO BE EVALUATED DAUX 4 
© x STATE VECTOR DAUX 5 
Cc RHS THE RIGHT HAND SIDE UF THE EQUATION S =F A DAUX 6 
Cc DAUX 7 
REAL KAR DAUX 8 
REAL IAs IT oMoK sMAgMASSoNCGoNL oN gMMAX DAUX 9 
INTEGER END »PTIME DAUX 10 
DIMENSIUN x (6) »RHS(6) oF (391) 9A(393) 9 INDEX (353) 9 DAUX 11 
Aj R(120) 0V(120) 9D(120) DAUX 12 
Cc DAUX 13 
COMMON /SHIP/ MASS oCINT »QA9CE 9CE2 9CE3 yDMUsEDMU,E2DMUsE30MU 9 BF 9BMMyDAUX 14 
NLoFLsIAgE (120) DAUX 15 
“COMMON /CUNST/ NCGsECG oP sDPR»RPD sGRAVTY »RHOsK »NUM»MA(120) »CDy9 TAs OAUX 16 
e B(120) sBETA HW (120) 9 TZ 9DRAGoWoXDoT oXPoMolTs DAUX 17 
0 DELTAS s TX9EST (120) »CoRO»KAR gMMAX (1-0) oe TEST (120) o DAUX 18 
N (120) »>PHALF DAUX 19 
COMMON /IN/ BM (120) ,81 (120) »VELIN DAUX 20 
COMMON/UUT /NPRINT »NPLOT sEND VAUX 21 
COMMON /SEAWAVE/ START »RISE ,RAMP DAUX 22 
COMMON /WAVE/ RoPT(120) oZMAgZWMA EMAS 5 ZZWMA 9 ZWEMA 9 ZOWMA 9ECMAZ 9 DAUX 23 
6 ZwOOT (120) DAUX 24 
Cc DAUX 25 
RAMP = RMP(TIME START RISE) DAUX 26 
PIH = PI/2, DAUX 27 
CT = C#TIME DAUX 28 
Cx6 = CUS(x(6)) DAUX 29 
SK6 = SIN(X(6)) DAUX 30 
CaeeeeeeSET VALUFS UF MA AND B DAUX 31 
00 75 I=1,NUM DAUX 32 
PT(L) = (X(4) E (1) #CX6eN(1) #SA66CT) HK DAUX 33 
R(I) = RU®COS(PT(I)) #RAMP DAUX 34 
Cc # * # # # # # # COMPUTE HW SUBMERGENCE OF A POINT AND R- THE WAVE DAUX 35 
Cc HW(T) IS IN THE FIXED COORUINATE SYSTEM DAUX 36 
Si 
