anchor.diff:            !
             temp0 = real(itemp,rk)*dtemp
             !
  -          write(tunit,"(f9.1,1x,es15.4)") temp0,pf(0,itemp)
  +          write(tunit,"(f8.1,1x,f15.4)") temp0,pf(0,itemp)
             !
           enddo
           !
  
positive.code:      !
       open(unit=tunit,file=trim(output)//".life",action='write',status='replace')
       !
       do ilevelf = 1,nlines
         !
         ! write to .life-file
         !
         if ( mod(nint(2.0_rk*jrot(ilevelf)),2)==0 ) then 
           write(my_fmt,'(a)')  '(1x,i11,1x,f12.6,1x,i6,1x,f7.1,1x,es12.4,3x)'
           write(tunit,my_fmt,advance="no") indices(ilevelf),energies(ilevelf),gtot(ilevelf),nint(jrot(ilevelf)),1.0_rk/Asum(ilevelf)
         else
           write(tunit,my_fmt,advance="no") indices(ilevelf),energies(ilevelf),gtot(ilevelf),jrot(ilevelf),1.0_rk/Asum(ilevelf)
         endif
         !
         do kitem = 1,maxitems
  
negative.code:           !
            temp0 = real(itemp,rk)*dtemp
            !
            write(tunit,"(f8.1,1x,f15.4)") temp0,pf(0,itemp)
            !
          enddo
          !
  