fabm_rest_def Subroutine

      subroutine fabm_rest_def(iou,it)
         !! define FABM restart variables in NetCDF context
      integer, intent(in) :: iou
         !! NetCDF file id
      integer, intent(in) :: it(:)
         !! time dimension

         integer :: n

         associate(VAR => model%interior_state_variables)
         do n=1,npelagic
            call defvar (trim(VAR(n)%name)//'1', iou, 4, it,
     &                 VAR(n)%minimum, VAR(n)%maximum, ' ', 'D',
     &                 trim(VAR(n)%long_name)//' at tau', ' ', 
     &                 trim(VAR(n)%units))
            call defvar (trim(VAR(n)%name)//'2', iou, 4, it,
     &                 VAR(n)%minimum, VAR(n)%maximum, ' ', 'D',
     &                 trim(VAR(n)%long_name)//' at tau+1', ' ', 
     &                 trim(VAR(n)%units))
         end do
         end associate

      end subroutine fabm_rest_def