diff --git a/bbb/bbb.v b/bbb/bbb.v index 26b92d29..27211be0 100755 --- a/bbb/bbb.v +++ b/bbb/bbb.v @@ -318,6 +318,7 @@ cfbgt real /0./ #Coef for the B x Grad(T) terms. cfjhf real /1./ #Coef for convective cur (fqp) heat flow jhswitch integer /0/ #Coef for the Joule-heating terms oldseec integer /1/ #Switch for Joule-heating bugfix +override integer /0/ #Switch to manually override checks on old models cf2ef real /0./ #Coef for ExB drift in 2-direction cfyef real /0./ #Coef for ExB drift in y-direction cftef real /0./ #Coef for ExB drift in toroidal direction diff --git a/bbb/odesetup.m b/bbb/odesetup.m index b811d195..5e894b1b 100755 --- a/bbb/odesetup.m +++ b/bbb/odesetup.m @@ -6548,6 +6548,11 @@ c_mpi call MPI_BARRIER(uedgeComm, ierr) write(*,*) "future versions of UEDGE" write(*,*) "Please set oldseec = 0 " write(*,*) "" + if (override .eq. 0) then + write(*,*) "To use the deprecated model oldseec = 1, manually" + write(*,*) "override the settings by using override=1" + call xerrab("Error: oldseec=1 used without override=1") + endif endif if (jhswitch .gt. 0) then write(*,*) " **** WARNING ****" @@ -6556,9 +6561,16 @@ c_mpi call MPI_BARRIER(uedgeComm, ierr) write(*,*) "from future versions of UEDGE." write(*,*) "Please set jhswitch = 0 " write(*,*) "" + if (override .eq. 0) then + write(*,*) "To use the deprecated model jhswitch > 0, manually" + write(*,*) "override the settings by using override=1" + call xerrab("Error: jhswitch>0 used without override=1") + endif endif + + c_mpi call MPI_BARRIER(uedgeComm, myfoo) imeth = inewton(igrid)