Merge pull request #419 from danieljprice/implicit-radiation

danieljprice · May 14, 2023 · 4f5d4eb · 4f5d4eb
2 parents 78e4948 + 4db8476
commit 4f5d4eb
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/main/deriv.F90 b/src/main/deriv.F90
@@ -69,7 +69,7 @@ subroutine derivs(icall,npart,nactive,xyzh,vxyzu,fxyzu,fext,divcurlv,divcurlB,&
  use derivutils,     only:do_timing
  use cons2prim,      only:cons2primall,cons2prim_everything,prim2consall
  use metric_tools,   only:init_metric
- use radiation_implicit, only:do_radiation_implicit
+ use radiation_implicit, only:do_radiation_implicit,ierr_failed_to_converge
  use options,        only:implicit_radiation,implicit_radiation_store_drad
  integer,      intent(in)    :: icall
  integer,      intent(inout) :: npart
@@ -170,7 +170,7 @@ subroutine derivs(icall,npart,nactive,xyzh,vxyzu,fxyzu,fext,divcurlv,divcurlB,&
  !
  if (do_radiation .and. implicit_radiation .and. dt > 0.) then
     call do_radiation_implicit(dt,npart,rad,xyzh,vxyzu,radprop,drad,ierr)
-    if (ierr /= 0) call fatal('radiation','Failed to converge')
+    if (ierr /= 0 .and. ierr /= ierr_failed_to_converge) call fatal('radiation','Failed in radiation')
  endif
 
 !

diff --git a/src/main/radiation_implicit.f90 b/src/main/radiation_implicit.f90
@@ -46,7 +46,7 @@ module radiation_implicit
  character(len=*), parameter :: label = 'radiation_implicit'
 
  private
- public :: do_radiation_implicit
+ public :: do_radiation_implicit,ierr_failed_to_converge
 
 contains