Initial values of variables in common blocks as specified in
the BLOCKDATA subroutine.
DOUBLE PRECISION
PRECIS (Set in INIT) Based on machine dependent
precision variables
RANGE 1.D300 Upper limit of double precision values
REAL
DFMIN 2. Controls number of degrees of freedom
SRMIN .01 flat spots in binding constraints
ALPST(2) 2*0. Controlling regularization
EXMAX The maximum single precision exponent (init)
GMNMX(2) 2*0. Range for gird points
PLEVEL(2,2) 4*.5 probabilities to use in selecting solutions
RSVMNX(2,2) 2*1., 2*0. Specifies alpha values for regularization
RUSER(551) 551*0. user defined contol variables to be used in any
user subroutine or to store results used in
later calculations
SRANGE 1.E35 Upper range for single precision values
INTEGER
IGRID 2 Quadrature Grid Specification:
spacing of solution grid
IQUAD 3 Quadrature Grid Specification:
type of calculation (direct, approximation)
IUNIT -1 scratch file unit
IWT 1 The type of weighted analysis
LINEPG 60 output lines/page
MIOERR 5 number of messages about input errors before
aborting
MPKMOM 5 the number of individual peaks for which
moments are to be computed
MQPITR 35 max. iterations in peak-constrained analysis
NEQ 0 Equality constraints or not
NERFIT 10 Number of residuals to compute
NG 31 Number of Grid points
NINTT 1 Whether Y's and T's are computed or input
NLINF 0 Second Sum or not
NORDER 2 Setting the regulizer
IAPACK(6) ' PCS-1' (This a label indicating version. We probably
should change it.)
ICRIT(2) 2*1 chosen solution to reject using probability.
IFORMT(70) '(5E15.6)' Format for reading Ts
IFORMW(70) '(5E15.6)' Format for reading Ws
IFORMY(70) '(5E15.6)' Format for reading Ys
IPLFIT(2) 2*2 controls when plots of the fit to data are made
IPLRES(2) 2*2 controls when weighted residuals will be plotted
IPRINT(2) 2*4 controls when plots of solution will be printed
ITITLE(80) title of run
IUSER(50) 9*0, 2, 7*0, 50, 32*0
user defined contol variables to be used in any
user subroutine or to store results used in
later calculations
IUSROU(2) 2*0 controls when USEROU is called
LSIGN(4,4) 16*0 starting grid point for peak-constrained analysis
MOMNMX(2) -1, 3 controls the moments that are calculated
NENDZ(2) 2, 2 Number of extra zeros to be placed at
ends of grid
NFLAT(4,2) 8*0 Max. number of Peak-constrained solutions
NGL (CONTIN, based on NG and NLINF)
total number of unknowns
NGLP1 (CONTIN, based on NGL + 1)
total number of unknowns +1
NIN 5 Input unit number for reading
NINEQ (CONTIN 0, SETNNG) Used for setting inequality
constraints
NNSGN(2) 2*0 Which Peak-Contrained analysis to be
performed
NOUT 6 Output unit number for writing
NQPROG(2) 6, 6 Number of alpha values in grid
NSGN(4) 4*0 Number of monotonic regions for peak-constrined
analysis
NY (WRITYT, SETWT) number of y's
LOGICAL
DOCHOS .TRUE. Whether to output the chosen solution
once again at the end of the analysis
DOMOM .TRUE. Whether to compute moment whenever
solution is printed
DOUSIN .TRUE. Whether to call USERIN
DOUSNQ .FALSE. Whether to call USERNQ
(inequality constraints)
LAST .TRUE. Whether this is the last input set
NEWPG1 .FALSE. Whether to start a new page at start of
new run
NONNEG .TRUE. Whether to constrain S(gm) to be
nonnegative
ONLY1 .TRUE. whether to compute second curve to be
plotted with solution
PRWT .TRUE. Whether to output suqre roots of
the least square weights
PRY .TRUE. Whether to output T and Y arrays
SIMULA .FALSE. Whether to compute simulated data or
to read in Y's
LUSER(30) 30*.FALSE.
user defined contol variables to be used in any
user subroutine or to store results used in
later calculations