thermal_conductivity_stainless_steel Function

  function thermal_conductivity_stainless_steel(c_pt,T,B) result(thermal_conductivity) bind(C)
      type(c_ptr), value :: c_pt
      real(dp), value :: T,B
      real(dp) :: thermal_conductivity    
      real(dp) :: tmin,tmax,TT,a0,a1,a2,a3,a4,a5
      type(stainless_steel_cfg_t), pointer :: stainless_steel_ptr
      call c_f_pointer(c_pt,stainless_steel_ptr)

      tmin=1.0d0
      tmax=1000.0d0
    
      TT=min(T,tmax)
      TT=max(TT,tmin)

      a0=-2.7522797_dp
      a1=1.47636_dp
      a2=-0.5109092_dp
      a3=0.2846192_dp
      a4=-0.063919_dp
      a5=0.0047142_dp

      thermal_conductivity = exp(a0 + a1*log(TT) + a2*log(TT)**2 + a3*log(TT)**3 + a4*log(TT)**4 + a5*log(TT)**5)
  end function thermal_conductivity_stainless_steel