If you have a literal real in a program and you leave out the 'd0', that number would be chopped/rounded to the precision of a 32-bit real, i.e., about seven significant decimal digits. You specified that you wanted to process numbers with 13+2 digits; for that, you need double precision.
How to write the code is dependent on where your numbers come from. A calculation, or from reading a data file? If the former, your variables being printed out would have to be double precision. If the latter, you have to choose an appropriate input format. Please read the formatted I/O chapter of a Fortran book to learn the details.