anchor.diff:          !
           !!!indexf_ = indexf ; indexi_ = indexi ; acoef_ = acoef
           !
  -        !cycle
  -        ! 20  exit
  +        !cycle loop_tran
  +        !!
  +        !333 print*,iostat_
  +        !!
  +        !stop
  +        !
        enddo loop_tran
        !
        close(tunit)
  
positive.code:            intband = intband + sum(abscoef_RAM(1:nswap_))
             !
             if (nswap_<1) then
                if (verbose>=5) write(out,"('... Finish swap')")
                exit loop_tran
             endif
             !
  
negative.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')
            !
            if (pressure<small_.and.(w(1:3)=='VOI'.or.w(1:3)=='PSE')) then
               ! for pressure= 0 Voigt is repalced by Doppler
               w = 'DOPPLER'
               !
               write(out,"('This is P=0 atm Voigt, which will be reated as Doppler')")
               !
            endif
  
            !
            proftype = trim(w)
            !
            if (trim(w(1:5))=="LOREN") ioffset = 500
  