*  ========================================= 
*  JOB: 
*

! { parameters passed }
set skip 100  ! skip size for the energy calculation
set i 0  ! counter
$SET$

! get parameters
!par and top files
$LIB$

! { psf }
open unit 1 read formatted name @path/scr/@pdb.dyna_solv.psf
read psf card unit 1
close unit 1

! { coordinates }
OPEN UNIT 1 READ FORMatted NAME @path/scr/@pdb.dyna_solv.crd
READ COORdinate CARD  UNIT 1
CLOSE UNIT 1

stream @path/inp/all-traj.dat

trajectory iread 21 nread @count skip @skip


LABEL conditions 

incr i by 1

traj read

$IF$
if @i LT @{conf_n} GOTO conditions
if @i EQ @{conf_n} GOTO stop


LABEL writing

NBONDS  cutnb 13.5 ctonnb 8.5 ctofnb 12.5 -
vswitch shift cdie eps 1 !-
SKIPe  ALL EXCLude ELECtrostatic VDW
inte sele segid @fc end sele segid @sc end

! this is also the total int. energy, for consistency check
set a = ?ENER
set b = ?VDW
set c = ?ELEC

open write formatted unit 10 name @path/scr/@pdb.conf_@i.pdb 
write coor pdb sele segid @fc .or. segid @sc end unit 10 
* pdb for @pdb
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 10

! Let's keep a record of the radii
scalar radius store 1
scalar wmain recall 1
! Let's write the radii in a crd file
open unit 11 write formatted name @path/scr/@pdb.conf_@i-rad.crd.his
write coor card sele segid @fc .or. segid @sc end unit 11
* crd with radii in WMAIN for @pdb
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 11
open unit 12  write formatted name @path/scr/@{fc_n}.conf_@i-rad.crd.his
write coor card sele segid @fc end unit 12
* crd with radii in WMAIN for @{fc_n}
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 12
open unit 13 write formatted name @path/scr/@{sc_n}.conf_@i-rad.crd.his
write coor card sele segid @sc end unit 13
* crd with radii in WMAIN for @{sc_n}
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 13

! Let's keep a record of the charges
scalar charges store 2
scalar wmain recall 2 
! Let's write the radii in a crd file
open unit 14 write formatted name @path/scr/@pdb.conf_@i-cha.crd.his
write coor card sele segid @fc .or. segid @sc end unit 14
* crd with radii in WMAIN for @pdb
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 14
open unit 15 write formatted name @path/scr/@{fc_n}.conf_@i-cha.crd.his
write coor card sele segid @fc end unit 15
* crd with radii in WMAIN for @{fc_n}
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 15
open unit 16 write formatted name @path/scr/@{sc_n}.conf_@i-cha.crd.his
write coor card sele segid @sc end unit 16
* crd with radii in WMAIN for @{sc_n}
* conformation  = @i TOTAL energy = @a
* Van der Waals   =  @b      Electrostatic = @c
*
close unit 16

GOTO conditions                                                                                                                  


LABEL stop
stop