anchor.diff:      if (minval(depth) .lt. 0) then
         write(17, "(A)") "Depths must be positive."
         call clean_stop(0, .FALSE.)
  -      stop 1
       end if
     end if
   
  
positive.code: 
    ! Finalize MPI
    call clean_stop(nTimeSteps, .TRUE.)
  
    stop
  end program aronnax
  
  ! ------------------------------ Primary routine ----------------------------
  
negative.code:        write(17, "(A, I0)") 'nProcY = ', nProcY
      end if
      call clean_stop(0, .FALSE.)
    end if
  
    ! myid starts at zero, so index these variables from zero.
  