anchor.diff:            !
             call readf(pressure)
             !
  +        case ("NQUADRATURE","NQUAD")
  +          !
  +          call readi(nquad)
  +          !
  +          if (nquad>nquadmax) call report ("nquad is > nquadmax, which then should be icreased ",.true.)
  +          !
           case ("PF","QSTAT")
             !
             call readf(partfunc)
  
positive.code:         x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
          !
          !   half width for Doppler profiling
          if (proftype(1:5)=='VOIGT'.or.proftype(1:3)=='PSE') then
            if (tranfreq<small_) cycle
            halfwidth0=dpwcoef*tranfreq
            x0 = sqrt(ln2)/halfwidth0*dfreq*0.5_rk
  
negative.code:      Asum = -1.0_rk
       !
       write(my_fmt,'(a)')  '(1x,i11,1x,f12.4,1x,f5.1,1x,es16.8,5x)'
       !
     endif
     !
     ! prepare the quadratures (Gauss-Hermite)
     !
     if (trim(proftype(1:5))=='VOI-Q') then 
       write(out,"(/'Number of quadrature poitns (Gauss-Hermite) = ',i7/)") nquad
       !
       allocate(weight(nquad),abciss(nquad),bnormq(nquad),stat=info)
       call ArrayStart('weight-abciss',info,size(weight),kind(weight))
       call ArrayStart('weight-abciss',info,size(abciss),kind(abciss))
       call ArrayStart('weight-abciss',info,size(bnormq),kind(bnormq))
       !
       call gauher(abciss,weight,nquad)
       !
     endif
     !
  