The phase behavior of suspensions of colloidal hard tetragonal parallelepipeds ("TPs") (also known as rectangular nanorods or nanobars) was studied by using Monte Carlo simulations to gain a detailed understanding of the effect of flat-faceted particles on inducing regular local packing and long range structural order. A TP particle has orthogonal sides with lengths a, b, and c, such that a=b and its aspect ratio is r=c/a. The phase diagram for such perfect TPs was mapped out for particle aspect ratios ranging from 0.125 to 5.0. Equation of state curves, order parameters, particle distribution functions, and snapshots were used to analyze the resulting phases. Given the athermal nature of the systems studied, it is the interplay of purely entropic forces that drives phase transitions amongst the structures observed that include crystal, columnar, smectic, parquet, and isotropic phases. In the parquet phase that occurs for 0.54