diff --git a/src/spicelib/devices/bsim4/b4reliability.c b/src/spicelib/devices/bsim4/b4reliability.c
index f873e4f03..0229d1e6b 100644
--- a/src/spicelib/devices/bsim4/b4reliability.c
+++ b/src/spicelib/devices/bsim4/b4reliability.c
@@ -372,7 +372,8 @@ BSIM4reliability_internal (BSIM4instance *here, CKTcircuit *ckt, unsigned int mo
                 }
             }
 
-            printf ("\n\nExtrapolation at %-.9g years:\n\t\t\t\tDeltaVth: %-.9gmV\n\n\n\n", ckt->CKTtargetFitting, here->relStruct->deltaVth * 1000) ;
+            fprintf (stdout, "\n\nExtrapolation at %-.9g seconds (%-.2g years):\n\t\t\t\tDeltaVth: %-.9gmV\n\n\n\n", 
+                ckt->CKTtargetFitting, ckt->CKTtargetFitting / 365 / 24 / 60 / 60, here->relStruct->deltaVth * 1000) ;
         }
 
         // Free the vgs average list
diff --git a/src/spicelib/devices/relmodel/relmodelcalcfitting.c b/src/spicelib/devices/relmodel/relmodelcalcfitting.c
index 2f51c1e45..28582b3b1 100644
--- a/src/spicelib/devices/relmodel/relmodelcalcfitting.c
+++ b/src/spicelib/devices/relmodel/relmodelcalcfitting.c
@@ -66,7 +66,8 @@ RELMODELcalculateFitting (unsigned int rows, unsigned int number_of_periods, dou
         f += gsl_vector_get (x, 2 * j) * sin (j * target * factor_for_2pi) ;
     }
 
-    fprintf (stderr, "\n\nExtrapolation at %-.9g seconds:\n\t\t\t\tDeltaVth: %-.9gmV\n\n\n\n", target, f * 1000) ;
+    fprintf (stderr, "\n\nExtrapolation at %-.9g seconds (%-.2g years):\n\t\t\t\tDeltaVth: %-.9gmV\n\n\n\n",
+        target, target / 31536000, f * 1000) ;
 
     *result = f ;