anchor.diff:      ! further storage can be tacked on at end of ipar for charged id's
       
       ! error codes
  -
       integer, parameter :: negative_photosphere_gas_pressure = -2
       integer, parameter :: negative_photosphere_density = -3
       integer, parameter :: bad_composition = -4
  +    
   contains
       
   	subroutine do_get_bc09_Teff(grav, Plight, Tb, Teff, flux)
  
positive.code:        		&   res,phase,chi)
         
         P = exp(res(i_lnP))
         rpar(iPph) = P
         eta = res(i_Theta) !1.0/TpT
         Gamma = res(i_Gamma)
         call get_thermal_conductivity(rho,Teff,chi, &
             & Gamma,eta,ionic,K,which_components=cond_exclude_sf) !cond_use_only_kap)
         kappa = 4.0*onethird*arad*clight*Teff**3/rho/K% kap
         rpar(iKph) = kappa
  	   dfdlnrho = 0.0
         
         photosphere = P - 2.0*onethird*gravity/kappa
      end function photosphere
  
  end module bc09
  
negative.code:     integer, parameter :: igrav = number_comp_rpar + 1
      integer, parameter :: itau = igrav + 1
      integer, parameter :: iTeff = itau + 1
      integer, parameter :: itemp = iTeff + 1
      integer, parameter :: ipres = itemp + 1
      integer, parameter :: iKph = ipres + 1
      integer, parameter :: iChi_rho = iKph + 1
      integer, parameter :: iChi_T = iChi_rho + 1
      integer, parameter :: del_ad = iChi_T + 1
  