00001 *
00002 *
00003 *
00004 *
00005 *
00006 *
00007 *
00008 *
00009 *
00010 *
00011 *
00012 *
00013 * Class members:
00014 *
00015 INTEGER jlim1,jlim2
00016 **** PARAMETER( jlim1=64, jlim2=1024) ! usualy enough
00017 **** PARAMETER( jlim1=256, jlim2=4096) ! denser grid
00018 PARAMETER( jlim1=512, jlim2=8192) ! even denser grid, standard
00019 **** PARAMETER( jlim1=1024, jlim2=16384) ! very dense grid but slow
00020
00021 DOUBLE PRECISION xx,yy,zint,zsum
00022 DOUBLE PRECISION swt,sswt
00023 INTEGER nevs,jmax
00024 COMMON /Vesk1/
00025 * xx( jlim2+1), ! x-grid
00026 * yy( jlim2+1), ! f(x) values
00027 * zint(jlim2+1), ! cumulative integral from x=xx(1)=0 to x=xx(j)
00028 * zsum, ! estimated integral (zint is normalized to one)
00029 * swt, ! sum of weights
00030 * sswt, ! sum of squares of weights
00031 * nevs, ! serial counter
00032 * jmax ! grid length
00033 SAVE
00034 *
00035 *
00036 * Class procedures:
00037 *
00038 * SUBROUTINE Vesk1_Initialize
00039 * SUBROUTINE Vesk1_Make
00040 * SUBROUTINE Vesk1_Finalize
00041 *
00042 *
00043 *
00044 *
00045 *