Added output of nest parameters

src/write_global_output.f90

@@ -17,7 +17,7 @@
 !------------------------------------------------------------------------------|
 
 subroutine write_global_output (params,istep,iter,current_time,osolve,ov, &
-                                vo,surface,cl,zi,outputtype)
+                                vo,surface,cl,zi,nest,outputtype)
 
 !------------------------------------------------------------------------------|
 !(((((((((((((((( Purpose of the routine  ))))))))))))))))))))))))))))))))))))))
@@ -69,7 +69,7 @@ character*5 outputtype
 integer i,j,is,k
 integer iproc,nproc,ierr
 character*4 cistep,citer
-character(len=7) :: outdir
+character(len=5) :: outdir
 
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
@@ -79,7 +79,7 @@ call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 
 if (iproc.eq.0) then
    if (params%nest) then
-     outdir='NESTOUT'
+     outdir='SSOUT'
    else
      outdir='OUT'
    endif
@@ -202,6 +202,12 @@ if (iproc.eq.0) then
                  j=1,2**params%levelmax_oct+1),i=1,2**params%levelmax_oct+1) 
    end if
 
+   if (params%nest) then
+     ! write nested model info - dwhipp 05/11
+     write(9) nest%sselemx,nest%sselemy,nest%sselemz,nest%xminls,nest%yminls,&
+              nest%zminls
+    endif
+
    close (9)
 
 !   do i=1,osolve%nnode