PHOS/shaker/shrndv.f

   1 *CMZ :          17/07/98  15.49.05  by  Federico Carminati
   2 *-- Author :
   3       SUBROUTINE SHRNDV
   4 c       =================
   5
   6 c       Prepare random generation vectors
   7
   8 *KEEP,SHPHYP.
   9       COMMON /SHPHYP/ JWEI,NDNDY,YLIM,PTLIM,JWEAK,JPI0,JETA,JPIC,JPRO,
  10      +                  JKAC,JKA0,JRHO,JOME,JPHI,JPSI,JDRY
  11 *KEEP,SHGENE.
  12       COMMON /SHGENE/ IEVT,NPI0,NETA,NPIC,NPRO,NKAC,NKA0,NRHO,NOME,
  13      +                  NPHI,NPSI,NDRY
  14 *KEEP,SHRAND.
  15       COMMON /SHRAND/ PISP1(100),PISP2(100),ETASP1(100),ETASP2(100),
  16      +                  PROSP(100),KAOSP(100)
  17 *KEEP,SHNORM.
  18       COMMON /SHNORM/ PINOR,PIRAT,ETANOR,ETARAT,RHONOR,OMENOR,PHINOR,
  19      +                  PSINOR,DRYNOR
  20 *KEEP,FUNINT.
  21       COMMON /FUNINT/ FINT
  22 *KEND.
  23
  24       EXTERNAL SHFPI,SHFETA,SHFPRO,SHFKAO,SHFRHO,SHFOME,SHFPHI
  25       EXTERNAL SHFPSI,SHFDRY
  26
  27       IF (JPI0.EQ.1.OR.JPIC.EQ.1) THEN
  28         CALL FUNPRE(SHFPI,PISP1,0.,2.)
  29         PINO1 = FINT
  30         CALL FUNPRE(SHFPI,PISP2,2.,PTLIM)
  31         PINO2 = FINT
  32         PINOR=PINO1+PINO2
  33         PIRAT=PINO1/PINOR
  34       ENDIF
  35
  36       IF (JETA.EQ.1) THEN
  37         CALL FUNPRE(SHFETA,ETASP1,0.,2.)
  38         ETANO1 = FINT
  39         CALL FUNPRE(SHFETA,ETASP2,2.,PTLIM)
  40         ETANO2 = FINT
  41         ETANOR=ETANO1+ETANO2
  42         ETARAT=ETANO1/ETANOR
  43       ENDIF
  44
  45       IF (JPRO.EQ.1) THEN
  46         CALL FUNPRE(SHFPRO,PROSP,0.,PTLIM)
  47       ENDIF
  48
  49       IF (JKAC.EQ.1.OR.JKA0.EQ.1) THEN
  50         CALL FUNPRE(SHFKAO,KAOSP,0.,PTLIM)
  51       ENDIF
  52
  53
  54       IF (JRHO.EQ.1) RHONOR = GAUSS(SHFRHO,0.,PTLIM,4.E-6)
  55       IF (JOME.EQ.1) OMENOR = GAUSS(SHFOME,0.,PTLIM,4.E-6)
  56       IF (JPHI.EQ.1) PHINOR = GAUSS(SHFPHI,0.,PTLIM,2.E-6)
  57       IF (JPSI.EQ.1) PSINOR = GAUSS(SHFPSI,0.,PTLIM,2.E-8)
  58       IF (JDRY.EQ.1) DRYNOR = GAUSS(SHFDRY,0.,PTLIM,1.E-7)
  59
  60       RETURN
  61       END