mfold-quik

来自「mfold」· 代码 · 共 23 行

#!/bin/sh
# This shell script folds all RNA or DNA sequences in a single file
# and creates simple output files.

# Abort subroutine
abort(){
    \rm -f mfold.log  fort.*
    if [ $# -gt 0 ] ; then
	    echo -e "$1"
    fi
    echo "Job Aborted"
    exit 1
	}

DATDIR=${MFOLDDAT:-$MFOLD/dat}
if [ $DATDIR = "/dat" ] ; then
    abort "The environment variable, MFOLDDAT, is not set."
fi

if [ $# = 0 ]; then
  echo -e "Usage is\nmfold SEQ='file_name' with optional parameters:
   [ NA=RNA (default) or DNA ] [ LC=sequence type (default = linear) ]
   [ T=temperature (default = 37