Tue Oct 11 08:35:56 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+sph.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 = amb.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:35:56 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: 192 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: LEU 6 has a net charge of -0.2480 !!! 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.01700 assigned positive charge : 64.94500 centred at (gu) : 91.28522 88.98808 86.72270 assigned negative charge : -96.96198 centred at (gu) : 91.51642 88.90321 86.52583 start vw surface at 0.1000000 fill in re-entrant regions at 0.1100000 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 = 1.9999996E-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.0000001E-02 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.1500000 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.3430022E-05 -7.4990385E-06 midg,1,midg; midg,igrid,midg -1.5819816E-05 -1.0894654E-05 1,midg,midg; igrid,midg,midg -9.7904531E-06 -1.7559214E-05 gauss-seidel spectral radius is 0.9936830 estimated iterations to convergence 101 setup time was (sec) 0.3500000 now iterating at: 08:35:56 rms-change max change #iterations 0.3150755 55.83022 at 10iterations 5.9296470E-02 9.125778 at 20iterations 1.1556951E-02 1.873032 at 30iterations 2.1728929E-03 0.4296341 at 40iterations 4.2590982E-04 7.5569153E-02 at 50iterations finished qdiffx linear iterations at : 08:35:59 total time elapsed so far: 2.760000 # loops : 50 mean,max change (kT/e) : 4.2590982E-04 7.5569153E-02 energy calculations done at 2.760000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 2.940000 DelPhi exited at 08:35:59 2.942u 0.148s 0:03.40 90.5% 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:35:59 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: 192 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: LEU 6 has a net charge of -0.2480 !!! 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.01700 assigned positive charge : 64.94500 centred at (gu) : 93.85563 86.96429 80.16816 assigned negative charge : -96.96198 centred at (gu) : 94.54925 86.70965 79.57732 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.4800000 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.4852760 -0.1748860 midg,1,midg; midg,igrid,midg -0.2974676 -0.2337275 1,midg,midg; igrid,midg,midg -0.3087755 -0.3979732 gauss-seidel spectral radius is 0.9991326 estimated iterations to convergence 268 setup time was (sec) 0.6200000 now iterating at: 08:36:00 rms-change max change #iterations 3.849104 260.6085 at 10iterations 2.099455 96.17169 at 20iterations 1.151501 45.76146 at 30iterations 0.6281291 30.61785 at 40iterations 0.3436731 12.48897 at 50iterations 0.1879320 6.732279 at 60iterations 0.1016805 4.005270 at 70iterations 5.5728786E-02 2.178051 at 80iterations 3.0307593E-02 1.091660 at 90iterations 1.6493469E-02 0.6695023 at 100iterations 8.9657363E-03 0.3072815 at 110iterations 4.8935665E-03 0.1849365 at 120iterations finished qdiffx linear iterations at : 08:36:06 total time elapsed so far: 6.490000 # loops : 120 mean,max change (kT/e) : 4.8935665E-03 0.1849365 energy calculations done at 6.490000 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:06 6.668u 0.188s 0:07.31 93.5% 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:06 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: 192 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: LEU 6 has a net charge of -0.2480 !!! 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.0000000E-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.01700 assigned positive charge : 64.94500 centred at (gu) : 95.78351 85.44641 75.25217 assigned negative charge : -96.96198 centred at (gu) : 96.82383 85.06450 74.36593 start vw surface at 0.1100000 fill in re-entrant regions at 0.5400000 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.5599999 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 1.190000 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.310000 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.618525 -1.039949 midg,1,midg; midg,igrid,midg -1.577899 -1.278753 1,midg,midg; igrid,midg,midg -1.820230 -2.012671 gauss-seidel spectral radius is 0.9994030 estimated iterations to convergence 323 setup time was (sec) 1.460000 now iterating at: 08:36:08 rms-change max change #iterations 7.839315 265.0367 at 10iterations 4.846381 145.5187 at 20iterations 2.981292 80.87311 at 30iterations 1.829712 41.91902 at 40iterations 1.119424 23.53503 at 50iterations 0.6823040 16.86317 at 60iterations 0.4148052 8.955307 at 70iterations 0.2519868 4.940487 at 80iterations 0.1541470 3.113848 at 90iterations 9.3304686E-02 1.866066 at 100iterations 5.6783315E-02 1.131256 at 110iterations 3.4469981E-02 0.6859932 at 120iterations 2.0967865E-02 0.4450302 at 130iterations 1.2761410E-02 0.2507401 at 140iterations 7.7665248E-03 0.1588135 at 150iterations 4.7272132E-03 0.1005554 at 160iterations 2.8698165E-03 5.4267883E-02 at 170iterations 1.7538343E-03 3.4744263E-02 at 180iterations finished qdiffx linear iterations at : 08:36:18 total time elapsed so far: 11.43000 # loops : 180 mean,max change (kT/e) : 1.7538343E-03 3.4744263E-02 energy calculations done at 11.43000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 11.65000 DelPhi exited at 08:36:19 11.655u 0.174s 0:12.28 96.2% 0+0k 0+0io 0pf+0w