~maddevelopers/mg5amcnlo/dim6_eft

« back to all changes in this revision

Viewing changes to Template/LO/Source/PDF/NNPDFDriver.f

  • Committer: olivier Mattelaer
  • Date: 2014-09-24 20:33:52 UTC
  • mfrom: (253.17.56 2.2.0)
  • Revision ID: olivier.mattelaer@uclouvain.be-20140924203352-bplvftjd8czlhagt
pass to 2.2.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
21
21
      integer nfl,nx,nq2,mem,rep
22
22
      double precision alphas
23
23
      double precision xgrid(100),logxgrid(100)
24
 
      double precision q2grid(50),logq2grid(50)
25
 
      double precision pdfgrid(0:100,14,100,50)
 
24
      double precision q2grid(60),logq2grid(60)
 
25
      double precision pdfgrid(0:100,14,100,60)
26
26
      logical hasphoton
27
27
      common /nnpdf/nfl,nx,nq2,mem,rep,hasphoton,alphas,xgrid,logxgrid,
28
28
     1     q2grid,logq2grid,pdfgrid
31
31
*
32
32
      nfl = 13
33
33
      nx = 100
34
 
      nq2 = 50
 
34
      nq2 = 60
35
35
      mem = 1
36
36
      rep = 0
37
37
      alphas = 0
54
54
      integer nfl,nx,nq2,mem,rep
55
55
      double precision alphas
56
56
      double precision xgrid(100),logxgrid(100)
57
 
      double precision q2grid(50),logq2grid(50)
58
 
      double precision pdfgrid(0:100,14,100,50)
 
57
      double precision q2grid(60),logq2grid(60)
 
58
      double precision pdfgrid(0:100,14,100,60)
59
59
      logical hasphoton
60
60
      common /nnpdf/nfl,nx,nq2,mem,rep,hasphoton,alphas,xgrid,logxgrid,
61
61
     1     q2grid,logq2grid,pdfgrid
78
78
      integer nfl,nx,nq2,mem,rep
79
79
      double precision alphas
80
80
      double precision xgrid(100),logxgrid(100)
81
 
      double precision q2grid(50),logq2grid(50)
82
 
      double precision pdfgrid(0:100,14,100,50)
 
81
      double precision q2grid(60),logq2grid(60)
 
82
      double precision pdfgrid(0:100,14,100,60)
83
83
      logical hasphoton
84
84
      common /nnpdf/nfl,nx,nq2,mem,rep,hasphoton,alphas,xgrid,logxgrid,
85
85
     1     q2grid,logq2grid,pdfgrid
161
161
      integer nfl,nx,nq2,mem,rep
162
162
      double precision alphas
163
163
      double precision xgrid(100),logxgrid(100)
164
 
      double precision q2grid(50),logq2grid(50)
165
 
      double precision pdfgrid(0:100,14,100,50)
 
164
      double precision q2grid(60),logq2grid(60)
 
165
      double precision pdfgrid(0:100,14,100,60)
166
166
      logical hasphoton
167
167
      common /nnpdf/nfl,nx,nq2,mem,rep,hasphoton,alphas,xgrid,logxgrid,
168
168
     1     q2grid,logq2grid,pdfgrid
171
171
      Q2 = Q*Q
172
172
*     check bounds
173
173
      if (x.lt.xmingrid.or.x.lt.xgrid(1).or.x.gt.xgrid(nx)) then
174
 
         write(6,*) "Parton interpolation: x out of range -- freezed"
 
174
c$$$         write(6,*) "Parton interpolation: x out of range -- freezed"
175
175
         if (x.lt.xgrid(1)) x = xgrid(1)
176
176
*         if (x.lt.xmingrid) x = xmingrid
177
177
         if (x.gt.xgrid(nx))x = xgrid(nx)
178
178
      endif
179
179
      if (Q2.lt.q2grid(1).or.Q2.gt.q2grid(nq2)) then
180
 
         write(6,*) "Parton interpolation: Q2 out of range -- freezed"
181
 
         write(6,*) "Q2 = ",Q2, " GeV2", q2grid(1)
 
180
c$$$         write(6,*) "Parton interpolation: Q2 out of range -- freezed"
 
181
c$$$         write(6,*) "Q2 = ",Q2, " GeV2", q2grid(1)
182
182
         if (Q2.lt.q2grid(1)) Q2 = q2grid(1)
183
183
         if (Q2.gt.q2grid(nq2)) Q2 = q2grid(nq2)
184
184
      endif