anchor.diff:          !
           Asum = -1.0_rk
           !
  -        write(my_fmt,'(a)')  '(1x,i11,1x,f12.6,1x,i6,1x,f7.1,1x,es12.4,3x)'
  +        write(my_fmt,'(a)')  '(1x,i11,1x,f12.6,1x,i6,1x,f7.1,1x,es16.8,3x)'
           !
         endif
         !
  
positive.code:      real(rk),intent(in) :: Jf,Ji
       real(rk) :: halfwidth,gamma_,n_,f
       integer(ik) :: ispecies,Jpp,Jp
       !
       ! if no species are specified in input the hwm value will be used for halfwidth instead 
       !if (Nspecies==0) then 
       !  f = hwhm
       !  return
       !endif
       !
       halfwidth = 0
       !
       do ispecies =1,Nspecies
         !
  
negative.code:   end subroutine gauher_ark
    !C  (C) Copr. 1986-92 Numerical Recipes Software
  
  
      !
    end module spectrum
  