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:     double precision, dimension(ns,1), intent(inout) :: spec
      double precision, dimension(1), intent(inout) :: mass,lbol
  
      integer :: zlo,zmet
  
      zlo = max(min(locate(log10(zlegend/0.0190),zpos),nz-1),1)
      do zmet=zlo,zlo+1
         if (has_ssp(zmet) .eq. 0) then
            call ssp(zmet)
         endif
      enddo
  
  
negative.code:   !f2py intent(hide) has_ssp
    integer, dimension(nz) :: has_ssp=0
  
    !f2py intent(hide) has_ssp_age
    integer, dimension(nz,nt) :: has_ssp_age=0
  
    
  contains
  
    subroutine setup(compute_vega_mags0)
  