[Pw_forum] problem in xspectra calculations

Niharika Joshi joshiniharika20 at yahoo.in
Thu Oct 13 09:19:56 CEST 2011


Thank you Matteo sir for your reply.
Replacing the 'if' statement didn't work.
I'm using a 64-bit machine with composerxe-2011.4.191 compiler. And yes, I'm using pseudo of the example for carbon.
-Niharika Joshi
(project student,IISER Pune)
--- On Tue, 11/10/11, Matteo Calandra <matteo.calandra at impmc.jussieu.fr> wrote:

From: Matteo Calandra <matteo.calandra at impmc.jussieu.fr>
Subject: [Pw_forum]  problem in xspectra calculations
To: pw_forum at pwscf.org
Date: Tuesday, 11 October, 2011, 12:25 PM


> Message: 4
> Date: Tue, 11 Oct 2011 15:01:33 +0530 (IST)
> From: Niharika Joshi<joshiniharika20 at yahoo.in>
> Subject: Re: [Pw_forum] problem in xspectra calculations
> To: PWSCF Forum<pw_forum at pwscf.org>
> Message-ID:
>     <1318325493.18035.YahooMailClassic at web95715.mail.in.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hello Paolo Sir,
> Thank you for your reply.
> I tried running the xanes calculation after making the change in xspectra.f90. But it is still giving me the same message 'Could not find the element? C in the table of K edge energies!' at the end and leaving the output file
 incomplete.
> For some reason, the condition for 'if-loop' in mygetK.f90 is not getting satisfied and thus it is not getting executed. I wrote the variables sym and seuilK_tab(6) . For both variables it writes 'C' but does not enter the if-loop.
> -Niharika Joshi
> (project student, IISER Pune)



Dear Niharita,

   which compiler/Machine are you using ? It seems it treats in a 
strange way string variables. Moreover, are you using the pseudo
of the example for Carbon or another one ?

You can try the following two options:
1) in ligne 66 in mygetK.f90 change line
if(sym.eq.seuilK_tab(i)%name) then
in
if(sym.eq.trim(adjustl(seuilK_tab(i)%name))) then

2) if it doesn't work, then replace in xspectra.f90
all the lines with (two lines) of
e_1s=mygetK(upf(xiabs)%psd)
with the actual value of the Carbon K-edge energy, namely
with
e_1s=284.2

All the best,
 M.

-- 
* * * *
Matteo Calandra, Directeur de Recherche (DR2)
Institut de Mineralogie et de Physique des Milieux Condenses de Paris
Universite Pierre et Marie Curie, tour 23, 3eme etage, case 115
4 Place Jussieu, 75252 Paris Cedex 05 France
Tel: +33-1-44 27 52 16       Fax: +33-1-44 27 37 85
http://www.impmc.jussieu.fr/~calandra
_______________________________________________
Pw_forum mailing list
Pw_forum at pwscf.org
http://www.democritos.it/mailman/listinfo/pw_forum
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20111013/93882d13/attachment.htm 


More information about the Pw_forum mailing list