[Pw_forum] vc-md only reach half of target temperature
Dear Forum, It try to heat vc-md calculations. But I can not reach the target temperature. When I set tempw = 100 , it will stable at about 38K When I set tempw=1000, it will stable at about 400K. I pasted part of my input bellow. Any one have some idea? Thank you very much. ''' calculation = 'vc-md' , prefix = 'alpo4-10' , restart_mode = 'from_scratch' , nstep = 3 , /iprint = 1 , dt = 20 , / pot_extrapolation='second-order' wfc_extrapolation='second-order' ion_temperature = 'berendsen' , tempw = 100 , / cell_dynamics = 'pr' , press = 0 , ''' WANG Riping 2012.10.12 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://pwscf.org/pipermail/pw_forum/attachments/20121011/132b03bf/attachment.html
[Pw_forum] vc-md could be stop by nstep after a slight change in code
Dear forum, The vc-md calculation in the present espresso-5.0 is unable to stop by setting nstep. I found the problem in the code file pw/src/vcsmd.f90. I modified the line '''ELSE IF ( calc == 'nm' .OR. calc == 'cm' ) THEN''' into '''ELSE IF ( calc == 'nm' .OR. calc == 'cm' .OR. calc == 'cd' ) THEN'''. Now it works for stoping calculations by nstep. Hope to make it corrected in next version. Yours. WANG Riping 2012.10.5 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://pwscf.org/pipermail/pw_forum/attachments/20121005/8cb4f8af/attachment.html
[Pw_forum] Pw_forum Digest, Vol 64, Issue 6
Dear Paolo When the link states that it is very easy to fix, but it did not tell me how to fix. Is there any more clear answer? Thank you. WANG On 3 October 2012 08:50, wrote: > SMessage: 4 > Date: Tue, 2 Oct 2012 21:30:15 +0200 > From: Paolo Giannozzi > Subject: Re: [Pw_forum] How to stop vc-md with given nstep? > To: PWSCF Forum > Message-ID: > Content-Type: text/plain; charset=US-ASCII; format=flowed > > > On Oct 2, 2012, at 11:09 , Riping Wang wrote: > > > > I am using vc-md with espress-5.0. > > I try to run nstep=2 equal to some integer. > > But It already run 428 stpes now, and it is still runing. > > it's a known problem: >http://www.democritos.it/pipermail/pw_forum/2012-July/024718.html > It is very easy to fix it, though. -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://pwscf.org/pipermail/pw_forum/attachments/20121003/40978876/attachment.html
[Pw_forum] How to stop vc-md with given nstep?
Dear developer, I am using vc-md with espress-5.0. I try to run nstep=2 equal to some integer. But It already run 428 stpes now, and it is still runing. Part of input file and output file are appended bellow. Could anyone help? Thank you. WANG Riping 2012.10.2 input: ''' calculation = 'vc-md' , prefix = 'alpo4-10' , restart_mode = 'from_scratch' , nstep = 2 , iprint = 1 , dt = 5 , tstress = .TRUE. , tprnfor = .TRUE. , disk_io = 'high' , wf_collect = .true. , pseudo_dir = /home/wang/espresso/pseudo , ''' output: ''' bravais-lattice index = 14 lattice parameter (alat) = 8.0856 a.u. unit-cell volume = 333.5567 (a.u.)^3 number of atoms/cell =6 number of atomic types=3 number of electrons =32.00 number of Kohn-Sham states= 16 kinetic-energy cutoff = 50. Ry charge density cutoff = 400. Ry convergence threshold = 1.0E-06 mixing beta = 0.7000 number of iterations used =8 plain mixing Exchange-correlation = SLA PW PBE PBE ( 1 4 3 4 0) EXX-fraction =0.00 nstep =2 ... Entering Dynamics; it = 428 time = 0.10333 pico-seconds ''' -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://pwscf.org/pipermail/pw_forum/attachments/20121002/084d5c77/attachment.html
[Pw_forum] oscillation of P T in cp
Dear developers, when I use vc-cp calculations in quantum espresso 5.0, I have a question: in this code, how much is normal oscillation of pressure and temperature that could be controlled? Thank you very much. WANG Riping 2012.9.6 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20120906/c26d1152/attachment.htm
[Pw_forum] symmetry not found in pw.x and ph.x
Dear pwforum, I am calculating scf for a ibrav=4 structure with espresso-4.3.2. But it can not find the symmetry. In the output file, it report "no symmetry found". I am also calculation with ph.x with the same. the phonon mode is not identified. I confirmed the atomic positions in the input file. It is consistent with 32 point symmetry (space group 154). I tried ibrav=2 and 12, they can find symmetry and phonon modes, even with atomic positions slightly different from point symmety formula.. So I am wondering: what kind of situation, the code can not find symmtery? Is there any suggestion to make it find symmetry? I attached the input file and output file of pw.x followiing. Thanks very much for you attention and helps. Best withes. WANG Riping 2011.8.11 calculation = scf , restart_mode = 'from_scratch' , prefix = alpo4-154.0_0 , pseudo_dir = '~/pseudo', outdir = outdir , nstep = 500 , max_seconds = 36 , etot_conv_thr = 5.0D-6 , forc_conv_thr = 1.0D10 , tprnfor = .true. / ibrav = 4 , celldm(1) = 9.471757419 , celldm(3) = 2.21251394 , nat = 18 , ntyp = 3 , ecutwfc = 50.0 , ecutrho = 400.0 , / conv_thr = 1.0d-10 , / ion_dynamics = none , / cell_dynamics = none , press = 0 , press_conv_thr = 0.05D0 , / ATOMIC_SPECIES O 15.9994 O.pbe-van_ak.UPF Al 26.981538 Al.pbe-n-van.UPF P 30.973761 P.pbe-van_ak.UPF ATOMIC_POSITIONS crystal Al 0.4692 0. 0.16671 0 0 Al 0. 0.4692 0.83330 1 0 Al 0.5308 0.5308 0.50001 1 0 P0.4694 0. 0.66671 0 0 P0. 0.4694 0.0 1 0 P0.5306 0.5306 0.1 1 0 O0.4160 0.2860 0.1011 O0.7140 0.1300 0.7678 O0.8700 0.5840 0.4344 O0.2860 0.4160 0.8989 O0.1300 0.7140 0.2322 O0.5840 0.8700 0.5655 O0.4150 0.2538 0.6140 O0.7462 0.1612 0.2807 O0.8388 0.5850 0.9473 O0.2538 0.4150 0.3860 O0.1612 0.7462 0.7193 O0.5850 0.8388 0.0526 K_POINTS automatic 4 4 3 0 0 0 Program PWSCF v.4.3.2 starts on 11Aug2011 at 11:58: 4 This program is part of the open-source Quantum ESPRESSO suite for quantum simulation of materials; please cite "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009); URL http://www.quantum-espresso.org;, in publications or presentations arising from this work. More details at http://www.quantum-espresso.org/wiki/index.php/Citing_Quantum-ESPRESSO Parallel version (MPI), running on 3 processors R & G space division: proc/pool =3 Current dimensions of program PWSCF are: Max number of different atomic species (ntypx) = 10 Max number of k-points (npk) = 4 Max angular momentum in pseudopotentials (lmaxx) = 3 Waiting for input... Reading input from stdin Warning: card ignored Warning: card ION_DYNAMICS = NONE , ignored Warning: card / ignored Warning: cardignored Warning: card CELL_DYNAMICS = NONE , ignored Warning: card PRESS = 0 , ignored Warning: card PRESS_CONV_THR = 0.05D0 , ignored Warning: card / ignored Subspace diagonalization in iterative solution of the eigenvalue problem: a serial algorithm will be used Parallelization info sticks: dense smooth PW G-vecs:dense smooth PW Min 826 41211673322259723822 Max 827 41311773331259863827 Sum24791237349 21997977939 11473 bravais-lattice index =4 lattice parameter (alat) = 9.4718 a.u. unit-cell volume =1628.2022 (a.u.)^3 number of atoms/cell = 18 number of atomic types=3 number of electrons =96.00 number of Kohn-Sham states= 48 kinetic-energy cutoff = 50. Ry charge density cutoff = 400. Ry convergence threshold = 1.0E-10 mixing beta = 0.7000 number of iterations used =8 plain mixing Exchange-correlation = SLA PW PBE PBE ( 1 4 3 4 0) EXX-fraction =0.00 celldm(1)= 9.471757 celldm(2)= 0.00 celldm(3)= 2.212514 celldm(4)= 0.00 celldm(5)= 0.00 celldm(6)= 0.00 crystal axes: (cart. coord. in units of alat) a(1) = ( 1.00 0.00 0.00 ) a(2) = ( -0.50 0.866025 0.00 ) a(3) = ( 0.00 0.00 2.212514 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( 1.00 0.577350 0.00 )
[Pw_forum] compiling QHA
Dear all, When I compile QHA in espresso-4.3, the following error occur: [wang at solei QHA]$ ./Compile ifort -O3 -c tetra.f ifort -O3 -c k_brillouin.f ifort -O3 -c generate_tetra.f ifort -O3 -c det3.f ifort -O3 -c det4.f ifort -static -o tetra.x tetra.o k_brillouin.o generate_tetra.o det3.o det4.o ld: cannot find -lm ld: cannot find -lpthread ld: cannot find -lc ld: cannot find -ldl ld: cannot find -lc make: *** [tetra] Error 1 ifort -O3 -c Debye.f90 ifort -O3 -c Debye_T.f ifort -O3 -c debye3.f ifort -O3 -c cheval.f ifort -O3 -c d1mach.f ifort -static -o Debye.x Debye.o Debye_T.o debye3.o cheval.o d1mach.o ld: cannot find -lm ld: cannot find -lpthread ld: cannot find -lc ld: cannot find -ldl ld: cannot find -lc make: *** [Debye_x] Error 1 ifort -FR -c Mean_square_displacement.f90 ifort -static -o Mean_square_displacement.x Mean_square_displacement.o ld: cannot find -lm ld: cannot find -lpthread ld: cannot find -lc ld: cannot find -ldl ld: cannot find -lc make: *** [MSD] Error 1 Could you please help me ? Thank you very much. WANG Riping 2011.6.22 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110622/2bbc6497/attachment-0001.htm
[Pw_forum] atomic positions in xsf from cppp.x
Dear Forum, when I use cppp.x to generate an xsf file for an vc-cp output file with nstep=1000, there will be cell vectors and atomic positions in xsf file. what is this cell vectors and atomic positions stand for? For average structure configuration of 1000 steps, or for average structure configuration of one period, just for the last configuration of the this calculation? Part of input file are following: calculation = 'vc-cp' , prefix = 'SiO2' , restart_mode = 'restart' , nstep = 1000 , ndr = 51 , ndw = 52 , Thank you very much. WANG Riping 2011.3.8 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110308/2c5abb9d/attachment.htm
[Pw_forum] ecfixed makes the volume different?
Dear Forum, Could anyone give me some hints about how to adjust "qcutz,q2sigma,ecfixed" in vc-cp? I have some difficulties. After wave function initialization, the input file are for zero pressure zero temperature are following. (After zero pressure zero temperature, I will elevate temperature with verlet algorithm.) But when I use different ecfixed with fixed qcutz and q2sigma, the volumes at equilibrium are different. some is close to that without "qcutz,q2sigma,ecfixed" used. some are 100 ( toal around 3000) different from that without "qcutz,q2sigma,ecfixed" used. (In my case, I am sure the fnosep, wmass, and emass is good enough.) I attached the voluem evolution figure. In the figure, from step 200-1200 is belong to the following input file. before 200 step is wave initialization, after 1200 step, I increased temperature. Thank you very much! WANG Riping 2011.3.7 calculation = 'vc-cp' , prefix = 'SiO2-mog' , restart_mode = 'restart' , nstep = 1000 , iprint = 1 , isave = 10 , dt = 5.0 , ndr = 51 , ndw = 52 , tstress = .TRUE. , tprnfor = .TRUE. , saverho = .TRUE. , disk_io = 'high' , /ekin_conv_thr = 1.0d-3 , /etot_conv_thr = 5.0d-3 , /forc_conv_thr = 1.0d-2 , pseudo_dir ='~/espresso/pseudo' , outdir = './' , / ibrav = 14 , celldm(1) = 16.021289268 , celldm(2) = 0.558448710 , celldm(3) = 1.245135500 , celldm(4) = -0.28830 , celldm(5) = -0.007432130 , celldm(6) = 0.64680 , nat = 36 , ntyp = 2 , ecutwfc = 30 , ecutrho = 240.0 , nr1b = 16 , nr2b = 16 , nr3b = 16 , qcutz = 150.0 , q2sigma = 2 , ecfixed = 16.8 , / electron_dynamics = 'sd' , emass = 400 , emass_cutoff = 3. , / ion_dynamics = 'sd' , / cell_dynamics = 'pr' , press = 0 , wmass = 300 , / ATOMIC_SPECIES O 16.00 O.pbe-van_ak.UPF Si 28.00 Si.pbe-n-van.UPF ATOMIC_POSITIONS (crystal) ... -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110307/172de12c/attachment.htm -- next part -- A non-text attachment was scrubbed... Name: volume.p Type: application/octet-stream Size: 3736 bytes Desc: not available Url : http://www.democritos.it/pipermail/pw_forum/attachments/20110307/172de12c/attachment.obj
[Pw_forum] atomic position in vc-relax CPMD
Dear Forum, How can I get average atomic positions at at given pressure and temperature with vc-relax QE-cpmd? I set the following control model: calculation = 'vc-cp' , prefix = 'SiO2' , restart_mode = 'restart', nstep = 2000 , iprint = 1 , isave = 1 , dt = 5.0 , ndr = 54 , ndw = 54 , tstress = .TRUE. , tprnfor = .TRUE. , saverho = .TRUE. , disk_io = 'high' , But it print atomic position information in output file every 10 steps. I hope it print c information every 1 step. and then I can make average position by matlab. I am not sure about this method good or not. Could any one give comment? Thank you WANG Riping 2011.3.4 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110304/a50e804b/attachment.htm
[Pw_forum] nstep for CPMD average physical quantities
To add information, I use vc-relax, and I want to get average quantities at certain pressure and temperature. WANG Riping 2011.3.4 On 4 March 2011 09:48, Riping WANG wrote: > Dear Forum, > > For QE-CPMD calcualtion, Could any told me how many nstep keeping is > necessary to get average physical quantities? > Thank you very much. > > WANG Riping > 20101.3.4 > > > -- > > ** > WANG Riping > Ph.D student, > Institute for Study of the Earth's Interior,Okayama University, > 827 Yamada, Misasa, Tottori-ken 682-0193, Japan > Tel: +81-858-43-3739(Office), 1215(Inst) > E-mail: wang.riping.81 at gmail.com > > ** > > > > -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110304/f5454c02/attachment.htm
[Pw_forum] nstep for CPMD average physical quantities
Dear Forum, For QE-CPMD calcualtion, Could any told me how many nstep keeping is necessary to get average physical quantities? Thank you very much. WANG Riping 20101.3.4 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110304/fe9c9d01/attachment.htm
[Pw_forum] temph and fnoseh in cp.x
Dear forum, Is it necessary to control temph in cp.x. At first, I did not control temph, and ignored fnoseh, I can keep the phonon temperature at target value very well. But later, I try to set fnoseh to control temph, it breaks during elevating phonon temperature. So, is anyone can give me some guide about temph and fnoseh? Is any good reference paper? Thanks very much. WANG Riping 2010.2.4 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110204/e717fa40/attachment-0001.htm
[Pw_forum] Pw_forum Digest, Vol 42, Issue 43
Der prof. Paolo, Thanks very much to your reply. 1. Re: ion_radius (Paolo Giannozzi) WANG Riping 2010.12.27 On 25 December 2010 12:36, wrote: > Send Pw_forum mailing list submissions to >pw_forum at pwscf.org > > To subscribe or unsubscribe via the World Wide Web, visit >http://www.democritos.it/mailman/listinfo/pw_forum > or, via email, send a message with subject or body 'help' to >pw_forum-request at pwscf.org > > You can reach the person managing the list at >pw_forum-owner at pwscf.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Pw_forum digest..." > > > Today's Topics: > > 1. Re: ion_radius (Paolo Giannozzi) > 2. how does the calculation scales with number of processors > (mohnish pandey) > 3. Re: how does the calculation scales with number ofprocessors > (Stefano de Gironcoli) > 4. Re: how does the calculation scales with number ofprocessors > (mohnish pandey) > 5. Re: how does the calculation scales with number ofprocessors > (Paolo Giannozzi) > 6. Merry Christmas (jiayudai) > 7. About the difference between the energy obtained from > geometry optimization and single energy calculation. (Hongsheng Zhao) > 8. About the relationship between enthalpy and pressure. > (Hongsheng Zhao) > 9. Re: About the relationship between enthalpy and pressure. > (Stefano de Gironcoli) > 10. Re: About the difference between the energy obtained from > geometry optimization and single energy calculation. (Prasenjit Ghosh) > > > -- > > Message: 1 > Date: Fri, 24 Dec 2010 09:56:00 +0100 > From: Paolo Giannozzi > Subject: Re: [Pw_forum] ion_radius > To: PWSCF Forum > Message-ID: <5B32AABF-A8DA-4C29-8850-9250578A7450 at democritos.it> > Content-Type: text/plain; charset=US-ASCII; delsp=yes; format=flowed > > > On Dec 24, 2010, at 7:56 , Riping WANG wrote: > > > In example18, for cp.x practice. I found the " ion_radius(1) > > and ion_radius(2) " in the input file. > > Could you please tell me something about this input parameter, and > > how to set? > > usuallly you do not need to set it: the default value should be fine. > ion_radius is the width of a gaussian distribution of charge replacing > the point charge of the (pseudo-)nucleus. It should be big enough > to ensure convergence of the total energy in G-space, small enough > to yield negligible overlap between gaussians. Valus around 0.8 are > usually fine. The total energy and forces should not change if > ion_radius changes. > > P. > --- > Paolo Giannozzi, Dept of Chemistry, Univ. Udine > via delle Scienze 208, 33100 Udine, Italy > Phone +39-0432-558216, fax +39-0432-558222 > > > > > > > -- > > Message: 2 > Date: Fri, 24 Dec 2010 15:12:19 +0530 > From: mohnish pandey > Subject: [Pw_forum] how does the calculation scales with number of >processors > To: PWSCF Forum > Message-ID: ><AANLkTimzcxdX+zRc4pHNYUQUU8Ph+qG7aiwggG=XyEnh at mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Dear QE users, > I am trying to see how the runtime scales with > number of processors. I did the same calculation using one node with eight > processors on one cluster and 8 nodes with total 56 processors on other > cluster but the time does not seem to scale linearly with the number of > processors. Can anybody give me an idea how does the time scale with number > of processors. > Thanks a lot in advance. > > -- > Regards, > MOHNISH, > - > Mohnish Pandey > Y6927262,5th Year dual degree student, > Department of Chemical Engineering, > IIT KANPUR, UP, INDIA > - > -- next part -- > An HTML attachment was scrubbed... > URL: > http://www.democritos.it/pipermail/pw_forum/attachments/20101224/094c0e6d/attachment-0001.htm > > -- > > Message: 3 > Date: Fri, 24 Dec 2010 11:20:34 +0100 > From: Stefano de Gironcoli > Subject: Re: [Pw_forum] how does the calculation scales with number of >processors > To: PWSCF Forum > Message-ID: <4D1473F2.9090302 at sissa.it> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > It strongly depends on the particular calculation and on the > parallelization strategy you choose. Linear scaling would be the ideal &g
[Pw_forum] ion_radius
Dear Forum, I am using espresso-4.2.1. In example18, for cp.x practice. I found the " ion_radius(1) and ion_radius(2) " in the input file. Could you please tell me something about this input parameter, and how to set? as following: ion_dynamics='verlet' , ion_temperature='nose' , tempw = .2 , fnosep = 250 , ion_radius(1) = 1.0 , ion_radius(2) = 1.0 , I can not find any explanation about ion_radius from espresso manual. Thanks very much. WANG Riping 2010.12.24 -- ** WANG Riping Ph.D student, Institute for Study of the Earth's Interior,Okayama University, 827 Yamada, Misasa, Tottori-ken 682-0193, Japan Tel: +81-858-43-3739(Office), 1215(Inst) E-mail: wang.riping.81 at gmail.com ** -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20101224/f278e2a2/attachment.htm