anchor.diff:            if (tranfreq<small_) cycle
             !
             halfwidth=dpwcoef*tranfreq
  +          !
           endif
           !
           if (offset<0) offset = ioffset*halfwidth
  
positive.code:           if (trim(w(1:5))=="LOREN") ioffset = 500
            if (trim(w(1:5))=="VOIGT") ioffset = 500
            if (trim(w(1:3))=="PSE") ioffset = 500
            !
         case ("HWHM","HALFWIDTH")
            !
  
negative.code:             !
              b = 1.0_rk/(xp-xm)
              !
              dnu_half = dfreq*0.5_rk
              !
              !$omp parallel do private(ipoint,dfreq_,xp,xm,de) shared(intens) schedule(dynamic)
              do ipoint=ib,ie
                 !
                 dfreq_=freq(ipoint)-tranfreq
                 !
                 xp = atan((dfreq_+dnu_half)/halfwidth)
                 xm = atan((dfreq_-dnu_half)/halfwidth)
                 !
                 de = (xp)-(xm)
                 !
  