Values output using list-directed I/O (`PRINT *, R, D') should be written with a field width, precision, and so on appropriate for the type (precision) of each value.
(Currently, no distinction is made between single-precision
and double-precision values
It is likely this item will require the
to be undertaken.
In the meantime, use of formatted I/O is recommended.
While it might be of little consolation,
g77 does support `FORMAT(F<WIDTH>.4)', for example,
as long as `WIDTH' is defined as a named constant
That at least allows some compile-time specification
of the precision of a data type,
perhaps controlled by preprocessing directives.