anchor.diff:  
     end subroutine
   
  +  subroutine interp_ssp(ns,zpos,tpos,spec,mass,lbol)
  +
  +    ! Return the SSPs interpolated to the target metallicity 
  +    !(zpos) and target age (tpos)
  +    
  +    implicit none
  +
  +    integer, intent(in) :: ns
  +    double precision, intent(in) :: zpos
  +    double precision, intent(in) :: tpos
  +
  +    double precision, dimension(ns,1), intent(inout) :: spec
  +    double precision, dimension(1), intent(inout) :: mass,lbol
  +
  +    double precision, dimension(nt) :: time
  +    
  +    integer :: zlo,zmet,tlo
  +
  +    zlo = max(min(locate(log10(zlegend/0.0190),zpos),nz-1),1)
  +    time = timestep_isoc(zlo,:)
  +    tlo = max(min(locate(time,tpos),nt-1),1)
  +
  +    do zmet=zlo,zlo+1
  +       if ((has_ssp_age(zmet,tlo) .eq. 0) .or. (has_ssp_age(zmet,tlo+1) .eq. 0)) then
  +          pset%ssp_gen_age = 0
  +          pset%ssp_gen_age(tlo:tlo+1) = 1
  +          call ssp(zmet)
  +          pset%ssp_gen_age = 1
  +       endif
  +    enddo
  +
  +    call ztinterp(zpos,spec,lbol,mass,tpos=tpos)
  +
  +    end subroutine
  +
  +  subroutine smooth_spectrum(ns,wave,spec,sigma_broad,minw,maxw)
  +    
  +    ! Smooth the spectrum by a gaussian of width sigma_broad
  +    
  +    implicit none
  +    integer, intent(in) :: ns
  +    double precision, intent(in) :: sigma_broad,minw,maxw
  +    double precision, dimension(ns), intent(in) :: wave
  +    double precision, dimension(ns), intent(inout) :: spec
  +    
  +    call smoothspec(wave,spec,sigma_broad,minw,maxw)
  +
  +  end subroutine
  +
     subroutine stellar_spectrum(ns,mact,logt,lbol,logg,phase,ffco,lmdot,wght,spec_out)
       
       ! Get a stellar spectrum for a given set of parameters
  
positive.code:     
    end subroutine 
    
    
    subroutine get_setup_vars(cvms)
  
  
negative.code:        call ztinterp(zpos,spec,lbol,mass,zpow=pset%pmetals)
         call compsp(0,1,outfile,mass,lbol,spec,pset,ocompsp)
      endif
  
    end subroutine
  
    subroutine get_spec(ns,n_age,spec_out)
  
      ! Get the grid of spectra for the computed CSP at all ages.
  