anchor.diff:            !
             if (mod(npoints,2)==0) npoints = npoints + 1
             !
  -        case ("ABSORPTION","EMISSION","LIFETIME","PHOENIX")
  +        case ("ABSORPTION","LIFETIME","OSCILLATOR")
             !
             specttype = trim(w)
             !
          case ("MEM","MEMORY")
             !
            call readf(memory_limit)
  +         !
  +       case ("OMP_PROCS","NPROCS","OMP_NUM_PROCS")
  +          !
  +         call readi(N_omp_procs)
             !
           case ("ISO","ISOTOPE")
             !
  
positive.code:           !
            call readi(verbose)
            !
          case ("PARTFUNC","PARTITION-FUNCTION","COOLING")
            !
            specttype = trim(w)
            !
  
negative.code:    !
     allocate(intens(npoints),stat=info)
     call ArrayStart('intens',info,size(intens),kind(intens))
     intens = 0
     !
     forall(ipoint=1:npoints) freq(ipoint)=freql+real(ipoint-1,rk)*dfreq
     !
  