Update flush implementation, fix unit number bug and check that

files are open at the end of the run before flushing.

configuration/driver/icedrv_MAIN.F90

@@ -35,6 +35,7 @@ program icedrv
       implicit none
 
       integer n
+      logical openflag
       character(len=*), parameter :: subname='(icedrv)'
 
       !-----------------------------------------------------------------
@@ -55,15 +56,24 @@ program icedrv
 
       write(ice_stdout, *) "ICEPACK COMPLETED SUCCESSFULLY "
 
-      call icedrv_system_flush(ice_stdout)
-      close (ice_stdout)
+      inquire(unit=ice_stdout,opened=openflag)
+      if (openflag) then
+         call icedrv_system_flush(ice_stdout)
+         close (ice_stdout)
+      endif
 
-      call icedrv_system_flush(nu_diag)
-      close (nu_diag)
+      inquire(unit=nu_diag,opened=openflag)
+      if (openflag) then
+         call icedrv_system_flush(nu_diag)
+         close (nu_diag)
+      endif
 
       do n = 1, nx
-         call icedrv_system_flush(nu_diag_out+n+1)
-         close (nu_diag_out+n-1)
+         inquire(unit=nu_diag_out+n-1,opened=openflag)
+         if (openflag) then
+            call icedrv_system_flush(nu_diag_out+n-1)
+            close (nu_diag_out+n-1)
+         endif
       enddo
 
       end program icedrv

configuration/driver/icedrv_diagnostics.F90

@@ -267,13 +267,12 @@ subroutine runtime_diags (dt)
              write(nu_diag_out+n-1,901) 'isotopic conc chg    = ',pdiso(n,k),k
           enddo
         endif
+        call icedrv_system_flush(nu_diag_out+n-1)
       end do
 899   format (43x,a24)
 900   format (a25,2x,f24.17)
 901   format (a25,2x,f24.17,i6)
 
-      call icedrv_system_flush(nu_diag_out+n-1)
-
       end subroutine runtime_diags
 
 !=======================================================================