Commit 4fcba5c5 authored by Daniel Wortmann's avatar Daniel Wortmann

Bugfixes for calls to new sorting routine

parent 97a2aae5
...@@ -61,7 +61,7 @@ CONTAINS ...@@ -61,7 +61,7 @@ CONTAINS
END DO END DO
gvacl(n2) = SQRT(REAL(gvac(1)**2+gvac(2)**2)) gvacl(n2) = SQRT(REAL(gvac(1)**2+gvac(2)**2))
ENDDO k_loop ENDDO k_loop
CALL sort(gindex,gvacl) CALL sort(gindex(:n2),gvacl)
DO j = 1,n2 DO j = 1,n2
! gvac1d, gvac2d are now ordered by increasing length ! gvac1d, gvac2d are now ordered by increasing length
gvac1d(j)=gvac1(gindex(j)) gvac1d(j)=gvac1(gindex(j))
......
...@@ -179,7 +179,7 @@ CONTAINS ...@@ -179,7 +179,7 @@ CONTAINS
END DO END DO
END DO END DO
CALL sort(index,e) CALL sort(index(:n),e)
! Check if no deep eigenvalue is found ! Check if no deep eigenvalue is found
IF (e_min-MINVAL(e(1:n))>1.0) THEN IF (e_min-MINVAL(e(1:n))>1.0) THEN
......
...@@ -233,7 +233,7 @@ CONTAINS ...@@ -233,7 +233,7 @@ CONTAINS
& (stars%mx2+stars%kv3(2,k))*(2*stars%mx1+1) +& & (stars%mx2+stars%kv3(2,k))*(2*stars%mx1+1) +&
& (stars%mx3+stars%kv3(3,k))*(2*stars%mx1+1)*(2*stars%mx2+1) & (stars%mx3+stars%kv3(3,k))*(2*stars%mx1+1)*(2*stars%mx2+1)
ENDDO ENDDO
CALL sort(index,stars%sk3,gsk3) CALL sort(index(:stars%ng3),stars%sk3,gsk3)
ALLOCATE (ig2p(stars%ng3)) ALLOCATE (ig2p(stars%ng3))
...@@ -634,7 +634,7 @@ CONTAINS ...@@ -634,7 +634,7 @@ CONTAINS
& (stars%mx2+stars%kv3(2,k))*(2*stars%mx1+1) +& & (stars%mx2+stars%kv3(2,k))*(2*stars%mx1+1) +&
& (stars%mx3+stars%kv3(3,k))*(2*stars%mx1+1)*(2*stars%mx2+1) & (stars%mx3+stars%kv3(3,k))*(2*stars%mx1+1)*(2*stars%mx2+1)
ENDDO ENDDO
CALL sort(index,stars%sk3,gsk3) CALL sort(index(:stars%ng3),stars%sk3,gsk3)
DO k = 1,stars%ng3 DO k = 1,stars%ng3
kv3rev(k,1) = stars%kv3(1,INDEX(k)) kv3rev(k,1) = stars%kv3(1,INDEX(k))
kv3rev(k,2) = stars%kv3(2,INDEX(k)) kv3rev(k,2) = stars%kv3(2,INDEX(k))
......
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