[cvs] / MG_ME / Models / sm / couplings.f Repository:
ViewVC logotype

Diff of /MG_ME/Models/sm/couplings.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.10.4.1, Fri Jul 18 14:28:02 2008 UTC revision 1.16, Wed Jun 2 00:18:38 2010 UTC
# Line 85  Line 85 
85        if(readlha) then        if(readlha) then
86           call lh_readin(param_name)           call lh_readin(param_name)
87           G = DSQRT(4d0*PI*ALFAS) ! use setting of the param_card.dat @ NLO           G = DSQRT(4d0*PI*ALFAS) ! use setting of the param_card.dat @ NLO
       endif  
 c  
 c     Strong coupling  
 c  
 c     As a rule we first check if a pdf has been chosen in the  
 c     run_card.dat (which has been already read at this stage).  
 c     If there pdfs in the initial state, then the alpha_s(MZ) used  
 c     is set to the corresponding value.  
   
       GG(1) = -G  
       GG(2) = -G  
88  c  c
89  c auxiliary local values  c auxiliary local values
90  c  c
# Line 200  Line 189 
189    
190        awidth = 0d0        awidth = 0d0
191    
192          endif
193    c
194    c     Strong coupling
195    c
196    c     As a rule we first check if a pdf has been chosen in the
197    c     run_card.dat (which has been already read at this stage).
198    c     If there pdfs in the initial state, then the alpha_s(MZ) used
199    c     is set to the corresponding value.
200    
201          GG(1) = -G
202          GG(2) = -G
203    
204  c----------------------------  c----------------------------
205  c end subroutine coupsm  c end subroutine coupsm
206  c----------------------------  c----------------------------
# Line 285  Line 286 
286           block_found=.true.           block_found=.true.
287             elseif(buff(1:5).eq.'decay') then             elseif(buff(1:5).eq.'decay') then
288                 n=n+1                 n=n+1
289                 l1=30                 l1=60
290                 if(index(buff,"#").ne.0) l1=index(buff,"#")-1 ! ignore comments                 if(index(buff,"#").ne.0) l1=index(buff,"#")-1 ! ignore comments
291                 read(buff(6:l1),*) ivalue(n),value(n)                 read(buff(6:l1),*) ivalue(n),value(n)
292                 name(n)="decay"                 name(n)="decay"
# Line 303  Line 304 
304            endif            endif
305              if(buff(1:1).ne.'#') then  !if it not a comment              if(buff(1:1).ne.'#') then  !if it not a comment
306                n=n+1                n=n+1
307                l1=30                l1=60
308                if(index(buff,"#").ne.0) l1=index(buff,"#")-1 ! ignore comments                if(index(buff,"#").ne.0) l1=index(buff,"#")-1 ! ignore comments
309  c  c
310  c  WARNING:... not all blocks have the same sintax!! You need to change it  c  WARNING:... not all blocks have the same sintax!! You need to change it
# Line 342  Line 343 
343         call set_it(n,ivalue,value,name,15,bn,mtaMS,1.777d0)         call set_it(n,ivalue,value,name,15,bn,mtaMS,1.777d0)
344         bn="mgckm"         bn="mgckm"
345         call set_it(n,ivalue,value,name,1,bn,vud,1d0)         call set_it(n,ivalue,value,name,1,bn,vud,1d0)
        bn="ldme"  
        call set_it(n,ivalue,value,name,1,bn,ldme_1S01,0.3867d0)  
        call set_it(n,ivalue,value,name,2,bn,ldme_1S08,0.35333d-2)  
        call set_it(n,ivalue,value,name,3,bn,ldme_3S11,1.16d0)  
        call set_it(n,ivalue,value,name,4,bn,ldme_3S18,1.06d-2)  
        call set_it(n,ivalue,value,name,5,bn,ldme_3P01,0d0)  
        call set_it(n,ivalue,value,name,6,bn,ldme_3P08,0d0)  
        call set_it(n,ivalue,value,name,7,bn,ldme_3P11,0d0)  
        call set_it(n,ivalue,value,name,8,bn,ldme_3P18,0d0)  
        call set_it(n,ivalue,value,name,9,bn,ldme_3P21,0d0)  
        call set_it(n,ivalue,value,name,10,bn,ldme_3P28,0d0)  
        call set_it(n,ivalue,value,name,11,bn,ldme_1P11,0d0)  
        call set_it(n,ivalue,value,name,12,bn,ldme_1P18,0d0)  
346         bn="mass"         bn="mass"
347         call set_it(n,ivalue,value,name,4,bn,cmass,1.4d0)         call set_it(n,ivalue,value,name,4,bn,cmass,1.4d0)
348         call set_it(n,ivalue,value,name,5,bn,bmass,4.7d0)         call set_it(n,ivalue,value,name,5,bn,bmass,4.7d0)
# Line 498  Line 486 
486        integer lun        integer lun
487        logical fopened        logical fopened
488        character*(*) filename        character*(*) filename
489        character*30  tempname        character*90  tempname
490        integer fine        integer fine
491        integer dirup        integer dirup,i
492    
493  c-----  c-----
494  c  Begin Code  c  Begin Code
# Line 508  Line 496 
496  c  c
497  c     first check that we will end in the main directory  c     first check that we will end in the main directory
498  c  c
       open(unit=lun,file="Source/makefile",status='old',err=20)  
       dirup=0  
       goto 100  
  20   close(lun)  
   
       open(unit=lun,file="../Source/makefile",status='old',err=30)  
       dirup=1  
       goto 100  
  30   close(lun)  
   
       open(unit=lun,file="../../Source/makefile",status='old',err=40)  
       dirup=2  
       goto 100  
  40   close(lun)  
   
       open(unit=lun,file="../../../Source/makefile",status='old',err=50)  
       dirup=3  
       goto 100  
  50   close(lun)  
   
  100  continue  
       close(lun)  
   
       fopened=.true.  
499        tempname=filename        tempname=filename
500        fine=index(tempname,' ')        fine=index(tempname,' ')
501        if(fine.eq.0) fine=len(tempname)        if(fine.eq.0) fine=len(tempname)
502          tempname=tempname(1:fine)
503  c  c
504  c         if I have to read a card  c         if I have to read a card
505  c  c
506            if(index(filename,"_card").gt.0) then            if(index(filename,"_card").gt.0) then
507               tempname='/Cards/'//tempname(1:fine)           tempname='./Cards/'//tempname
              fine=fine+7  
508        endif        endif
509    
       if(dirup.eq.0) open(unit=lun,file=tempname(1:fine),status='old',err=110)  
       if(dirup.eq.1) open(unit=lun,file='../'//tempname(1:fine),status='old',err=110)  
       if(dirup.eq.2) open(unit=lun,file='../../'//tempname(1:fine),status='old',err=110)  
       if(dirup.eq.3) open(unit=lun,file='../../../'//tempname(1:fine),status='old',err=110)  
       return  
510    
511   110  fopened=.false.        fopened=.false.
512        close(lun)        do i=0,5
513        write (*,*) 'Warning: file ',tempname(1:fine),' is not in the main directory'           open(unit=lun,file=tempname,status='old',ERR=30)
514             fopened=.true.
515             write(*,*) 'read model file',tempname
516             exit
517     30      tempname='../'//tempname
518             if (i.eq.5)then
519                write(*,*) 'Warning: file ',tempname,' is not correct'
520                stop
521             endif
522          enddo
523    
524    
525        return        return
526        end        end
527    
528    
529    

Legend:
Removed from v.1.10.4.1  
changed lines
  Added in v.1.16

Mail to administrator
ViewVC Help
Powered by ViewVC 1.0.0