Hi
Below is script for interaction energy calculation. I call this script from bash script.
Code:
!================================================================
! Crystal definition and cutoffs
!================================================================
set ci  10.   ! inner switching cutoff
set rc  12.   ! cutoff
set ctl 14.   ! list-cutoff
set xo  14.   ! cutoff for crystal

!================================================================
!non bonded cut-offs
!================================================================
 NBOND CUTNB @ctl CUTIM @ctl CTOFnb @rc CTONnb @ci -
       atom vatom vdistance                        -
       VSWITCH  SHIFT  CDIE eps 1.0  e14fac 1.0  wmin 1.5

set res1 1

LABEL innerloop
!============================================================
! Write data to file
!============================================================
OPEN UNIT 28 WRITE FORM NAME ./solv_inter_ene/@i/inter_ene_solv_@res1.@i.dat

!============================================================
! Read TRAJ
!============================================================
OPEN READ UNIT 51 FILE NAME ./recen/@i/unfold_1le0.joball.@i.recen.dcd
! OPEN READ UNIT 51 FILE NAME ./extract/unfold_1le0.joball.@i.extract.dcd
TRAJECTORY QUERY UNIT 51
traj iread 51 nskip 1 begin 0 stop ?NSTEP skip 

SET frame 0
LABEL TimeFrame
TRAJ READ

!============================================================
! Example inter Select resid <residue_number> END Select resid <residue_number> show end
! interaction each residue w.r.t solution
!============================================================
inter SELECT resid @res1 END SELECT .NOT. segid PRT END

!==================================================
! writing output
!==================================================
write title unit 28
* @frame ?ener ?vdw ?elec
*

increment frame
!if frame .LT. 15 GOTO TimeFrame         !!!! small run to check that script runs correctly
if frame .LT. ?NFILE GOTO TimeFrame

close unit 28

increment res1
if res1 .LT. 13 GOTO innerloop
STOP


This calculation takes 20 mins for 1 Amino Acid with all solvent. I have 12 amino acids so it will take 240 mins for 1 trajectory. I have 96 trajectories.

Is there any better way (faster way) of doing this analysis. am i doing something wrong.