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