| 3c5d1739 |
1 | subroutine numberPDF(noe) |
| 2 | implicit none |
| 3 | integer nset,noe |
| 4 | nset = 1 |
| 5 | call numberPDFM(nset,noe) |
| 6 | return |
| 7 | end |
| 8 | c |
| 9 | subroutine listPDF(nset,imem,parm) |
| 10 | implicit none |
| 11 | include 'parmsetup.inc' |
| 12 | character*16 name(nmxset) |
| 13 | integer nmem(nmxset),ndef(nmxset) |
| 14 | common/NAME/name,nmem,ndef,mem |
| 15 | character*16 s1 |
| 16 | integer i,j,mem,imem,noe,nop,listN(nmxset),listP(nmxset),type |
| 17 | real*8 parmL(nmxset,0:noemax,nopmax),parm(nopmax) |
| 18 | integer nset |
| 19 | save listN,listP,parmL |
| 20 | * |
| 21 | mem=imem |
| 22 | if (mem.gt.listN(nset)) then |
| 23 | write(*,*) 'Maximum number of PDFs in list exceeded: ', |
| 24 | . mem,' > ',listN |
| 25 | write(*,*) 'Returning most likely PDF' |
| 26 | mem=0 |
| 27 | endif |
| 28 | if (mem.lt.0) then |
| 29 | write(*,*) 'Negative PDF member requested: ',mem |
| 30 | write(*,*) 'Returning most likely PDF' |
| 31 | mem=0 |
| 32 | endif |
| 33 | do i=1,listP(nset) |
| 34 | parm(i)=parmL(nset,mem,i) |
| 35 | enddo |
| 36 | return |
| 37 | * |
| 38 | entry nopPDF(nset,nop) |
| 39 | nop=listP(nset) |
| 40 | return |
| 41 | * |
| 42 | entry numberPDFM(nset,noe) |
| 43 | if(NAME(nset).eq.'MRSTgrid' |
| 44 | +.or.NAME(nset).eq.'MRST98grid' |
| 45 | +.or.NAME(nset).eq.'A02' |
| 46 | +.or.NAME(nset).eq.'A02M' |
| 47 | +.or.NAME(nset).eq.'CTEQ5grid' |
| 48 | +.or.NAME(nset).eq.'CTEQ6grid' |
| 49 | +.or.NAME(nset).eq.'CTEQ6ABgrid' !*** added *** |
| 50 | +.or.NAME(nset).eq.'SASG' |
| 51 | +.or.NAME(nset).eq.'GRVG' |
| 52 | +.or.NAME(nset).eq.'DOG' |
| 53 | +.or.NAME(nset).eq.'DGG' |
| 54 | +.or.NAME(nset).eq.'LACG' |
| 55 | +.or.NAME(nset).eq.'GSG' |
| 56 | +.or.NAME(nset).eq.'GSG96' |
| 57 | +.or.NAME(nset).eq.'ACFGP' |
| 58 | +.or.NAME(nset).eq.'WHITG' |
| 59 | +.or.NAME(nset).eq.'OWP' |
| 60 | +.or.NAME(nset).eq.'SMRSP' |
| 61 | +.or.NAME(nset).eq.'GRVP' |
| 62 | +.or.NAME(nset).eq.'ABFKWP' |
| 63 | c +.or.NAME.eq.'H12000' |
| 64 | c +.or.NAME.eq.'GRV' |
| 65 | +) then |
| 66 | noe = nmem(nset) |
| 67 | else |
| 68 | noe=listN(nset) |
| 69 | endif |
| 70 | |
| 71 | return |
| 72 | * |
| 73 | entry InitListPDF(nset) |
| 74 | type=-1 |
| 75 | read(1,*) s1,listN(nset),listP(nset) |
| 76 | c print *,s1,listN,listP |
| 77 | if (index(s1,'list').eq.1) then |
| 78 | type=1 |
| 79 | do i=0,listN(nset) |
| 80 | read(1,*) (parmL(nset,i,j),j=1,listP(nset)) |
| 81 | enddo |
| 82 | c print *,parmL(0,1) |
| 83 | endif |
| 84 | if (type.lt.0) then |
| 85 | write(*,*) 'File description error:' |
| 86 | write(*,*) 'Unknown parameter list type ',s1 |
| 87 | stop |
| 88 | endif |
| 89 | return |
| 90 | * |
| 91 | end |
| 92 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame