anchor.diff:      
   contains
   	subroutine check_okay(msg,ierr)
  -		use iso_fortran_env, only : error_unit
   		character(len=*), intent(in) :: msg
   		integer, intent(inout) :: ierr
   		if (ierr /= 0) then
  
positive.code: program test_NScool
      use iso_fortran_env, only : output_unit
      use NScool_def
      use NScool_lib
  
  
negative.code:     ! (end of the epochs that were specified in the inlist).
      call get_NScool_info_ptr(NScool_id,s,ierr)    
      do i = 1, s% number_epochs
          write(output_unit,'(2(a," = ",es11.4))') 't/d',s% t_monitor(i),'; obs. Teff/K', s% Teff_monitor(i)
      end do
      
      call NScool_shutdown    
  