+ else if (velnanoconst.ne.0 ) then
+ if (refstr) then
+ call rms_nac_nnc(rms,frac,frac_nn,co,.false.)
+ write (line1,'(i10,f15.2,3f12.3,f7.2,2f6.3,f12.3,f10.1,2f8.2, &
+ f9.3,i5,$)') &
+ itime,totT,EK,potE+potEcomp(51),totE+potEcomp(51),&
+ rms,frac,frac_nn,kinetic_T,t_bath,gyrate(),&
+ vecsim,vectrue,me
+ format1="a133"
+!C print *,"CHUJOWO"
+ else
+!C print *,'A CHUJ',potEcomp(23)
+ write (line1,'(i10,f15.2,8f12.3,i5,$)')&
+ itime,totT,EK,potE+potEcomp(51),totE+potEcomp(51), &
+ kinetic_T,t_bath,gyrate(),&
+ vecsim,vectrue,me
+ format1="a114"
+ endif