From c19853c7b78f21b5dcefe6e048e0c63d0bbac68a Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Tue, 24 May 2011 17:21:46 +0000
Subject: [PATCH] Added output of nest parameters

---
 src/write_global_output.f90 | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/src/write_global_output.f90 b/src/write_global_output.f90
index b9dd73ab..eef6a5c2 100644
--- a/src/write_global_output.f90
+++ b/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
-- 
GitLab