00001
00002 c
00003 c Copyright (c) 1986,1987,1988,1989,1990,1991,1992,1993,
00004 c by Steve McMillan, Drexel University, Philadelphia, PA.
00005 c
00006 c All rights reserved.
00007 c
00008 c Redistribution and use in source and binary forms are permitted
00009 c provided that the above copyright notice and this paragraph are
00010 c duplicated in all such forms and that any documentation,
00011 c advertising materials, and other materials related to such
00012 c distribution and use acknowledge that the software was developed
00013 c by the author named above.
00014 c
00015 c THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR
00016 c IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED
00017 c WARRANTIES OF MERCHANTIBILITY AND FITNESS FOR A PARTICULAR PURPOSE.
00018 c
00019
00020 subroutine compoz(fpn,fge1,lpow)
00021 save
00022 c
00023 c decomposes a floating point number fpn into
00024 c fge1*10**(lpow) where: abs(fge1) is .ge. 1
00025 c
00026 if (fpn.ne.0.) then
00027 fabs = abs(fpn)
00028 lpow = alog10(fabs)
00029 fge1 = fpn*10.**(-lpow)
00030 if (abs(fge1).gt.0.99999) return
00031 fge1 = 10.*fge1
00032 lpow = lpow-1
00033 else
00034 fge1 = 0.
00035 lpow = 0
00036 end if
00037 c
00038 end
