Commit 5a5622c8 authored by Matthias Redies's avatar Matthias Redies

use modern loops

parent f1323dfc
......@@ -203,14 +203,14 @@ CONTAINS
ENDDO
ENDIF
CALL intgr1(dens,rnot,h,n,a)
DO 110 i = 1,n
DO i = 1,n
rh(i) = dens(i)/rad(i)
110 CONTINUE
ENDDO
CALL intgr1(rh,rnot,h,n,b)
fisr = b(n)
DO 120 i = 1,n
DO i = 1,n
vr(i) = (a(i)+rad(i)* (fisr-b(i))-z)
120 CONTINUE
ENDDO
!+ta
DO ispin = 1, input%jspins
DO i = 1,n
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment