function rsign(a,b)
      if (b.ge.0.) then
         rsign=abs(a)
      else
         rsign=-abs(a)
      endif
      return
      end