#!/bin/bash
set -x
trap read debug
x3dna_utils cp_std ADNA
##############################################
##  ADNA from fiber
##############################################
rebuild -base_p adna.par adna.pdb
rebuild adna.par adna.alc
rebuild -atomic adna.par adna.pdb
find_pair adna.pdb stdout | analyze stdin
frame_mol -max -2,12 ref_frames.dat adna.alc adnab.alc
alc2img -al -s=20 adnab.alc adnab.ps
alc2img -r adnab.alc adnab.r3d
echo "by rotation x 90" | rotate_mol -r=stdin adnab.alc adnac.alc
alc2img -al -s=20 adnac.alc adnac.ps
rm Atomic*.pdb