/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/readallparam.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/readallparam.f

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

revision 1.2 by pam-fi, Tue May 30 16:30:37 2006 UTC revision 1.5 by pam-fi, Fri Aug 17 16:47:16 2007 UTC
# Line 149  c         if(C2F_DEBUG.eq.1)print * Line 149  c         if(C2F_DEBUG.eq.1)print *
149                
150        c2f_error = 0        c2f_error = 0
151                        
152  ******retrieve ANGULAR BINNING info  *     ======================================================
153    *     read ANGULAR BINNING info
154    *     ======================================================
155        fname_binning='binning.dat'        fname_binning='binning.dat'
156  c      if(C2F_DEBUG.eq.1)  c      if(C2F_DEBUG.eq.1)
157        if(VERBOSE)        if(VERBOSE)
# Line 188  c      print*,nangbin Line 190  c      print*,nangbin
190                
191    
192    
193    *     ======================================================
194    *     read eta2-eta3-eta4 parameters
195    *     ======================================================
196        do ieta=2,4               !loop on eta 2,3,4                do ieta=2,4               !loop on eta 2,3,4        
 ******retrieve correction parameters  
197   200     format(' Opening eta',i1,' files...')   200     format(' Opening eta',i1,' files...')
198  c         if(C2F_DEBUG.eq.1)write(*,200)ieta  c         if(C2F_DEBUG.eq.1)write(*,200)ieta
199           if(VERBOSE)write(*,200)ieta           if(VERBOSE)write(*,200)ieta
# Line 214  c                  if(C2F_DEBUG.eq.1) Line 218  c                  if(C2F_DEBUG.eq.1)
218       $                 ,c2f_path(1:c2f_pathlen)       $                 ,c2f_path(1:c2f_pathlen)
219       $                 //fname_param(1:LNBLNK(fname_param))       $                 //fname_param(1:LNBLNK(fname_param))
220                    c2f_error = 1                    c2f_error = 1
221                    return  c                  return
222                      goto 2001
223                 endif                 endif
224                 netaval=0                 netaval=0
225                 do ival=1,netavalmax                 do ival=1,netavalmax
# Line 248  c                        if(C2F_DEBUG.eq Line 253  c                        if(C2F_DEBUG.eq
253                    netaval=netaval+1                    netaval=netaval+1
254                 enddo                 enddo
255   2000          close(10)   2000          close(10)
256     2001          continue
257  *               print*,'... done'  *               print*,'... done'
258              enddo              enddo
259           enddo           enddo
260    
261        enddo                     !end loop on eta 2,3,4        enddo                     !end loop on eta 2,3,4
262    
263    *     ======================================================
264    *     read landi correction
265    *     ======================================================
266          if(VERBOSE)print*,' Opening landi-correction files...'
267    
268     301  format('lcorr-l',i1,'.dat')
269          do ilad=1,nladders_view
270            
271             write(fname_param,301)ilad
272             open(10,
273         $        FILE=c2f_path(1:c2f_pathlen)
274         $        //fname_param(1:LNBLNK(fname_param))
275         $        ,STATUS='OLD'
276         $        ,IOSTAT=iostat
277         $        )
278             if(iostat.ne.0)then
279                if(.true.)
280         $           print*,'READETAPARAM: ** Error opening file ** '
281         $           ,c2f_path(1:c2f_pathlen)
282         $           //fname_param(1:LNBLNK(fname_param))
283    c            c2f_error = 1 !tolto solo temporaneamente
284    c     return
285                goto 3001
286             endif
287            
288             do iang=1,nangbin
289                
290                read(10,*,IOSTAT=iostat)
291         $           ii,aal,aar,(fcorr(iv,ilad,iang),iv=1,nviews)
292                
293                write(*,*)
294         $           ii,aal,aar,(fcorr(iv,ilad,iang),iv=1,nviews)
295    
296                if(ii.ne.iang)print*,'parametri eta incasinati'
297                if(aal.ne.angl(iang))print*,'parametri eta incasinati'
298                if(aar.ne.angr(iang))print*,'parametri eta incasinati'
299                
300                if(iostat.ne.0)then
301                   if(.true.)
302         $              print*,'READETAPARAM: '
303         $              //'*** Error reading file *** '
304         $              ,c2f_path(1:c2f_pathlen)
305         $              //fname_param(1:LNBLNK(fname_param))
306                   c2f_error = 1
307                   goto 3000        
308                endif
309            
310             enddo                  !end loop on angular bins
311     3000    close(10)
312     3001    continue
313            
314          enddo                     ! end loop on ladders
315    
316        return        return
317        end        end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23