aa0 
IF (EST (NUM) .LT23075) STOP 3 
COMPUTE NO AND BM ARRAYS 
DO 32 I=1,NUM 
IF(EST(I) .GE.0.75) GO TO 30 
NO (T) =-0046875# (100-SQRT (EST (1) /00375-(EST (1) 0075) ##2,0) ) 
BM (1)=.375#SQRT (1.0-(EST(1) 7. 75-10) ##200) 
GO TO 32 
30 NO(I)=0.0 
BM(I) = 0.375 
32 CONTINUE 
Ceeeteee2COMPUTE CONSTANTS AND INITIALIZE ARRAYS 
ANANNQANAANANANANANAN 
oO 
Cc 
IN 
M=W/GRAVTY 
RHO=1.99 
IT=M#RG*RG 
K = 2.%PT/LAMBDA 
C=SQRT (GRAVTY/K) 
NU=RO#K 
PHALF = (PI/2.)#RHO 
BETA = BETA#RPD 
CO = COS(BFTA) 
TA = TAN(BETA) 
DO 60 I=1.NUM 
E(I) = ECG-EST(I) 
N(I) = NCG*NO(I) 
MMAX(1) = KAR®PHALF®BM(I) #BM(1) 
TEST(I) = (20%BM(1)#TA)/PI1 
CONTINUE 
END=END*1 
RETURN 
END 
SUBROUTINE PLUTER(FX»XA»HMAX »LAMBUA, IB yNWAVE 9 IPT) 
PUT: 
FX A TwO DIMENSIONAL ARRAY CONTAINING PITCH AND 
HEAVE VALUES AT EACH TIME STEP 
KA INITIAL TIME 
HMA x TIME INTERVAL» PTIME*HMAX = INTERVAL BETWEEN 
FX VALUES 
LAMBDA WAVELENGTH USED IN CALCULATING PITCH AND 
HEAVE RATIOES 
IB NUMBER OF FX VALUES 
NWAVE START OF VALUES AFTER WAVE IS COMPLETELY ON 
REAL ITsK »LAMBDA 9MoMA oMMAX gNoNCG 
INTEGER END 
DIMENSION FX (29400) »FMIN(2) »FMAX (2) sNVAR (2) 
COMMON /CONST/ NCGsECGoPI sDPR»RPD»GRAVTY »RHO 9K »NUMo9MA (120) sCDo TAs 
e B(120) sBETAsHW(120) »TZsDRAGoWoXDoToXPoMoITs 
DELTAS 9 TX oEST (120) sCoROoKAgMMAX (120) sTEST (120) 5 
N(120) »>PHALF 
@ 
COMMON/OUT /NPRINT sNPLOT »END 
C# # # # ¢ # # # SET UP VALUES FOR PLOT AND CREATE PLOT 
66 
INPUT109 
INPUT110 
INPUT1L11 
INPUT112 
INPUT113 
INPUT114 
INPUT115 
INPUT116 
INPUT117 
INPUT118 
INPUT119 
INPUT120 
INPUT121 
INPUT122 
INPUT123 
INPUT124 
INPUT125 
INPUT126 
INPUT127 
INPUT128 
INPUT129 
INPUT130 
INPUT131 
INPUT132 
INPUT133 
INPUT134 
INPUT135 
INPUT136 
INPUT137 
INPUT138 
INPUT139 
INPUT140 
INPUT141 
PLOTER 
PLOTER 
PLOTER 
PLOTER 
PLOTER 
PLOTER 
PLOTER 
PLOTER 
PLOTERIO 
PLOTER1L1 
PLOTER12 
PLOTER13 
PLOTERI14 
PLOTERIS 
PLOTERI6 
PLOTERI7 
PLOTERI8 
PLOTERIO 
PLOTER20 
PLOTER21 
PLOTER22 
PLOTER23 
PLOTER24 
PLOTER2S 
PLOTER26 
PLOTER27 
ODINONL WHR 
