1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
5 C Calculates the partial decay widths of
6 C the Higgs bosons into sfermions.
7 C calculated by X. Tata
10 C 10/23/93: modified by H. Baer, 10/8/96
11 C Intra-flavor sfermion mixing is neglected
12 C for all flavors EXCEPT for stops, sbottoms and staus.
16 C It is assumed that the A-terms are real.
17 C In addition, all coefficients of the sfermion
18 C trilinear terms from the superpotential
19 C EXCEPT the stop (AAT), sbottom (AAB) and stau (AAL)
20 C coefficients are set to zero.
22 C ===> Code for the general case removing all these
23 C artificial restrictions is present below.
24 C The preceeding restrictions are specified
25 C by giving special values to some variables
26 C This is discussed in two sections beginning
27 C with the symbols (*@&*) in the code below.
29 C-----------------------------------------------------------------------
30 #if defined(CERNLIB_IMPNONE)
33 #include "isajet/sslun.inc"
34 #include "isajet/sssm.inc"
35 #include "isajet/sspar.inc"
36 #include "isajet/sstype.inc"
39 REAL SR2,PI,GG,TW2,BETA,DSA,DCA,DSB,DCB,MH
40 REAL EP,TANB,COTB,ATERM,MSFMIX,THETSF,SIN2B
41 REAL TEMP,TEMP1,TEMP2,YA1,YA2
42 REAL SINA,COSA,SINA2,COSA2,M1,M2,M12,LAMB
43 REAL SINAU,COSAU,SINAD,COSAD
44 REAL A11,A22,A12,B11,B22,B12,C11,C12,C21,C22
46 REAL DWSFL,DWSFH,DWSFP,DWSFC,SSXLAM
47 REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
48 DOUBLE PRECISION SSMQCD
49 DIMENSION ATERM(12),MSFMIX(12,2),THETSF(12)
50 DIMENSION ASQ(10,3),BSQ(9),CSQ(6,4)
51 DIMENSION DWSF(12,4),DWSFL(12,4),DWSFH(12,4)
52 DIMENSION DWSFP(12,4),DWSFC(6,4)
53 INTEGER II,IJ,JJ,IC,IJU,IJD,NUMH
58 TW2=SN2THW/(1.0-SN2THW)
59 GG=SQRT(4.0*PI*ALFAEM/SN2THW)
72 C (all A-terms are assumed to be real)
73 C The A-terms are loaded into the array ATERM(12)
74 C in the following way:
75 C ATERM(1)=selectron A-term
76 C ATERM(2)=smuon A-term
77 C ATERM(3)=stau A-term
78 C ATERM(4)=up squark A-term
79 C ATERM(5)=charm squark A-term
80 C ATERM(6)=down squark A-term
81 C ATERM(7)=strange squark A-term
82 C ATERM(8)=sbottom A-term
83 C ATERM(9)=stop A-term
84 C ATERM(10)=selectronic sneutrino A-term
85 C ATERM(11)=smuonic sneutrino A-term
86 C ATERM(12)=stauonic sneutrino A-term
98 C Set mixing parameters.
99 C The intra-flavor-mixed sfermion masses are loaded into
100 C the array MSFMIX(12,2) where (#,1) is the lighter
101 C mixed sfermion mass of a given flavor and (#,2) is the
102 C heavier sfermion mass. The sfermionic mixing angles are
103 C loaded into the array THETSF(12). The identities of the
104 C elements of these arrays are given below:
105 C MSFMIX(1,*)=mixed selectron masses
106 C THETSF(1)=selectron mixing angle
107 C MSFMIX(2,*)=mixed smuon masses
108 C THETSF(2)=smuon mixing angle
109 C MSFMIX(3,*)=mixed stau masses
110 C THETSF(3)=stau mixing angle
111 C MSFMIX(4,*)=mixed up squark masses
112 C THETSF(4)=up squark mixing angle
113 C MSFMIX(5,*)=mixed charm squark masses
114 C THETSF(5)=charm squark mixing angle
115 C MSFMIX(6,*)=mixed down squark masses
116 C THETSF(6)=down squark mixing angle
117 C MSFMIX(7,*)=mixed strange squark masses
118 C THETSF(7)=strange squark mixing angle
119 C MSFMIX(8,*)=mixed sbottom masses
120 C THETSF(8)=sbottom mixing angle
121 C MSFMIX(9,*)=mixed stop masses
122 C THETSF(9)=stop mixing angle
123 C For sneuterinos MSFMIX(#,2)=0.0, THETSF(#)=0.0 ; #=10-12
124 C Yukawa contributions from D-terms to the sneutrino masses
125 C are supposed to be added in here.
126 C MSFMIX(10,1)= selectronic sneutrino mass with D-terms
127 C MSFMIX(11,1)= smuonic sneutrino mass with D-terms
128 C MSFMIX(12,1)= stauonic sneutrino mass with D-terms
136 C (*@&*) 10/24/93 - Special conditions used ---
137 C set all mixing angles EXCEPT stop, sbottom, stau to zero.
138 C For all EXCEPT st, sb and stau, set mixed sfermion masses
139 C to bare sfermion masses:
140 C MSFMIX(#,1) = Left sfermion mass
141 C MSFMIX(#,2) = Right sfermion mass ; # = 1-8
143 C MSFMIX(9,1) = AMT1SS
144 C MSFMIX(9,2) = AMT2SS , etc.
146 C (The choice of which to call Left and which to call
147 C Right is based on the definition of the sfermion
148 C mixing angle theta_sf :
149 C sfermion_1 = cos(theta_sf) * sfermion_L
150 C - sin(theta_sf) * sfermion_R
151 C sfermion_2 = sin(theta_sf) * sfermion_L
152 C + cos(theta_sf) * sfermion_R
153 C Thus if we set theta_sf = 0, then
154 C sfermion_1 = sfermion_L
155 C and sfermion_2 = sfermion_R . )
188 ELSE IF(NUMH.EQ.2) THEN
190 ELSE IF(NUMH.EQ.3) THEN
193 ELSE IF(NUMH.EQ.4) THEN
197 ASMB=SUALFS(AMBT**2,.36,AMTP,3)
198 MBMB=AMBT*(1.-4*ASMB/3./PI)
199 MBQ=SSMQCD(DBLE(MBMB),DBLE(MH))
200 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
201 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
203 MTQ=SSMQCD(DBLE(MTMT),DBLE(MH))
206 C Scalar neutral Higgses --> sfermions
207 C partial decay widths
210 TEMP=GG*AMW*SIN(BETA-ALFAH)/2.0
213 ELSE IF(NUMH.EQ.2) THEN
214 TEMP=-GG*AMW*COS(BETA-ALFAH)/2.0
219 TEMP1=TEMP*(1.0-TW2/3.0)
220 TEMP2=GG*YA1/(AMW*DSB)
221 ASQ(4,1)=TEMP1-TEMP2*AMUP**2
222 ASQ(5,1)=TEMP1-TEMP2*AMCH**2
223 ASQ(9,1)=TEMP1-TEMP2*MTQ**2
225 TEMP1=-TEMP*(1.0+TW2/3.0)
226 TEMP2=GG*YA2/(AMW*DCB)
227 ASQ(6,1)=-TEMP1-TEMP2*AMDN**2
228 ASQ(7,1)=-TEMP1-TEMP2*AMST**2
229 ASQ(8,1)=-TEMP1-TEMP2*MBQ**2
231 ASQ(10,1)=TEMP*(1.0+TW2)
233 TEMP2=GG*YA2/(AMW*DCB)
234 ASQ(1,1)=TEMP1-TEMP2*AME**2
235 ASQ(2,1)=TEMP1-TEMP2*AMMU**2
236 ASQ(3,1)=TEMP1-TEMP2*AMTAU**2
238 TEMP1=4.0*TEMP*TW2/3.0
239 TEMP2=GG*YA1/(AMW*DSB)
240 ASQ(4,2)=TEMP1-TEMP2*AMUP**2
241 ASQ(5,2)=TEMP1-TEMP2*AMCH**2
242 ASQ(9,2)=TEMP1-TEMP2*MTQ**2
244 TEMP1=-2.0*TEMP*TW2/3.0
245 TEMP2=GG*YA2/(AMW*DCB)
246 ASQ(6,2)=TEMP1-TEMP2*AMDN**2
247 ASQ(7,2)=TEMP1-TEMP2*AMST**2
248 ASQ(8,2)=TEMP1-TEMP2*MBQ**2
252 TEMP2=GG*YA2/(AMW*DCB)
253 ASQ(1,2)=TEMP1-TEMP2*AME**2
254 ASQ(2,2)=TEMP1-TEMP2*AMMU**2
255 ASQ(3,2)=TEMP1-TEMP2*AMTAU**2
257 TEMP1=GG/(2.0*AMW*DSB)
258 ASQ(4,3)=(EP*YA2 + ATERM(4)*YA1)*TEMP1*AMUP
259 ASQ(5,3)=(EP*YA2 + ATERM(5)*YA1)*TEMP1*AMCH
260 ASQ(9,3)=(EP*YA2 + ATERM(9)*YA1)*TEMP1*MTQ
262 TEMP1=GG/(2.0*AMW*DCB)
263 ASQ(6,3)=(ATERM(6)*YA2 + EP*YA1)*TEMP1*AMDN
264 ASQ(7,3)=(ATERM(7)*YA2 + EP*YA1)*TEMP1*AMST
265 ASQ(8,3)=(ATERM(8)*YA2 + EP*YA1)*TEMP1*MBQ
268 ASQ(1,3)=(ATERM(1)*YA2 + EP*YA1)*TEMP1*AME
269 ASQ(2,3)=(ATERM(2)*YA2 + EP*YA1)*TEMP1*AMMU
270 ASQ(3,3)=(ATERM(3)*YA2 + EP*YA1)*TEMP1*AMTAU
275 TEMP1=1.0/(16.0*PI*MH**3)
277 TEMP1=3.0/(16.0*PI*MH**3)
287 A11=ASQ(IJ,1)*COSA2+ASQ(IJ,2)*SINA2
288 $ -2.0*ASQ(IJ,3)*SINA*COSA
289 LAMB=SSXLAM(MH**2,M1**2,M2**2)
290 DWSF(IJ,1)=TEMP1*SQRT(LAMB)*A11**2
291 ELSE IF(MH.LE.M12) THEN
299 A22=ASQ(IJ,1)*SINA2+ASQ(IJ,2)*COSA2
300 $ +2.0*ASQ(IJ,3)*SINA*COSA
301 LAMB=SSXLAM(MH**2,M1**2,M2**2)
302 DWSF(IJ,2)=TEMP1*SQRT(LAMB)*A22**2
303 ELSE IF(MH.LE.M12) THEN
311 A12=(ASQ(IJ,1)-ASQ(IJ,2))*SINA*COSA
312 $ +ASQ(IJ,3)*(COSA2-SINA2)
313 LAMB=SSXLAM(MH**2,M1**2,M2**2)
314 DWSF(IJ,3)=TEMP1*SQRT(LAMB)*A12**2
315 ELSE IF(MH.LE.M12) THEN
319 DWSF(IJ,4)=DWSF(IJ,3)
323 DWSFL(IJ,JJ)=DWSF(IJ,JJ)
325 ELSE IF(NUMH.EQ.2) THEN
327 DWSFH(IJ,JJ)=DWSF(IJ,JJ)
333 C Now take care of sneutrinos.
340 LAMB=SSXLAM(MH**2,M1**2,M2**2)
341 DWSF(IJ,1)=SQRT(LAMB)*(ASQ(10,1))**2
343 ELSE IF(MH.LE.M12) THEN
351 DWSFL(IJ,JJ)=DWSF(IJ,JJ)
353 ELSE IF(NUMH.EQ.2) THEN
355 DWSFH(IJ,JJ)=DWSF(IJ,JJ)
363 C Pseudocalar neutral Higgses --> sfermions
364 C partial decay widths
366 233 TEMP1=GG/(2.0*AMW)
367 BSQ(1)=TEMP1*AME*(EP-TANB*ATERM(1))
368 BSQ(2)=TEMP1*AMMU*(EP-TANB*ATERM(2))
369 BSQ(3)=TEMP1*AMTAU*(EP-TANB*ATERM(3))
370 BSQ(4)=TEMP1*AMUP*(EP-COTB*ATERM(4))
371 BSQ(5)=TEMP1*AMCH*(EP-COTB*ATERM(5))
372 BSQ(6)=TEMP1*AMDN*(EP-TANB*ATERM(6))
373 BSQ(7)=TEMP1*AMST*(EP-TANB*ATERM(7))
374 BSQ(8)=TEMP1*MBQ*(EP-TANB*ATERM(8))
375 BSQ(9)=TEMP1*MTQ*(EP-COTB*ATERM(9))
379 TEMP1=1.0/(16.0*PI*MH**3)
381 TEMP1=3.0/(16.0*PI*MH**3)
391 B11=-2.0*COSA*SINA*BSQ(IJ)
392 LAMB=SSXLAM(MH**2,M1**2,M2**2)
393 DWSFP(IJ,1)=TEMP1*SQRT(LAMB)*B11**2
394 ELSE IF(MH.LE.M12) THEN
403 LAMB=SSXLAM(MH**2,M1**2,M2**2)
404 DWSFP(IJ,2)=TEMP1*SQRT(LAMB)*B22**2
405 ELSE IF(MH.LE.M12) THEN
412 B12=(COSA2-SINA2)*BSQ(IJ)
413 LAMB=SSXLAM(MH**2,M1**2,M2**2)
414 DWSFP(IJ,3)=TEMP1*SQRT(LAMB)*B12**2
415 ELSE IF(MH.LE.M12) THEN
418 DWSFP(IJ,4)=DWSFP(IJ,3)
427 C Charged Higgses --> sfermions
428 C partial decay widths
431 CSQ(1,1)=GG*(TEMP1+(TANB*AMDN**2 + COTB*AMUP**2)/AMW)/SR2
432 CSQ(2,1)=GG*(TEMP1+(TANB*AMST**2 + COTB*AMCH**2)/AMW)/SR2
433 CSQ(3,1)=GG*(TEMP1+(TANB*MBQ**2 + COTB*MTQ**2)/AMW)/SR2
434 CSQ(4,1)=GG*(TEMP1 + (TANB*AME**2)/AMW)/SR2
435 CSQ(5,1)=GG*(TEMP1 + (TANB*AMMU**2)/AMW)/SR2
436 CSQ(6,1)=GG*(TEMP1 + (TANB*AMTAU**2)/AMW)/SR2
438 TEMP1=GG*(COTB+TANB)/(SR2*AMW)
439 CSQ(1,2)=TEMP1*AMUP*AMDN
440 CSQ(2,2)=TEMP1*AMCH*AMST
441 CSQ(3,2)=TEMP1*MTQ*MBQ
447 CSQ(1,3)=TEMP1*AMUP*(EP-COTB*ATERM(4))
448 CSQ(2,3)=TEMP1*AMCH*(EP-COTB*ATERM(5))
449 CSQ(3,3)=TEMP1*MTQ*(EP-COTB*ATERM(9))
454 CSQ(1,4)=TEMP1* AMDN*(EP-TANB*ATERM(6))
455 CSQ(2,4)=TEMP1* AMST*(EP-TANB*ATERM(7))
456 CSQ(3,4)=TEMP1* MBQ*(EP-TANB*ATERM(8))
457 CSQ(4,4)=TEMP1* AME*(EP-TANB*ATERM(1))
458 CSQ(5,4)=TEMP1* AMMU*(EP-TANB*ATERM(2))
459 CSQ(6,4)=TEMP1* AMTAU*(EP-TANB*ATERM(3))
462 TEMP1=3.0/(16.0*PI*MH**3)
466 ELSE IF(IC.EQ.2) THEN
469 ELSE IF(IC.EQ.3) THEN
473 SINAU=SIN(THETSF(IJU))
474 COSAU=COS(THETSF(IJU))
475 SINAD=SIN(THETSF(IJD))
476 COSAD=COS(THETSF(IJD))
482 C11=COSAU*COSAD*CSQ(IC,1)
483 $ + SINAU*SINAD*CSQ(IC,2)
484 $ - SINAU*COSAD*CSQ(IC,3)
485 $ - COSAU*SINAD*CSQ(IC,4)
486 LAMB=SSXLAM(MH**2,M1**2,M2**2)
487 DWSFC(IC,1)=TEMP1*SQRT(LAMB)*C11**2
488 ELSE IF(MH.LE.M12) THEN
496 C12=COSAU*SINAD*CSQ(IC,1)
497 $ - SINAU*COSAD*CSQ(IC,2)
498 $ - SINAU*SINAD*CSQ(IC,3)
499 $ + COSAU*COSAD*CSQ(IC,4)
500 LAMB=SSXLAM(MH**2,M1**2,M2**2)
501 DWSFC(IC,2)=TEMP1*SQRT(LAMB)*C12**2
502 ELSE IF(MH.LE.M12) THEN
510 C21=SINAU*COSAD*CSQ(IC,1)
511 $ - COSAU*SINAD*CSQ(IC,2)
512 $ + COSAU*COSAD*CSQ(IC,3)
513 $ - SINAU*SINAD*CSQ(IC,4)
514 LAMB=SSXLAM(MH**2,M1**2,M2**2)
515 DWSFC(IC,3)=TEMP1*SQRT(LAMB)*C21**2
516 ELSE IF(MH.LE.M12) THEN
524 C22=SINAU*SINAD*CSQ(IC,1)
525 $ + COSAU*COSAD*CSQ(IC,2)
526 $ + COSAU*SINAD*CSQ(IC,3)
527 $ - SINAU*COSAD*CSQ(IC,4)
528 LAMB=SSXLAM(MH**2,M1**2,M2**2)
529 DWSFC(IC,4)=TEMP1*SQRT(LAMB)*C22**2
530 ELSE IF(MH.LE.M12) THEN
537 C Now calculate the sleptonic
538 C partial decay widths of the
542 TEMP1=1.0/(16.0*PI*MH**3)
546 ELSE IF(IC.EQ.5) THEN
549 ELSE IF(IC.EQ.6) THEN
553 SINAD=SIN(THETSF(IJD))
554 COSAD=COS(THETSF(IJD))
560 C11=COSAD*CSQ(IC,1)-SINAD*CSQ(IC,4)
561 LAMB=SSXLAM(MH**2,M1**2,M2**2)
562 DWSFC(IC,1)=TEMP1*SQRT(LAMB)*C11**2
563 ELSE IF(MH.LE.M12) THEN
571 C12=SINAD*CSQ(IC,1)+COSAD*CSQ(IC,4)
572 LAMB=SSXLAM(MH**2,M1**2,M2**2)
573 DWSFC(IC,2)=TEMP1*SQRT(LAMB)*C12**2
574 ELSE IF(MH.LE.M12) THEN
582 CALL SSSAVE(ISHL,DWSFL(1,1),ISEL,-ISEL,0,0,0)
583 CALL SSSAVE(ISHL,DWSFL(1,2),ISER,-ISER,0,0,0)
584 CALL SSSAVE(ISHL,DWSFL(2,1),ISMUL,-ISMUL,0,0,0)
585 CALL SSSAVE(ISHL,DWSFL(2,2),ISMUR,-ISMUR,0,0,0)
586 CALL SSSAVE(ISHL,DWSFL(3,1),ISTAU1,-ISTAU1,0,0,0)
587 CALL SSSAVE(ISHL,DWSFL(3,2),ISTAU2,-ISTAU2,0,0,0)
588 CALL SSSAVE(ISHL,DWSFL(3,3),ISTAU1,-ISTAU2,0,0,0)
589 CALL SSSAVE(ISHL,DWSFL(3,4),ISTAU2,-ISTAU1,0,0,0)
590 CALL SSSAVE(ISHL,DWSFL(4,1),ISUPL,-ISUPL,0,0,0)
591 CALL SSSAVE(ISHL,DWSFL(4,2),ISUPR,-ISUPR,0,0,0)
592 CALL SSSAVE(ISHL,DWSFL(5,1),ISCHL,-ISCHL,0,0,0)
593 CALL SSSAVE(ISHL,DWSFL(5,2),ISCHR,-ISCHR,0,0,0)
594 CALL SSSAVE(ISHL,DWSFL(6,1),ISDNL,-ISDNL,0,0,0)
595 CALL SSSAVE(ISHL,DWSFL(6,2),ISDNR,-ISDNR,0,0,0)
596 CALL SSSAVE(ISHL,DWSFL(7,1),ISSTL,-ISSTL,0,0,0)
597 CALL SSSAVE(ISHL,DWSFL(7,2),ISSTR,-ISSTR,0,0,0)
598 CALL SSSAVE(ISHL,DWSFL(8,1),ISBT1,-ISBT1,0,0,0)
599 CALL SSSAVE(ISHL,DWSFL(8,2),ISBT2,-ISBT2,0,0,0)
600 CALL SSSAVE(ISHL,DWSFL(8,3),ISBT1,-ISBT2,0,0,0)
601 CALL SSSAVE(ISHL,DWSFL(8,4),ISBT2,-ISBT1,0,0,0)
602 CALL SSSAVE(ISHL,DWSFL(9,1),ISTP1,-ISTP1,0,0,0)
603 CALL SSSAVE(ISHL,DWSFL(9,2),ISTP2,-ISTP2,0,0,0)
604 CALL SSSAVE(ISHL,DWSFL(9,3),ISTP1,-ISTP2,0,0,0)
605 CALL SSSAVE(ISHL,DWSFL(9,4),ISTP2,-ISTP1,0,0,0)
606 CALL SSSAVE(ISHL,DWSFL(10,1),ISNEL,-ISNEL,0,0,0)
607 CALL SSSAVE(ISHL,DWSFL(11,1),ISNML,-ISNML,0,0,0)
608 CALL SSSAVE(ISHL,DWSFL(12,1),ISNTL,-ISNTL,0,0,0)
610 CALL SSSAVE(ISHH,DWSFH(1,1),ISEL,-ISEL,0,0,0)
611 CALL SSSAVE(ISHH,DWSFH(1,2),ISER,-ISER,0,0,0)
612 CALL SSSAVE(ISHH,DWSFH(2,1),ISMUL,-ISMUL,0,0,0)
613 CALL SSSAVE(ISHH,DWSFH(2,2),ISMUR,-ISMUR,0,0,0)
614 CALL SSSAVE(ISHH,DWSFH(3,1),ISTAU1,-ISTAU1,0,0,0)
615 CALL SSSAVE(ISHH,DWSFH(3,2),ISTAU2,-ISTAU2,0,0,0)
616 CALL SSSAVE(ISHH,DWSFH(3,3),ISTAU1,-ISTAU2,0,0,0)
617 CALL SSSAVE(ISHH,DWSFH(3,4),ISTAU2,-ISTAU1,0,0,0)
618 CALL SSSAVE(ISHH,DWSFH(4,1),ISUPL,-ISUPL,0,0,0)
619 CALL SSSAVE(ISHH,DWSFH(4,2),ISUPR,-ISUPR,0,0,0)
620 CALL SSSAVE(ISHH,DWSFH(5,1),ISCHL,-ISCHL,0,0,0)
621 CALL SSSAVE(ISHH,DWSFH(5,2),ISCHR,-ISCHR,0,0,0)
622 CALL SSSAVE(ISHH,DWSFH(6,1),ISDNL,-ISDNL,0,0,0)
623 CALL SSSAVE(ISHH,DWSFH(6,2),ISDNR,-ISDNR,0,0,0)
624 CALL SSSAVE(ISHH,DWSFH(7,1),ISSTL,-ISSTL,0,0,0)
625 CALL SSSAVE(ISHH,DWSFH(7,2),ISSTR,-ISSTR,0,0,0)
626 CALL SSSAVE(ISHH,DWSFH(8,1),ISBT1,-ISBT1,0,0,0)
627 CALL SSSAVE(ISHH,DWSFH(8,2),ISBT2,-ISBT2,0,0,0)
628 CALL SSSAVE(ISHH,DWSFH(8,3),ISBT1,-ISBT2,0,0,0)
629 CALL SSSAVE(ISHH,DWSFH(8,4),ISBT2,-ISBT1,0,0,0)
630 CALL SSSAVE(ISHH,DWSFH(9,1),ISTP1,-ISTP1,0,0,0)
631 CALL SSSAVE(ISHH,DWSFH(9,2),ISTP2,-ISTP2,0,0,0)
632 CALL SSSAVE(ISHH,DWSFH(9,3),ISTP1,-ISTP2,0,0,0)
633 CALL SSSAVE(ISHH,DWSFH(9,4),ISTP2,-ISTP1,0,0,0)
634 CALL SSSAVE(ISHH,DWSFH(10,1),ISNEL,-ISNEL,0,0,0)
635 CALL SSSAVE(ISHH,DWSFH(11,1),ISNML,-ISNML,0,0,0)
636 CALL SSSAVE(ISHH,DWSFH(12,1),ISNTL,-ISNTL,0,0,0)
638 CALL SSSAVE(ISHA,DWSFP(1,3),ISEL,-ISER,0,0,0)
639 CALL SSSAVE(ISHA,DWSFP(1,4),ISER,-ISEL,0,0,0)
640 CALL SSSAVE(ISHA,DWSFP(2,3),ISMUL,-ISMUR,0,0,0)
641 CALL SSSAVE(ISHA,DWSFP(2,4),ISMUR,-ISMUL,0,0,0)
642 CALL SSSAVE(ISHA,DWSFP(3,1),ISTAU1,-ISTAU1,0,0,0)
643 CALL SSSAVE(ISHA,DWSFP(3,2),ISTAU2,-ISTAU2,0,0,0)
644 CALL SSSAVE(ISHA,DWSFP(3,3),ISTAU1,-ISTAU2,0,0,0)
645 CALL SSSAVE(ISHA,DWSFP(3,4),ISTAU2,-ISTAU1,0,0,0)
646 CALL SSSAVE(ISHA,DWSFP(4,3),ISUPL,-ISUPR,0,0,0)
647 CALL SSSAVE(ISHA,DWSFP(4,4),ISUPR,-ISUPL,0,0,0)
648 CALL SSSAVE(ISHA,DWSFP(5,3),ISCHL,-ISCHR,0,0,0)
649 CALL SSSAVE(ISHA,DWSFP(5,4),ISCHR,-ISCHL,0,0,0)
650 CALL SSSAVE(ISHA,DWSFP(6,3),ISDNL,-ISDNR,0,0,0)
651 CALL SSSAVE(ISHA,DWSFP(6,4),ISDNR,-ISDNL,0,0,0)
652 CALL SSSAVE(ISHA,DWSFP(7,3),ISSTL,-ISSTR,0,0,0)
653 CALL SSSAVE(ISHA,DWSFP(7,4),ISSTR,-ISSTL,0,0,0)
654 CALL SSSAVE(ISHA,DWSFP(8,1),ISBT1,-ISBT1,0,0,0)
655 CALL SSSAVE(ISHA,DWSFP(8,2),ISBT2,-ISBT2,0,0,0)
656 CALL SSSAVE(ISHA,DWSFP(8,3),ISBT1,-ISBT2,0,0,0)
657 CALL SSSAVE(ISHA,DWSFP(8,4),ISBT2,-ISBT1,0,0,0)
658 CALL SSSAVE(ISHA,DWSFP(9,1),ISTP1,-ISTP1,0,0,0)
659 CALL SSSAVE(ISHA,DWSFP(9,2),ISTP2,-ISTP2,0,0,0)
660 CALL SSSAVE(ISHA,DWSFP(9,3),ISTP1,-ISTP2,0,0,0)
661 CALL SSSAVE(ISHA,DWSFP(9,4),ISTP2,-ISTP1,0,0,0)
663 CALL SSSAVE(ISHC,DWSFC(1,1),ISUPL,-ISDNL,0,0,0)
664 CALL SSSAVE(ISHC,DWSFC(1,2),ISUPR,-ISDNR,0,0,0)
665 CALL SSSAVE(ISHC,DWSFC(2,1),ISCHL,-ISSTL,0,0,0)
666 CALL SSSAVE(ISHC,DWSFC(2,2),ISCHR,-ISSTR,0,0,0)
667 CALL SSSAVE(ISHC,DWSFC(3,1),ISTP1,-ISBT1,0,0,0)
668 CALL SSSAVE(ISHC,DWSFC(3,2),ISTP1,-ISBT2,0,0,0)
669 CALL SSSAVE(ISHC,DWSFC(3,3),ISTP2,-ISBT1,0,0,0)
670 CALL SSSAVE(ISHC,DWSFC(3,4),ISTP2,-ISBT2,0,0,0)
671 CALL SSSAVE(ISHC,DWSFC(4,1),-ISEL,ISNEL,0,0,0)
672 CALL SSSAVE(ISHC,DWSFC(5,1),-ISMUL,ISNML,0,0,0)
673 CALL SSSAVE(ISHC,DWSFC(6,1),-ISTAU1,ISNTL,0,0,0)
674 CALL SSSAVE(ISHC,DWSFC(6,2),-ISTAU2,ISNTL,0,0,0)