diff --git a/dshr/dshr_mod.F90 b/dshr/dshr_mod.F90
index 579296a9a..f253adc05 100644
--- a/dshr/dshr_mod.F90
+++ b/dshr/dshr_mod.F90
@@ -1346,7 +1346,7 @@ subroutine dshr_orbital_update(Time, logunit,  maintask, eccen, obliqr, lambm0,
     integer           :: orb_year ! orbital year for current orbital computation
     character(len=CL) :: msgstr   ! temporary
     logical           :: lprint
-    logical           :: first_time = .true.
+    logical, save     :: first_time = .true.
     character(len=*) , parameter :: subname = "(dshr_orbital_update)"
     !-------------------------------------------