function ML18_to_3D_dTeff(Teff,logg)
      real*8 Teff, logg, a(8), Shift, Teff0, logg0, ML18_to_3D_dTeff
c 
c     IN: Teff, effective temperature (K)
c         logg, surface gravity (g in cm/s2)
c
c     OUT: Teff correction (K)
c      ; Fri Jun 14 20:16:39 2013 B.F.
      
      A(1)=1.0947335E-03
      A(2)=-1.8716231E-01
      A(3)=1.9350009E-02
      A(4)=6.4821613E-01
      A(5)=-2.2863187E-01
      A(6)=5.8699232E-01
      A(7)=-1.0729871E-01
      A(8)=1.1009070E-01

      Teff0=(Teff-10000.0)/1000.00
      logg0=(logg-8.00000)/1.00000

      Shift=A(1)+(A(2)+A(7)*Teff0+A(8)*logg0)*exp(-(A(3)+A(5)*Teff0+
     *A(6)*logg0)**2*((Teff0-A(4))**2))

      ML18_to_3D_dTeff = Shift*1000.00+0.00000
      return
      end 

      function ML18_to_3D_dlogg(Teff,logg)
      real*8 Teff, logg, a(12), Shift,Teff0,logg0,ML18_to_3D_dlogg
c 
c     IN: Teff, effective temperature (K)
c         logg, surface gravity (g in cm/s2)
c
c     OUT: log g correction (g in cm/s2)
c    ; Fri Jun 14 21:09:56 2013 B.F.

      A(1)=7.5209868E-04
      A(2)=-9.2086619E-01
      A(3)=3.1253746E-01
      A(4)=-1.0348176E+01
      A(5)=6.5854716E-01
      A(6)=4.2849862E-01
      A(7)=-8.8982873E-02
      A(8)=1.0199718E+01
      A(9)=4.9277883E-02
      A(10)=-8.6543477E-01
      A(11)=3.6232756E-03
      A(12)=-5.8729354E-02

      Teff0=(Teff-10000.0)/1000.00
      logg0=(logg-8.00000)/1.00000
      Shift=(A(1)+A(5)*exp(-A(6)*((Teff0-A(7))**2)))+A(2)*
     *exp(-A(3)*((Teff0-(A(4)+A(8)*exp(-(A(9)+A(11)*
     *Teff0+A(12)*logg0)**2*((Teff0-A(10))**2))))**2))

      ML18_to_3D_dlogg = Shift
      return
      end 

      function ML17_to_3D_dTeff(Teff,logg)
      real*8 Teff, logg, a(8), Shift, Teff0, logg0, ML17_to_3D_dTeff
c 
c     IN: Teff, effective temperature (K)
c         logg, surface gravity (g in cm/s2)
c
c     OUT: Teff correction (K)
c    ; Tue Jun 18 18:58:49 2013 B.F.

      A(1)=-1.0461690E-03
      A(2)=-2.6846737E-01
      A(3)=3.0654611E-01
      A(4)=1.8025848E+00
      A(5)=1.5006909E-01
      A(6)=1.0125295E-01
      A(7)=-5.2933335E-02
      A(8)=-1.3414353E-01

      Teff0=(Teff-10000.0)/1000.00
      logg0=(logg-8.00000)/1.00000
      Shift=A(1)+(A(2)+A(5)*Teff0+(A(6)+A(7)*Teff0+A(8)*logg0)*logg0)*
     *exp(-A(3)*((Teff0-A(4))**2))

      ML17_to_3D_dTeff = Shift*1000.00+0.00000
      return
      end 

      function ML17_to_3D_dlogg(Teff,logg)
      real*8 Teff, logg, a(11), Shift, Teff0, logg0, ML17_to_3D_dlogg
c 
c     IN: Teff, effective temperature (K)
c         logg, surface gravity (g in cm/s2)
c
c     OUT: log g correction (g in cm/s2)
c    ; Tue Jun 18 18:56:55 2013 B.F.

      A(1)= 1.1922481E-03
      A(2)=-2.7230889E-01
      A(3)=-6.7437328E-02
      A(4)=-8.7753624E-01
      A(5)= 1.4936511E-01
      A(6)=-1.9749393E-01
      A(7)= 4.1687626E-01
      A(8)= 3.8195432E-01
      A(9)=-1.4141054E-01
      A(10)=-2.9439950E-02
      A(11)=1.1908339E-01

      Teff0=(Teff-10000.0)/1000.00
      logg0=(logg-8.00000)/1.00000
      Shift=A(1)+A(2)*exp(-(A(3)+(A(5)+A(7)*exp(-A(8)*
     *((Teff0-A(9))**2)))*Teff0+(A(6)+A(10)*Teff0+A(11)*logg0)*
     *logg0)**2*((Teff0-A(4))**2))
      
      ML17_to_3D_dlogg = Shift
      return
      end