Fixes for icooling=8

src/main/eos_stamatellos.f90

@@ -20,8 +20,9 @@ module eos_stamatellos
  implicit none
  real,allocatable,public :: optable(:,:,:)
  real,allocatable,public :: Gpot_cool(:),duFLD(:),gradP_cool(:),lambda_FLD(:),urad_FLD(:) !gradP_cool=gradP/rho
+ real,allocatable,public :: ttherm_store(:),teqi_store(:)
  character(len=25), public :: eos_file= 'myeos.dat' !default name of tabulated EOS file
- logical,public :: doFLD = .True.
+ logical,public :: doFLD = .True., floor_energy = .False.
  integer,public :: iunitst=19
  integer,save :: nx,ny ! dimensions of optable read in
 
@@ -38,6 +39,8 @@ subroutine init_S07cool()
  allocate(duFLD(npart))
  allocate(lambda_fld(npart))
  allocate(urad_FLD(npart))    
+ allocate(ttherm_store(npart))    
+ allocate(teqi_store(npart))    
  urad_FLD(:) = 0d0
  open (unit=iunitst,file='EOSinfo.dat',status='replace')    
  if (doFLD) then
@@ -54,6 +57,8 @@ subroutine finish_S07cool()
  if (allocated(duFLD)) deallocate(duFLD)
  if (allocated(lambda_fld)) deallocate(lambda_fld)
  if (allocated(urad_FLD)) deallocate(urad_FLD)
+ if (allocated(ttherm_store)) deallocate(ttherm_store)
+ if (allocated(teqi_store)) deallocate(teqi_store)
  close(iunitst)
 end subroutine finish_S07cool