AMBER Archive (2008)

Subject: Re: AMBER: increasing temperature smoothly

From: Germain Vallverdu (germain.vallverdu_at_lcp.u-psud.fr)
Date: Mon Sep 22 2008 - 06:31:56 CDT


Hi

Thank you for your previous help, my input now work correctly, but I am
still not sure of what it will do.
This is my input

***********************************************************************
equilibration a volume constant thermalisation avec contrainte prot
&cntrl
 nstlim=50000, dt=0.002,
 irest=0, ntx=1, t=0.0, ig=190908,
 ntpr=1000, ntwr=1000, ntwx=1000, iwrap=1,
 tempi=100.0, temp0=300.0, ntt=3, gamma_ln=1.0,
 ntb=1, ntp=0,
 ntc=2, ntf=2,
 cut=10.0
 ntr=1, restraint_wt=5.0, restraintmask=":1-227",
 nmropt=1,
/
&wt
 type='TEMP0',
 istep1=1,
 istep2=20000,
 istep2=0,
 imult=0,
 iinc=0,
 value1=100,
 value2=300,
/
&wt
 type='END'
/
****************************************

And I got that on the output :

           Begin reading energy term weight changes/NMR restraints
 WEIGHT CHANGES:
 TEMP0 1 0 100.000000 300.000000 0 0
 Warning: ISTEP1 >= ISTEP2

 RESTRAINTS:
                          ** No restraint defined **

                  Done reading weight changes/NMR restraints

It looks like sander does not correctly read the istep1 and istep2 or I
do not understand the warning. I put a the end of this mail the mdout file, I think TEMP0 is
always 100K and do not smoothly increase after 20000 steps TEMP is still
around 100K.

An other thing, I change my input in order to use new flags for
restraint with ntr=1 (restraint_wt and restraintmask). I got the correct
number of atoms but the restraint_wt value never appears in the output
file (only at the beginning in my input), is it ok and is it a way to
check that the value is correctly read ?

Thanks

Germain

mdout :

          -------------------------------------------------------
          Amber 9 SANDER 2006
          -------------------------------------------------------

| Run on 09/22/2008 at 05:59:33
  [-O]verwriting output

File Assignments:
| MDIN: eq_v.in
| MDOUT: eq_v.out
|INPCRD: mintot.restrt
| PARM: ../../ceru_pH8.top
|RESTRT: eq_v.restrt
| REFC: mintot.restrt
| MDVEL: mdvel
| MDEN: mden
| MDCRD: eq_v.crd
|MDINFO: mdinfo
|INPDIP: inpdip
|RSTDIP: rstdip

|INPTRA: inptraj
|

 Here is the input file:

equilibration a volume constant thermalisation avec contrainte prot
&cntrl
 nstlim=50000, dt=0.002,
 irest=0, ntx=1, t=0.0, ig=190908,
 ntpr=1000, ntwr=1000, ntwx=1000, iwrap=1,
 tempi=100.0, temp0=300.0, ntt=3, gamma_ln=1.0,
 ntb=1, ntp=0,
 ntc=2, ntf=2,
 cut=10.0
 ntr=1, restraint_wt=5.0, restraintmask=":1-227",
 nmropt=1,
/
&wt
 type='TEMP0',
 istep1=1,
 istep2=20000,
 istep2=0,
 imult=0,
 iinc=0,
 value1=100,
 value2=300,
/
&wt
 type='END'
/
                                                                               

--------------------------------------------------------------------------------
   1. RESOURCE USE:
--------------------------------------------------------------------------------

| Flags: MPI
 getting new box info from bottom of inpcrd
| INFO: Old style inpcrd file read

| peek_ewald_inpcrd: Box info found
|Largest sphere to fit in unit cell has radius = 27.854
| New format PARM file being parsed.
| Version = 1.000 Date = 07/24/08 Time = 18:04:28
 NATOM = 27219 NTYPES = 17 NBONH = 25405 MBONA = 1851
 NTHETH = 4047 MTHETA = 2506 NPHIH = 7734 MPHIA = 4855
 NHPARM = 0 NPARM = 0 NNB = 51324 NRES = 8110
 NBONA = 1851 NTHETA = 2506 NPHIA = 4855 NUMBND = 58
 NUMANG = 123 NPTRA = 44 NATYP = 38 NPHB = 1
 IFBOX = 1 NMXRS = 40 IFCAP = 0 NEXTRA = 0
 NCOPY = 0

| Memory Use Allocated
| Real 1581266
| Hollerith 171426
| Integer 863206
| Max Pairs 1959768
| nblistReal 326628
| nblist Int 848268
| Total 29915 kbytes
| Duplicated 0 dihedrals
| Duplicated 0 dihedrals

     BOX TYPE: RECTILINEAR

--------------------------------------------------------------------------------
   2. CONTROL DATA FOR THE RUN
--------------------------------------------------------------------------------

                                                                                

General flags:
     imin = 0, nmropt = 1

Nature and format of input:
     ntx = 1, irest = 0, ntrx = 1

Nature and format of output:
     ntxo = 1, ntpr = 1000, ntrx = 1, ntwr = 1000
     iwrap = 1, ntwx = 1000, ntwv = 0, ntwe = 0
     ioutfm = 0, ntwprt = 0, idecomp = 0, rbornstat= 0

Potential function:
     ntf = 2, ntb = 1, igb = 0, nsnb = 25
     ipol = 0, gbsa = 0, iesp = 0
     dielc = 1.00000, cut = 10.00000, intdiel = 1.00000
     scnb = 2.00000, scee = 1.20000

Frozen or restrained atoms:
     ibelly = 0, ntr = 1

Molecular dynamics:
     nstlim = 50000, nscm = 1000, nrespa = 1
     t = 0.00000, dt = 0.00200, vlimit = 20.00000

Langevin dynamics temperature regulation:
     ig = 190908
     temp0 = 300.00000, tempi = 100.00000, gamma_ln= 1.00000

SHAKE:
     ntc = 2, jfastw = 0
     tol = 0.00001

NMR refinement options:
     iscale = 0, noeskp = 1, ipnlty = 1, mxsub = 1
     scalm = 100.00000, pencut = 0.10000, tausw = 0.10000

Ewald parameters:
     verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1
     vdwmeth = 1, eedmeth = 1, netfrc = 1
     Box X = 55.707 Box Y = 67.721 Box Z = 71.744
     Alpha = 90.000 Beta = 90.000 Gamma = 90.000
     NFFT1 = 60 NFFT2 = 72 NFFT3 = 72
     Cutoff= 10.000 Tol =0.100E-04
     Ewald Coefficient = 0.27511
     Interpolation order = 4

    LOADING THE CONSTRAINED ATOMS AS GROUPS

   5. REFERENCE ATOM COORDINATES

      
     Mask :1-227; matches 3580 atoms

--------------------------------------------------------------------------------
   3. ATOMIC COORDINATES AND VELOCITIES
--------------------------------------------------------------------------------

                                                                                
 begin time read from input coords = 0.000 ps

           Begin reading energy term weight changes/NMR restraints
 WEIGHT CHANGES:
 TEMP0 1 0 100.000000 300.000000 0 0
 Warning: ISTEP1 >= ISTEP2

 RESTRAINTS:
                          ** No restraint defined **

                  Done reading weight changes/NMR restraints

 Number of triangulated 3-point waters found: 7878
| Atom division among processors:
| 0 3414 6804 10206 13608 17010 20412 23814
| 27219

     Sum of charges from parm topology file = -0.00000010
     Forcing neutrality...
| Running AMBER/MPI version on 8 nodes

--------------------------------------------------------------------------------
   4. RESULTS
--------------------------------------------------------------------------------

 ---------------------------------------------------
 APPROXIMATING switch and d/dx switch using CUBIC SPLINE INTERPOLATION
 using 5000.0 points per unit in tabled values
 TESTING RELATIVE ERROR over r ranging from 0.0 to cutoff
| CHECK switch(x): max rel err = 0.2738E-14 at 2.422500
| CHECK d/dx switch(x): max rel err = 0.8350E-11 at 2.736960
 ---------------------------------------------------
| Local SIZE OF NONBOND LIST = 1285982
| TOTAL SIZE OF NONBOND LIST = 9860526

 NSTEP = 0 TIME(PS) = 0.000 TEMP(K) = 144.71 PRESS = 0.0
 Etot = -110687.6567 EKtot = 8088.1491 EPtot = -118775.8058
 BOND = 118.2113 ANGLE = 543.0451 DIHED = 1870.1948
 1-4 NB = 708.8992 1-4 EEL = 10994.4530 VDWAALS = 16431.8065
 EELEC = -149442.4157 EHBOND = 0.0000 RESTRAINT = 0.0000
 Ewald error estimate: 0.4026E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 1000 TIME(PS) = 2.000 TEMP(K) = 98.34 PRESS = 0.0
 Etot = -93779.4603 EKtot = 5496.3685 EPtot = -99275.8287
 BOND = 283.2210 ANGLE = 906.7699 DIHED = 1944.4381
 1-4 NB = 736.5023 1-4 EEL = 10999.8051 VDWAALS = 13401.2728
 EELEC = -127725.5356 EHBOND = 0.0000 RESTRAINT = 177.6977
 EAMBER (non-restraint) = -99453.5264
 Ewald error estimate: 0.6712E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 2000 TIME(PS) = 4.000 TEMP(K) = 101.68 PRESS = 0.0
 Etot = -93806.1331 EKtot = 5683.2197 EPtot = -99489.3529
 BOND = 301.4622 ANGLE = 981.5110 DIHED = 1952.6719
 1-4 NB = 732.8871 1-4 EEL = 10957.1901 VDWAALS = 13806.9959
 EELEC = -128411.9235 EHBOND = 0.0000 RESTRAINT = 189.8524
 EAMBER (non-restraint) = -99679.2053
 Ewald error estimate: 0.1988E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 3000 TIME(PS) = 6.000 TEMP(K) = 101.30 PRESS = 0.0
 Etot = -94196.6017 EKtot = 5661.5856 EPtot = -99858.1872
 BOND = 300.3195 ANGLE = 928.1096 DIHED = 1943.6064
 1-4 NB = 739.8021 1-4 EEL = 10991.5404 VDWAALS = 14146.4001
 EELEC = -129104.9920 EHBOND = 0.0000 RESTRAINT = 197.0268
 EAMBER (non-restraint) = -100055.2141
 Ewald error estimate: 0.4545E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 4000 TIME(PS) = 8.000 TEMP(K) = 100.95 PRESS = 0.0
 Etot = -94522.2647 EKtot = 5642.5383 EPtot = -100164.8030
 BOND = 307.4636 ANGLE = 943.1072 DIHED = 1960.5755
 1-4 NB = 733.9056 1-4 EEL = 10998.7824 VDWAALS = 14211.2785
 EELEC = -129506.2367 EHBOND = 0.0000 RESTRAINT = 186.3209
 EAMBER (non-restraint) = -100351.1239
 Ewald error estimate: 0.8982E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 5000 TIME(PS) = 10.000 TEMP(K) = 100.06 PRESS = 0.0
 Etot = -94708.5195 EKtot = 5592.5980 EPtot = -100301.1175
 BOND = 293.3545 ANGLE = 950.9084 DIHED = 1944.1587
 1-4 NB = 744.2482 1-4 EEL = 10991.8788 VDWAALS = 14249.4491
 EELEC = -129676.2080 EHBOND = 0.0000 RESTRAINT = 201.0928
 EAMBER (non-restraint) = -100502.2103
 Ewald error estimate: 0.2167E-03
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 6000 TIME(PS) = 12.000 TEMP(K) = 100.03 PRESS = 0.0
 Etot = -94809.5298 EKtot = 5591.0975 EPtot = -100400.6273
 BOND = 310.6405 ANGLE = 956.0766 DIHED = 1951.8952
 1-4 NB = 727.4308 1-4 EEL = 10993.9714 VDWAALS = 14416.8062
 EELEC = -129953.0676 EHBOND = 0.0000 RESTRAINT = 195.6197
 EAMBER (non-restraint) = -100596.2470
 Ewald error estimate: 0.1127E-03
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 7000 TIME(PS) = 14.000 TEMP(K) = 100.77 PRESS = 0.0
 Etot = -94881.9311 EKtot = 5632.2554 EPtot = -100514.1865
 BOND = 315.5354 ANGLE = 917.0407 DIHED = 1943.1465
 1-4 NB = 734.1075 1-4 EEL = 11004.2035 VDWAALS = 14646.7410
 EELEC = -130263.7498 EHBOND = 0.0000 RESTRAINT = 188.7887
 EAMBER (non-restraint) = -100702.9752
 Ewald error estimate: 0.7061E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 8000 TIME(PS) = 16.000 TEMP(K) = 100.36 PRESS = 0.0
 Etot = -94978.5408 EKtot = 5609.3902 EPtot = -100587.9311
 BOND = 274.8586 ANGLE = 953.0856 DIHED = 1948.5728
 1-4 NB = 739.5875 1-4 EEL = 10998.4873 VDWAALS = 14691.2717
 EELEC = -130372.2568 EHBOND = 0.0000 RESTRAINT = 178.4621
 EAMBER (non-restraint) = -100766.3932
 Ewald error estimate: 0.1038E-03
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 9000 TIME(PS) = 18.000 TEMP(K) = 101.15 PRESS = 0.0
 Etot = -94878.1360 EKtot = 5653.7020 EPtot = -100531.8380
 BOND = 294.9783 ANGLE = 952.5742 DIHED = 1948.6110
 1-4 NB = 727.0705 1-4 EEL = 10964.4690 VDWAALS = 14758.9956
 EELEC = -130370.5142 EHBOND = 0.0000 RESTRAINT = 191.9776
 EAMBER (non-restraint) = -100723.8155
 Ewald error estimate: 0.1021E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 10000 TIME(PS) = 20.000 TEMP(K) = 99.85 PRESS = 0.0
 Etot = -95064.1909 EKtot = 5580.6555 EPtot = -100644.8465
 BOND = 302.4455 ANGLE = 945.5122 DIHED = 1950.4055
 1-4 NB = 727.4089 1-4 EEL = 10990.9253 VDWAALS = 14723.5195
 EELEC = -130468.2563 EHBOND = 0.0000 RESTRAINT = 183.1929
 EAMBER (non-restraint) = -100828.0394
 Ewald error estimate: 0.7719E-05
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 11000 TIME(PS) = 22.000 TEMP(K) = 100.12 PRESS = 0.0
 Etot = -95176.9619 EKtot = 5595.9174 EPtot = -100772.8793
 BOND = 290.6507 ANGLE = 965.1600 DIHED = 1939.1947
 1-4 NB = 735.4330 1-4 EEL = 10997.1517 VDWAALS = 14862.2462
 EELEC = -130750.9726 EHBOND = 0.0000 RESTRAINT = 188.2570
 EAMBER (non-restraint) = -100961.1363
 Ewald error estimate: 0.2384E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 12000 TIME(PS) = 24.000 TEMP(K) = 100.79 PRESS = 0.0
 Etot = -95090.5945 EKtot = 5633.5191 EPtot = -100724.1135
 BOND = 290.0453 ANGLE = 956.8040 DIHED = 1944.7565
 1-4 NB = 734.8787 1-4 EEL = 11010.1155 VDWAALS = 14655.9865
 EELEC = -130510.3061 EHBOND = 0.0000 RESTRAINT = 193.6059
 EAMBER (non-restraint) = -100917.7195
 Ewald error estimate: 0.1392E-03
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 13000 TIME(PS) = 26.000 TEMP(K) = 100.97 PRESS = 0.0
 Etot = -95214.9848 EKtot = 5643.1427 EPtot = -100858.1276
 BOND = 295.0800 ANGLE = 924.4669 DIHED = 1948.8010
 1-4 NB = 732.6476 1-4 EEL = 10998.3931 VDWAALS = 14675.5983
 EELEC = -130618.0450 EHBOND = 0.0000 RESTRAINT = 184.9305
 EAMBER (non-restraint) = -101043.0581
 Ewald error estimate: 0.5939E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 14000 TIME(PS) = 28.000 TEMP(K) = 101.27 PRESS = 0.0
 Etot = -95171.5861 EKtot = 5660.1735 EPtot = -100831.7597
 BOND = 286.9686 ANGLE = 945.2925 DIHED = 1948.9398
 1-4 NB = 745.5049 1-4 EEL = 11028.5898 VDWAALS = 14745.0702
 EELEC = -130733.7912 EHBOND = 0.0000 RESTRAINT = 201.6656
 EAMBER (non-restraint) = -101033.4253
 Ewald error estimate: 0.3847E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 15000 TIME(PS) = 30.000 TEMP(K) = 99.83 PRESS = 0.0
 Etot = -95236.4020 EKtot = 5579.6272 EPtot = -100816.0292
 BOND = 291.2661 ANGLE = 962.5928 DIHED = 1947.1570
 1-4 NB = 733.8280 1-4 EEL = 10974.0593 VDWAALS = 14744.8398
 EELEC = -130658.0259 EHBOND = 0.0000 RESTRAINT = 188.2538
 EAMBER (non-restraint) = -101004.2829
 Ewald error estimate: 0.1587E-03
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 16000 TIME(PS) = 32.000 TEMP(K) = 100.85 PRESS = 0.0
 Etot = -95325.3572 EKtot = 5636.5220 EPtot = -100961.8791
 BOND = 287.7508 ANGLE = 955.6392 DIHED = 1945.3710
 1-4 NB = 732.5010 1-4 EEL = 11003.3367 VDWAALS = 14782.1138
 EELEC = -130854.7970 EHBOND = 0.0000 RESTRAINT = 186.2053
 EAMBER (non-restraint) = -101148.0845
 Ewald error estimate: 0.7738E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 17000 TIME(PS) = 34.000 TEMP(K) = 99.79 PRESS = 0.0
 Etot = -95302.3383 EKtot = 5577.6139 EPtot = -100879.9521
 BOND = 303.1840 ANGLE = 950.9092 DIHED = 1942.3307
 1-4 NB = 736.1861 1-4 EEL = 10988.0747 VDWAALS = 14728.4117
 EELEC = -130719.4911 EHBOND = 0.0000 RESTRAINT = 190.4426
 EAMBER (non-restraint) = -101070.3947
 Ewald error estimate: 0.6943E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 18000 TIME(PS) = 36.000 TEMP(K) = 100.94 PRESS = 0.0
 Etot = -95287.1372 EKtot = 5641.9601 EPtot = -100929.0973
 BOND = 293.4712 ANGLE = 948.4529 DIHED = 1946.1384
 1-4 NB = 737.6260 1-4 EEL = 11001.3173 VDWAALS = 14694.1809
 EELEC = -130743.7090 EHBOND = 0.0000 RESTRAINT = 193.4249
 EAMBER (non-restraint) = -101122.5222
 Ewald error estimate: 0.6305E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 19000 TIME(PS) = 38.000 TEMP(K) = 99.67 PRESS = 0.0
 Etot = -95371.8561 EKtot = 5570.7119 EPtot = -100942.5680
 BOND = 295.8473 ANGLE = 936.0317 DIHED = 1949.5244
 1-4 NB = 739.2216 1-4 EEL = 10973.2399 VDWAALS = 14708.5423
 EELEC = -130746.2158 EHBOND = 0.0000 RESTRAINT = 201.2407
 EAMBER (non-restraint) = -101143.8087
 Ewald error estimate: 0.6972E-05
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

 NSTEP = 20000 TIME(PS) = 40.000 TEMP(K) = 100.34 PRESS = 0.0
 Etot = -95411.6817 EKtot = 5608.3597 EPtot = -101020.0414
 BOND = 306.2103 ANGLE = 936.1500 DIHED = 1942.5709
 1-4 NB = 734.2413 1-4 EEL = 10958.8019 VDWAALS = 14704.4575
 EELEC = -130796.4196 EHBOND = 0.0000 RESTRAINT = 193.9463
 EAMBER (non-restraint) = -101213.9877
 Ewald error estimate: 0.7502E-04
 ------------------------------------------------------------------------------

 NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 0.000
===============================================================================

******************************* end

Ross Walker a écrit :
> Hi Germain,
>
> Firstly you don't have nmropt=1 in you cntrl namelist. I think the
> restraints also need to come last, after the &wt namelists.
>
> All the best
> Ross
>
>
>> -----Original Message-----
>> From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu] On Behalf
>> Of Germain Vallverdu
>> Sent: Friday, September 19, 2008 8:18 AM
>> To: amber_at_scripps.edu
>> Subject: AMBER: increasing temperature smoothly
>>
>> Dear amber users
>>
>> I would like to know if the following input will do what I want. I would
>> like to start a simulation with temp0=100K and in the first 20000 steps,
>> I would like temp0 increase linearly from 100K to 300K and then keep the
>> value of 300K for the remainder of the run.
>>
>> equilibration a volume constant thermalisation avec contrainte prot
>> &cntrl
>> nstlim=50000, dt=0.002,
>> irest=0, ntx=1, t=0.0, ig=300708,
>> ntpr=1000, ntwr=1000, ntwx=1000,
>> tempi=100.0, temp0=300.0, ntt=3,
>> ntb=1, ntp=0,
>> ntc=2, ntf=2,
>> cut=10.0
>> ntr=1,
>> /
>> Keep protein fixed with weak restraints
>> 10.0
>> RES 1 227
>> END
>> END
>> &wt
>> type='TEMP0',
>> istep1=1,
>> istep2=20000,
>> istep2=0,
>> imult=0,
>> iinc=0,
>> value1=100,
>> value2=300,
>> /
>> &wt
>> type='END'
>> /
>>
>> I got the following error message
>>
>> PGFIO-F-252/formatted read/unit=5/operation attempted after end of file.
>> File name = eq_v.in formatted, sequential access record = 34
>> In source file _rgroup.f, at line number 175
>>
>>
>> Thanks
>>
>> Germain
>>
>> --
>> Germain Vallverdu
>> Laboratoire de Chimie Physique
>> Université Paris Sud 11
>> germain.vallverdu_at_lcp.u-psud.fr
>> ma page perso <http://pagesperso.lcp.u-psud.fr/vallverdu/>
>> 01 69 15 30 38 / 06 88 59 08 87
>> -----------------------------------------------------------------------
>> The AMBER Mail Reflector
>> To post, send mail to amber_at_scripps.edu
>> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
>> to majordomo_at_scripps.edu
>>
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
> to majordomo_at_scripps.edu
>

-- 
Germain Vallverdu
Laboratoire de Chimie Physique
Université Paris Sud 11
germain.vallverdu_at_lcp.u-psud.fr
ma page perso <http://pagesperso.lcp.u-psud.fr/vallverdu/>
01 69 15 30 38 / 06 88 59 08 87

-- Germain Vallverdu Laboratoire de Chimie Physique Université Paris Sud 11 germain.vallverdu_at_lcp.u-psud.fr ma page perso <http://pagesperso.lcp.u-psud.fr/vallverdu/> 01 69 15 30 38 / 06 88 59 08 87 ----------------------------------------------------------------------- The AMBER Mail Reflector To post, send mail to amber_at_scripps.edu To unsubscribe, send "unsubscribe amber" (in the *body* of the email) to majordomo_at_scripps.edu