simulator('spectre)

PERIOD = <your_clock_period_for_NRZ_data>
TSTART = <start_time_for_eye_diagram_analysis>
TSTART = <stop_time_for_eye_diagram_analysis>
vthreshold = <your desired threshold voltage>

simulator_results = strcat("<path_to_your_results_directory")

eye_waveform = eyeDiagram(v("vout" ?result "tran-tran" ?resultsDir simulator_results) TSTART TSTOP 2*PERIOD ?autoCenter t )

sprintf(fname_csv "eye_height_measurements.csv")
fp = outfile(fname_csv "w")
fprintf(fp "Time (sec),Eye height (V)\n")

delta_t = PERIOD/10.0
time = 0.50*PERIOD
max_eye_height = -1.0
max_eye_height_time = 0.0

while(time < (1.50*PERIOD + 0.10*delta_t)

eye_height = eyeHeightAtXY(eye_waveform time vthreshold ?output "total")
fprintf(fp "%1.6e,%1.6e\n" time eye_height)

if( (eye_height > max_eye_height) then
	max_eye_height = eye_height
	max_eye_height_time = time
	)
time = time + delta_t) ; end while loop

fprintf(fp "Maximum eye height of %1.6e V at %1.6e s\n" max_eye_height max_eye_height_time)
close(fp)

exit()