Skip to content

Commit

Permalink
add new var nbat which dictates batch behavior
Browse files Browse the repository at this point in the history
  • Loading branch information
kent0 committed Mar 24, 2024
1 parent cd4f92a commit aa94c75
Show file tree
Hide file tree
Showing 4 changed files with 8 additions and 3 deletions.
2 changes: 1 addition & 1 deletion code/MOR
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
4 changes: 4 additions & 0 deletions code/mpar.f
Original file line number Diff line number Diff line change
Expand Up @@ -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)
Expand Down Expand Up @@ -534,6 +537,7 @@ subroutine bcastmpar
call bcast(tbarrseq,isize)
call bcast(nintp,isize)
call bcast(iaug,isize)
call bcast(nbat,isize)

! reals

Expand Down
2 changes: 1 addition & 1 deletion code/pod.f
Original file line number Diff line number Diff line change
Expand Up @@ -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'

Expand Down
3 changes: 2 additions & 1 deletion code/rom.f
Original file line number Diff line number Diff line change
Expand Up @@ -717,6 +717,7 @@ subroutine mor_init_params
gz=0.

nintp=0
nbat=1

if (nio.eq.0) write (6,*) 'exiting mor_init_params'

Expand Down Expand Up @@ -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)
Expand Down

0 comments on commit aa94c75

Please sign in to comment.