pro colors_adjust,filin,filout ; this routine reads a colors save file and adjusts various colors ; in the following ways: ; for 3500K <= teff <= 8500K ; r-i is increased by +.006*(teff-4800.)/1000. independent of logg, logz ; i-z is increased by -.021*(teff-5000.)/1000. independent of logg, logz ; J-H is increased by -.036*(teff-5000.)/1000. independent of logg, logz ; H-K is increased by +.011*(teff-5000.)/1000. independent of logg, logz ; for teff > 7500K, the increases applied are those appropriate to 7500K ; g-d51 is increasec by -0.23*(teff-4300)/1000. for teff <= 4300K, indep of ; logg, logz ; it then saves the result to filout ; constants restore,filin ; find models cooler and hotter than 7500K sh=where(teff gt 7500.,nsh) sc=where(teff le 7500.,nsc) ; r-i c75=+.006*(7500.-4800.)/1000. colors(2,sh)=colors(2,sh)+c75 colors(2,sc)=colors(2,sc)+.006*(teff(sc)-4800.)/1000. ; i-z c75=-.021*(7500.-5000.)/1000. colors(3,sh)=colors(3,sh)+c75 colors(3,sc)=colors(3,sc)-.021*(teff(sc)-5000.)/1000. ; J-H c75=-.036*(7500.-5000.)/1000. colors(5,sh)=colors(5,sh)+c75 colors(5,sc)=colors(5,sc)-.036*(teff(sc)-5000.)/1000. ; J-K c75=+.011*(7500.-5000.)/1000. colors(6,sh)=colors(6,sh)+c75 colors(6,sc)=colors(6,sc)+.011*(teff(sc)-5000.)/1000. ; g-d51 sc1=where(teff le 4300. and logg ge 3.5,nsc1) sc2=where(teff le 4300. and logg lt 3.5,nsc2) colors(8,sc1)=colors(8,sc1)-.23*(teff(sc1)-4300.)/1000. colors(8,sc2)=colors(8,sc2)-(1.-(3.5-logg(sc2))*(0.8/3.5))*(.23)*(teff(sc2)-4300)/1000. save,colors,names,teff,logg,logz,file=filout end