function log10(number)
*              =====
*     Number of decimals
*     ------------------
      n=iabs(number)
      do 10 k=1,25
      n=n/10
      log10=k
      if(number.lt.0) log10=k+1
      if(n.eq.0) return
10    continue
      stop 'log10: looping'
      end