double precision IP,ISC(ntyp+1),mp,msc(ntyp+1), 
     & d_t_work(MAXRES6),d_t_work_new(MAXRES6),d_t(3,0:MAXRES2),
     & d_t_new(3,0:MAXRES2),vtot(maxres2),
     & d_af_work(MAXRES6),d_as_work(MAXRES6),
     & d_t_old(3,0:MAXRES2),d_a_old(3,0:MAXRES2),d_a_short(3,0:MAXRES2),
     & Gmat(MAXRES2,MAXRES2),Ginv(MAXRES2,MAXRES2),A(MAXRES2,MAXRES2),
     & d_a(3,0:MAXRES2),d_a_work(6*MAXRES),
     & Gsqrp(MAXRES2,MAXRES2),Gsqrm(MAXRES2,MAXRES2),
     & Gvec(maxres2,maxres2),Geigen(maxres2)
      integer dimen,dimen1,dimen3
      common /inertia/ IP,ISC,mp,MSC
      common /lagrange/ d_t,d_t_old,d_t_new,d_t_work,d_t_work_new,d_a,
     & d_a_old,d_a_work,d_af_work,d_as_work,d_a_short,
     & A,Ginv,Gmat,Gvec,Geigen,Gsqrp,Gsqrm,
     & vtot,dimen,dimen1,dimen3