anchor.diff:                 intens(ipoint)=intens(ipoint)+voigt_humlicek(tranfreq_i,tranfreq,halfwidth,halfwidth0)*abscoef
                  !
               enddo
  -            !$omp end parallel do             !
  +            !$omp end parallel do      
  +            !
           case ('MAX');
             !
             ipoint =  max(nint( ( tranfreq-freql)/dfreq )+1,1)
  
positive.code:           !
            call readi(ioffset)
            !
         case ("MASS")
            !
            call readf(meanmass)
  
negative.code:    !
     maxitems = maxitems-4
     !
     if (verbose>=3) print*,"nlevels: (total)",nlevels," (selected)",nlines
     !
     allocate(energies(nlines),Jrot(nlines),gtot(nlines),indices(nlevels),stat=info)
     call ArrayStart('energies',info,size(energies),kind(energies))
     call ArrayStart('Jrot',info,size(Jrot),kind(Jrot))
     call ArrayStart('gtot',info,size(gtot),kind(gtot))
     call ArrayStart('indices',info,size(indices),kind(indices))
     !
     if (trim(specttype)=='LIFETIME') THEN 
       !
       ! allocate the matrix for the sum of the A-coeffs 
       !
       allocate(Asum(nlines),stat=info)
       call ArrayStart('Asum',info,size(Asum),kind(Asum))
       !
       Asum = -1.0_rk
       !
       write(my_fmt,'(a)')  '(1x,i11,1x,f12.4,1x,f5.1,1x,es16.8,5x)'
       !
     endif
     !
     allocate(quantum_numbers(maxitems,iline),stat=info)
     call ArrayStart('quantum_numbers',info,size(quantum_numbers),kind(quantum_numbers))
     !
     if (verbose>=3) call MemoryReport
     !
     quantum_numbers = ""
     !
     indices = 0
  