3 C*********************************************************************
5 SUBROUTINE PYWIDT_HIJING(KFLR,RMAS,WDTP,WDTE)
7 C...Calculates full and partial widths of resonances.
8 #include "ludat1_hijing.inc"
9 #include "ludat2_hijing.inc"
10 #include "ludat3_hijing.inc"
11 #include "pypars_hijing.inc"
12 #include "pyint1_hijing.inc"
13 #include "pyint4_hijing.inc"
14 DIMENSION WDTP(0:40),WDTE(0:40,0:5)
16 C...Some common constants.
24 C...Reset width information.
34 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
35 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
36 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 110
39 WDTP(I)=(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
42 WDTP(0)=WDTP(0)+WDTP(I)
43 IF(MDME(IDC,1).GT.0) THEN
44 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
45 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
46 WDTE(I,0)=WDTE(I,MDME(IDC,1))
47 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)
51 ELSEIF(KFLA.EQ.23) THEN
53 IF(MINT(61).EQ.1) THEN
54 EI=KCHG(IABS(MINT(15)),1)/3.
58 GZMZ=PMAS(23,2)*PMAS(23,1)
60 GZI=EI*VI/(8.*XW*(1.-XW))*SQM*(SQM-SQMZ)/
61 & ((SQM-SQMZ)**2+GZMZ**2)
62 ZZI=(VI**2+AI**2)/(16.*XW*(1.-XW))**2*SQM**2/
63 & ((SQM-SQMZ)**2+GZMZ**2)
64 IF(MSTP(43).EQ.1) THEN
65 C...Only gamma* production included
68 ELSEIF(MSTP(43).EQ.2) THEN
69 C...Only Z0 production included
73 ELSEIF(MINT(61).EQ.2) THEN
80 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
81 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
82 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 120
88 IF(MINT(61).EQ.0) THEN
89 WDTP(I)=3.*(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
90 & SQRT(MAX(0.,1.-4.*RM1))*RADC
91 ELSEIF(MINT(61).EQ.1) THEN
92 WDTP(I)=3.*((GGI*EF**2+GZI*EF*VF+ZZI*VF**2)*
93 & (1.+2.*RM1)+ZZI*AF**2*(1.-4.*RM1))*
94 & SQRT(MAX(0.,1.-4.*RM1))*RADC
95 ELSEIF(MINT(61).EQ.2) THEN
96 GGF=3.*EF**2*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))*RADC
97 GZF=3.*EF*VF*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))*RADC
98 ZZF=3.*(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
99 & SQRT(MAX(0.,1.-4.*RM1))*RADC
103 C...Z0 -> l+ + l-, nu + nub
107 WDTP(I)=(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
108 & SQRT(MAX(0.,1.-4.*RM1))
109 IF(MINT(61).EQ.0) THEN
110 WDTP(I)=(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
111 & SQRT(MAX(0.,1.-4.*RM1))
112 ELSEIF(MINT(61).EQ.1) THEN
113 WDTP(I)=((GGI*EF**2+GZI*EF*VF+ZZI*VF**2)*
114 & (1.+2.*RM1)+ZZI*AF**2*(1.-4.*RM1))*
115 & SQRT(MAX(0.,1.-4.*RM1))
116 ELSEIF(MINT(61).EQ.2) THEN
117 GGF=EF**2*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
118 GZF=EF*VF*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
119 ZZF=(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
120 & SQRT(MAX(0.,1.-4.*RM1))
126 IF(MINT(61).EQ.0) THEN
127 WDTP(I)=0.25*CF**2*(1.-4.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
128 ELSEIF(MINT(61).EQ.1) THEN
129 WDTP(I)=0.25*(GGI+GZI*CF+ZZI*CF**2)*(1.-4.*RM1)*
130 & SQRT(MAX(0.,1.-4.*RM1))
131 ELSEIF(MINT(61).EQ.2) THEN
132 GGF=0.25*(1.-4.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
133 GZF=0.25*CF*(1.-4.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
134 ZZF=0.25*CF**2*(1.-4.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
138 WDTP(0)=WDTP(0)+WDTP(I)
139 IF(MDME(IDC,1).GT.0) THEN
140 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
141 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
142 WDTE(I,0)=WDTE(I,MDME(IDC,1))
143 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)
144 VINT(111)=VINT(111)+GGF*WID2
145 VINT(112)=VINT(112)+GZF*WID2
146 VINT(114)=VINT(114)+ZZF*WID2
149 IF(MSTP(43).EQ.1) THEN
150 C...Only gamma* production included
153 ELSEIF(MSTP(43).EQ.2) THEN
154 C...Only Z0 production included
159 ELSEIF(KFLA.EQ.24) THEN
161 DO 130 I=1,MDCY(24,3)
163 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
164 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
165 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 130
168 WDTP(I)=3.*(2.-RM1-RM2-(RM1-RM2)**2)*
169 & SQRT(MAX(0.,(1.-RM1-RM2)**2-4.*RM1*RM2))*
170 & VCKM((I-1)/4+1,MOD(I-1,4)+1)*RADC
173 C...W+/- -> l+/- + nu
174 WDTP(I)=(2.-RM1-RM2-(RM1-RM2)**2)*
175 & SQRT(MAX(0.,(1.-RM1-RM2)**2-4.*RM1*RM2))
178 WDTP(0)=WDTP(0)+WDTP(I)
179 IF(MDME(IDC,1).GT.0) THEN
180 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
181 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
182 WDTE(I,0)=WDTE(I,MDME(IDC,1))
183 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)
187 ELSEIF(KFLA.EQ.25) THEN
189 DO 170 I=1,MDCY(25,3)
191 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
192 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
193 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 170
196 WDTP(I)=3.*RM1*(1.-4.*RM1)*SQRT(MAX(0.,1.-4.*RM1))*RADC
200 WDTP(I)=RM1*(1.-4.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
203 C...H0 -> g + g; quark loop contribution only
207 EPS=(2.*PMAS(J,1)/RMAS)**2
209 IF(EPS.GT.1.E-4) THEN
211 RLN=LOG((1.+ROOT)/(1.-ROOT))
215 PHIRE=0.25*(RLN**2-PARU(1)**2)
216 PHIIM=0.5*PARU(1)*RLN
218 PHIRE=-(ASIN(1./SQRT(EPS)))**2
221 ETARE=ETARE+0.5*EPS*(1.+(EPS-1.)*PHIRE)
222 ETAIM=ETAIM+0.5*EPS*(EPS-1.)*PHIIM
224 ETA2=ETARE**2+ETAIM**2
225 WDTP(I)=(AS/PARU(1))**2*ETA2
228 C...H0 -> gamma + gamma; quark, charged lepton and W loop contributions
231 DO 150 J=1,3*MSTP(1)+1
232 IF(J.LE.2*MSTP(1)) THEN
234 EPS=(2.*PMAS(J,1)/RMAS)**2
235 ELSEIF(J.LE.3*MSTP(1)) THEN
238 EPS=(2.*PMAS(10+JL,1)/RMAS)**2
240 EPS=(2.*PMAS(24,1)/RMAS)**2
243 IF(EPS.GT.1.E-4) THEN
245 RLN=LOG((1.+ROOT)/(1.-ROOT))
249 PHIRE=0.25*(RLN**2-PARU(1)**2)
250 PHIIM=0.5*PARU(1)*RLN
252 PHIRE=-(ASIN(1./SQRT(EPS)))**2
255 IF(J.LE.2*MSTP(1)) THEN
256 ETARE=ETARE+0.5*3.*EJ**2*EPS*(1.+(EPS-1.)*PHIRE)
257 ETAIM=ETAIM+0.5*3.*EJ**2*EPS*(EPS-1.)*PHIIM
258 ELSEIF(J.LE.3*MSTP(1)) THEN
259 ETARE=ETARE+0.5*EJ**2*EPS*(1.+(EPS-1.)*PHIRE)
260 ETAIM=ETAIM+0.5*EJ**2*EPS*(EPS-1.)*PHIIM
262 ETARE=ETARE-0.5-0.75*EPS*(1.+(EPS-2.)*PHIRE)
263 ETAIM=ETAIM+0.75*EPS*(EPS-2.)*PHIIM
266 ETA2=ETARE**2+ETAIM**2
267 WDTP(I)=(AEM/PARU(1))**2*0.5*ETA2
270 C...H0 -> gamma + Z0; quark, charged lepton and W loop contributions
273 DO 160 J=1,3*MSTP(1)+1
274 IF(J.LE.2*MSTP(1)) THEN
278 EPS=(2.*PMAS(J,1)/RMAS)**2
279 EPSP=(2.*PMAS(J,1)/PMAS(23,1))**2
280 ELSEIF(J.LE.3*MSTP(1)) THEN
285 EPS=(2.*PMAS(10+JL,1)/RMAS)**2
286 EPSP=(2.*PMAS(10+JL,1)/PMAS(23,1))**2
288 EPS=(2.*PMAS(24,1)/RMAS)**2
289 EPSP=(2.*PMAS(24,1)/PMAS(23,1))**2
293 IF(EPS.GT.1.E-4) THEN
294 RLN=LOG((1.+ROOT)/(1.-ROOT))
298 PHIRE=0.25*(RLN**2-PARU(1)**2)
299 PHIIM=0.5*PARU(1)*RLN
300 PSIRE=-(1.+0.5*ROOT*RLN)
301 PSIIM=0.5*PARU(1)*ROOT
303 PHIRE=-(ASIN(1./SQRT(EPS)))**2
305 PSIRE=-(1.+SQRT(EPS-1.)*ASIN(1./SQRT(EPS)))
310 IF(EPSP.GT.1.E-4) THEN
311 RLN=LOG((1.+ROOT)/(1.-ROOT))
315 PHIREP=0.25*(RLN**2-PARU(1)**2)
316 PHIIMP=0.5*PARU(1)*RLN
317 PSIREP=-(1.+0.5*ROOT*RLN)
318 PSIIMP=0.5*PARU(1)*ROOT
320 PHIREP=-(ASIN(1./SQRT(EPSP)))**2
322 PSIREP=-(1.+SQRT(EPSP-1.)*ASIN(1./SQRT(EPSP)))
325 FXYRE=EPS*EPSP/(8.*(EPS-EPSP))*(1.-EPS*EPSP/(EPS-EPSP)*(PHIRE-
326 & PHIREP)+2.*EPS/(EPS-EPSP)*(PSIRE-PSIREP))
327 FXYIM=EPS*EPSP/(8.*(EPS-EPSP))*(-EPS*EPSP/(EPS-EPSP)*(PHIIM-
328 & PHIIMP)+2.*EPS/(EPS-EPSP)*(PSIIM-PSIIMP))
329 F1RE=EPS*EPSP/(2.*(EPS-EPSP))*(PHIRE-PHIREP)
330 F1IM=EPS*EPSP/(2.*(EPS-EPSP))*(PHIIM-PHIIMP)
331 IF(J.LE.2*MSTP(1)) THEN
332 ETARE=ETARE-3.*EJ*VJ*(FXYRE-0.25*F1RE)
333 ETAIM=ETAIM-3.*EJ*VJ*(FXYIM-0.25*F1IM)
334 ELSEIF(J.LE.3*MSTP(1)) THEN
335 ETARE=ETARE-EJ*VJ*(FXYRE-0.25*F1RE)
336 ETAIM=ETAIM-EJ*VJ*(FXYIM-0.25*F1IM)
338 ETARE=ETARE-SQRT(1.-XW)*(((1.+2./EPS)*XW/SQRT(1.-XW)-
339 & (5.+2./EPS))*FXYRE+(3.-XW/SQRT(1.-XW))*F1RE)
340 ETAIM=ETAIM-SQRT(1.-XW)*(((1.+2./EPS)*XW/SQRT(1.-XW)-
341 & (5.+2./EPS))*FXYIM+(3.-XW/SQRT(1.-XW))*F1IM)
344 ETA2=ETARE**2+ETAIM**2
345 WDTP(I)=(AEM/PARU(1))**2*(1.-(PMAS(23,1)/RMAS)**2)**3/XW*ETA2
348 C...H0 -> Z0 + Z0, W+ + W-
349 WDTP(I)=(1.-4.*RM1+12.*RM1**2)*SQRT(MAX(0.,1.-4.*RM1))/
353 WDTP(0)=WDTP(0)+WDTP(I)
354 IF(MDME(IDC,1).GT.0) THEN
355 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
356 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
357 WDTE(I,0)=WDTE(I,MDME(IDC,1))
358 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)
362 ELSEIF(KFLA.EQ.32) THEN
364 C.... this line added by A.M.
367 IF(MINT(61).EQ.1) THEN
368 EI=KCHG(IABS(MINT(15)),1)/3.
372 GZMZ=PMAS(23,2)*PMAS(23,1)
376 GZPMZP=PMAS(32,2)*PMAS(32,1)
378 GZI=EI*VI/(8.*XW*(1.-XW))*SQM*(SQM-SQMZ)/
379 & ((SQM-SQMZ)**2+GZMZ**2)
380 GZPI=EI*VPI/(8.*XW*(1.-XW))*SQM*(SQM-SQMZP)/
381 & ((SQM-SQMZP)**2+GZPMZP**2)
382 ZZI=(VI**2+AI**2)/(16.*XW*(1.-XW))**2*SQM**2/
383 & ((SQM-SQMZ)**2+GZMZ**2)
384 ZZPI=2.*(VI*VPI+AI*API)/(16.*XW*(1.-XW))**2*
385 & SQM**2*((SQM-SQMZ)*(SQM-SQMZP)+GZMZ*GZPMZP)/
386 & (((SQM-SQMZ)**2+GZMZ**2)*((SQM-SQMZP)**2+GZPMZP**2))
387 ZPZPI=(VPI**2+API**2)/(16.*XW*(1.-XW))**2*SQM**2/
388 & ((SQM-SQMZP)**2+GZPMZP**2)
389 IF(MSTP(44).EQ.1) THEN
390 C...Only gamma* production included
396 ELSEIF(MSTP(44).EQ.2) THEN
397 C...Only Z0 production included
403 ELSEIF(MSTP(44).EQ.3) THEN
404 C...Only Z'0 production included
410 ELSEIF(MSTP(44).EQ.4) THEN
411 C...Only gamma*/Z0 production included
415 ELSEIF(MSTP(44).EQ.5) THEN
416 C...Only gamma*/Z'0 production included
420 ELSEIF(MSTP(44).EQ.6) THEN
421 C...Only Z0/Z'0 production included
426 ELSEIF(MINT(61).EQ.2) THEN
434 DO 180 I=1,MDCY(32,3)
436 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
437 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
438 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 180
446 IF(MINT(61).EQ.0) THEN
447 WDTP(I)=3.*(VPF**2*(1.+2.*RM1)+APF**2*(1.-4.*RM1))*
448 & SQRT(MAX(0.,1.-4.*RM1))*RADC
449 ELSEIF(MINT(61).EQ.1) THEN
450 WDTP(I)=3.*((GGI*EF**2+GZI*EF*VF+GZPI*EF*VPF+ZZI*VF**2+
451 & ZZPI*VF*VPF+ZPZPI*VPF**2)*(1.+2.*RM1)+(ZZI*AF**2+
452 & ZZPI*AF*APF+ZPZPI*APF**2)*(1.-4.*RM1))*
453 & SQRT(MAX(0.,1.-4.*RM1))*RADC
454 ELSEIF(MINT(61).EQ.2) THEN
455 GGF=3.*EF**2*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))*RADC
456 GZF=3.*EF*VF*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))*RADC
457 GZPF=3.*EF*VPF*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))*RADC
458 ZZF=3.*(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
459 & SQRT(MAX(0.,1.-4.*RM1))*RADC
460 ZZPF=3.*(VF*VPF*(1.+2.*RM1)+AF*APF*(1.-4.*RM1))*
461 & SQRT(MAX(0.,1.-4.*RM1))*RADC
462 ZPZPF=3.*(VPF**2*(1.+2.*RM1)+APF**2*(1.-4.*RM1))*
463 & SQRT(MAX(0.,1.-4.*RM1))*RADC
467 C...Z'0 -> l+ + l-, nu + nub
473 IF(MINT(61).EQ.0) THEN
474 WDTP(I)=(VPF**2*(1.+2.*RM1)+APF**2*(1.-4.*RM1))*
475 & SQRT(MAX(0.,1.-4.*RM1))
476 ELSEIF(MINT(61).EQ.1) THEN
477 WDTP(I)=((GGI*EF**2+GZI*EF*VF+GZPI*EF*VPF+ZZI*VF**2+
478 & ZZPI*VF*VPF+ZPZPI*VPF**2)*(1.+2.*RM1)+(ZZI*AF**2+
479 & ZZPI*AF*APF+ZPZPI*APF**2)*(1.-4.*RM1))*
480 & SQRT(MAX(0.,1.-4.*RM1))
481 ELSEIF(MINT(61).EQ.2) THEN
482 GGF=EF**2*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
483 GZF=EF*VF*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
484 GZPF=EF*VPF*(1.+2.*RM1)*SQRT(MAX(0.,1.-4.*RM1))
485 ZZF=(VF**2*(1.+2.*RM1)+AF**2*(1.-4.*RM1))*
486 & SQRT(MAX(0.,1.-4.*RM1))
487 ZZPF=(VF*VPF*(1.+2.*RM1)+AF*APF*(1.-4.*RM1))*
488 & SQRT(MAX(0.,1.-4.*RM1))
489 ZPZPF=(VPF**2*(1.+2.*RM1)+APF**2*(1.-4.*RM1))*
490 & SQRT(MAX(0.,1.-4.*RM1))
494 WDTP(0)=WDTP(0)+WDTP(I)
495 IF(MDME(IDC,1).GT.0) THEN
496 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
497 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
498 WDTE(I,0)=WDTE(I,MDME(IDC,1))
499 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)
500 VINT(111)=VINT(111)+GGF
501 VINT(112)=VINT(112)+GZF
502 VINT(113)=VINT(113)+GZPF
503 VINT(114)=VINT(114)+ZZF
504 VINT(115)=VINT(115)+ZZPF
505 VINT(116)=VINT(116)+ZPZPF
508 IF(MSTP(44).EQ.1) THEN
509 C...Only gamma* production included
515 ELSEIF(MSTP(44).EQ.2) THEN
516 C...Only Z0 production included
522 ELSEIF(MSTP(44).EQ.3) THEN
523 C...Only Z'0 production included
529 ELSEIF(MSTP(44).EQ.4) THEN
530 C...Only gamma*/Z0 production included
534 ELSEIF(MSTP(44).EQ.5) THEN
535 C...Only gamma*/Z'0 production included
539 ELSEIF(MSTP(44).EQ.6) THEN
540 C...Only Z0/Z'0 production included
546 ELSEIF(KFLA.EQ.37) THEN
548 DO 190 I=1,MDCY(37,3)
550 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
551 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
552 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 190
555 WDTP(I)=3.*((RM1*PARU(121)+RM2/PARU(121))*
556 & (1.-RM1-RM2)-4.*RM1*RM2)*
557 & SQRT(MAX(0.,(1.-RM1-RM2)**2-4.*RM1*RM2))*RADC
560 C...H+/- -> l+/- + nu
561 WDTP(I)=((RM1*PARU(121)+RM2/PARU(121))*
562 & (1.-RM1-RM2)-4.*RM1*RM2)*
563 & SQRT(MAX(0.,(1.-RM1-RM2)**2-4.*RM1*RM2))
566 WDTP(0)=WDTP(0)+WDTP(I)
567 IF(MDME(IDC,1).GT.0) THEN
568 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
569 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
570 WDTE(I,0)=WDTE(I,MDME(IDC,1))
571 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)
575 ELSEIF(KFLA.EQ.40) THEN
577 DO 200 I=1,MDCY(40,3)
579 RM1=(PMAS(IABS(KFDP(IDC,1)),1)/RMAS)**2
580 RM2=(PMAS(IABS(KFDP(IDC,2)),1)/RMAS)**2
581 IF(SQRT(RM1)+SQRT(RM2).GT.1..OR.MDME(IDC,1).LT.0) GOTO 200
591 WDTP(0)=WDTP(0)+WDTP(I)
592 IF(MDME(IDC,1).GT.0) THEN
593 WDTE(I,MDME(IDC,1))=WDTP(I)*WID2
594 WDTE(0,MDME(IDC,1))=WDTE(0,MDME(IDC,1))+WDTE(I,MDME(IDC,1))
595 WDTE(I,0)=WDTE(I,MDME(IDC,1))
596 WDTE(0,0)=WDTE(0,0)+WDTE(I,0)