Added variable to define where output is written for normal versus nested models

src/write_global_output.f90

@@ -69,6 +69,7 @@ character*5 outputtype
 integer i,j,is,k
 integer iproc,nproc,ierr
 character*4 cistep,citer
+character(len=7) :: outdir
 
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
@@ -79,11 +80,16 @@ call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 if (iproc.eq.0) then
    call int_to_char (cistep,4,istep)
    select case (outputtype)
+      if (params%nest) then
+        outdir='NESTOUT'
+      else
+        outdir='OUT'
+      endif
       case ('debug')
          call int_to_char (citer,4,iter)
-         open (9,file='OUT/time_'//cistep//'_'//citer//'.bin',status='unknown',form='unformatted')
+         open (9,file=trim(outdir)//'/time_'//cistep//'_'//citer//'.bin',status='unknown',form='unformatted')
       case ('final')
-         open (9,file='OUT/time_'//cistep//'.bin',status='unknown',form='unformatted')
+         open (9,file=trim(outdir)//'/time_'//cistep//'.bin',status='unknown',form='unformatted')
 !         open (1234,file='OUT/time_'//cistep//'.dat',status='replace')
       case default
          call stop_run ('pb with argument in write_global_output$')