anchor.diff:     logical :: eof
      character(len=cl) :: w
      !
  -   integer(ik) :: imol
  +   integer(ik) :: imol,iostat_
      real(rk)    :: gf,gi,mem_t
      character(55) ch_q,ch_broad
      !
  
positive.code:        case('GAUSSIAN','GAUSS','DOPPL','DOPPLER','RECT','BOX','BIN','STICKS','STICK','GAUS0','DOPP0',&
              'LOREN','LORENTZIAN','LORENTZ','MAX','VOIGT','PSEUDO','PSE-ROCCO','PSE-LIU','VOI-QUAD','PHOENIX')
            !
            proftype = trim(w)
            !
            if (trim(w(1:5))=="LOREN") ioffset = 500
  
negative.code:      !
     else
       !
       if (verbose>=2) print('(/"Total intensity = ",es16.8,", T = ",f18.4)'), intband,temp
       !
     endif
     !
  