if( type==1) then     ! Use Starr formulation
      TT = TI - ZEROC
      if    (TT < TSTARR1                   ) then
         EX = (TT*(TT*(TT*(TT*(TT*(TT*S16+S15)+S14)+S13)+S12)+S11)+S10)
      elseif(TT >= TSTARR1 .and. TT < TSTARR2) then
         W = (TSTARR2 - TT)/(TSTARR2-TSTARR1)
         EX =       W *(TT*(TT*(TT*(TT*(TT*(TT*S16+S15)+S14)+S13)+S12)+S11)+S10) &
              + (1.-W)*(TT*(TT*(TT*(TT*(TT*(TT*S26+S25)+S24)+S23)+S22)+S21)+S20)
      elseif(TT >= TSTARR2 .and. TT < TSTARR3) then
         EX = (TT*(TT*(TT*(TT*(TT*(TT*S26+S25)+S24)+S23)+S22)+S21)+S20)
      elseif(TT >= TSTARR3 .and. TT < TSTARR4) then
         W = (TSTARR4 - TT)/(TSTARR4-TSTARR3)
         EX =       W *(TT*(TT*(TT*(TT*(TT*(TT*S26+S25)+S24)+S23)+S22)+S21)+S20) &
              + (1.-W)*(TT*(TT*(TT*(TT*(TT*(TT*BI6+BI5)+BI4)+BI3)+BI2)+BI1)+BI0)
      else
         EX = (TT*(TT*(TT*(TT*(TT*(TT*BI6+BI5)+BI4)+BI3)+BI2)+BI1)+BI0)
      endif
   elseif(type==2) then !  Fit used in CAM.
      TT = ZEROC/TI
      EX = DI(0) * exp( -(DI(1)/TT + DI(2)*log(TT) + DI(3)*TT) )
   elseif(type==3) then   !  Murphy and Koop (2005, QJRMS)
      EX = exp( CI(0)+ CI(1)/TI + CI(2)*alog(TI) + CI(3)*TI )
   endif