+c write (iout,*) "Calling rmsnat"
+ rms_closest(i) = rmsnat(i)
+ call TMscore_sub(rmsd,gdt_ts_closest(i),gdt_ha_closest(i),
+ & tmscore_closest(i),cfname,.true.)
+c write (iout,*) "Family",i," rmsd",rmsd,"gdt_ts",
+c & gdt_ts_closest(i)," gdt_ha",gdt_ha_closest(i),
+c & "tmscore",tmscore_closest(i)
+ call center
+ call pdbout(totfree_gr(i)/beta_h(ib),rms_closest(i),titel)