diff --git a/src/netcdf/netcdf.f90 b/src/netcdf/netcdf.f90
index a4ee55c8c..44ba0e58d 100644
--- a/src/netcdf/netcdf.f90
+++ b/src/netcdf/netcdf.f90
@@ -68,7 +68,7 @@ module function netcdf_get_old_t_final_system(self, param) result(old_t_final)
       real(DP), dimension(1)                    :: val
       real(DP) :: KE_orb_orig, KE_spin_orig, PE_orig, Ltmp
 
-      call param%nciu%open(param, readonly=.true.)
+      call param%nciu%open(param)
       call check( nf90_inquire_dimension(param%nciu%ncid, param%nciu%time_dimid, len=itmax) )
       call check( nf90_inquire_dimension(param%nciu%ncid, param%nciu%id_dimid, len=idmax) )
       allocate(vals(idmax))
@@ -110,7 +110,6 @@ module function netcdf_get_old_t_final_system(self, param) result(old_t_final)
       end if
 
       deallocate(vals)
-      call param%nciu%close()
       
       return
    end function netcdf_get_old_t_final_system