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

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

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

revision 1.37 by pam-fi, Fri Dec 5 08:26:47 2008 UTC revision 1.39 by pam-fi, Mon Jan 26 14:01:41 2009 UTC
# Line 376  c         if(DEBUG.EQ.1)iprint=1 Line 376  c         if(DEBUG.EQ.1)iprint=1
376           if(VERBOSE.EQ.1)iprint=1           if(VERBOSE.EQ.1)iprint=1
377           if(DEBUG.EQ.1)iprint=2           if(DEBUG.EQ.1)iprint=2
378           call mini2(jstep,ifail,iprint)           call mini2(jstep,ifail,iprint)
379           if(ifail.ne.0) then  cc         if(ifail.ne.0) then
380             if(ifail.ne.0.or.CHI2.ne.CHI2) then !new
381                if(CHI2.ne.CHI2)CHI2=-9999. !new
382              if(VERBOSE.EQ.1)then              if(VERBOSE.EQ.1)then
383                 print *,                 print *,
384       $              '*** MINIMIZATION FAILURE *** (after refinement) '       $              '*** MINIMIZATION FAILURE *** (after refinement) '
385       $              ,iev       $              ,iev
   
386              endif              endif
387           endif           endif
388                    
# Line 3011  c                                 chi2=- Line 3012  c                                 chi2=-
3012  *     **********************************************************  *     **********************************************************
3013    
3014                                if(chi2.le.0.)goto 666                                              if(chi2.le.0.)goto 666              
3015                                  if(chi2.ge.1.e08)goto 666 !OPTIMIZATION
3016                                  if(chi2.ne.chi2)goto 666  !OPTIMIZATION
3017    
3018  *     --------------------------  *     --------------------------
3019  *     STORE candidate TRACK INFO  *     STORE candidate TRACK INFO
# Line 4102  c                  call xyz_PAM(0,icl,is Line 4105  c                  call xyz_PAM(0,icl,is
4105                 alfayz2_av_nt(iyz)=alfayz2_av(iyz)                 alfayz2_av_nt(iyz)=alfayz2_av(iyz)
4106                 nnn=nnn+ptcloud_yz(iyz)                 nnn=nnn+ptcloud_yz(iyz)
4107              enddo              enddo
4108              do ipt=1,nnn              do ipt=1,min(ndblt_max_nt,nnn)
4109                 db_cloud_nt(ipt)=db_cloud(ipt)                 db_cloud_nt(ipt)=db_cloud(ipt)
4110               enddo               enddo
4111           endif           endif
# Line 4115  c                  call xyz_PAM(0,icl,is Line 4118  c                  call xyz_PAM(0,icl,is
4118                 alfaxz3_av_nt(ixz)=alfaxz3_av(ixz)                 alfaxz3_av_nt(ixz)=alfaxz3_av(ixz)
4119                 nnn=nnn+ptcloud_xz(ixz)                               nnn=nnn+ptcloud_xz(ixz)              
4120              enddo              enddo
4121              do ipt=1,nnn              do ipt=1,min(ntrpt_max_nt,nnn)
4122                tr_cloud_nt(ipt)=tr_cloud(ipt)                tr_cloud_nt(ipt)=tr_cloud(ipt)
4123               enddo               enddo
4124           endif           endif
Legend:
Removed from v.1.37  
changed lines
  Added in v.1.39
  ViewVC Help
Powered by ViewVC 1.1.23