anchor.diff:           if (partfunc_do) then
              if (j/=j0) then 
                 !
  -              if (specttype(1:4)=='PART') then 
  +              if (proftype(1:4)=='PART') then 
                   if (verbose>=1.and.npoints<1000) print('("!",f8.1,3(1x,'//npoints_fmt//'es20.8))'),jrot(i),pf(ipartf,1:npoints)
                 else 
                   if (verbose>=4) print('("|",f8.1,1x,es16.8)'),jrot(i),partfunc
  
positive.code:   integer(ik) :: N_omp_procs=1
    !
    integer(ik)   :: GNS=1,npoints=1001,nchar=1,nfiles=1,ipartf=0,verbose=3,ioffset = 10,iso=1
    real(rk)      :: temp=298.0,partfunc=-1.0,freql=-small_,freqr= 20000.0,thresh=1.0d-90,halfwidth=1e-2,meanmass=1.0,maxtemp=10000.0
    real(rk)      :: voigt_gamma = 0.05, voigt_n = 0.44, offset = -25.0, pressure = 1.0_rk 
    real(rk)      :: enermax = 1e6, abscoef_thresh = 1.0d-50, abundance = 1.0d0
    real(rk)      :: S_crit = 1e-29      ! cm/molecule, HITRAN cut-off paramater 
  
negative.code:       ! In case of specttype = partfunc the partition function will be computed for a series of temperatures. 
        ! npoints stands for the number of temperature steps, and the frequency limits as temperature limits
        !
        if (proftype(1:4)=='PART'.or.proftype(1:4)=='COOL') then
          !
          if (trim(enrfilename)=="none") then 
             !
  