Either Forcite, GULP or Discover molecular dynamics tasks seems not have any option for cooling/heating cycles, Could anyone help me how to do this? also, how about change ensembles during one MD simulation at various stages? I know LAMMPS can easily do this, but I am trying to switch to MS for its wonderful interface. Would this be achievable using a perl script? If so, could anyone give me some example script? Thanks  very much!

HY