diff --git a/code/MOR b/code/MOR
index 13d20abc..dd857130 100644
--- a/code/MOR
+++ b/code/MOR
@@ -121,7 +121,7 @@ c     parameters and common blocks for MOR
 
       integer ad_nsteps,ad_iostep,ad_step,ad_qstep
       common /moripar/ ad_nsteps,ad_iostep,ad_step,ad_qstep,inus,
-     $                 navg_step,mb,nb,nbo,nplay,nintp,iaug
+     $                 navg_step,mb,nb,nbo,nplay,nintp,iaug,nbat
 
       common /morivars/ i0,j0,k0,i1,j1,k1,nns,ns,nskip,navg,ncloc,npr,
      $                  npart,isolve,ic1,ic2,jc1,jc2,kc1,kc2,idirf
diff --git a/code/mpar.f b/code/mpar.f
index 083c8c1d..441d441e 100644
--- a/code/mpar.f
+++ b/code/mpar.f
@@ -157,6 +157,9 @@ subroutine mpar_read(ierr)
       call finiparser_getbool(i_out,'general:setbases',ifnd)
       if (ifnd.eq.1) ifsetbases=i_out.eq.1
 
+      call finiparser_getdbl(d_out,'general:nbat',ifnd)
+      if (ifnd.eq.1) nbat=nint(d_out)
+
       ibuoy=0
 
       call finiparser_getdbl(d_out,'buoyancy:magnitude',ifnd)
@@ -534,6 +537,7 @@ subroutine bcastmpar
       call bcast(tbarrseq,isize)
       call bcast(nintp,isize)
       call bcast(iaug,isize)
+      call bcast(nbat,isize)
 
       ! reals
 
diff --git a/code/pod.f b/code/pod.f
index 65c697f0..f9352c92 100644
--- a/code/pod.f
+++ b/code/pod.f
@@ -824,7 +824,7 @@ subroutine gengraml2(gram,s,ms,mdim)
          if (nio.eq.0) write (6,1) j,gram(1,j)
       enddo
 
-      call breduce(gram,ms*ms,1)
+      call breduce(gram,ms*ms,nbat)
 
       if (nio.eq.0) write (6,*) 'exiting gengraml2'
 
diff --git a/code/rom.f b/code/rom.f
index 31317467..771db656 100644
--- a/code/rom.f
+++ b/code/rom.f
@@ -717,6 +717,7 @@ subroutine mor_init_params
       gz=0.
 
       nintp=0
+      nbat=1
 
       if (nio.eq.0) write (6,*) 'exiting mor_init_params'
 
@@ -1343,7 +1344,7 @@ subroutine setc(cl,fname)
                      rtmp1(i,1)=vlsc2(tb(1,i,1),cu,n)
                   endif
                enddo
-               call breduce(rtmp1,nb,1)
+               call breduce(rtmp1,nb,nbat)
                do i=1,nb
                   call setc_local(cl,rtmp1(i,1),
      $               ic1,ic2,jc1,jc2,kc1,kc2,i,j,k)