Cc 
65 
Cc 
Cc 
Cc 
Cc 
Cc 
70 
75 
74 
76 
77 
78 
79 
80 
si 
82 
85 
Cc 
Cc # # 
Cc 
15 
17 
Cc # * 
18 
C # # 
HW(I) = X(5)-E (1) #SX6eN(1) #CX6-R (1) 
IF (HW(1) GT.) GO TO 65 
CRAFT IS NOT SUBMERGED 
MA(I) = 0. 
Bl(1I)=0. 
B(I) = 0. 
GO TO 75 
V(I) -RO#K*SIN(PT (1) ) *RAMP 
0 (1) Hw(I)7(CA6—-V (1) #SX6) 
D(I) 1S IN THE BODY AAIS SYSTEM AND IS THE SUBMERGENCE 
IF (D(1) GE.TEST(I)) GO TU 70 
CRAFT IS PARTLY SUBMERGED 
B(I) = O(1)#(1./TA)#PIH 
BL(T> D(T)#(1-/TA) #PI 
MA(I) KAR#PHALF #B (I) #B(1) 
GO TO 75 
CHINE IS IMMERSED 
Bl ARRAY IS USED FOR THE INTEGRALS OVER THE PORTION 
OF THE HULL FOR WHICH THE CHINE IS NOT IMMERSED 
MA(I)=MMAX(1) 
B(I)=BM(1) 
B1l(I)=06 
CONT INU 
IF (NPRINT.LTe%) GO TO 85 
WRITE (6974) TIME 
FORMAT(" TIME = "yF10.4) 
WRITE (6976) (X(T) » T=1 96) 
WRITE (6977) (R(I) » I=] NUM) 
WRITE (6978) (HW(I) »I=1 NUM) 
WRITE (6079) ( B(I) »I=l NUM) 
WRITE (658%) (V(1I) o I= »9NUM) 
WRITE (6581) (O(1) » l=] NUM) 
WRITE (6982) (MA(I) » T=) »NUM) 
FORMAT(" X(I) "“96(2X0E120c6)) 
FORMAT (* R(T) %910F 10.4) 
FORMAT ( HW(I)%510F10.4) 
FORMAT (" 8(1)%510F10.4) 
FORMAT (* v(I)%510F10.4) 
FORMAT (" (I) %910F 10.4) 
FORMAT (" MA(I) “slOF1064) 
CONTINUE 
* @ # # # » COMPUTES NL AND FL ANU THE ASSOCIATED INTERGALS 
CALL FUNCT (X) 
IF (NPRINT.LT2%)GO TO 17 
WRITE (6915) TXoFLoDRAGoTZoWgNboXDoT oXP 
FORMAT(" ",10E12.6) 
CONTINUE 
# # # # # COMPUTE THE F VECTOR 
F (lol) = \TX#FL#SX6=DRAG#CX6 
F(151)=0.0 
F(20l) = T2*FL#CX6*DRAG*SX6eW 
F (391) =NL=DRAG#XD+T#xXP 
IF (NPRINT.LT.3)G0 TO 18 
WRITE (6910) (F (Tol) »I=193) 
CONTINUE 
# &# # # # COMPUTE THE A MATRIX 
A(lol) = M+MASS#SX6#SX6 
58 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUXK 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
VAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
OAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
DAUX 
