From c613654ee5514a09058c5a0f0ad9e1a533892d25 Mon Sep 17 00:00:00 2001 From: Andreas Holm <60451789+holm10@users.noreply.github.com> Date: Thu, 11 Jan 2024 15:41:03 -0800 Subject: [PATCH] Makes isfniycbozero variable --- bbb/bbb.v | 2 +- bbb/oderhs.m | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/bbb/bbb.v b/bbb/bbb.v index ff53ebd5..76d27562 100755 --- 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 index 8b28674f..dfa7392c 100755 --- 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