#############################################################
## Hacer los respectivos Restrain (CA)                     ##
#############################################################

# Put harmonic restraints on the CA atoms for heating and unit cell
# equilibration.  Use the coordinates from the last simulation frame.

set all [atomselect top all]
$all set beta 0
set sel [atomselect top "resname C and z<-8 and z>-15"]
$sel set beta 1
$all writepdb rest_bPh.pdb
exit
