+ else if (dist1cut) then
+ call rms_template(1,rmsd1)
+ call rms_template(2,rmsd2)
+ if (refstr) then
+ call rms_nac_nnc(rms,frac,frac_nn,co,.false.)
+ write (line1,
+ & '(i10,f15.2,3f12.3,f7.2,4f6.3,3f12.3,i5,2f7.2,$)')
+ & itime,totT,EK,potE+potEcomp(27),totE+potEcomp(27),
+ & rms,frac,frac_nn,co,amax,kinetic_T,t_bath,gyrate(),me,
+ & rmsd1,rmsd2
+ format1="a147"
+ else
+ write (line1,'(i10,f15.2,7f12.3,i5,2f7.2,$)')
+ & itime,totT,EK,potE+potEcomp(27),totE+potEcomp(27),
+ & amax,kinetic_T,t_bath,gyrate(),me,rmsd1,rmsd2
+ format1="a128"
+ endif
+