-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathwjc.f
More file actions
118 lines (106 loc) · 3.87 KB
/
wjc.f
File metadata and controls
118 lines (106 loc) · 3.87 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
C **********************************************************************
SUBROUTINE WJC1 (q,u,w,vt,vs)
C
C Deuteron wavefunction from WJC-1 (Gross et al.) NN potential model.
C Gross & Stadler, Phys. Rev. C78, 014005 (2008).
C
C Note: q in 1/fm, wfns in fm^3/2.
C
C Wave function normalization
C \int dq q^2 (u^2+w^2+vt^2+vs^2) + V' term = 1
C so that wave functions themselves normalized to ~105%
C => renormalize to 1 for structure functions
C
C **********************************************************************
IMPLICIT NONE
INTEGER id,nq
PARAMETER (nq=60)
REAL*8 q,u,w,vt,vs,DQDVAL,
& qgrid(nq),ugrid(nq),wgrid(nq),vtgrid(nq),vsgrid(nq),dum
LOGICAL init /.FALSE./
REAL*8 pi,hcM,hcG
SAVE
pi = 4*DATAN(1.D0)
hcM = 197.327D0 ! GeV.fm conversion factor
hcG = 0.197327D0 ! MeV.fm conversion factor
IF (init) GO TO 999 ! Data already read
C...Read data from file
OPEN (10, FORM='FORMATTED',
& FILE='/u/home/wmelnitc/Work/EMC/D/Wfn/wjc-1.dat',
& STATUS='OLD')
C...Momentum space [qgrid in MeV, ugrid in GeV^-3/2]
DO id=1,nq
READ (10,*) qgrid(id),
& ugrid(id), wgrid(id), vtgrid(id), vsgrid(id)
qgrid(id) = qgrid(id) / hcM ! MeV => 1/fm
ugrid(id) = ugrid(id) * hcG**1.5D0 ! GeV^-3/2 => fm^3/2
wgrid(id) = wgrid(id) * hcG**1.5D0
vtgrid(id) = vtgrid(id) * hcG**1.5D0
vsgrid(id) = vsgrid(id) * hcG**1.5D0
ENDDO
PRINT *, '... WJC-1 model read...'
init = .TRUE.
C...Evaluate wavefunction
c 999 u = DQDVAL (q,nq,qgrid,ugrid,.FALSE.)
c w = DQDVAL (q,nq,qgrid,wgrid,.FALSE.)
c vt = DQDVAL (q,nq,qgrid,vtgrid,.FALSE.)
c vs = DQDVAL (q,nq,qgrid,vsgrid,.FALSE.)
999 CALL Pinterp (qgrid,ugrid,nq,q,u,dum,1)
CALL Pinterp (qgrid,wgrid,nq,q,w,dum,1)
CALL Pinterp (qgrid,vtgrid,nq,q,vt,dum,1)
CALL Pinterp (qgrid,vsgrid,nq,q,vs,dum,1)
RETURN
END
C **********************************************************************
SUBROUTINE WJC2 (q,u,w,vt,vs)
C
C Deuteron wavefunction from WJC-2 (Gross et al.) NN potential model.
C Gross & Stadler, Phys. Rev. C78, 014005 (2008).
C
C Note: q in 1/fm, wfns in fm^3/2.
C
C Wave function normalization
C \int dq q^2 (u^2+w^2+vt^2+vs^2) + V' term = 1
C so that wave functions themselves normalized to ~102%
C => renormalize to 1 for structure functions
C
C **********************************************************************
IMPLICIT NONE
INTEGER id,nq
PARAMETER (nq=60)
REAL*8 q,u,w,vt,vs,
& qgrid(nq),ugrid(nq),wgrid(nq),vtgrid(nq),vsgrid(nq),dum
LOGICAL init /.FALSE./
REAL*8 pi,hcM,hcG
SAVE
pi = 4*DATAN(1.D0)
hcM = 197.327D0 ! GeV.fm conversion factor
hcG = 0.197327D0 ! MeV.fm conversion factor
IF (init) GO TO 999 ! Data already read
C...Read data from file
OPEN (10, FORM='FORMATTED',
& FILE='/u/home/wmelnitc/Work/EMC/D/Wfn/wjc-2.dat',
& STATUS='OLD')
C...Momentum space [qgrid in MeV, ugrid in GeV^-3/2]
DO id=1,nq
READ (10,*) qgrid(id),
& ugrid(id), wgrid(id), vtgrid(id), vsgrid(id)
qgrid(id) = qgrid(id) / hcM ! MeV => 1/fm
ugrid(id) = ugrid(id) * hcG**1.5D0 ! GeV^-3/2 => fm^3/2
wgrid(id) = wgrid(id) * hcG**1.5D0
vtgrid(id) = vtgrid(id) * hcG**1.5D0
vsgrid(id) = vsgrid(id) * hcG**1.5D0
ENDDO
PRINT *, '... WJC-2 model read...'
init = .TRUE.
C...Evaluate wavefunction
c 999 u = DQDVAL (q,nq,qgrid,ugrid,.FALSE.)
c w = DQDVAL (q,nq,qgrid,wgrid,.FALSE.)
c vt = DQDVAL (q,nq,qgrid,vtgrid,.FALSE.)
c vs = DQDVAL (q,nq,qgrid,vsgrid,.FALSE.)
999 CALL Pinterp (qgrid,ugrid,nq,q,u,dum,1)
CALL Pinterp (qgrid,wgrid,nq,q,w,dum,1)
CALL Pinterp (qgrid,vtgrid,nq,q,vt,dum,1)
CALL Pinterp (qgrid,vsgrid,nq,q,vs,dum,1)
RETURN
END