Makes UEDGE naggier when using old models

LLNL · Oct 27, 2023 · f7f6042 · f7f6042
1 parent 4b6df21
commit f7f6042
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 0 deletions.
diff --git 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
@@ -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)