Difference between revisions of "Makerestart"
Jump to navigation
Jump to search
import>Csw34 |
import>Csw34 |
||
| Line 17: | Line 17: | ||
# VARIABLE ASSIGNMENT |
# VARIABLE ASSIGNMENT |
||
| + | |||
# Prints the line in the data file containing 'STEPS' and assigns $oldsteps to the second column |
# Prints the line in the data file containing 'STEPS' and assigns $oldsteps to the second column |
||
oldsteps=$(sed -n '/STEPS/p' data | awk '{print $2'}) |
oldsteps=$(sed -n '/STEPS/p' data | awk '{print $2'}) |
||
| Line 24: | Line 25: | ||
# SOME PRINTING |
# SOME PRINTING |
||
| + | |||
echo 'Setting up restart run in the current directory' |
echo 'Setting up restart run in the current directory' |
||
echo 'Number of steps completed='$oldsteps |
echo 'Number of steps completed='$oldsteps |
||
| Line 30: | Line 32: | ||
# CREATING THE NEW DATA FILE |
# CREATING THE NEW DATA FILE |
||
| + | |||
# NOTE that the last sed command ensures that if the data file has been used in a restart run before, |
# NOTE that the last sed command ensures that if the data file has been used in a restart run before, |
||
# the 'RESTORE' line is removed to ensure it only appears once in the new file! |
# the 'RESTORE' line is removed to ensure it only appears once in the new file! |
||
| Line 41: | Line 44: | ||
# FILE HANDLING |
# FILE HANDLING |
||
| + | |||
mkdir restart |
mkdir restart |
||
cp GMIN.dump input.crd restart/ |
cp GMIN.dump input.crd restart/ |
||
Revision as of 12:01, 9 May 2008
This script can be run in a directory containing the output from a CGMIN run to automatically setup a restart. The directory you run it in must include the following files:
- input.crd
- data
- GMIN.dump - if your run crashed and you only have GMIN.dump.save, copy it to GMIN.dump and then run the script
This script works for CHARMM with GMIN only but it is very simple to modify it for use with AMBER! All you need to change is the files that are copied across in the file handling section!
#!/bin/bash
########################################################################### # Created by Chris Whittleston (csw34@cam.ac.uk) on 9th May 2008 ## # Script to setup a restart run for GMIN ## # Usage: './makerestart <number of steps>' ## # e.g. if you want to do another 1000 steps then use './makerestart 1000'## ###########################################################################
# VARIABLE ASSIGNMENT
# Prints the line in the data file containing 'STEPS' and assigns $oldsteps to the second column
oldsteps=$(sed -n '/STEPS/p' data | awk '{print $2'})
steps=$oldsteps
# This increments the $steps variable by the number used in the first arguement
let "steps += $1"
# SOME PRINTING
echo 'Setting up restart run in the current directory' echo 'Number of steps completed='$oldsteps echo 'Number of additional steps requested='$1 echo 'Total number of steps after restart='$steps
# CREATING THE NEW DATA FILE
# NOTE that the last sed command ensures that if the data file has been used in a restart run before,
# the 'RESTORE' line is removed to ensure it only appears once in the new file!
# Prints the top of the data file up to (but not including) the line containing 'STEPS' to a new file, 'data.top'
awk 'NR==1,/STEPS/{print x};{x=$0}' data | sed -n '/RESTORE/!p' > data.top
# Adds back in the 'STEPS' keyword and also the 'RESTORE' command to specify a restart
echo 'STEPS '$steps' 1.0' >> data.top
echo 'RESTORE GMIN.dump' >> data.top
# Adds the bottom bit of the data file back on (NOT including the 'STEPS' line from the old one)
awk '/STEPS/,0' data | sed -n '2,$p' | sed -n '/RESTORE/!p' >> data.top
# FILE HANDLING
mkdir restart cp GMIN.dump input.crd restart/ sed -n '2,$p' data.top > restart/data rm data.top