diff --git a/src/netcdf/netcdf.f90 b/src/netcdf/netcdf.f90 index 5e05c68e5..ce4959005 100644 --- a/src/netcdf/netcdf.f90 +++ b/src/netcdf/netcdf.f90 @@ -63,20 +63,14 @@ module function netcdf_get_old_t_final_system(self, param) result(old_t_final) ! Result real(DP) :: old_t_final ! Internals - !class(swiftest_nbody_system), allocatable :: tmpsys - !class(swiftest_parameters), allocatable :: tmpparam - !integer(I4B) :: ierr integer(I4B) :: itmax real(DP), dimension(:), allocatable :: tvals - - call check( nf90_open(param%outfile, NF90_NOWRITE, param%nciu%ncid) ) + call param%nciu%open(param) call check( nf90_inquire_dimension(param%nciu%ncid, param%nciu%time_dimid, len=itmax) ) allocate(tvals(itmax)) - call check( nf90_inq_varid(param%nciu%ncid, TIME_DIMNAME, param%nciu%time_varid) ) call check( nf90_get_var(param%nciu%ncid, param%nciu%time_varid, tvals(:)) ) - call check( nf90_close(param%nciu%ncid) ) old_t_final = tvals(itmax)