FX A two-dimensional array, containing time histories to be plotted 
XA Initial time 
HMAX Time-interval increment; time interval between values in FX is given by 
HMA X*PTIME 
LAMBDA Wavelength 
IB Number of values to be plotted 
NWAVE Position in FX at which wave is completely turned on 
Function RMP (T, START, RISE) 
The RMP is a function that calculates a value between O and | corresponding to time T, 
based on a straight line from time START with a value of 0 to time START plus RISE with 
a value of 1. It is used to lower the initial wave amplitude to avoid large transients at start 
of the computations. 
The arguments are: 
Wy Actual time 
START Time at which to begin the ramp from 0 to 1 
RISE Duration of rise from O to | 
The function reaches the value 1 at time START plus RISE, if the rise is 0.0, RMP will 
return a value of 0.5. 
Subroutine TRAP (F, DX, NPTS, ANS) 
This routine performs the evaluation of an integral using a trapezoidal approximation. 
The argument variables are defined as follows: 
F Array of integrand values 
DX Increments at which F is evaluated 
NPTS Number of values in F 
ANS Result, which is equal to 
NPTS 
DX bs F(i) - 0.5 [F(1) + rovers) 
i=l 
PROGRAM PLTHSP 
This program uses a data file created by program MAIN to create CALCOMP plots. The 
data are read from logical Unit 9 and are rewritten on Unit 7 for CALCOMP input. Program 
PLTHSP sets the tape output unit equal to 7 and calls SUBROUTINE CALPHI to execute the 
plot procedures. 
48 
