Bug fixes

src/main/cooling_radapprox.f90

@@ -112,7 +112,7 @@ subroutine radcool_update_du(i,xi,yi,zi,rhoi,ui,duhydro,Tfloor)
  use physcon,  only:steboltz,pi,solarl,kb_on_mh,piontwo,rpiontwo
  use units,    only:umass,udist,unit_density,unit_ergg,utime,unit_pressure
  use eos_stamatellos, only:getopac_opdep,getintenerg_opdep,gradP_cool,&
-          ttherm_store,teqi_store,opac_store
+          ttherm_store,ueqi_store,opac_store
  use part,       only:xyzmh_ptmass,igas,eos_vars,iTemp
  integer,intent(in) :: i
  real,intent(in) :: xi,yi,zi,rhoi

src/main/eos_stamatellos.f90

@@ -20,7 +20,7 @@ module eos_stamatellos
  implicit none
  real,allocatable,public :: optable(:,:,:)
  real,allocatable,public :: gradP_cool(:)!gradP_cool=gradP/rho
- real,allocatable,public :: ttherm_store(:),teqi_store(:),opac_store(:),duSPH(:)
+ real,allocatable,public :: ttherm_store(:),ueqi_store(:),opac_store(:),duSPH(:)
  character(len=25), public :: eos_file= 'eos_lom.dat' !default name of tabulated EOS file
  logical,public :: floor_energy = .False.
  integer,public :: iunitst=19

src/main/step_leapfrog.F90

@@ -456,7 +456,7 @@ subroutine step(npart,nactive,t,dtsph,dtextforce,dtnew)
 !$omp shared(dustevol,ddustevol,use_dustfrac) &
 !$omp shared(dustprop,ddustprop,dustproppred) &
 !$omp shared(xyzmh_ptmass,vxyz_ptmass,fxyz_ptmass,nptmass,massoftype) &
-!$omp shared(dtsph,ieos,ufloor,icooling) &
+!$omp shared(dtsph,ieos,ufloor,icooling,Tfloor) &
 !$omp shared(ibin,ibin_old,ibin_sts,twas,timei,use_sts,dtsph_next,ibin_wake,sts_it_n) &
 !$omp shared(ibin_dts,nbinmax) &
 !$omp private(dti,hdti) &