stop --> mesa_error

Author: pmocz

num/private/mod_newton.f90

@@ -639,11 +639,11 @@ subroutine do_newton(nz, nvar, x1, xold1, AF1, ldAF, neq, matrix_type, mljac, mu
             exit
          end if
          if (is_bad(residual_norm)) then
-            call oops('residual_norm is a a bad number (NaN or Infinity)')
+            call oops('residual_norm is a bad number (NaN or Infinity)')
             exit
          end if
          if (is_bad(max_residual)) then
-            call oops('max_residual is a a bad number (NaN or Infinity)')
+            call oops('max_residual is a bad number (NaN or Infinity)')
             exit
          end if
          residual_norm = abs(residual_norm)

star/private/conv_premix.f90

@@ -731,7 +731,7 @@ subroutine create_zone_info_ (s, kc_t, kc_b, zi)
        if (in_conv .AND. s%mlt_mixing_type(kf) /= convective_mixing) then
 
           ! Transitioning out of a convective zone; complete
-          ! definition of new zone info and add to to the table
+          ! definition of new zone info and add to the table
 
           zi_new%kc_t = kf
           zi_new%vc_t = s%mlt_vc(kf+1)

star/private/diffusion_procs.f90

@@ -1447,7 +1447,7 @@ subroutine set1_g_rad_mombarg( &
             return
          end if
 
-         if (ierr /= 0) stop 'set1_g_rad'  !return
+         if (ierr /= 0) call mesa_error(__FILE__,__LINE__,'set1_g_rad')  !return
 
          do ii = 1, 17  !kk
             do i = 1, nc

star/private/evolve_support.f90

@@ -197,7 +197,7 @@ subroutine set_current_to_old(s)
             call set_qs(s, s% nz, s% q, s% dq, ierr)
             if (ierr /= 0) then
                write(*,*) 'set_current_to_old failed in set_qs'
-               stop
+               call mesa_error(__FILE__,__LINE__,'set_current_to_old')
             end if
             call set_m_and_dm(s)
             call set_dm_bar(s, s% nz, s% dm, s% dm_bar)
@@ -212,7 +212,7 @@ subroutine set_current_to_old(s)
                   if (is_bad_num(s% omega(k)) .or. abs(s% omega(k)) > 1d50) then
                      if (s% stop_for_bad_nums) then
                         write(*,2) 's% omega(k)', k, s% omega(k)
-                        stop 'set_current_to_old'
+                        call mesa_error(__FILE__,__LINE__,'set_current_to_old')
                      end if
                   end if
                end do

star/private/mesh_adjust.f90

@@ -344,14 +344,14 @@ subroutine do_mesh_adjust( &
                density_new, energy_new, op_err)
             if (op_err /= 0) then
                write(*,2) 'failed for do1_lnT', k
-               stop
                write(message,*) 'do1_lnT for k', k
                ierr = op_err
+               call mesa_error(__FILE__,__LINE__,'mesh adjust: do1_lnT')
             end if
             if (is_bad(energy_new(k)) .or. is_bad(dq(k))) then
                write(*,2) 'energy_new', k, energy_new(k)
                write(*,2) 'dq', k, dq(k)
-               stop ''
+               call mesa_error(__FILE__,__LINE__,'mesh adjust: bad energy_new or dq')
             end if
 
          end do

star/private/rsp_def.f90

@@ -101,7 +101,7 @@ module rsp_def
       real(dp), pointer :: xa(:)
       real(dp) :: X, Z, Y, abar, zbar, z53bar, XC, XN, XO, Xne
 
-      ! these for for rsp.f90 period and work calculations
+      ! these are for rsp.f90 period and work calculations
       real(dp) :: ETOT, EGRV, ETHE, EKIN, EDE_start, ECON, &
          TE, ELSTA, TEFF, E0, TT1, TE_start, T0, UN, ULL, &
          VMAX, RMAX, LMAX, LMIN, EKMAX, EKMIN, EKMAXL, EKDEL, &

star/private/star_solver.f90

@@ -465,7 +465,7 @@ subroutine do_solver_work
                end if
 
                if (is_bad_num(residual_norm)) then
-                  call oops('residual_norm is a a bad number (NaN or Infinity)')
+                  call oops('residual_norm is a bad number (NaN or Infinity)')
                   if (s% stop_for_bad_nums) then
                      write(*,1) 'residual_norm', residual_norm
                      call mesa_error(__FILE__,__LINE__,'solver')
@@ -474,7 +474,7 @@ subroutine do_solver_work
                end if
 
                if (is_bad_num(max_residual)) then
-                  call oops('max_residual is a a bad number (NaN or Infinity)')
+                  call oops('max_residual is a bad number (NaN or Infinity)')
                   if (s% stop_for_bad_nums) then
                      write(*,1) 'max_residual', max_residual
                      call mesa_error(__FILE__,__LINE__,'solver')