/[PAMELA software]/yoda/techmodel/forroutines/calorimeter/calpulse.for
ViewVC logotype

Diff of /yoda/techmodel/forroutines/calorimeter/calpulse.for

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4.0 by kusanagi, Sun Mar 6 04:33:02 2005 UTC revision 4.1 by kusanagi, Sat May 28 08:40:40 2005 UTC
# Line 1  Line 1 
1  C  C
2  C  Written by Mirko Boezio and Emiliano Mocchiutti  C  Written by Mirko Boezio and Emiliano Mocchiutti
3  C  C
4  C     * Version: 3.1.3 *  C     * Version: 3.1.4 *
5    C
6    C     3.1.3 - 3.1.4: (2005-05-26) the decoding is wrong!!! we are reading all strip in one plane per time and not one strip over all planes per time!!!!
7  C  C
8  C     3.1.2 - 3.1.3: (2005-02-23) added some printout in debugging mode.  C     3.1.2 - 3.1.3: (2005-02-23) added some printout in debugging mode.
9  C  C
# Line 210  c     if we run out of vector give an er Line 212  c     if we run out of vector give an er
212  c      c    
213           if (ic.gt.(lung-1)) then           if (ic.gt.(lung-1)) then
214              me = 1              me = 1
215    c            print *,'1 me = 1'
216              call clearall              call clearall
217              do i = 1, 4              do i = 1, 4
218                 error(i) = 129                             error(i) = 129            
# Line 245  C Line 248  C
248           enddo           enddo
249           if (check.ne.vect(length)) then           if (check.ne.vect(length)) then
250              ERROR(contr) = 132              ERROR(contr) = 132
251    c            print *,'2 me = 1'
252              me = 1              me = 1
253              goto 200              goto 200
254           endif                   endif        
# Line 271  C Line 275  C
275  c  c
276   49      continue   49      continue
277  c  c
278           do j = 1,96           do i = 1,11
279              do i = 1,11              do j = 1,96
280                 icb = icb + 1                 icb = icb + 1
281                 if (k.eq.1) then                 if (k.eq.1) then
282                    calpuls(k,i,97-j) = vect(icb)                    calpuls(k,i,97-j) = vect(icb)
# Line 281  c Line 285  c
285                 endif                 endif
286             enddo             enddo
287           enddo           enddo
288    c         print *,'3 me = 0'
289           me = 0           me = 0
290  c  c
291   50      continue   50      continue
Legend:
Removed from v.4.0  
changed lines
  Added in v.4.1
  ViewVC Help
Powered by ViewVC 1.1.23