aboutsummaryrefslogtreecommitdiff
path: root/orrs/src/ncorr.f
blob: e3a7b3233a4bbb364e2abbd8f28f4a94c15819e4 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
      PROGRAM NCORR
      PARAMETER (NH=13)
      REAL H(NH), RT(NH),DN(NH)
C
      DATA H /  2.3,    2.4,    2.5,    2.6,    2.8,    3.0,    3.4,
     &          4.0,    5.0,    7.0,   10.0,   15.0,   20.0 /
C
      DATA RT / 4000., 1820.,   700.,   270.,   100.,   76.,    52., 
     &          34.,    26.,     14.,    9.0,   6.4,    5.0 /
C
      DATA DN / 0.0067, 0.0064, 0.0076, 0.0104, 0.0200, 0.0323, 0.0529,
     &          0.0727, 0.110 , 0.169 , 0.253 , 0.390 , 0.526  /
C
      REAL X1(NH), Y1(NH), F1(NH), G1(NH)
      REAL X2(NH), Y2(NH), F2(NH), G2(NH)
C
      IDEV = 6
      SIZE = 8.0
      CH = 0.02
C
      HMAX = 20.0
      DH = 2.0
C
      FMAX = 0.5
      DF = 0.1
C
      FMAX = 0.05
      DF = 0.005
C
      PAR = 0.75
C
      HWT = 1.0/HMAX
      FWT = PAR/FMAX
C     
      N = NH
C
      DO 20 I=1, N
        HB = 1.0/(H(I)-1.0)
C
        X1(I) = 4.0*H(I)
        Y1(I) = DN(I)
        G1(I) = 0.13 * (0.215/HB)
        F1(I) = 0.13 * (0.215/HB)
     &        - 0.0345 * EXP(-15.0*(HB-0.65)**2)
        F1(I) = 0.028*(H(I)-1.0)
     &        - 0.0345 * EXP(-(3.87/(H(I)-1.0) - 2.52)**2)
C
        X2(I) = HMAX* HB
        Y2(I) = DN(I)
        G2(I) = 0.13 * (0.215/HB)
        F2(I) = 0.13 * (0.215/HB)
     &        - 0.0345 * EXP(-15.0*(HB-0.65)**2)
        F2(I) = 0.028*(H(I)-1.0)
     &        - 0.0345 * EXP(-(3.87/(H(I)-1.0) - 2.52)**2)
C
 20   CONTINUE
C
C
      CALL PLOTS(0,-999,IDEV)
      CALL FACTOR(SIZE)
C
      CALL PLOT(8.0*CH,8.0*CH,-3)
C
      CALL PLOTON
C
      CALL XAXIS(0.0,-PAR,1.0,0.1,0.0,0.1,CH,1)
C
      CALL XAXIS(0.0,0.0,1.0,DH*HWT,0.0,DH,CH,1)
      CALL YAXIS(0.0,0.0,PAR,DF*FWT,0.0,DF,CH,3)
C
      CALL XYPLOT(N,X1,Y1,0.0,HWT,0.0,FWT,1,0.3*CH,+1)
      CALL XYPLOT(N,X1,F1,0.0,HWT,0.0,FWT,1,0.3*CH, 0)
      CALL XYPLOT(N,X1,G1,0.0,HWT,0.0,FWT,5,0.3*CH, 0)
C
      CALL XYPLOT(N,X2,Y2,0.0,HWT,0.0,FWT,3,0.3*CH,+5)
      CALL XYPLOT(N,X2,F2,0.0,HWT,0.0,FWT,2,0.3*CH, 0)
      CALL XYPLOT(N,X2,G2,0.0,HWT,0.0,FWT,3,0.3*CH, 0)
C
      CALL PLOTOF
C
      WRITE(*,*) 'Hit <cr>'
      READ(*,1000) ANS
 1000 FORMAT(A4)
C
      CALL PLOT(0.0,0.0,+999)
      STOP
      END