**************************************
C     PROGRAM SCRNYMOD.f
C*************************************
      COMMON/EN/EN(251),P(251),B(251),DE(251)
      COMMON/PFLUX/FDIF(251),FINT(251)
      COMMON/SUNSP/W(600),T(600),NW 
                                                                                                                                                                                                      
CCCCC DATABASE FOR P AND He FLUENCES AND FLUXES
      COMMON/PARflnP/Cpn(7,10),Gpn(7,10),AGpn(7,10),BGpn(7,10)
      COMMON/PARflnHe/CHEn(7,10),GHEn(7,10),AGHEn(7,10),BGHEn(7,10)                                                                                                                                            
      COMMON/PARflxP/Cpx(7,10),Gpx(7,10),AGpx(7,10),BGpx(7,10)
      COMMON/PARflxHe/CHEx(7,10),GHEx(7,10),AGHEx(7,10),BGHEx(7,10)                                                                                                                                                                                                             COMMON/EN/EN(251),P(251),B(251),DE(251)

      COMMON/BASE/E0(13),P0(7)
      COMMON/PROB/PB(801),FXPRB(13,801),FPB(13,7),FF(801)
      COMMON/ZZZ/PROB(7),ANSEP
      COMMON/DISTFUN/FX30DST(801)
      COMMON/PAR/NE,NPB,NEN,NF,NEI(13),DM(12)
      COMMON/NNN/NOM(13),SS(13)
      COMMON/FLUX/MNX
      COMMON/PARFIN/C(7),G(7),AG(7),BG(7)
      COMMON/ZRZ/NZ,RZ(28),AQ(28)  
      COMMON/ENI/ENI(13)
      COMMON/NUMB/BN(10),PRB(7)
      DIMENSION PP(4,7),CC(7,10)
      DATA BN/1.,2.,4.,8.,16.,32.,64.,128.,256.,512./
      DATA PRB/0.9,0.7,0.5,0.3,0.1,0.032,0.01/

CCCCC    ELEMENT  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCC is possible NZ=1,2,   6,7,8, 10,11,12,13,14,15,16,17,18,19,20,
CCCCC                22,  24,  26,  28     INDICATOR BOARD - ID
	READ *,NZ
CCCCC CHOISE, YOU NEED IN FLUENCE (MNX=1) OR PEAK FLUX (MNX=2) CALCULATION
CCCCCCCCCCCCCCCCCCCCCCCCCCCC -  IB - CCCCCCCCCCCC
	READ *,MNX


CCCCC CHOISE, THE PROBABILITY WITH WHAT THE GIVEN SPECTRA IS EXCEEDED
CCCCC IF PRB=0.9 NPB=1;  IF PRB=0.7 NPB=2; IF PRB=0.5 NPB=3
CCCCC IF PRB=0.3 NPB=4;  IF PRB=0.1 NPB=5; IF PRB=0.032 NPB=9
CCCCC IF PRB=0.01 NPB=7;
	READ *,NPB
CCCCC CCCCCCCCCCCC  -  IB -  CCCCCCCCCCCCCCCCCCCCCCCCC  
CCCCC TIME INTERVAL (BEGIN 1 - END 2) (YEAR, MONTH) 
	READ *,TY1
	READ *,M1
	READ *,TY2
	READ *,M2
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  -  IB  - CCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

CCCCCC FIXATION OF THE CHOISE
      PRINT 840,NZ
  840 FORMAT(12H Particle Z=,I2)
      IF(MNX.EQ.1)PRINT 511
  511 FORMAT(34H Fluence in [part/(cm*2*MeV/nucl)] ) 
      IF(MNX.EQ.2)PRINT 512
  512 FORMAT(42H Peak Flux in [part/(cm*2*s*sr*MeV/nucl)] ) 
 
      PRINT 841,PRB(NPB)
  841 FORMAT(13H Probability=,F5.3)

CCCCC      DETERMINATION OF THE MEAN ACCEPTED SEP NUMBER AS MODEL PARAMETER
CCCCC      DETEMINATION of THE MODEL PARAMETER ANSEP (NUMBER EVENTS F>10**5)
      CALL NSEPW(TY1,M1,TY2,M2,AN)     


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCC   DETERMINATION OF THE ARBITARY SCALES
      CALL ARBITARY(AN)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCC FOR PROBABILITIES: 0.9, 0.7, 0.5, 0.25, 0.1, 0.032 AND 0.01
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

CCCCC Bazis (p or He) output parameters for exact input parameter "an" are
CCCCC here determined 
      DO 333 I=1,10
      IF(BN(I).GE.AN)GOTO 334
  333 CONTINUE
  334 NN=I-1
     
CCCCC Bazis output parameters for exact input parameter "an" is determined 
      NC=0
CCCCC CC - IT IS COMMON WORKING MASSIVE FOR EVERY PARAMETER 
CCCCC THERE IS 4 OUTPUT PARAMETERS AT ALL (C,G,AG,BG)
  444 NC=NC+1 

CCCCC NPB - PROBABILITY LEVEL
      DO 1 I=1,10
      IF(NZ.GE.2)GOTO 31
      IF(NC.EQ.1.AND.MNX.EQ.1)CC(NPB,I)=CPn(NPB,I)     
      IF(NC.EQ.2.AND.MNX.EQ.1)CC(NPB,I)=GPn(NPB,I)
      IF(NC.EQ.3.AND.MNX.EQ.1)CC(NPB,I)=AGPn(NPB,I)      
      IF(NC.EQ.4.AND.MNX.EQ.1)CC(NPB,I)=BGPn(NPB,I)
      IF(NC.EQ.1.AND.MNX.EQ.2)CC(NPB,I)=CPX(NPB,I)     
      IF(NC.EQ.2.AND.MNX.EQ.2)CC(NPB,I)=GPX(NPB,I) 
      IF(NC.EQ.3.AND.MNX.EQ.2)CC(NPB,I)=AGPX(NPB,I)      
      IF(NC.EQ.4.AND.MNX.EQ.2)CC(NPB,I)=BGPX(NPB,I)
      GOTO 1
   31 CONTINUE
      IF(NC.EQ.1.AND.MNX.EQ.1)CC(NPB,I)=CHen(NPB,I)     
      IF(NC.EQ.2.AND.MNX.EQ.1)CC(NPB,I)=GHeN(NPB,I)
      IF(NC.EQ.3.AND.MNX.EQ.1)CC(NPB,I)=AGHeN(NPB,I)      
      IF(NC.EQ.4.AND.MNX.EQ.1)CC(NPB,I)=BGHeN(NPB,I)
      IF(NC.EQ.1.AND.MNX.EQ.2)CC(NPB,I)=CHeX(NPB,I)     
      IF(NC.EQ.2.AND.MNX.EQ.2)CC(NPB,I)=GHeX(NPB,I)
      IF(NC.EQ.3.AND.MNX.EQ.2)CC(NPB,I)=AGHeX(NPB,I)      
      IF(NC.EQ.4.AND.MNX.EQ.2)CC(NPB,I)=BGHeX(NPB,I)
    1 CONTINUE
CCCCC 
      CN=BN(NN)
      CNN=BN(NN+1)
      IF(CC(NPB,NN+1).LT.CC(NPB,NN))GOTO 111
      PP(NC,NPB)=CC(NPB,NN)+(CC(NPB,NN+1)-CC(NPB,NN))*
     *(ALOG10(AN)-ALOG10(CN))/(ALOG10(CNN)-ALOG10(CN))
      GOTO 222
  111 CONTINUE
      PP(NC,NPB)=CC(NPB,NN+1)+(CC(NPB,NN)-CC(NPB,NN+1))*(ALOG10(CNN)-
     *ALOG10(AN))/(ALOG10(CNN)-ALOG10(CN))
  222 CONTINUE

      IF(NC.LT.4)GOTO 444
           
      PP(1,NPB)=10.**PP(1,NPB)
     
CCCCC PARAMETERS FOR P AND He: FLUENCES - FN,  AND PEAK FLUXES - FX  
      IF(NZ.EQ.1.AND.MNX.EQ.1)OPEN(UNIT=22,FILE='PAR_FNNZ.DAT')
      IF(NZ.EQ.1.AND.MNX.EQ.2)OPEN(UNIT=22,FILE='PAR_FXNZ.DAT')
      IF(NZ.GE.2.AND.MNX.EQ.1)OPEN(UNIT=22,FILE='PAR_FNNZ.DAT')
      IF(NZ.GE.2.AND.MNX.EQ.2)OPEN(UNIT=22,FILE='PAR_FXNZ.DAT')

      WRITE(22,341)PROB(NPB),PP(1,NPB),PP(2,NPB),PP(3,NPB),PP(4,NPB)
      CLOSE(22)
      
      IF(NZ.GT.2)GOTO 998
      
      A1=PP(1,NPB)     
      A2=PP(2,NPB)
      A3=PP(3,NPB)          
      A4=PP(4,NPB)          
      CALL SPECTERC(A1,A2,A3,A4)      
CCCCC OUTPUT OF INTEGRAL AND DIFFERENTIAL ENERFY SPECTRA
      
      IF(NPB.EQ.1.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_9NZ.DAT')
      IF(NPB.EQ.2.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_7NZ.DAT')
      IF(NPB.EQ.3.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_5NZ.DAT')
      IF(NPB.EQ.4.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_3NZ.DAT')
      IF(NPB.EQ.5.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_1NZ.DAT')
      IF(NPB.EQ.6.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_032NZ.DAT')
      IF(NPB.EQ.7.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_01NZ.DAT')

      IF(NPB.EQ.1.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_9NZ.DAT')
      IF(NPB.EQ.2.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_7NZ.DAT')
      IF(NPB.EQ.3.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_5NZ.DAT')
      IF(NPB.EQ.4.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_3NZ.DAT')
      IF(NPB.EQ.5.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_1NZ.DAT')
      IF(NPB.EQ.6.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_032NZ.DAT')
      IF(NPB.EQ.7.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_01NZ.DAT')

	  WRITE(11,545)(EN(I),FDIF(I),FINT(I),I=1,NEN,5)
      CLOSE(11)
C  838 CONTINUE
CCCC  THERE WERE ENERGY SPECTRA - DIFF AND INTEGR: 
CCC   CYCLE BY PROBABILITY FOR P AND HE

  545 FORMAT(3E12.5)
  341 FORMAT(2X,7HPROBAB=,F7.4,2X,3H D=,E11.4,2X,3H G=,F7.4,
     *2X,4H AG=,F6.3,2X,4H BG=,F6.3)
      IF(NZ.LE.2) GOTO 999

  998 CONTINUE

CCCCCCCCCCC THE PARAMETERS FOR NZ>2 ARE CALCULATED
      IF(NZ.GE.30)NZ=30
      GOTO(999,999,997,997,997,996,996,996,997,996,
     *     996,996,996,996,996,996,996,996,996,996,
     *     997,996,996,996,997,996,997,996,997,997),NZ
  996 CONTINUE
CCCCCC  CORRECTIONS FOR THE SPECTRAL PARAMETERS OF HEAVY YONS

CCCCCC  Correction of spectral coefficient D 
      Z=NZ
         
CCCCC IF FLUENCE

      IF(MNX.EQ.2)GOTO 881
      IF(NZ.GE.6.AND.NZ.LE.8)S1=1./PRB(NPB)**(1.35E-06*Z**4.57)
      IF(NZ.GE.6.AND.NZ.LE.8)S2=9.39E-07*Z**4.6/PRB(NPB)**
     *(0.026+0.0121*Z)
      IF(NZ.GE.10.AND.NZ.LE.17)S1=1./PRB(NPB)**(0.00297*Z**1.18)
      IF(NZ.GE.10.AND.NZ.LE.17)S2=0.00244*Z**1.122/PRB(NPB)**(0.1584
     *-0.0013*Z)
      IF(NZ.GE.18.AND.NZ.LE.28)S1=1./PRB(NPB)**(0.04*Z**0.588)
      IF(NZ.GE.18.AND.NZ.LE.28)S2=0.02854*Z**0.572/PRB(NPB)**(0.1215
     *-0.00129*Z)
      GOTO 771
CCCCC IN CASE OF PEAK FLUX
  881 CONTINUE
      IF(NZ.LT.19)S1=0.986/PRB(NPB)**(0.00628*Z-0.0302)
      IF(NZ.GE.19)S1=0.928/PRB(NPB)**(0.0086*Z+0.0214)
      IF(NZ.LT.19)S2=(0.00652*Z-0.030)/PRB(NPB)**0.171
      IF(NZ.GE.10)S2=(0.0085*Z+0.0186)/PRB(NPB)**0.102
  771 PP(1,NPB)=PP(1,NPB)*RZ(NZ)/RZ(2)*S1*AN**S2
     
CCCCCC  Correction of spectral INDEX G - NO CORRECTION 

 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCCC  Correction of spectral INDEX AG
      IF(MNX.EQ.2)GOTO 882
CCCCC IN CASE OF  FLUENCE
      IF(NZ.LE.17)DA=0.
      IF(NZ.GE.18)DA=(-0.0045-0.0054*ALOG(PRB(NPB)))*ALOG(AN)
      GOTO 772
  882 CONTINUE
      IF(NZ.LE.17)DA=0.
      IF(NZ.GE.18)DA=(0.0034-(0.002+0.00013*Z)*ALOG(PRB(NPB)))*ALOG(AN)
  772 PP(3,NPB)=PP(3,NPB)+DA
      
CCCCCC  Correction of spectral coefficient BG
      IF(MNX.EQ.2)GOTO 883
CCCCC For fluences
      IF(NZ.GE.10.AND.NZ.LE.17)G1=1.2E-04*Z-3.10E-04
      IF(NZ.GE.10.AND.NZ.LE.17)G2=9.4E-05*Z-3.72E-04
      IF(NZ.GE.18)G1=6.9E-05*Z+0.00233
      IF(NZ.GE.18)G2=8.26E-05*Z+9.05E-04
      IF(NZ.GE.10)PP(4,NPB)=PP(4,NPB)*AN**(G1-G2*ALOG(PRB(NPB)))
      GOTO 773
CCCCC For peak Fluxes
  883 CONTINUE
      IF(NZ.GE.10.AND.NZ.LE.18)G1=1.67E-04*z-8.9E-04
      IF(NZ.GE.10.AND.NZ.LE.18)G2=0.465
      IF(NZ.GE.19)G1=1.12E-04*Z+3.24E-03
      IF(NZ.GE.19)G2=0.383
      IF(NZ.GE.10)PP(4,NPB)=PP(4,NPB)*AN**(G1/PRB(NPB)**G2)
  773 CONTINUE

      
      A1=PP(1,NPB)
      A2=PP(2,NPB)
      A3=PP(3,NPB)          
      A4=PP(4,NPB)          
      CALL SPECTERC(A1,A2,A3,A4)      
CCCCC OUTPUT OF INTEGRAL AND DIFFERENTIAL ENERFY SPECTRA

CCCCC IN FOLLOWING FILES NZ NEEDED TO REPLACE BY NUMBER (6...OR...10...OR 24)      

      IF(NPB.EQ.1.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_9NZ.DAT')
      IF(NPB.EQ.2.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_7NZ.DAT')
      IF(NPB.EQ.3.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_5NZ.DAT')
      IF(NPB.EQ.4.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_3NZ.DAT')
      IF(NPB.EQ.5.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_1NZ.DAT')
      IF(NPB.EQ.6.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_032NZ.DAT')
      IF(NPB.EQ.7.AND.MNX.EQ.1)OPEN(UNIT=11,FILE='FN0_01NZ.DAT')

      IF(NPB.EQ.1.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_9NZ.DAT')
      IF(NPB.EQ.2.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_7NZ.DAT')
      IF(NPB.EQ.3.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_5NZ.DAT')
      IF(NPB.EQ.4.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_3NZ.DAT')
      IF(NPB.EQ.5.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_1NZ.DAT')
      IF(NPB.EQ.6.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_032NZ.DAT')
      IF(NPB.EQ.7.AND.MNX.EQ.2)OPEN(UNIT=11,FILE='FX0_01NZ.DAT')


      WRITE(11,545)(EN(I),FDIF(I),FINT(I),I=1,NEN,5)
      CLOSE(11)
  738 CONTINUE      
      IF(MNX.EQ.1)OPEN(UNIT=41,FILE='PAR_FNNZ.DAT')
      IF(MNX.EQ.2)OPEN(UNIT=41,FILE='PAR_FXNZ.DAT')
      
      WRITE(41,341)PROB(NPB),PP(1,NPB),PP(2,NPB),PP(3,NPB),PP(4,NPB)     
      CLOSE(41)

      GOTO 999

  997 PRINT 991,NZ
  991 FORMAT(2X,25HTHE FLUXES OF ELEMENT NZ=,I3,15H ARE EQUAL TO 0)


  999 END

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC      
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC      
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCC DETERMINATION OF THE MODEL PARAMETER ANSEP
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      SUBROUTINE NSEPW(TY1,M1,TY2,M2,ANSEP)
      COMMON/PAR/NE,NPB,NEN,NF,NEI(13),DM(12)
      COMMON/SUNSP/WS(600),T(600),NW    
      DO 11 I=1,NW
      TT=I
   11 T(I)=1968.958333+0.083333*TT
      T1=TY1
      DO 1 I=1,M1
      IF(I.LT.M1)T1=T1+DM(I)/365.
      IF(I.EQ.M1)T1=T1+DM(I)/730.
    1 CONTINUE
      T2=TY2
      DO 2 I=1,M2
      IF(I.LT.M2)T2=T2+DM(I)/365.
      IF(I.EQ.M2)T2=T2+DM(I)/730.
    2 CONTINUE
      DO 3 I=1,NW
      IF(T(I).LT.T1.OR.T(I).GT.T2)GOTO 3
      WWM=WWM+WS(I)
    3 CONTINUE

      PRINT 541,T1,T2,WWM
  541 FORMAT(15H Time Interval=,F7.2,1H-,F7.2/
     *32H Sum of month sunspot numbers W=,F7.1)

CCCC  ANSEP - NO EVENTS F>10**5; BNSEP - NO EVENTS F>10**6 (>30 MeV)
      ANSEP=0.01354*WWM
      BNSEP=ANSEP/2.
      WWY=WWM/12.
      PRINT 542,ANSEP,BNSEP
  542 FORMAT(20H Model Parameter NN=F6.2/
     *22H Mean SEP Number NSEP=,F6.2)  
C      OPEN(UNIT=11,FILE='NSEP.DAT')
C      WRITE(11,4)ANSEP,BSEP,WWM
C    4 FORMAT(
     
      RETURN
      END
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      SUBROUTINE SPECTERC(CP,GAMMA,ALPHA,AK)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      COMMON/EN/EN(251),P(251),B(251),DE(251)
      COMMON/PFLUX/FDIF(251),FINT(251)
      COMMON/PAR/NE,NPB,NEN,NPN,NEI(13),DM(12)
      COMMON/E30/EN0,P0,EN1,P1
      SSS=CP/(GAMMA*AK-1.)/(P(NEN)/P1)**(GAMMA*AK-1.)*P1/B(NEN)
      FINT(NEN)=SSS
      FDIF(NEN)=CP/(P(NEN)/P1)**(GAMMA*AK)/B(NEN)
      DO 22 I=1,NEN-1
      J=NEN-I
      GAM=GAMMA*AK
      IF(EN(J).LT.EN1.AND.EN(J).GT.EN0)GAM=GAMMA
      IF(EN(J).LT.EN0)GAM=GAMMA*(EN(J)/EN0)**ALPHA
           
      FDIF(J)=CP/(P(J)/P1)**GAM/B(J)
      SSS=CP/(P(J)/P1)**GAM*DE(J)/B(J)+SSS
      FINT(J)=SSS
   22 CONTINUE
  333 CONTINUE
      RETURN
      END
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

      SUBROUTINE ARBITARY(ANSEP)
      COMMON/FLUX/MNX
      COMMON/EN/EN(251),P(251),B(251),DE(251)
      COMMON/PROB/PB(801),FXPRB(13,801),FPB(13,7),FF(801)
      COMMON/PAR/NE,NPB,NEN,NF,NEI(13),DM(12)
      COMMON/E30/EN0,P0,EN1,P1
      COMMON/ZRZ/NZ,RZ(28),AQ(28)

      DO 1 I=1,NEN
      T=I
      EN(I)=3.16228*10.**((T-1.)/100.)
      DE(I)=3.16228*(10.**((T-0.5)/100.)-10.**((T-1.5)/100.))
      P(I)=SQRT(EN(I)*(EN(I)+2.*938.))
      B(I)=P(I)/SQRT(P(I)**2+938.**2)
    1 CONTINUE

      EN0=30.0
      P0=SQRT(EN0*(EN0+2.*938.))
      EN1=79.4
      P1=SQRT(EN1*(EN1+2.*938.))
C      PRINT 456, EN0,P0,EN1,P1
C  456 FORMAT(4F10.3)
C      PAUSE
      IF(NZ.EQ.1)CZ=1.
      IF(NZ.EQ.2)CZ=0.01
      IF(NZ.GT.2)CZ=0.00001
      IF(MNX.EQ.1)CNM=1.
      IF(MNX.EQ.2)CNM=1.0E-06
      IF(ANSEP.LE.2.)CC=0.01*CNM*CZ
      IF(ANSEP.GT.2.AND.ANSEP.LE.8.)CC=1.*CNM*CZ
      IF(ANSEP.GT.8.AND.ANSEP.LE.32.)CC=10.*CNM*CZ
      IF(ANSEP.GT.32.AND.ANSEP.LT.128.)CC=100.*CNM*CZ
      IF(ANSEP.GE.128.)CC=1000.*CNM*CZ

      DO 2 I=1,NF
      T=I-1
      PB(I)=CC*10.**(T*.02)
      FF(I)=1.0E-04*10.**(T*0.02)*CNM
    2 CONTINUE
      RETURN
      END
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      BLOCK DATA
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCC DATABASE FOR P AND He FLUENCES AND FLUXES
      COMMON/PARflnP/Cpn(7,10),Gpn(7,10),AGpn(7,10),BGpn(7,10)
      COMMON/PARflnHe/CHen(7,10),GHen(7,10),AGHen(7,10),BGHen(7,10)                                                                                                                                            
      COMMON/PARflxP/Cpx(7,10),Gpx(7,10),AGpx(7,10),BGpx(7,10)
      COMMON/PARflxHe/CHex(7,10),GHex(7,10),AGHex(7,10),BGHex(7,10)                                                                                                                                                                                                             COMMON/EN/EN(251),P(251),B(251),DE(251)

      COMMON/SUNSP/W(600),T(600),NW                                                                                                                                                                                                                                                                                                                                                                                                              
      DIMENSION W69(12),W70(120),W80(120),W90(120),
     *W00(120),W10(108)
      EQUIVALENCE(W(1),W69(1))
      EQUIVALENCE(W(13),W70(1))
      EQUIVALENCE(W(133),W80(1))
      EQUIVALENCE(W(253),W90(1))
      EQUIVALENCE(W(373),W00(1))
      EQUIVALENCE(W(493),W10(1))     
      COMMON/DISTF/GF,FPD,FN0,FX0
      COMMON/DISTG/G0,DG
      COMMON/DISTA/CA,GAF,GAG,DA
      COMMON/ZZZ/PROB(7),ANSEP
      COMMON/ENI/ENI(13)
      COMMON/PAR/NE,NPB,NEN,NF,NEI(13),DM(12)
      COMMON/ZRZ/NZ,RZ(28),AQ(28)
CCCRZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ RZ
      DATA RZ/1.,2.8E-03,0.,0.,0.,3.18E-05,7.56E-06,6.29E-05,0.,9.8E-06,
     *6.62E-07,1.22E-05,8.7E-07,8.63E-06,4.0E-08,1.65E-06,1.25E-08,
     *1.6E-07,2.6E-08,4.6E-07,0.,1.6E-08,0.,8.6E-07,0.,4.59E-06,0.,
     *2.48E-07/
CCCAQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ AQ
      DATA AQ/0.,2.,2.,2.,2.,2.03,2.12,2.13,0.,2.30,2.44,2.38,2.55,2.39,
     *2.44,2.48,2.60,2.92,2.96,2.99,3.0,3.0,3.0,3.0,3.0,3.035,3.08,3.12/
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      DATA GF,FPD,FN0,FX0/1.32,0.12,8.0E09,9.3E03/
CCC   Distribution of the spectral indexes
      DATA G0,DG/0.76,0.15/
      DATA CA,GAF,GAG,DA/1.16,0.059,0.143,0.0777/
CCC   ENERGY THRESOLDS
C     NE - NUMBER OF THE ENERGY THRESHOLDS
C     NE=NF - NUMBER OF ENERGY SCALE INTERVALS
C     NPB - NUMBER OF CONFIDENCE LEVELS
C     NF - NUMBER OF THE FLUENCE INTERVALS

      DATA NE,NPB,NEN,NF,NW/13,7,251,801,600/
C     PROB - CONFIDENCE LEVELS (0.001=99.9%, 0.1=90%, 0.9=10%)
      DATA PROB/0.9,0.70,0.5,0.25,0.1,0.0316,0.01/
CCCCC DAYS OF MONTH DM(12)
      DATA DM/31.,28.,31.,30.,31.,30.,31.,31.,30.,31.,30.,31./
      DATA W69/
     *110.0,109.6,108.0,106.4,106.2,106.1,105.9,106.5,105.4,104.1,104.6, 1969
     *104.9/
      DATA W70/
     *105.6,106.0,106.2,106.1,105.8,105.3,103.8,101.0,97.2,93.9,89.4,    1970
     *84.1,
     *80.4,77.8,74.4,70.9,68.1,66.7,65.4,64.6,65.8,66.2,66.8,69.4,       1971
     *70.8,71.2,72.4,73.4,72.9,70.5,68.2,65.5,62.2,60.6,58.7,55.1,       1972
     *50.9,46.5,44.2,42.7,40.7,39.1,37.5,36.1,34.4,32.6,31.8,31.5,       1973
     *32.7,34.4,34.0,33.9,34.6,34.5,34.0,33.1,32.1,30.2,27.5,25.2,       1974
     *23.0,22.1,21.2,18.6,16.8,16.0,15.0,14.3,14.4,15.4,16.1,16.3,       1975
     *15.2,13.2,12.2,12.6,12.5,12.2,12.9,14.0,14.3,13.5,13.5,14.8,       1976
     *16.7,18.1,20.0,22.2,24.2,26.3,29.0,33.4,39.1,45.6,51.9,56.9,       1977
     *61.3,64.5,69.6,76.9,83.2,89.3,97.4,104.0,108.4,111.1,113.3,117.7,  1978
     *123.7,130.9,136.5,141.1,147.3,153.0,155.0,155.4,155.7,157.8,162.3, 1979
     *164.5/
      DATA W80/
     *163.9,162.6,160.9,158.7,156.3,154.7,152.8,150.3,150.1,150.2,147.7, 1980
     *142.7,
     *140.3,141.5,143.0,143.4,142.9,141.5,140.3,141.1,142.8,142.2,138.8, 1981
     *137.8,
     *137.0,133.3,129.2,124.3,119.9,117.3,115.2,109.4,101.0,95.7,94.7,   1982
     *94.6,
     *92.8,90.3,85.9,81.5,77.1,70.5,65.5,65.8,67.9,68.2,66.8,64.0,       1983
     *60.2,56.4,53.0,49.8,47.5,46.5,44.1,39.6,33.9,29.0,24.8,21.8,       1984
     *20.6,19.8,18.7,18.4,18.4,18.2,17.5,17.3,17.5,17.4,16.9,15.3,       1985
     *13.9,13.2,13.0,13.7,14.4,13.9,13.8,13.2,12.3,13.2,14.9,16.3,       1986
     *17.6,19.6,22.1,24.4,26.5,29.6,33.7,37.3,41.4,46.0,49.1,53.8,       1987
     *60.2,67.0,73.7,80.0,86.1,94.8,104.3,113.7,121.3,125.3,130.4,137.6, 1988
     *142.0,145.0,149.7,153.5,156.9,158.4,158.5,157.7,156.6,157.4,157.5, 1989
     *153.5/
      DATA W90/
     *150.6,152.9,152.0,149.3,147.0,143.8,140.6,140.5,142.1,142.1,141.7, 1990
     *143.9,
     *147.6,147.6,146.6,146.5,145.5,145.2,146.3,146.6,144.9,141.7,138.1, 1991
     *131.7,
     *123.7,115.4,108.2,103.3,100.3,97.1,90.7,84.0,79.5,76.4,74.4,73.2,  1992
     *71.4,69.3,66.6,63.6,59.9,56.1,54.7,52.3,48.4,44.9,41.2,38.4,       1993
     *36.6,34.8,34.1,33.7,32.5,30.8,28.5,26.8,26.6,26.5,26.2,25.9,       1994
     *24.2,23.0,22.1,20.6,19.2,18.2,17.0,15.4,13.4,12.1,11.3,10.8,       1995
     *10.4,10.1, 9.7, 8.4, 8.0, 8.5, 8.4, 8.3, 8.4, 8.8, 9.8,10.4,       1996
     *10.5,11.0,13.5,16.5,18.4,20.4,22.7,25.1,28.4,31.9,35.0,39.0,       1997
     *43.8,49.0,53.5,56.6,59.4,62.5,65.5,67.8,69.5,70.5,73.0,77.9,       1998
     *82.6,84.6,83.8,85.5,90.5,93.1,94.3,97.5,102.3,107.8,111.0,         1999
     *111.1/

      DATA W00/
     *112.9,116.8,119.9,120.8,119.0,118.7,119.8,118.6,116.3,114.5,112.7, 2000
     *112.0,
     *108.7,104.0,104.8,107.5,108.6,109.8,111.7,113.6,114.1,114.0,115.5, 2001
     *114.6,
     *113.5,114.6,113.3,110.5,108.8,106.2,102.7,98.7,94.6,90.5,85.2,     2002
     *82.0,
     *80.8,78.5,74.0,70.1,67.6,65.0,61.8,60.1,59.6,58.2,56.8,54.8,       2003
     *52.0,49.4,47.2,45.6,43.9,41.7,40.2,39.3,37.6,35.9,35.4,35.3,       2004
     *34.7,34.0,33.6,31.7,29.0,28.9,29.2,27.5,25.9,25.5,24.9,23.0,       2005
     *20.8,18.7,17.4,17.1,17.3,16.3,15.3,15.6,15.6,14.2,12.7,12.1,       2006
cccccccccccccc  w at NOAA  cccccccccccccccccccccccccccccccccccccccccccc
     *12.4,12.3,12.4,12.3,12.1,11.8,11.7,11.5,11.4,11.3,11.,5.,          2007                      
     *10.,11.3,12.5,5.,5.,5.5,6.,6.5,7.,9.5,11.5,13.5,                   2008
     *15.,17.,20.5,23.,27.,31.,36.,40.,43.,47.3,51.8,55./                2009
      DATA W10/                                                          
     *59.3,63.3,67.2,71.,74.8,78.4,81.9,85.2,88.4,91.4,94.2,96.9,        2010
     *99.4,101.6,103.6,105.5,107.1,108.6,109.8,110.8,111.7,112.3,112.8,  2011
     *113.0,                                                                 
     *113.1,112.9,112.6,112.,111.4,110.6,109.6,108.5,107.2,105.8,        2012
     *104.4,102.7,
     *101.0,98.2,97.2,95.2,93.1,91.,88.8,86.6,84.3,82.,79.7,77.3,        2013
     *74.9,72.6,70.2,67.9,65.4,63.2,60.9,58.6,56.4,54.1,52.,49.8,        2014
     *47.7,45.6,43.6,41.7,39.9,38.,36.,34.5,32.8,31.2,29.7,28.1,         2015
     *30.,28.8,27.5,26.3,25.,24.4,23.8,23.1,22.5,21.4,20.3,19.1,         2016
     *18.,17.6,17.3,16.9,16.5,15.6,14.8,13.9,13.,12.6,12.3,11.9,         2017
     *11.5,11.6,11.8,11.9,12.,13.1,14.3,15.4,16.5,17.6,18.75,19.9/       2018

      DATA Cpn/2.646, 2.646, 2.646, 3.631, 4.709, 5.683, 6.304,
     *2.860, 2.860, 3.471, 4.417, 5.335, 6.095, 6.565,
     *3.001, 3.707, 4.283, 5.106, 5.830, 6.411, 6.772,
     *3.860, 4.517, 5.021, 5.681, 6.232, 6.663, 6.954,
     *4.694, 5.240, 5.639, 6.141, 6.551, 6.877, 7.106,
     *5.414, 5.847, 6.149, 6.519, 6.817, 7.064, 7.251,
     *6.042, 6.361, 6.579, 6.845, 7.059, 7.252, 7.396,
     *6.573, 6.800, 6.952, 7.140, 7.299, 7.444, 7.554,
     *7.021, 7.182, 7.291, 7.419, 7.536, 7.639, 7.726,
     *7.423, 7.530, 7.608, 7.699, 7.780, 7.855, 7.914/

      DATA Gpn/5.754,  5.754,  5.754,  6.327,  6.448,  6.580,  6.649,
     *5.892,  5.892,  6.255,  6.414,  6.525,  6.642,  6.626,
     *5.935,  6.291,  6.385,  6.481,  6.616,  6.605,  6.602,
     *6.266,  6.384,  6.436,  6.556,  6.611,  6.620,  6.547,
     *6.356,  6.437,  6.508,  6.594,  6.607,  6.608,  6.524,
     *6.423,  6.492,  6.562,  6.609,  6.627,  6.618,  6.501,
     *6.494,  6.542,  6.577,  6.607,  6.627,  6.587,  6.480,
     *6.542,  6.566,  6.585,  6.607,  6.610,  6.564,  6.484,
     *6.546,  6.582,  6.596,  6.629,  6.622,  6.591,  6.506,
     *6.573,  6.609,  6.610,  6.633,  6.622,  6.604,  6.567/

      DATA AGpn/.0689,.0689,.0689,.1085,.1871,0.2721,.3159,
     *.0721,   .0721,   .0955,   .1573,   .2377,   .3030,   .3169,
     *.0752,   .1038,   .1399,   .2076,   .2780,   .3133,   .3175,
     *.1027,   .1446,   .1858,   .2521,   .2983,   .3161,   .3009,
     *.1413,   .1885,   .2297,   .2796,   .3066,   .3113,   .2900,
     *.1825,   .2255,   .2590,   .2924,   .3054,   .3004,   .2732,
     *.2215,   .2543,   .2756,   .2957,   .3009,   .2880,   .2564,
     *.2491,   .2708,   .2826,   .2933,   .2925,   .2780,   .2476,
     *.2623,   .2780,   .2855,   .2908,   .2875,   .2705,   .2422,
     *.2737,   .2829,   .2860,   .2874,   .2821,   .2666,   .2418/

      DATA BGpn/1.075,  1.075,  1.075,   .803,   .772,   .722,   .665,
     *.987,   .987,   .826,   .781,   .747,   .687,   .635,
     *.974,   .817,   .788,   .760,   .711,   .662,   .611,
     *.826,   .789,   .772,   .732,   .684,   .629,   .591,
     *.793,   .769,   .746,   .702,   .656,   .605,   .573,
     *.772,   .747,   .718,   .676,   .630,   .586,   .560,
     *.748,   .720,   .693,   .653,   .609,   .575,   .551,
     *.722,   .695,   .670,   .634,   .600,   .569,   .546,
     *.701,   .673,   .653,   .619,   .591,   .563,   .548,
     *.680,   .655,   .639,   .612,   .589,   .566,   .549/




      DATA CHen/
     *-.119,   -.119,   -.119,    .934,   1.980,   2.975,   3.668,
     * .134,    .134,    .780,   1.696,   2.615,   3.440,   3.987,
     * .266,   1.013,   1.571,   2.390,   3.155,   3.810,   4.268,
     *1.162,   1.801,   2.295,   2.982,   3.590,   4.115,   4.495,
     *1.976,   2.517,   2.931,   3.476,   3.957,   4.378,   4.697,
     *2.701,   3.146,   3.472,   3.895,   4.267,   4.609,   4.886,
     *3.335,   3.683,   3.929,   4.255,   4.549,   4.827,   5.058,
     *3.885,   4.148,   4.334,   4.579,   4.807,   5.030,   5.216,
     *4.366,   4.559,   4.697,   4.881,   5.057,   5.231,   5.379,
     *4.789,   4.934,   5.039,   5.171,   5.303,   5.436,   5.553/


      DATA GHen/
     *6.827,  6.827,  6.827,  7.353,  7.621,  7.582,  7.489,
     *6.873,  6.873,  7.260,  7.564,  7.655,  7.543,  7.472,
     *6.955,  7.313,  7.513,  7.630,  7.579,  7.462,  7.401,
     *7.311,  7.543,  7.631,  7.621,  7.529,  7.427,  7.371,
     *7.507,  7.621,  7.631,  7.575,  7.464,  7.407,  7.376,
     *7.602,  7.620,  7.587,  7.508,  7.432,  7.375,  7.362,
     *7.621,  7.581,  7.551,  7.473,  7.399,  7.388,  7.374,
     *7.601,  7.556,  7.510,  7.465,  7.421,  7.402,  7.446,
     *7.556,  7.524,  7.501,  7.470,  7.424,  7.443,  7.508,
     *7.550,  7.507,  7.482,  7.477,  7.442,  7.451,  7.503/


      DATA AGHen/
     *.0922,   .0922,   .0922,   .0735,   .0939,   .1027,   .0917,
     *.0808,   .0808,   .0739,   .0856,   .1038,   .1013,   .0814,
     *.0825,   .0751,   .0827,   .0984,   .1047,   .0908,   .0663,
     *.0761,   .0859,   .0951,   .1039,   .0986,   .0778,   .0487,
     *.0838,   .0959,   .1019,   .1020,   .0876,   .0623,   .0309,
     *.0938,   .1016,   .1024,   .0945,   .0756,   .0451,   .0144,
     *.1004,   .1006,   .0977,   .0849,   .0617,   .0328,  -.0005,
     *.1017,   .0977,   .0902,   .0755,   .0518,   .0192,  -.0112,
     *.0976,   .0903,   .0827,   .0661,   .0408,   .0107,  -.0198,
     *.0936,   .0838,   .0749,   .0587,   .0327,  -.0002,  -.0323/

      DATA BGHen/
     *1.094,  1.094,  1.094,   .750,   .672,   .636,   .588,
     *1.004,  1.004,   .789,   .687,   .648,   .609,   .557,
     *.977,   .772,   .704,   .660,  .629,     .582,   .534,
     *.777,   .698,   .668,   .639,  .602,    .555,   .511,
     *.705,   .667,   .648,   .616,   .577,   .529,   .482,
     *.671,   .646,   .627,   .594,   .552,   .505,   .464,
     *.648,   .626,   .604,   .570,   .530,   .483,   .444,
     *.627,   .604,   .583,   .548,   .509,   .465,   .424,
     *.607,   .584,   .562,   .529,   .493,   .450,   .409,
     *.586,   .565,   .546,   .513,   .481,   .441,   .402/

      Data Cpx/
     *-3.026, -3.026, -3.026, -2.056,  -.948,   .037,   .667,
     *-2.875, -2.875, -2.268, -1.273,  -.331,   .447,   .917,
     *-2.829, -2.105, -1.478,  -.595,   .156,   .746,  1.116,
     *-2.071, -1.325,  -.762,  -.048,   .536,   .986,  1.293,
     *-1.285, -.636,  -.185,   .382,   .825,  1.186,  1.432,
     *-.609,  -.088,   .275,   .711,  1.042,  1.340,  1.554,
     *-.050,   .352,   .626,   .958,  1.225,  1.472,  1.648,
     *.380,    .684,   .892,  1.153,  1.379,  1.581,  1.736,
     *.709,    .941,  1.104,  1.313,  1.504,  1.678,  1.816,
     *.954,   1.141,  1.272,  1.451,  1.613,  1.771,  1.885/



      Data Gpx/
     *5.726,  5.726,  5.726,  6.375,  6.474,  6.599,  6.608,
     *5.885,  5.885,  6.266,  6.407,  6.535,  6.597,  6.623,
     *5.966,  6.312,  6.408,  6.491,  6.617,  6.651,  6.585,
     *6.311,  6.426,  6.478,  6.595,  6.633,  6.629,  6.439,
     *6.414,  6.457,  6.573,  6.615,  6.616,  6.547,  6.473,
     *6.519,  6.629,  6.622,  6.598,  6.633,  6.466,  6.336,
     *6.571,  6.629,  6.640,  6.590,  6.596,  6.456,  6.358,
     *6.631,  6.656,  6.644,  6.570,  6.492,  6.394,  6.260,
     *6.619,  6.619,  6.587,  6.536,  6.439,  6.322,  6.156,
     *6.611,  6.582,  6.547,  6.485,  6.383,  6.190,  6.077/

      Data Agpx/
     *.0629,   .0629,   .0629,   .1118,   .1933,   .2773,   .3130,
     *.0638,   .0638,   .0948,   .1626,   .2440,   .3032,   .3210,
     *.0680,   .1045,   .1479,   .2195,   .2851,   .3214,   .3149,
     *.1062,   .1584,   .2042,   .2701,   .3121,   .3248,   .2946,
     *.1620,   .2122,   .2581,   .3014,   .3189,   .3107,   .2891,
     *.2186,   .2690,   .2970,   .3192,   .3224,   .2922,   .2541,
     *.2644,   .3004,   .3151,   .3179,   .3113,   .2789,   .2391,
     *.3003,   .3197,   .3238,   .3123,   .2924,   .2557,   .2073,
     *.3177,   .3210,   .3158,   .3000,   .2717,   .2282,   .1755,
     *.3185,   .3165,   .3026,   .2854,   .2521,   .1953,   .1420/

      Data Bgpx/
     *1.074,   1.074,  1.074,   .790,   .767,   .716,   .667,
     *.974,   .974,   .814,   .780,   .743,   .691,   .634,
     *.941,   .802,   .781,   .758,   .707,   .651,   .607,
     *.800,   .778,   .765,   .722,   .675,   .622,   .600,
     *.777,   .762,   .732,   .693,   .648,   .606,   .572,
     *.751,   .719,   .700,   .664,   .613,   .585,   .559,
     *.725,   .693,   .668,   .633,   .590,   .564,   .543,
     *.692,   .660,   .636,   .606,   .578,   .551,   .539,
     *.662,   .632,   .612,   .583,   .561,   .542,   .535,
     *.631,   .607,   .589,   .567,   .548,   .542,   .531/

      Data CHex/
     *-5.783, -5.783, -5.783, -4.743, -3.674, -2.665, -1.964,
     *-5.577, -5.577, -4.939, -3.997, -3.043, -2.204, -1.646,
     *-5.522, -4.786, -4.190, -3.318, -2.528, -1.840, -1.387,
     *-4.751, -4.054, -3.494, -2.748, -2.098, -1.553, -1.166,
     *-3.997, -3.354, -2.877, -2.272, -1.757, -1.308,  -.964,
     *-3.311, -2.777, -2.388, -1.897, -1.475, -1.088,  -.795,
     *-2.735, -2.302, -1.992, -1.592, -1.243,  -.905,  -.647,
     *-2.275, -1.922, -1.671, -1.341, -1.037,  -.748,  -.511,
     *-1.900, -1.611, -1.403, -1.122,  -.858,  -.601,  -.398,
     *-1.592, -1.354, -1.173,  -.931,  -.695,  -.470,  -.285/

      Data GHex/
     *6.820,  6.820,  6.820,  7.357,  7.599,  7.586,  7.466,
     *6.828,  6.828,  7.230,  7.600,  7.652,  7.502,  7.436,
     *6.893,  7.325,  7.537,  7.640,  7.639,  7.454,  7.422,
     *7.316,  7.608,  7.651,  7.659,  7.518,  7.381,  7.392,
     *7.595,  7.638,  7.600,  7.542,  7.441,  7.402,  7.369,
     *7.644,  7.643,  7.584,  7.466,  7.393,  7.347,  7.365,
     *7.608,  7.563,  7.466,  7.386,  7.365,  7.327,  7.391,
     *7.547,  7.465,  7.413,  7.369,  7.343,  7.383,  7.446,
     *7.486,  7.401,  7.372,  7.323,  7.335,  7.436,  7.541,
     *7.379,  7.347,  7.318,  7.314,  7.367,  7.476,  7.611/

      Data AgHex/
     *.0822,   .0822,   .0882,   .0727,   .0909,   .1029,   .0897,
     *.0726,   .0726,   .0681,   .0879,   .1048,   .0989,   .0804,
     *.0734,   .0713,   .0822,   .0988,   .1085,   .0896,   .0654,
     *.0712,   .0872,   .0969,   .1074,   .0968,   .0706,   .0454,
     *.0876,   .0983,   .1040,   .1020,   .0844,   .0584,   .0288,
     *.1000,   .1062,   .1052,   .0908,   .0685,   .0378,   .0093,
     *.1058,   .1028,   .0924,   .0746,   .0517,   .0181,  -.0074,
     *.1017,   .0914,   .0796,   .0595,   .0325,   .0037,  -.0212,
     *.0921,   .0769,   .0640,   .0391,   .0152,  -.0096,  -.0339,
     *.0746,   .0589,   .0448,   .0219,  -.0003,  -.0247,  -.0458/

      Data BgHex/
     *1.091,  1.091,  1.091,  .745,   .673,   .637,   .588,
     *1.002,  1.002,   .785,  .680,   .647,   .612,   .558,
     *.966,   .755,   .694,   .657,   .620,   .582,   .529,
     *.752,   .679,   .660,   .632,   .601,   .555,   .506,
     *.679,   .658,   .645,   .613,   .571,   .522,   .481,
     *.656,   .634,   .617,   .586,   .544,   .498,   .458,
     *.636,   .613,   .595,   .561,   .517,   .476,   .434,
     *.612,   .588,   .566,   .530,   .491,   .448,   .412,
     *.584,   .560,   .539,   .505,   .468,   .425,   .387,
     *.559,   .533,   .513,   .481,   .445,   .405,   .370/
                                                  


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      END
