Use diag print level for output in ssids analyse

Assigning these informational prints to the "warnings and errors" print
level was not quite right. Now using the "basic diagnostics" level
instead. Therefore also using unit_diagnostics now, like elsewhere in
SSIDS' analyse phase, e.g. subroutines analyse_double and
ssids_analyse_coord_double.

src/ssids/anal.f90

@@ -364,8 +364,9 @@ subroutine find_subtree_partition(nnodes, sptr, sparent, rptr, options, &
        if (st .ne. 0) return
     end do
 
-    if (options%print_level >= 0) then
-       print *, "[find_subtree_partition] load_balance = ", best_load_balance
+    if ((options%print_level .ge. 1) .and. (options%unit_diagnostics .ge. 0)) then
+       write (options%unit_diagnostics,*) &
+            "[find_subtree_partition] load_balance = ", best_load_balance
     end if
     ! Consolidate adjacent non-children nodes into same part and regen exec_alloc
     !print *
@@ -1014,12 +1015,14 @@ subroutine analyse_phase(n, ptr, row, ptr2, row2, order, invp, &
     if (st .ne. 0) go to 100
 
     ! Sort out subtrees
-    if (options%print_level >= 0) then
-       print *, "Input topology"
+    if ((options%print_level .ge. 1) .and. (options%unit_diagnostics .ge. 0)) then
+       write (options%unit_diagnostics,*) "Input topology"
        do i = 1, size(akeep%topology)
-         print *, "Region ", i, " with ", akeep%topology(i)%nproc, " cores"
+         write (options%unit_diagnostics,*) &
+              "Region ", i, " with ", akeep%topology(i)%nproc, " cores"
          if(size(akeep%topology(i)%gpus).gt.0) &
-            print *, "---> gpus ", akeep%topology(i)%gpus
+            write (options%unit_diagnostics,*) &
+                 "---> gpus ", akeep%topology(i)%gpus
        end do
     end if
     call find_subtree_partition(akeep%nnodes, akeep%sptr, akeep%sparent,           &