LHAPDF/lhapdf5.3.1/wrapevolve.f

   1       subroutine evolvePDF(x,Q,f)
   2       implicit none
   3       integer nset,imem,Eorder,IP2
   4       real*8 x,Q,P2,Q2fit,f(-6:6),alfas,a,photon
   5       nset = 1
   6       call evolvePDFM(nset,x,Q,f)
   7       return
   8 c
   9       entry evolvePDFp(x,Q,P2,IP2,f)
  10       nset = 1
  11       call evolvePDFpM(nset,x,Q,P2,IP2,f)
  12       return
  13 c
  14       entry evolvePDFa(x,Q,a,f)
  15       nset = 1
  16       call evolvePDFaM(nset,x,Q,a,f)
  17       return
  18 c
  19       entry evolvePDFphoton(x,Q,f,photon)
  20       nset = 1
  21       call evolvePDFphotonM(nset,x,Q,f,photon)
  22       return
  23 c
  24       entry initPDF(imem)
  25       nset = 1
  26       call initPDFM(nset,imem)
  27
  28       return
  29       end
  30 c
  31       subroutine evolvePDFaM(nset,xin,Qin,a,f)
  32       implicit none
  33       real*8 x,Q,a,f(-6:6)
  34       real*8 ruv,rdv,ru,rd,rs,rc,rb,rt,rg
  35       integer nset,iimem,j
  36       real*8 xin,qin,q2in
  37       character*20 lparm
  38       real*8 xmin,xmax,q2min,q2max
  39
  40       call getlhaparm(18,lparm)
  41       if(lparm.ne.'EXTRAPOLATE') then
  42         call getnmem(nset,iimem)
  43         call getminmaxm(nset,iimem,xmin,xmax,q2min,q2max)
  44         x = max(xmin,min(xmax,xin))
  45         q2in = qin**2
  46         q = sqrt(max(0d0,q2min,min(q2max,q2in)))
  47       else
  48         x = xin
  49         q = qin
  50       endif
  51
  52       call eks98(x,q,a,ruv,rdv,ru,rd,rs,rc,rb,rt,rg)
  53       call evolvePDFM(nset,x,Q,f)
  54       f(0) = f(0)*rg
  55       f(1) = f(1)*rdv-f(-1)*(rdv-rd)
  56       f(2) = f(2)*ruv-f(-2)*(ruv-ru)
  57       f(3) = f(3)*rs
  58
  59       f(4) = f(4)*rc
  60       f(5) = f(5)*rb
  61       f(6) = f(6)*rt
  62       f(-6) = f(-6)*rt
  63       f(-5) = f(-5)*rb
  64       f(-4) = f(-4)*rc
  65       f(-3) = f(-3)*rs
  66       f(-2) = f(-2)*ru
  67       f(-1) = f(-1)*rd
  68
  69       return
  70       end
  71 c
  72       subroutine evolvePDFM(nset,xin,Qin,f)
  73       implicit none
  74       include 'parmsetup.inc'
  75       integer Eorder,index,imem
  76       character*16 name(nmxset)
  77       integer nmem(nmxset),ndef(nmxset),mem,ip2
  78       common/NAME/name,nmem,ndef,mem
  79       integer iset,iimem
  80       common/SET/iset,iimem
  81       integer nset,j
  82       real*8 x,xin,Q,Qin,Q2fit,alfas,p2,q2in
  83       real*8 f(-6:6),photon
  84       character*20 lparm
  85       real*8 xmin,xmax,q2min,q2max
  86       save
  87 *
  88       call setnset(nset)
  89 c
  90 c      print *,'this is evolvePDFM, name=',nset,name(nset)
  91 c   set all f's to 0.0d0 at start
  92 c      do j = -6,6
  93 c        f(j) = 0.0d0
  94 c      enddo
  95
  96       call getlhaparm(18,lparm)
  97       if(lparm.ne.'EXTRAPOLATE') then
  98         call getnmem(nset,iimem)
  99         call getminmaxm(nset,iimem,xmin,xmax,q2min,q2max)
 100         x = max(xmin,min(xmax,xin))
 101         q2in = qin**2
 102         q = sqrt(max(0d0,q2min,min(q2max,q2in)))
 103       else
 104         x = xin
 105         q = qin
 106       endif
 107
 108       if (name(nset).eq.'QCDNUM') call QCDNUMevolve(x,Q,f)
 109       if (name(nset).eq.'QCDNUM_MRST') call QCDNUMevolve(x,Q,f)
 110       if (name(nset).eq.'QCDNUM_MRST3') call QCDNUM3evolve(x,Q,f)
 111       if (name(nset).eq.'QCDNUM_MRST4') call QCDNUM4evolve(x,Q,f)
 112       if (name(nset)(1:11).eq.'QCDNUM_ZEUS') call ZEUSevolve(x,Q,f)
 113       if (name(nset).eq.'CTEQ5grid') call CTEQ5evolve(x,Q,f)
 114       if (name(nset).eq.'CTEQ6grid') call CTEQ6evolve(x,Q,f)
 115       if (name(nset).eq.'CTEQ65grid') call CTEQ65evolve(x,Q,f)
 116       if (name(nset).eq.'CTEQ65cgrid') call CTEQ65cevolve(x,Q,f)
 117       if (name(nset).eq.'CTEQ6ABgrid') call CTEQ6evolve(x,Q,f)
 118       if (name(nset).eq.'EVLCTEQ') call EVLCTEQevolve(x,Q,f)
 119       if (name(nset).eq.'MRSTgrid') call MRSTevolve(x,Q,f)
 120       if (name(nset).eq.'MRST3grid') call MRSTevolve(x,Q,f)
 121       if (name(nset).eq.'MRST4grid') call MRSTevolve(x,Q,f)
 122       if (name(nset).eq.'MRST98grid') call MRST98evolve(x,Q,f)
 123       if (name(nset).eq.'MRST2006grid') call MRST2006evolve(x,Q,f)
 124       if (name(nset).eq.'MRSTpdf') call QCDNUMevolve(x,Q,f)
 125       if (name(nset).eq.'MRST') call QCDNUMevolve(x,Q,f)
 126       if (name(nset).eq.'A02') call A02evolve(x,Q,f)
 127       if (name(nset).eq.'A02M') call A02Mevolve(x,Q,f)
 128       if (name(nset).eq.'H12000') call H1evolve(x,Q,f)
 129       if (name(nset).eq.'GRV') call GRVevolve(x,Q,f)
 130       if (name(nset).eq.'OWP') call OWPevolve(x,Q,f)
 131       if (name(nset).eq.'SMRSP') call SMRSPevolve(x,Q,f)
 132       if (name(nset).eq.'GRVP0') call GRVP0evolve(x,Q,f)
 133       if (name(nset).eq.'GRVP1') call GRVP1evolve(x,Q,f)
 134       if (name(nset).eq.'ABFKWP') call ABFKWPevolve(x,Q,f)
 135       if (name(nset).eq.'USER') call USERevolve(x,Q,f)
 136       return
 137 *
 138       entry evolvePDFpM(nset,xin,Qin,P2,IP2,f)
 139 c
 140       call setnset(nset)
 141 c
 142       call getlhaparm(18,lparm)
 143       if(lparm.ne.'EXTRAPOLATE') then
 144         call getnmem(nset,iimem)
 145         call getminmaxm(nset,iimem,xmin,xmax,q2min,q2max)
 146         x = max(xmin,min(xmax,xin))
 147         q2in = qin**2
 148         q = sqrt(max(0d0,q2min,min(q2max,q2in)))
 149       else
 150         x = xin
 151         q = qin
 152       endif
 153 c
 154       if(name(nset).eq.'SASG') call SASGevolvep(x,Q,P2,IP2,f)
 155       if(name(nset).eq.'GRVG0') call GRVGevolvep0(x,Q,P2,IP2,f)
 156       if(name(nset).eq.'GRVG1') call GRVGevolvep1(x,Q,P2,IP2,f)
 157       if (name(nset).eq.'DOG0') call DOGevolvep0(x,Q,P2,IP2,f)
 158       if (name(nset).eq.'DOG1') call DOGevolvep1(x,Q,P2,IP2,f)
 159       if (name(nset).eq.'DGG') call DGGevolvep(x,Q,P2,IP2,f)
 160       if (name(nset).eq.'LACG') call LACGevolvep(x,Q,P2,IP2,f)
 161       if (name(nset).eq.'GSG0') call GSGevolvep0(x,Q,P2,IP2,f)
 162       if (name(nset).eq.'GSG1') call GSGevolvep1(x,Q,P2,IP2,f)
 163       if (name(nset).eq.'GSG960') call GSG96evolvep0(x,Q,P2,IP2,f)
 164       if (name(nset).eq.'GSG961') call GSG96evolvep1(x,Q,P2,IP2,f)
 165       if (name(nset).eq.'ACFGP') call ACFGPevolvep(x,Q,P2,IP2,f)
 166       if (name(nset).eq.'WHITG') call WHITGevolvep(x,Q,P2,IP2,f)
 167       return
 168 c
 169       entry evolvePDFphotonM(nset,xin,qin,f,photon)
 170
 171 c
 172       call setnset(nset)
 173 c
 174       call getlhaparm(18,lparm)
 175       if(lparm.ne.'EXTRAPOLATE') then
 176         call getnmem(nset,iimem)
 177         call getminmaxm(nset,iimem,xmin,xmax,q2min,q2max)
 178         x = max(xmin,min(xmax,xin))
 179         q2in = qin**2
 180         q = sqrt(max(0d0,q2min,min(q2max,q2in)))
 181       else
 182         x = xin
 183         q = qin
 184       endif
 185 c
 186       if(name(nset).eq.'MRST4qed') then
 187         call MRSTqedevolve(x,Q,f,photon)
 188       else
 189         photon = 0.0d0
 190       endif
 191
 192       return
 193
 194
 195       entry readevolve(nset)
 196 *
 197       read(1,*) name(nset)
 198 c      print *, 'this is readevolve', name(nset)
 199 *
 200       call setnset(nset)
 201 c
 202       if (name(nset).eq.'QCDNUM') call QCDNUMread(nset)
 203       if (name(nset).eq.'QCDNUM_MRST') call QCDNUMread(nset)
 204       if (name(nset).eq.'QCDNUM_MRST3') call QCDNUM3read(nset)
 205       if (name(nset).eq.'QCDNUM_MRST4') call QCDNUM4read(nset)
 206       if (name(nset)(1:11).eq.'QCDNUM_ZEUS') call ZEUSread(nset)
 207       if (name(nset).eq.'CTEQ5grid') call CTEQ5read(nset)
 208       if (name(nset).eq.'CTEQ6grid') call CTEQ6read(nset)
 209       if (name(nset).eq.'CTEQ65grid') call CTEQ65read(nset)
 210       if (name(nset).eq.'CTEQ65cgrid') call CTEQ65cread(nset)
 211       if (name(nset).eq.'CTEQ6ABgrid') call CTEQ6read(nset)
 212       if (name(nset).eq.'EVLCTEQ') call EVLCTEQread(nset)
 213       if (name(nset).eq.'MRSTgrid') call MRSTread(nset)
 214       if (name(nset).eq.'MRST3grid') call MRSTread(nset)
 215       if (name(nset).eq.'MRST4grid') call MRSTread(nset)
 216       if (name(nset).eq.'MRST4qed') call MRSTqedread(nset)
 217       if (name(nset).eq.'MRST98grid') call MRST98read(nset)
 218       if (name(nset).eq.'MRST2006grid') call MRST2006read(nset)
 219       if (name(nset).eq.'MRSTpdf') call QCDNUMread(nset)
 220       if (name(nset).eq.'MRST') call QCDNUMread(nset)
 221       if (name(nset).eq.'A02') call A02read(nset)
 222       if (name(nset).eq.'A02M') call A02Mread(nset)
 223       if (name(nset).eq.'H12000') call H1read(nset)
 224       if (name(nset).eq.'GRV') call GRVread(nset)
 225       if (name(nset).eq.'SASG') call SASGread(nset)
 226       if (name(nset).eq.'GRVG0' .OR.
 227      +    name(nset).eq.'GRVG1') call GRVGread(nset)
 228       if (name(nset).eq.'DOG0' .OR.
 229      +    name(nset).eq.'DOG1') call DOGread(nset)
 230       if (name(nset).eq.'DGG') call DGGread(nset)
 231       if (name(nset).eq.'LACG') call LACGread(nset)
 232       if (name(nset).eq.'GSG0' .OR.
 233      +    name(nset).eq.'GSG1') call GSGread(nset)
 234       if (name(nset).eq.'GSG960' .OR.
 235      +    name(nset).eq.'GSG961') call GSG96read(nset)
 236       if (name(nset).eq.'ACFGP') call ACFGPread(nset)
 237       if (name(nset).eq.'WHITG') call WHITGread(nset)
 238       if (name(nset).eq.'OWP') call OWPread(nset)
 239       if (name(nset).eq.'SMRSP') call SMRSPread(nset)
 240       if (name(nset).eq.'GRVP0' .OR.
 241      +    name(nset).eq.'GRVP1') call GRVPread(nset)
 242       if (name(nset).eq.'ABFKWP') call ABFKWPread(nset)
 243       if (name(nset).eq.'USER') call USERread(nset)
 244       return
 245 *
 246       entry alfasevolve(nset,alfas,Qin)
 247 c
 248       call setnset(nset)
 249 c      print *,'from alpfasevolveM',nset,Q,name(nset)
 250
 251       call getlhaparm(18,lparm)
 252       if(lparm.ne.'EXTRAPOLATE') then
 253         call getnmem(nset,iimem)
 254         call getq2minm(nset,iimem,q2min)
 255         call getq2maxm(nset,iimem,q2max)
 256         q2in = qin**2
 257         q = sqrt(max(0d0,q2min,min(q2max,q2in)))
 258       else
 259         q = qin
 260       endif
 261 c
 262 c      print *,'calling alfasevolve',alfas,q
 263       if (name(nset).eq.'QCDNUM') call QCDNUMalfa(alfas,Q)
 264       if (name(nset).eq.'QCDNUM_MRST') call QCDNUMalfa(alfas,Q)
 265       if (name(nset).eq.'QCDNUM_MRST3') call QCDNUM3alfa(alfas,Q)
 266       if (name(nset).eq.'QCDNUM_MRST4') call QCDNUM4alfa(alfas,Q)
 267       if (name(nset)(1:11).eq.'QCDNUM_ZEUS') call ZEUSalfa(alfas,Q)
 268       if (name(nset).eq.'CTEQ5grid') call CTEQ5alfa(alfas,Q)
 269       if (name(nset).eq.'CTEQ6grid') call CTEQ6alfa(alfas,Q)
 270       if (name(nset).eq.'CTEQ65grid') call CTEQ65alfa(alfas,Q)
 271       if (name(nset).eq.'CTEQ65cgrid') call CTEQ65alfa(alfas,Q)
 272       if (name(nset).eq.'EVLCTEQ') call EVLCTEQalfa(alfas,Q)
 273       if (name(nset).eq.'MRSTgrid') call MRSTalfa(5,alfas,Q)
 274       if (name(nset).eq.'MRST3grid') call MRSTalfa(3,alfas,Q)
 275       if (name(nset).eq.'MRST4grid') call MRSTalfa(4,alfas,Q)
 276       if (name(nset).eq.'MRST4qed') call MRSTalfa(4,alfas,Q)
 277       if (name(nset).eq.'MRST98grid') call MRST98alfa(alfas,Q)
 278       if (name(nset).eq.'MRST2006grid') call MRST2006alfa(5,alfas,Q)
 279       if (name(nset).eq.'MRSTpdf') call QCDNUMalfa(alfas,Q)
 280       if (name(nset).eq.'MRST') call QCDNUMalfa(alfas,Q)
 281       if (name(nset).eq.'A02') call A02alfa(alfas,Q)
 282       if (name(nset).eq.'A02M') call A02Malfa(alfas,Q)
 283       if (name(nset).eq.'H12000') call H1alfa(alfas,Q)
 284       if (name(nset).eq.'GRV') call GRValfa(alfas,Q)
 285       if (name(nset).eq.'SASG') call SASGalfa(alfas,Q)
 286       if (name(nset).eq.'GRVG0' .OR.
 287      +    name(nset).eq.'GRVG1') call GRVGalfa(alfas,Q)
 288       if (name(nset).eq.'DOG0' .OR.
 289      +    name(nset).eq.'DOG1') call DOGalfa(alfas,Q)
 290       if (name(nset).eq.'DGG') call DGGalfa(alfas,Q)
 291       if (name(nset).eq.'LACG') call LACGalfa(alfas,Q)
 292       if (name(nset).eq.'GSG0' .OR.
 293      +    name(nset).eq.'GSG1') call GSGalfa(alfas,Q)
 294       if (name(nset).eq.'GSG960' .OR.
 295      +    name(nset).eq.'GSG961') call GSG96alfa(alfas,Q)
 296       if (name(nset).eq.'ACFGP') call ACFGPalfa(alfas,Q)
 297       if (name(nset).eq.'WHITG') call WHITGalfa(alfas,Q)
 298       if (name(nset).eq.'OWP') call OWPalfa(alfas,Q)
 299       if (name(nset).eq.'SMRSP') call SMRSPalfa(alfas,Q)
 300       if (name(nset).eq.'GRVP0' .OR.
 301      +    name(nset).eq.'GRVP1') call GRVPalfa(alfas,Q)
 302       if (name(nset).eq.'ABFKWP') call ABFKWPalfa(alfas,Q)
 303       if (name(nset).eq.'USER') call USERalfa(alfas,Q)
 304       return
 305 *
 306       entry initevolution(nset,Eorder,Q2fit)
 307 c
 308       call setnset(nset)
 309 c
 310       if (name(nset).eq.'QCDNUM') call QCDNUMinit(nset,Eorder,Q2fit)
 311       if (name(nset)(1:11).eq.'QCDNUM_ZEUS')
 312      + call ZEUSinit(nset,Eorder,Q2fit)
 313       if (name(nset).eq.'CTEQ5grid') call CTEQ5init(Eorder,Q2fit)
 314       if (name(nset).eq.'CTEQ6grid') call CTEQ6init(Eorder,Q2fit)
 315       if (name(nset).eq.'CTEQ65grid') call CTEQ65init(Eorder,Q2fit)
 316       if (name(nset).eq.'CTEQ65cgrid') call CTEQ65init(Eorder,Q2fit)
 317       if (name(nset).eq.'EVLCTEQ') call EVLCTEQinit(nset,Eorder,Q2fit)
 318       if (name(nset).eq.'MRSTgrid') call MRSTinit(Eorder,Q2fit)
 319       if (name(nset).eq.'MRST3grid') call MRSTinit(Eorder,Q2fit)
 320       if (name(nset).eq.'MRST4grid') call MRSTinit(Eorder,Q2fit)
 321       if (name(nset).eq.'MRST4qed') call MRSTinit(Eorder,Q2fit)
 322       if (name(nset).eq.'MRST98grid') call MRST98init(Eorder,Q2fit)
 323       if (name(nset).eq.'MRST2006grid') call MRST2006init(Eorder,Q2fit)
 324       if (name(nset).eq.'QCDNUM_MRST') then
 325         call QCDNUMinit(nset,Eorder,Q2fit)
 326         call QNLSET('BMARK',.TRUE.)
 327       endif
 328       if (name(nset).eq.'QCDNUM_MRST3') then
 329         call QCDNUM3init(nset,Eorder,Q2fit)
 330         call QNLSET('BMARK',.TRUE.)
 331       endif
 332       if (name(nset).eq.'QCDNUM_MRST4') then
 333         call QCDNUM4init(nset,Eorder,Q2fit)
 334         call QNLSET('BMARK',.TRUE.)
 335       endif
 336       if (name(nset).eq.'A02') call A02init
 337       if (name(nset).eq.'A02M') call A02Minit
 338       if (name(nset).eq.'H12000') call H1init(Eorder,Q2fit)
 339       if (name(nset).eq.'GRV') call GRVinit(Eorder,Q2fit)
 340       if (name(nset).eq.'SASG') call SASGinit(Eorder,Q2fit)
 341       if (name(nset).eq.'GRVG0' .OR.
 342      +    name(nset).eq.'GRVG1') call GRVGinit(Eorder,Q2fit)
 343       if (name(nset).eq.'DOG0' .OR.
 344      +    name(nset).eq.'DOG1') call DOGinit(Eorder,Q2fit)
 345       if (name(nset).eq.'DGG') call DGGinit(Eorder,Q2fit)
 346       if (name(nset).eq.'LACG') call LACGinit(Eorder,Q2fit)
 347       if (name(nset).eq.'GSG0' .OR.
 348      +    name(nset).eq.'GSG1') call GSGinit(Eorder,Q2fit)
 349       if (name(nset).eq.'GSG960' .OR.
 350      +    name(nset).eq.'GSG961') call GSG96init(Eorder,Q2fit)
 351       if (name(nset).eq.'ACFGP') call ACFGPinit(Eorder,Q2fit)
 352       if (name(nset).eq.'WHITG') call WHITGinit(Eorder,Q2fit)
 353       if (name(nset).eq.'OWP') call OWPinit(Eorder,Q2fit)
 354       if (name(nset).eq.'SMRSP') call SMRSPinit(Eorder,Q2fit)
 355       if (name(nset).eq.'GRVP0' .OR.
 356      +    name(nset).eq.'GRVP1') call GRVPinit(Eorder,Q2fit)
 357       if (name(nset).eq.'ABFKWP') call ABFKWPinit(Eorder,Q2fit)
 358       if (name(nset).eq.'USER') call USERinit(nset,Eorder,Q2fit)
 359       return
 360 *
 361       entry initPDFM(nset,imem)
 362 c
 363       call setnset(nset)
 364       call setnmem(nset,imem)
 365 c
 366 c      print *,'entered initPDFM,',nset,imem,name(nset)
 367       iimem = imem
 368       if (name(nset).eq.'QCDNUM') then
 369          call InitEvolvePDF(nset,imem)
 370          call QCDNUMpdf(nset)
 371       endif
 372       if (name(nset).eq.'QCDNUM_MRST') then
 373          call InitEvolvePDF(nset,imem)
 374          call QCDNUMpdf(nset)
 375       endif
 376       if (name(nset).eq.'QCDNUM_MRST3') then
 377          call InitEvolvePDF(nset,imem)
 378          call QCDNUM3pdf(nset)
 379       endif
 380       if (name(nset).eq.'QCDNUM_MRST4') then
 381          call InitEvolvePDF(nset,imem)
 382          call QCDNUM4pdf(nset)
 383       endif
 384       if (name(nset)(1:11).eq.'QCDNUM_ZEUS') then
 385          call InitEvolvePDF(nset,imem)
 386          call ZEUSpdf(nset)
 387       endif
 388       if (name(nset).eq.'MRST') then
 389          call InitEvolvePDF(nset,imem)
 390          call QCDNUMpdf(nset)
 391       endif
 392       if (name(nset).eq.'MRSTpdf') then
 393          call InitEvolvePDF(nset,imem)
 394          call QCDNUMpdf(nset)
 395       endif
 396       if (name(nset).eq.'EVLCTEQ') then
 397          call InitEvolvePDF(nset,imem)
 398          call EVLCTEQpdf(nset)
 399 c          call EVLCTEQpdf(nset,imem)
 400       endif
 401       if (name(nset).eq.'CTEQ65grid') then
 402          call CTEQ6NewAlpha(nset,imem)
 403       endif
 404       if (name(nset).eq.'CTEQ65cgrid') then
 405          call CTEQ6NewAlpha(nset,imem)
 406       endif
 407       if (name(nset).eq.'CTEQ6ABgrid') then
 408          call CTEQ6NewAlpha(nset,imem)
 409 c         call CTEQ6pdf(nset)
 410       endif
 411       if (name(nset).eq.'H12000') then
 412          call InitEvolvePDF(nset,imem)
 413          call H1pdf(imem)
 414       endif
 415       if (name(nset).eq.'USER') then
 416          call InitEvolvePDF(nset,imem)
 417          call USERpdf(nset)
 418       endif
 419       if (name(nset).eq.'CTEQ5grid') call CTEQ5pdf(imem)
 420       if (name(nset).eq.'CTEQ6grid') call CTEQ6pdf(imem)
 421       if (name(nset).eq.'CTEQ65grid') call CTEQ65pdf(imem)
 422       if (name(nset).eq.'CTEQ65cgrid') call CTEQ65pdf(imem)
 423       if (name(nset).eq.'CTEQ6ABgrid') call CTEQ6pdf(imem)
 424       if (name(nset).eq.'MRSTgrid') call MRSTpdf(imem)
 425       if (name(nset).eq.'MRST3grid') call MRSTpdf(imem)
 426       if (name(nset).eq.'MRST4grid') call MRSTpdf(imem)
 427       if (name(nset).eq.'MRST4qed') call MRSTpdf(imem)
 428       if (name(nset).eq.'MRST98grid') call MRST98pdf(imem)
 429       if (name(nset).eq.'MRST2006grid') call MRST2006pdf(imem)
 430       if (name(nset).eq.'A02') call A02pdf(imem)
 431       if (name(nset).eq.'A02M') call A02Mpdf(imem)
 432       if (name(nset).eq.'GRV0' .OR.
 433      +    name(nset).eq.'GRV1') call GRVpdf(imem)
 434       if (name(nset).eq.'SASG') call SASGpdf(imem)
 435       if (name(nset).eq.'GRVG') call GRVGpdf(imem)
 436       if (name(nset).eq.'DOG0' .OR.
 437      +    name(nset).eq.'DOG1') call DOGpdf(imem)
 438       if (name(nset).eq.'DGG') call DGGpdf(imem)
 439       if (name(nset).eq.'LACG') call LACGpdf(imem)
 440       if (name(nset).eq.'GSG0' .OR.
 441      +    name(nset).eq.'GSG1') call GSGpdf(imem)
 442       if (name(nset).eq.'GSG960' .OR.
 443      +    name(nset).eq.'GSG961') call GSG96pdf(imem)
 444       if (name(nset).eq.'ACFGP') call ACFGPpdf(imem)
 445       if (name(nset).eq.'WHITG') call WHITGpdf(imem)
 446       if (name(nset).eq.'OWP') call OWPpdf(imem)
 447       if (name(nset).eq.'SMRSP') call SMRSPpdf(imem)
 448       if (name(nset).eq.'GRVP0' .OR.
 449      +    name(nset).eq.'GRVP1') call GRVPpdf(imem)
 450       if (name(nset).eq.'ABFKWP') call ABFKWPpdf(imem)
 451       return
 452 *
 453       end