--- DarthVader/TrackerLevel2/src/F77/tricircle.f	2006/11/21 17:12:59	1.2
+++ DarthVader/TrackerLevel2/src/F77/tricircle.f	2014/01/16 15:30:00	1.5
@@ -31,12 +31,16 @@
       integer npoints           !fit number of points
       real dep(npoints),indep(npoints) !dependent and independent variables
 
-      real angle(npoints)       !angle between the tangent line in the input points and 
+c      real angle(npoints)       !angle between the tangent line in the input points and 
+c                                ! the independent variable axis
+c      real residual(npoints)    !residuals
+c      real chi                  !sum of squared residuals
+c      real xc,zc,radius         !circle parameters   
+      double precision angle(npoints)       !angle between the tangent line in the input points and 
                                 ! the independent variable axis
-
-      real residual(npoints)    !residuals
-      real chi                  !sum of squared residuals
-      real xc,zc,radius         !circle parameters
+      double precision residual(npoints)    !residuals EM GCC4.7
+      double precision chi              !sum of squared residuals EM GCC4.7
+      double precision xc,zc,radius         !circle parameters   EM GCC4.7
       
       integer eflag             !error flag =1 if the procedure fails
 
@@ -51,8 +55,10 @@
       integer ifail             !=-1 if singular matrix error, =0 if not singular
       integer jfail             !=0 if determinant can be evaluated, =-1 if determinat is probably too small, =+1 if too large
 
-      parameter (big=1.e4)      !just a number greater than C(npoints,3)
-      double precision xxc(big),zzc(big),rrr(big) !centres and radii to be averaged
+      integer ibig ! EM GCC4.7
+c      parameter (ibig=1.e4)      !just a number greater than C(npoints,3)
+      parameter (ibig=10000)      !just a number greater than C(npoints,3)   EM GCC 4.7
+      double precision xxc(ibig),zzc(ibig),rrr(ibig) !centres and radii to be averaged
 
       double precision tmp1,tmp2,tmp(npoints)    !temp variables
 
@@ -174,20 +180,20 @@
             call DFACT(3,d,3,ir,ifail,detd,jfail)
             if(ifail.eq.-1) then
               if(DEBUG)then
-              print*,'tricircle: ERROR: singular matrix D:'
-              do i=1,3
-                print*,(d(i,j),j=1,3)
-              enddo
+                 print*,'tricircle: ERROR: singular matrix D:'
+                 do i=1,3
+                    print*,(d(i,j),j=1,3)
+                 enddo
               endif
               eflag=1
-            endif
-            if(jfail.eq.-1) then
+           endif
+           if(jfail.eq.-1) then
               if(DEBUG)then
-              print*
-     $             ,'tricircle: ERROR: matrix D: determinant too small?'
-              do i=1,3
-                print*,(d(i,j),j=1,3)
-              enddo
+                 print*
+     $      ,'tricircle: ERROR: matrix D: determinant too small?'
+                 do i=1,3
+                    print*,(d(i,j),j=1,3)
+                 enddo
               endif
               eflag=1
             elseif(jfail.eq.1) then
@@ -212,7 +218,7 @@
               enddo
               endif
               eflag=1
-            endif
+           endif
             if(jfail.eq.-1) then
               if(DEBUG)then
               print*