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(3,maxres2),d_as(3,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),
     & d_a(3,0:MAXRES2),d_a_work(6*MAXRES),
     & DM(MAXRES2),DU1(MAXRES2),DU2(MAXRES2),DMorig(MAXRES2),
     & DU1orig(MAXRES2),DU2orig(MAXRES2)
      integer dimen,dimen1,dimen3,dimenp,dimen_chain,iposd_chain
      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,d_as,d_af_work,d_as_work,d_a_short,
     & DM,DU1,DU2,DMorig,DU1orig,DU2orig,
     & vtot,dimen,dimen1,dimen3,dimenp,dimen_chain(maxchain),
     & iposd_chain(maxchain)