Tue Oct 11 08:36:19 PDT 2011 #!/bin/csh -f # input: # $1 receptor coordinates file; # $2 grid size; # $3 # of iterations for the first run; # $4 # of iterations for the second run; # $5 # of iterations for the third run; # output: # delphi map: $1.phi # log file: delphi.out set RECCRG = rec+sph2.crg set GRID = 179 set n1 = 50 set n2 = 120 set n3 = 180 set DELPHI = $DOCK_BASE/private/delphi set VDW = vdw.siz #set VDW = vdw.hydren set CRG = amb2.crg.oxt set PRM = genric.prm # three step focussing calc. if (-e ARCDAT) /bin/rm ARCDAT cp $VDW fort.11 cp $CRG fort.12 cp $RECCRG fort.13 echo "gsize=$GRID, perfil=20, bndcon=2, linit=$n1\n" >! fort.10 cat $PRM >> fort.10 time $DELPHI if ($status != 0) exit(100) echo "gsize=$GRID, perfil=60, bndcon=3, linit=$n2\n" >! fort.10 cat $PRM >> fort.10 mv fort.14 fort.18 time $DELPHI if ($status != 0) exit(100) echo "gsize=$GRID, perfil=90, bndcon=3, linit=$n3\n" >! fort.10 cat $PRM >> fort.10 mv fort.14 fort.18 time $DELPHI if ($status != 0) exit(100) /bin/rm fort.1[0-3] fort.18 ARCDAT mv fort.14 $RECCRG:r.phi ___________________DelPhi II____________________ / | A program to solve the PB equation | | in 3D, using non-linear form, incorporating | | 2 dielectric regions, ionic strength, periodic | | and focussing boundary conditions, utilizing | | stripped optimum successive over-relaxation | | and an improved algorithm for mapping the | | Mol. Surface to the finite-Difference grid | __________________ ____________________/ DelPhi II program started on Tue Oct 11 2011 at 08:36:19 opening parameter file fort.10 atom radii read from file fort.11 ! default extended atom radii based loosely ! on mike connolly's MS program- note H's are 0 # of radius parameter records: 18 atomic charges read from file fort.12 # of charge parameter records: 201 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 582 !!! WARNING: ALA 1 has a net charge of -0.2480 !!! WARNING: ASN 2 has a net charge of -0.9360 !!! WARNING: LYS 3 has a net charge of -0.1810 !!! WARNING: GLY 4 has a net charge of -0.2480 !!! WARNING: PHE 5 has a net charge of -0.2480 !!! WARNING: LEV 6 has a net charge of -0.6480 !!! WARNING: LEU 7 has a net charge of -0.2480 !!! WARNING: VAL 9 has a net charge of -0.2480 !!! WARNING: ARG 10 has a net charge of -0.9880 !!! WARNING: LYS 11 has a net charge of -0.1810 !!! WARNING: GLY 12 has a net charge of -0.2480 !!! WARNING: ASN 13 has a net charge of -0.9360 !!! WARNING: LEU 14 has a net charge of -0.2480 !!! WARNING: ARG 16 has a net charge of -0.9880 !!! WARNING: CYS 18 has a net charge of -0.3830 !!! WARNING: LEU 19 has a net charge of -0.2480 !!! WARNING: CYS 23 has a net charge of -0.3830 !!! WARNING: SER 24 has a net charge of -0.5580 !!! WARNING: ARG 25 has a net charge of -0.9880 !!! WARNING: ARG 26 has a net charge of -0.9880 !!! WARNING: ARG 27 has a net charge of -0.9880 !!! WARNING: ALA 28 has a net charge of -0.2480 !!! WARNING: PHE 29 has a net charge of -0.2480 !!! WARNING: PHE 30 has a net charge of -0.2480 !!! WARNING: ALA 31 has a net charge of -0.2480 !!! WARNING: LEU 32 has a net charge of -0.2480 !!! WARNING: LEU 33 has a net charge of -0.2480 !!! WARNING: SER 34 has a net charge of -0.5580 !!! WARNING: LEU 35 has a net charge of -0.2480 !!! WARNING: SER 36 has a net charge of -0.5580 !!! WARNING: ALA 37 has a net charge of -0.2480 !!! WARNING: THR 38 has a net charge of -0.5580 !!! WARNING: ASP 39 has a net charge of -1.2480 !!! WARNING: ALA 40 has a net charge of -0.2480 !!! WARNING: PHE 41 has a net charge of -0.2480 !!! WARNING: TRP 42 has a net charge of -0.5420 !!! WARNING: ALA 43 has a net charge of -0.2480 !!! WARNING: LYS 44 has a net charge of -0.1810 !!! WARNING: TYR 45 has a net charge of -0.5870 !!! WARNING: THR 46 has a net charge of -0.5580 !!! WARNING: ALA 47 has a net charge of -0.2480 !!! WARNING: CYS 48 has a net charge of -0.3830 !!! WARNING: GLU 49 has a net charge of -1.2480 !!! WARNING: SER 50 has a net charge of -0.5580 !!! WARNING: ALA 51 has a net charge of -0.2480 !!! WARNING: ARG 52 has a net charge of -0.9880 !!! WARNING: ASN 53 has a net charge of -0.9360 !!! WARNING: ARG 55 has a net charge of -0.9880 !!! WARNING: GLU 56 has a net charge of -1.2480 !!! WARNING: LYS 57 has a net charge of -0.1810 !!! WARNING: LEU 58 has a net charge of -0.2480 !!! WARNING: ASN 59 has a net charge of -0.9360 !!! WARNING: GLU 60 has a net charge of -1.2480 !!! WARNING: CYS 61 has a net charge of -0.3830 !!! WARNING: LEU 62 has a net charge of -0.2480 !!! WARNING: GLU 63 has a net charge of -1.2480 !!! WARNING: GLY 64 has a net charge of -0.2480 !!! WARNING: GLU 68 has a net charge of -1.2480 !!! WARNING: ASN 73 has a net charge of -0.9360 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 2.9999999E-02 grid size : 179 scale,in grids/A, set to be: 0.9202533 object centred at (gu) : 0.0000000E+00 0.0000000E+00 0.0000000E+00 inner,outer dielectrics : 2.060000 78.50000 ionic strength (M) : 0.1450000 debye length (A) : 8.001814 ion exclusion radius (A) : 2.000000 probe radius (A) : 1.400000 boundary conditions : dipolar x,y,z periodic bc. flags : F F F # of linear iterations : 50 # of non-linear iterations : 0 concentration map output : F spherical charge distbn. : F INSIGHT format output : F site potential output : F modified atom file output : T map file label : qdiffxas: qdiffxs4 with an improved surfacing routine convergence graph turned off potential listings turned off box fill (%): 20.00000 xmin,xmax (A): -12.05200 26.63300 ymin,ymax (A): 21.55400 54.12100 zmin,zma (A): 104.9920 142.1620 x,y,z range (A): 38.68500 32.56700 37.16998 scale (grids/A): 0.9202533 object centre (A): 7.290501 37.83750 123.5770 number of atom coordinates read : 582 total number of charged atoms : 477 net assigned charge : -32.41700 assigned positive charge : 64.94500 centred at (gu) : 91.28522 88.98808 86.72270 assigned negative charge : -97.36198 centred at (gu) : 91.47092 88.90482 86.58081 start vw surface at 0.1000000 fill in re-entrant regions at 0.1000000 boundary points on vw surface= 4376 nv = 96 ne = 84 npr = 12243 no. pairs analyzed = 12243 no. exposed pairs = 1233 no. arc points = 8959 no. surface atoms = 411 nbur = 171 mkacc time = 2.9999994E-02 writing accessible surface arcs data to file ARCDAT grid for indexing accessible points = 1.400000 m= 417 mr = 1807 m= 8 mr = 87 m= 0 mr = 0 time to grow re-entrant surface = 3.9999999E-02 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.1400000 number of dielectric boundary points 2907 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 117 iepsmp to db, and charging done at 0.2500000 number of grid points assigned charge 2785 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid -2.3170596E-05 -7.7979730E-06 midg,1,midg; midg,igrid,midg -1.6022319E-05 -1.1057202E-05 1,midg,midg; igrid,midg,midg -1.0147353E-05 -1.7446851E-05 gauss-seidel spectral radius is 0.9936830 estimated iterations to convergence 101 setup time was (sec) 0.3700000 now iterating at: 08:36:19 rms-change max change #iterations 0.3163691 55.83022 at 10iterations 5.9570365E-02 9.136459 at 20iterations 1.1605638E-02 1.861755 at 30iterations 2.1807298E-03 0.4302330 at 40iterations 4.2744607E-04 7.4752808E-02 at 50iterations finished qdiffx linear iterations at : 08:36:22 total time elapsed so far: 2.780000 # loops : 50 mean,max change (kT/e) : 4.2744607E-04 7.4752808E-02 energy calculations done at 2.780000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 2.960000 DelPhi exited at 08:36:23 2.959u 0.129s 0:03.54 86.7% 0+0k 0+0io 0pf+0w ___________________DelPhi II____________________ / | A program to solve the PB equation | | in 3D, using non-linear form, incorporating | | 2 dielectric regions, ionic strength, periodic | | and focussing boundary conditions, utilizing | | stripped optimum successive over-relaxation | | and an improved algorithm for mapping the | | Mol. Surface to the finite-Difference grid | __________________ ____________________/ DelPhi II program started on Tue Oct 11 2011 at 08:36:23 opening parameter file fort.10 atom radii read from file fort.11 ! default extended atom radii based loosely ! on mike connolly's MS program- note H's are 0 # of radius parameter records: 18 atomic charges read from file fort.12 # of charge parameter records: 201 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 582 !!! WARNING: ALA 1 has a net charge of -0.2480 !!! WARNING: ASN 2 has a net charge of -0.9360 !!! WARNING: LYS 3 has a net charge of -0.1810 !!! WARNING: GLY 4 has a net charge of -0.2480 !!! WARNING: PHE 5 has a net charge of -0.2480 !!! WARNING: LEV 6 has a net charge of -0.6480 !!! WARNING: LEU 7 has a net charge of -0.2480 !!! WARNING: VAL 9 has a net charge of -0.2480 !!! WARNING: ARG 10 has a net charge of -0.9880 !!! WARNING: LYS 11 has a net charge of -0.1810 !!! WARNING: GLY 12 has a net charge of -0.2480 !!! WARNING: ASN 13 has a net charge of -0.9360 !!! WARNING: LEU 14 has a net charge of -0.2480 !!! WARNING: ARG 16 has a net charge of -0.9880 !!! WARNING: CYS 18 has a net charge of -0.3830 !!! WARNING: LEU 19 has a net charge of -0.2480 !!! WARNING: CYS 23 has a net charge of -0.3830 !!! WARNING: SER 24 has a net charge of -0.5580 !!! WARNING: ARG 25 has a net charge of -0.9880 !!! WARNING: ARG 26 has a net charge of -0.9880 !!! WARNING: ARG 27 has a net charge of -0.9880 !!! WARNING: ALA 28 has a net charge of -0.2480 !!! WARNING: PHE 29 has a net charge of -0.2480 !!! WARNING: PHE 30 has a net charge of -0.2480 !!! WARNING: ALA 31 has a net charge of -0.2480 !!! WARNING: LEU 32 has a net charge of -0.2480 !!! WARNING: LEU 33 has a net charge of -0.2480 !!! WARNING: SER 34 has a net charge of -0.5580 !!! WARNING: LEU 35 has a net charge of -0.2480 !!! WARNING: SER 36 has a net charge of -0.5580 !!! WARNING: ALA 37 has a net charge of -0.2480 !!! WARNING: THR 38 has a net charge of -0.5580 !!! WARNING: ASP 39 has a net charge of -1.2480 !!! WARNING: ALA 40 has a net charge of -0.2480 !!! WARNING: PHE 41 has a net charge of -0.2480 !!! WARNING: TRP 42 has a net charge of -0.5420 !!! WARNING: ALA 43 has a net charge of -0.2480 !!! WARNING: LYS 44 has a net charge of -0.1810 !!! WARNING: TYR 45 has a net charge of -0.5870 !!! WARNING: THR 46 has a net charge of -0.5580 !!! WARNING: ALA 47 has a net charge of -0.2480 !!! WARNING: CYS 48 has a net charge of -0.3830 !!! WARNING: GLU 49 has a net charge of -1.2480 !!! WARNING: SER 50 has a net charge of -0.5580 !!! WARNING: ALA 51 has a net charge of -0.2480 !!! WARNING: ARG 52 has a net charge of -0.9880 !!! WARNING: ASN 53 has a net charge of -0.9360 !!! WARNING: ARG 55 has a net charge of -0.9880 !!! WARNING: GLU 56 has a net charge of -1.2480 !!! WARNING: LYS 57 has a net charge of -0.1810 !!! WARNING: LEU 58 has a net charge of -0.2480 !!! WARNING: ASN 59 has a net charge of -0.9360 !!! WARNING: GLU 60 has a net charge of -1.2480 !!! WARNING: CYS 61 has a net charge of -0.3830 !!! WARNING: LEU 62 has a net charge of -0.2480 !!! WARNING: GLU 63 has a net charge of -1.2480 !!! WARNING: GLY 64 has a net charge of -0.2480 !!! WARNING: GLU 68 has a net charge of -1.2480 !!! WARNING: ASN 73 has a net charge of -0.9360 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 2.9999999E-02 grid size : 179 scale,in grids/A, set to be: 2.760760 object centred at (gu) : 0.0000000E+00 0.0000000E+00 0.0000000E+00 inner,outer dielectrics : 2.060000 78.50000 ionic strength (M) : 0.1450000 debye length (A) : 8.001814 ion exclusion radius (A) : 2.000000 probe radius (A) : 1.400000 boundary conditions : focussing x,y,z periodic bc. flags : F F F # of linear iterations : 120 # of non-linear iterations : 0 concentration map output : F spherical charge distbn. : F INSIGHT format output : F site potential output : F modified atom file output : T map file label : qdiffxas: qdiffxs4 with an improved surfacing routine convergence graph turned off potential listings turned off box fill (%): 60.00000 xmin,xmax (A): -12.05200 26.63300 ymin,ymax (A): 21.55400 54.12100 zmin,zma (A): 104.9920 142.1620 x,y,z range (A): 38.68500 32.56700 37.16998 scale (grids/A): 2.760760 object centre (A): 7.290501 37.83750 123.5770 number of atom coordinates read : 582 total number of charged atoms : 477 net assigned charge : -32.41700 assigned positive charge : 64.94500 centred at (gu) : 93.85563 86.96429 80.16816 assigned negative charge : -97.36198 centred at (gu) : 94.41276 86.71446 79.74225 start vw surface at 0.1100000 fill in re-entrant regions at 0.1800000 boundary points on vw surface= 43341 reading accessible surface arcs data from file ARCDAT no. of arc points read = 8959 grid for indexing accessible points = 1.400000 m= 16401 mr = 24423 m= 5064 mr = 11053 m= 1065 mr = 3291 m= 142 mr = 653 m= 13 mr = 79 m= 0 mr = 2 time to grow re-entrant surface = 0.1700000 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.3800000 number of dielectric boundary points 26525 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 0 iepsmp to db, and charging done at 0.4900000 number of grid points assigned charge 3816 setting boundary conditions focussing boundary condition read from file fort.18 focussing potential map: qdiffxas: qdiffxs4 with an improved surfacing routine original scale (grids/A) : 0.9202533 object centre at (A) : 7.290501 37.83750 123.5770 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -0.4857395 -0.1836505 midg,1,midg; midg,igrid,midg -0.2991049 -0.2354224 1,midg,midg; igrid,midg,midg -0.3214226 -0.3983139 gauss-seidel spectral radius is 0.9991326 estimated iterations to convergence 268 setup time was (sec) 0.6300000 now iterating at: 08:36:23 rms-change max change #iterations 3.864420 260.6083 at 10iterations 2.109643 96.17169 at 20iterations 1.157132 46.49427 at 30iterations 0.6314718 30.62744 at 40iterations 0.3456742 12.45815 at 50iterations 0.1890996 6.733944 at 60iterations 0.1022767 3.998983 at 70iterations 5.6097683E-02 2.188744 at 80iterations 3.0515503E-02 1.081696 at 90iterations 1.6607700E-02 0.6911163 at 100iterations 9.0246154E-03 0.3086243 at 110iterations 4.9256915E-03 0.1872101 at 120iterations finished qdiffx linear iterations at : 08:36:29 total time elapsed so far: 6.500000 # loops : 120 mean,max change (kT/e) : 4.9256915E-03 0.1872101 energy calculations done at 6.510000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 6.670000 DelPhi exited at 08:36:30 6.677u 0.161s 0:07.12 95.9% 0+0k 0+0io 0pf+0w ___________________DelPhi II____________________ / | A program to solve the PB equation | | in 3D, using non-linear form, incorporating | | 2 dielectric regions, ionic strength, periodic | | and focussing boundary conditions, utilizing | | stripped optimum successive over-relaxation | | and an improved algorithm for mapping the | | Mol. Surface to the finite-Difference grid | __________________ ____________________/ DelPhi II program started on Tue Oct 11 2011 at 08:36:30 opening parameter file fort.10 atom radii read from file fort.11 ! default extended atom radii based loosely ! on mike connolly's MS program- note H's are 0 # of radius parameter records: 18 atomic charges read from file fort.12 # of charge parameter records: 201 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 582 !!! WARNING: ALA 1 has a net charge of -0.2480 !!! WARNING: ASN 2 has a net charge of -0.9360 !!! WARNING: LYS 3 has a net charge of -0.1810 !!! WARNING: GLY 4 has a net charge of -0.2480 !!! WARNING: PHE 5 has a net charge of -0.2480 !!! WARNING: LEV 6 has a net charge of -0.6480 !!! WARNING: LEU 7 has a net charge of -0.2480 !!! WARNING: VAL 9 has a net charge of -0.2480 !!! WARNING: ARG 10 has a net charge of -0.9880 !!! WARNING: LYS 11 has a net charge of -0.1810 !!! WARNING: GLY 12 has a net charge of -0.2480 !!! WARNING: ASN 13 has a net charge of -0.9360 !!! WARNING: LEU 14 has a net charge of -0.2480 !!! WARNING: ARG 16 has a net charge of -0.9880 !!! WARNING: CYS 18 has a net charge of -0.3830 !!! WARNING: LEU 19 has a net charge of -0.2480 !!! WARNING: CYS 23 has a net charge of -0.3830 !!! WARNING: SER 24 has a net charge of -0.5580 !!! WARNING: ARG 25 has a net charge of -0.9880 !!! WARNING: ARG 26 has a net charge of -0.9880 !!! WARNING: ARG 27 has a net charge of -0.9880 !!! WARNING: ALA 28 has a net charge of -0.2480 !!! WARNING: PHE 29 has a net charge of -0.2480 !!! WARNING: PHE 30 has a net charge of -0.2480 !!! WARNING: ALA 31 has a net charge of -0.2480 !!! WARNING: LEU 32 has a net charge of -0.2480 !!! WARNING: LEU 33 has a net charge of -0.2480 !!! WARNING: SER 34 has a net charge of -0.5580 !!! WARNING: LEU 35 has a net charge of -0.2480 !!! WARNING: SER 36 has a net charge of -0.5580 !!! WARNING: ALA 37 has a net charge of -0.2480 !!! WARNING: THR 38 has a net charge of -0.5580 !!! WARNING: ASP 39 has a net charge of -1.2480 !!! WARNING: ALA 40 has a net charge of -0.2480 !!! WARNING: PHE 41 has a net charge of -0.2480 !!! WARNING: TRP 42 has a net charge of -0.5420 !!! WARNING: ALA 43 has a net charge of -0.2480 !!! WARNING: LYS 44 has a net charge of -0.1810 !!! WARNING: TYR 45 has a net charge of -0.5870 !!! WARNING: THR 46 has a net charge of -0.5580 !!! WARNING: ALA 47 has a net charge of -0.2480 !!! WARNING: CYS 48 has a net charge of -0.3830 !!! WARNING: GLU 49 has a net charge of -1.2480 !!! WARNING: SER 50 has a net charge of -0.5580 !!! WARNING: ALA 51 has a net charge of -0.2480 !!! WARNING: ARG 52 has a net charge of -0.9880 !!! WARNING: ASN 53 has a net charge of -0.9360 !!! WARNING: ARG 55 has a net charge of -0.9880 !!! WARNING: GLU 56 has a net charge of -1.2480 !!! WARNING: LYS 57 has a net charge of -0.1810 !!! WARNING: LEU 58 has a net charge of -0.2480 !!! WARNING: ASN 59 has a net charge of -0.9360 !!! WARNING: GLU 60 has a net charge of -1.2480 !!! WARNING: CYS 61 has a net charge of -0.3830 !!! WARNING: LEU 62 has a net charge of -0.2480 !!! WARNING: GLU 63 has a net charge of -1.2480 !!! WARNING: GLY 64 has a net charge of -0.2480 !!! WARNING: GLU 68 has a net charge of -1.2480 !!! WARNING: ASN 73 has a net charge of -0.9360 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 3.9999999E-02 grid size : 179 scale,in grids/A, set to be: 4.141140 object centred at (gu) : 0.0000000E+00 0.0000000E+00 0.0000000E+00 inner,outer dielectrics : 2.060000 78.50000 ionic strength (M) : 0.1450000 debye length (A) : 8.001814 ion exclusion radius (A) : 2.000000 probe radius (A) : 1.400000 boundary conditions : focussing x,y,z periodic bc. flags : F F F # of linear iterations : 180 # of non-linear iterations : 0 concentration map output : F spherical charge distbn. : F INSIGHT format output : F site potential output : F modified atom file output : T map file label : qdiffxas: qdiffxs4 with an improved surfacing routine convergence graph turned off potential listings turned off box fill (%): 90.00000 xmin,xmax (A): -12.05200 26.63300 ymin,ymax (A): 21.55400 54.12100 zmin,zma (A): 104.9920 142.1620 x,y,z range (A): 38.68500 32.56700 37.16998 scale (grids/A): 4.141140 object centre (A): 7.290501 37.83750 123.5770 number of atom coordinates read : 582 total number of charged atoms : 477 net assigned charge : -32.41700 assigned positive charge : 64.94500 centred at (gu) : 95.78351 85.44641 75.25217 assigned negative charge : -97.36198 centred at (gu) : 96.61909 85.07173 74.61333 start vw surface at 0.1100000 fill in re-entrant regions at 0.5300000 boundary points on vw surface= 99156 reading accessible surface arcs data from file ARCDAT no. of arc points read = 8959 grid for indexing accessible points = 1.400000 m= 47340 mr = 59755 m= 23100 mr = 35789 m= 10033 mr = 17969 m= 3566 mr = 7497 m= 1002 mr = 2668 m= 215 mr = 731 m= 39 mr = 163 m= 3 mr = 28 m= 0 mr = 0 time to grow re-entrant surface = 0.5200000 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 1.140000 number of dielectric boundary points 59854 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 0 iepsmp to db, and charging done at 1.260000 number of grid points assigned charge 3816 setting boundary conditions focussing boundary condition read from file fort.18 focussing potential map: qdiffxas: qdiffxs4 with an improved surfacing routine original scale (grids/A) : 2.760760 object centre at (A) : 7.290501 37.83750 123.5770 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -2.620437 -1.082290 midg,1,midg; midg,igrid,midg -1.585276 -1.286280 1,midg,midg; igrid,midg,midg -1.889358 -2.014197 gauss-seidel spectral radius is 0.9994030 estimated iterations to convergence 323 setup time was (sec) 1.410000 now iterating at: 08:36:31 rms-change max change #iterations 7.882366 265.0370 at 10iterations 4.877199 145.5187 at 20iterations 3.002034 80.87122 at 30iterations 1.843804 41.57292 at 40iterations 1.129555 24.22997 at 50iterations 0.6887297 17.72733 at 60iterations 0.4189128 8.997574 at 70iterations 0.2546009 4.981106 at 80iterations 0.1558531 3.119100 at 90iterations 9.4380289E-02 1.868408 at 100iterations 5.7440255E-02 1.151321 at 110iterations 3.4875192E-02 0.6804199 at 120iterations 2.1226138E-02 0.4532013 at 130iterations 1.2913283E-02 0.2502651 at 140iterations 7.8567378E-03 0.1572351 at 150iterations 4.7829319E-03 0.1028900 at 160iterations 2.9038303E-03 5.3611755E-02 at 170iterations 1.7734886E-03 3.4690857E-02 at 180iterations finished qdiffx linear iterations at : 08:36:41 total time elapsed so far: 10.64000 # loops : 180 mean,max change (kT/e) : 1.7734886E-03 3.4690857E-02 energy calculations done at 10.64000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 10.81000 DelPhi exited at 08:36:41 10.813u 0.183s 0:11.31 97.1% 0+0k 0+0io 0pf+0w