anchor.diff:    integer(ik)   :: nquad = 20      ! Number of quadrature points
     integer(hik)   :: N_to_RAM = -1000 ! Lines to keep in RAM
     !
  -  character(len=cl) :: specttype="absorption",enrfilename="none",intfilename(nfiles_max),proftype="DOPPL",output="output"
  +  character(len=cl) :: specttype="ABSORPTION",enrfilename="NONE",intfilename(nfiles_max),proftype="DOPPL",output="output"
     integer(ik)   :: intJvalue(nfiles_max)
     character(4) a_fmt
     character(9) b_fmt
  
positive.code:        !
         if (stick_hitran) then
           !
           write(out,"(i3,f12.6,e10.3,e10.3,f5.4,f5.3,f10.4,f4.2,f8.6)",advance="no") &
                       iso,tranfreq,abscoef,acoef,&
                       species(1)%gamma,species(2)%gamma,&
                       energyi,species(1)%n,species(1)%delta
  
negative.code:      !
       open(unit=tunit,file=trim(intfilename(i)),action='read',status='old',iostat=iostat_)
       !
       eof = .false. ! this will be used to control if the end of file is reached
       !
       if (verbose>=3) print('("  Processing ",a,"...")'), trim(intfilename(i))
  