MINICERN/packlib/kernlib/kerngen/tcgen/pscale.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/02/15 17:49:49  mclareni
   6 * Kernlib
   7 *
   8 *
   9 #include "kerngen/pilot.h"
  10       FUNCTION PSCALE (NSC,NMAX,A,NST)
  11 C
  12 C CERN PROGLIB# M215    PSCALE          .VERSION KERNFOR  1.0   681105
  13 C ORIG. 05/11/68
  14 C
  15 C--                FIND THE LARGEST POWER PSCALE=10.**NSC
  16 C-                 SUCH THAT A*PSCALE HAS AT MOST NST DIGITS ON THE
  17 C-                                    LEFT OF THE DECIMAL POINT
  18 C-                 NSC LIMITED TO NSC.LE.NMAX
  19 C
  20 C
  21       NSC=NMAX
  22       IF (A.EQ.0.) GO TO 9
  23       N=100.+ LOG10 (ABS (A))+.0001
  24       N=NST+99-N
  25       IF (N.LT.NSC) NSC=N
  26     9 PSCALE=10.**NSC
  27       RETURN
  28       END