1 # Amber results for this point

2 #   NSTEP       ENERGY          RMS            GMAX         NAME    NUMBER

3 #      1      -3.8830E-03     3.9148E-01     1.0210E+00     C2          2

4 #

5 # BOND    =        0.0008  ANGLE   =        0.0008  DIHED      =        0.0048

6 # VDWAALS =        0.0000  EEL     =        0.0000  HBOND      =        0.0000

7 # 1-4 VDW =       -0.0211  1-4 EEL =        0.0109  RESTRAINT  =        0.0000

8

9 load_amber_coords inpcrd=ethane.inpcrd prmtop=ethane.prmtop coords=ethane.c

10

11 matrix e

12 energy energy=e coords=ethane.c theory=dl_poly  : [ list \

13                                                use_amber_prmtop=yes  \

14                                                amber_prmtop_file=ethane.prmtop \

15                                                scale14 = [ list [ expr 1 / 1.2 ] 0.5  ] \

16                                                save_dl_poly_files = yes \

17                                                list_option=full ]

18

19 set final_energy  [ get_matrix_element matrix=e indices= {0 0 } ]

20 delete_object e

21

22 # This test result is the ChemShell result at implementation

23 # Amber result is -3.8830E-03 kcal / mol which using the ChemShell

24 # conversion factor is -6.1879E-06

25

26 set etest  -0.00000619

27 if { abs( $final_energy - $etest) > 1.0e-8 } {

28    chemerr "Failed .. E= $final_energy, should be $etest"

29 } else {

30    return "OK"

31 }