set PROJECT enk
begin
open enkrms a
for {$i = 1} {$i < 26} {incr i} {
readFile coordinate filename = struc11.cor
subset define Subset "*:*:min" "C*"
molGeom get coords prodrug "*:*:min"
for {$k = 1} {$k < 5} {incr k} {
readFile coordinate filename = struc$i$k.cor
set energy [readFile coordinate filename = struc$i$k.cor]
subset define Subset "*:*:all$k" "C*"
molGeom get coords superimpose "*:*:all$k" 
geometry RMS rms $prodrug $superimpose  
set rms [vector RMS]
puts file4 "$rms"
unset rms
unset RMS
unset superimpose
unset energy
}
unset prodrug 
}