1 SUBROUTINE ch3piset(JJ)
11 SUBROUTINE ch3piget(JJ)
18 IF (j3pi.EQ.1.OR.j3pi.EQ.2)
THEN
21 write(*,*)
'FROM value_parameter.f CH3PIGET, wrong J3PI=',j3pi
26 SUBROUTINE getff2pirho(JJ)
28 include
'../parameter.inc'
33 SUBROUTINE olachnl(SIGN)
35 COMMON / jaki / jak1,jak2,jakp,jakm,ktom
36 INTEGER JAK1,JAK2,JAKP,JAKM,KTOM
40 IF (ktom.EQ.1.OR.ktom.EQ.-1)
THEN
42 ELSEIF (ktom.EQ.2)
THEN
45 print *,
'STOP IN OLACHNL: KTOM=',ktom
56 DATA pi /3.141592653589793238462643/
104 coefr(1,4)= 1.0/sqrt(2.)/3.0
105 coefr(2,4)=-1.0/sqrt(2.)/3.0
110 coefr(1,5)=-sqrt(2.)/3.0
111 coefr(2,5)= sqrt(2.)/3.0
126 coefr(5,7)=-sqrt(2.0/3.0)
132 subroutine rchl_parameters(KAK)
151 include
'../parameter.inc'
156 COMMON /chanopt/ j3pi
173 ELSE IF(kak.EQ.5)
THEN
201 gamma_kst = 0.06672d0
202 gamma_kstpr = 0.240d0
204 ELSE IF(kak.EQ.71)
THEN
210 gamma_kstpr = 0.206d0
211 gamma_rcht = -0.043d0
216 mkst = (mksp +mks0)/2.
232 beta_rho = -0.318551d0
240 fk_rpt = fpi_rpt*1.198d0
241 gv_rpt = fpi_rpt*fpi_rpt/fv_rpt
262 ELSE IF (j3pi.EQ.2)
THEN
263 alpsig = 1.139486d0*0.63d0
264 betasig = 1.139486d0*0.63d0
265 gamsig = 0.889769d0*0.63d0
266 delsig = 0.889769d0*0.63d0
277 ELSE IF(j3pi.EQ.2)
THEN
286 call rchl_reparam(0,iwarm,kak)
287 IF (iwarm.EQ.1)
RETURN
318 mmpi_av = (mpiz+2.*mpic)/3.d0
332 g2 = mro/(192.*pi*pi*sqrt(2.)*fv_rpt)*3.
337 c1256 = -3/96./pi**2*fv_rpt*mro/sqrt(2.)/fpi_rpt**2
341 d3 = -mro**2/(64.*pi*pi*fpi_rpt**2)
345 ht0 = -1.2400398216503017d-2
349 c1_kpi = lap_kpi/mpic**2
350 c2_kpi = (lapp_kpi - lap_kpi**2)/2.d0/mpic**4
359 b00_3piscal = 0.268/mmpi_av**2
360 c00_3piscal = -0.0139/mmpi_av**4
361 d00_3piscal = -0.00139/mmpi_av**6
362 x00_3piscal = 36.77*mmpi_av**2
363 a02_3piscal = -0.0444
364 b02_3piscal = -0.0857/mmpi_av**2
365 c02_3piscal = -0.00221/mmpi_av**4
366 d02_3piscal = -0.000129/mmpi_av**6
367 x02_3piscal = -21.62*mmpi_av**2
433 call rchl_reparam(1,iwarm,kak)
437 subroutine rchl_reparam(IMODE,IWARM,KAK)
438 include
'../parameter.inc'
439 common / params / p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15
441 DOUBLE PRECISION P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12,P13,P14,P15
446 IF (imode.EQ.-1)
THEN
483 fk_rpt = fpi_rpt*1.198d0
484 gv_rpt = fpi_rpt*fpi_rpt/fv_rpt