Added passing of nest to write_global_output subroutine

src/DOUAR.f90

@@ -785,7 +785,7 @@ do while (istep.le.params%nstep)
          !--------------------------------------------------------------------------------------
          if (params%debug>=2) then
             call show_time (total,step,inc,1,'write global output$')
-            call write_global_output(params,istep,iter,current_time,osolve,ov,vo,surface,cl,zi,'debug')
+            call write_global_output(params,istep,iter,current_time,osolve,ov,vo,surface,cl,zi,nest,'debug')
             call mpi_barrier (mpi_comm_world,ierr)
 !            call slices (params,osolve,ov,vo,sections,istep,iter,iter_nl)
          end if
@@ -1420,7 +1420,7 @@ do while (istep.le.params%nstep)
    !--------------------------------------------------------------------------------------
    !--------------------------------------------------------------------------------------
    call show_time (total,step,inc,1,'Write output$')
-   call write_global_output (params,istep,iter-1,current_time,osolve,ov,vo,surface,cl,zi,'final')
+   call write_global_output (params,istep,iter-1,current_time,osolve,ov,vo,surface,cl,zi,nest,'final')
    call mpi_barrier (mpi_comm_world,ierr)
 
    !--------------------------------------------------------------------------------------