← Back to branch summary

~ubuntu-branches/ubuntu/oneiric/ifeffit/oneiric

« back to all changes in this revision

Viewing changes to src/feff6/M/totale.f

  • Committer: Bazaar Package Importer
  • Author(s): Alessio Treglia
  • Date: 2009-03-08 23:15:54 UTC
  • mfrom: (5.1.1 sid)
  • Revision ID: james.westby@ubuntu.com-20090308231554-pubym5d5by779fpn
Tags: 2:1.2.11d-1ubuntu1
* Merge from debian unstable, remaining changes:
  - Use --install-layout=deb option to build with Python 2.6.
* Update maintainer as per spec.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/feff6/M/totale.f
 
1
      subroutine totale (dval)
 
2
      implicit double precision (a-h,o-z)
 
3
      save
 
4
      include 'atom.h'
 
5
 
 
6
       external exchee
 
7
      dc(1)=1
 
8
      do 10 i=1,np
 
9
   10 dp(i)=d(i)/dr(i)
 
10
      if (nuc.le.0) go to 30
 
11
      do 20 i=1,nuc
 
12
   20 dp(i)=d(i)*(3-dr(i)*dr(i)/(dr(nuc)*dr(nuc)))/(dr(nuc)+dr(nuc))
 
13
      dc(1)=4
 
14
   30 call somm (dr,dp,dq,dpas,dc(1),0,np)
 
15
      dc(1)=-z*dc(1)
 
16
      do 40 i=1,np
 
17
      dp(i)=d(i)*dvf(i)
 
18
      dvn(i)=d(i)*dvn(i)
 
19
   40 d(i)=d(i)*exchee(d(i),dr(i))
 
20
      dc(2)=2
 
21
      dc(3)=1
 
22
      dc(5)=2
 
23
      if (nuc.ne.0) dc(3)=4
 
24
      call somm (dr,dp,dq,dpas,dc(3),0,np)
 
25
      call somm (dr,dvn,dq,dpas,dc(5),0,np)
 
26
      call somm (dr,d,dq,dpas,dc(2),0,np)
 
27
      dc(4)=dval-dc(3)
 
28
      dval=dval-half*dc(5)-dc(2)
 
29
      dc(2)=dc(3)-dc(1)-dc(5)-dc(2)
 
30
      dc(3)=half*dc(5)
 
31
      if (iprint .ge. 5)  write(16,50) dval,dc(4),dc(3),dc(2),dc(1)
 
32
   50 format (1h0,5x,'et=',1pe14.7,5x,'ec=',1pe14.7,5x,'ee=',1pe14.7,5x,
 
33
     1 'ex=',1pe14.7,5x,'en=',1pe14.7)
 
34
      return
 
35
      end