anchor.diff:  
          rho = p/e/(gamma - 1.0_dp_t)
   
  -       if (present(T)) T = mu*m_amu*e/(1.5_dp_t*k)
  +       if (present(T)) T = mu*m_amu*e*(gamma - 1.0_dp_t)/k
          if (present(s)) s = c_v*log(p/rho**gamma) 
   
       else if (input == eos_input_ps) then
  
positive.code: 
         e = p / (rho * (gamma - 1.0_dp_t))
  
         if (present(T)) T = mu*m_amu*e/(1.5_dp_t*k)
         if (present(s)) s = c_v*log(p/rho**gamma) 
  
      else if (input == eos_input_e) then
  
negative.code: 
         e = p / (rho * (gamma - 1.0_dp_t))
  
         if (present(T)) T = mu*m_amu*e*(gamma - 1.0_dp_t)/k
  
      else
         print *, 'ERROR: invalid input in EOS'
  