Makes isfniycbozero variable

LLNL · Jan 11, 2024 · c613654 · c613654
1 parent d716325
commit c613654
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/bbb/bbb.v b/bbb/bbb.v
@@ -451,7 +451,7 @@ isnicore(nispmx)  integer   /1,30*0/      +input #switch for ion-density core B.
 				    #=4, use impur. source terms (impur only)
 				    #=5, set d(ni)/dy=-ni/lynicore at midp &
                                     #    ni constant poloidally
-isfniycbozero   integer /0/ +input # Switch for divergence-free fluxes on core boundary
+isfniycbozero   real /0./ +input # Switch for divergence-free fluxes on core boundary
                     #=0, allows divergence-free fluxes to modify net core flux
                     #=1, redistributes fluxes due to divergence-free term 
                     #    without affecting the net core boundary flux

diff --git a/bbb/oderhs.m b/bbb/oderhs.m
@@ -3346,7 +3346,7 @@ ccc          iysptrx is the last closed flux surface (see S.R. nphygeo)
             end do
             fniycboave = fniycboave / (ixpt2(1) - ixpt1(1))
             do ix = ixpt1(1)+1, ixpt2(1)
-              fniycbo(ix, ifld) = fniycbo(ix, ifld) - fniycboave
+              fniycbo(ix, ifld) = fniycbo(ix, ifld) - isfniycbozero*fniycboave
             end do
           end do
       else if (isfniycbozero .lt. 0) then