anchor.diff:      pset%evtype=evtype
   
       has_ssp(:) = 0
  -
  +    has_ssp_age(:,:) = 0
  +    
     end subroutine
   
     subroutine set_csp_params(smooth_velocity0,vactoair_flag0,redshift_colors0,&
  
positive.code:     integer, intent(in) :: zi
      pset%zmet = zi
      call ssp_gen(pset, mass_ssp_zz(:,zi),lbol_ssp_zz(:,zi),&
                   spec_ssp_zz(:,:,zi))
      has_ssp(zi) = 1
  
    end subroutine
  
    subroutine get_ssp_spec(ns,n_age,n_z,ssp_spec_out,ssp_mass_out,ssp_lbol_out)
  
negative.code:     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
  
  