| 1 |
C |
| 2 |
C Written by Mirko Boezio and Emiliano Mocchiutti |
| 3 |
C |
| 4 |
C * Version: 2.17.1 * |
| 5 |
C |
| 6 |
C----------------------------------------------------- |
| 7 |
SUBROUTINE RIEMPI(m,lunga,lleng,lung,pari,vect,vecta) |
| 8 |
C----------------------------------------------------- |
| 9 |
|
| 10 |
IMPLICIT NONE |
| 11 |
C |
| 12 |
integer m, lunga ,lleng, lung, i, bi, bit,pari |
| 13 |
INTEGER*2 VECT(30000) |
| 14 |
INTEGER*1 VECTA(lung) |
| 15 |
C |
| 16 |
INTEGER iev,iev2 |
| 17 |
C |
| 18 |
INTEGER NPLA, NCHA, LENSEV |
| 19 |
PARAMETER (NPLA=22,NCHA=96,LENSEV=NPLA*NCHA) |
| 20 |
INTEGER merror(4) |
| 21 |
integer*2 e2(4) |
| 22 |
INTEGER contr |
| 23 |
REAL DEXY(2,NPLA,NCHA), dexyc(2,npla,ncha), base(2,npla,6) |
| 24 |
|
| 25 |
|
| 26 |
real calselftrig(4,7), calIItrig(4), calstripshit(4), |
| 27 |
& calDSPtaberr(4), calevnum(4) |
| 28 |
|
| 29 |
|
| 30 |
REAL calped(4,11,96), calgood(4,11,96), calthr(4,11,6) |
| 31 |
REAL calrms(4,11,96), calbase(4,11,6), calvar(4,11,6) |
| 32 |
REAL calpuls(4,11,96) |
| 33 |
real perror(4) |
| 34 |
integer stwerr(4),dump |
| 35 |
|
| 36 |
COMMON / evento / IEV, stwerr,perror, |
| 37 |
& dexy,dexyc,base, |
| 38 |
& calselftrig,calIItrig, |
| 39 |
& calstripshit,calDSPtaberr,calevnum |
| 40 |
|
| 41 |
save / evento / |
| 42 |
|
| 43 |
COMMON / calib / IEV2, calped, calgood, calthr, calrms, |
| 44 |
& calbase, |
| 45 |
& calvar, |
| 46 |
& calpuls |
| 47 |
|
| 48 |
save / calib / |
| 49 |
c |
| 50 |
COMMON / VARIE / dump, CONTR |
| 51 |
SAVE / VARIE / |
| 52 |
C |
| 53 |
do i = 1, (lunga - lleng) |
| 54 |
vect(i) = 0 |
| 55 |
if (m.le.lung) then |
| 56 |
do bit=0, 7 |
| 57 |
if (pari.eq.1.and.m.eq.lung) then |
| 58 |
bi = ibits(vecta(m+1),bit,1) |
| 59 |
if (bi.eq.1) vect(i) = ibset(vect(i),bit) |
| 60 |
endif |
| 61 |
if (m.lt.lung) then |
| 62 |
bi = ibits(vecta(m+1),bit,1) |
| 63 |
if (bi.eq.1) vect(i) = ibset(vect(i),bit) |
| 64 |
endif |
| 65 |
bi = ibits(vecta(m),bit,1) |
| 66 |
if (bi.eq.1) vect(i) = ibset(vect(i),bit+8) |
| 67 |
enddo |
| 68 |
endif |
| 69 |
m = m + 2 |
| 70 |
enddo |
| 71 |
c |
| 72 |
17 FORMAT(2X,'Elemento:',2X,I4,2X,' word:',2X,Z8) |
| 73 |
return |
| 74 |
end |