← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to Template/NLO/SubProcesses/leshouche_inc_chooser.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/NLO/SubProcesses/leshouche_inc_chooser.f
10
10
      integer maxflow
11
11
      parameter (maxflow=999)
12
12
      integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc),
13
 
     &     icolup(2,nexternal,maxflow)
14
 
      common /c_leshouche_inc/idup,mothup,icolup
 
13
     &     icolup(2,nexternal,maxflow),niprocs
 
14
      common /c_leshouche_inc/idup,mothup,icolup,niprocs
15
15
      logical firsttime
16
16
      data firsttime /.true./
17
17
      include 'leshouche_decl.inc'
18
18
      common/c_leshouche_idup_d/ idup_d
19
 
      save mothup_d, icolup_d
 
19
      save mothup_d, icolup_d, niprocs_d
20
20
      
21
21
c
22
22
      if (maxproc_used.gt.maxproc) then
31
31
      endif
32
32
 
33
33
      if (firsttime) then
34
 
        call read_leshouche_info(idup_d,mothup_d,icolup_d)
 
34
        call read_leshouche_info(idup_d,mothup_d,icolup_d,niprocs_d)
35
35
        firsttime = .false.
36
36
      endif
37
37
 
38
 
      do j=1,maxproc_used
 
38
      niprocs=niprocs_d(nFKSprocess)
 
39
      do j=1,niprocs
39
40
         do i=1,nexternal
40
41
            IDUP(i,j)=IDUP_D(nFKSprocess,i,j)
41
42
            MOTHUP(1,i,j)=MOTHUP_D(nFKSprocess,1,i,j)
54
55
      end
55
56
 
56
57
 
57
 
      subroutine read_leshouche_info(idup_d,mothup_d,icolup_d)
 
58
      subroutine read_leshouche_info(idup_d,mothup_d,icolup_d,niprocs_d)
58
59
C read the various information from the configs_and_props_info.dat file
59
60
      implicit none
60
61
      include "nexternal.inc"
91
92
               icolup_d(1,1,nexternal,j)=-99999 ! should not be used
92
93
               icolup_d(1,2,nexternal,j)=-99999
93
94
            enddo
 
95
            niprocs_d(1)=maxproc_used
94
96
            return
95
97
         endif
96
98
      endif
106
108
          do k=1,nexternal
107
109
            idup_d(i,k,j)=itmp_array(k)
108
110
          enddo
 
111
          niprocs_d(i)=j
109
112
        else if (buff(:1).eq.'M') then
110
113
        ! idup
111
114
        ! I  i   j   l   id1 ..idn -> MOTHUP_D(i,j,k,l)=idk