Skip to content
Snippets Groups Projects
Commit fd33d5b7 authored by Wiarda, Dorothea's avatar Wiarda, Dorothea
Browse files

Since Maxkda is no longer used, don't compare against it

parent 24ac1b53
No related branches found
No related tags found
No related merge requests found
......@@ -327,23 +327,14 @@ C *** Determine the size of the adjustment to be Made
IF (X1(I).LT.0.1667D0) A = X21(I)*2.0D0
IF (X1(I).GE.0.1667D0) A = X21(I)
call reallocate_real_data(Energb, Ii, 100)
Energb(Ii) = X2(I+1)
IF (Ii.GT.Maxkda) STOP '[#6 in Adjust_Auxil in dat/mdat5]'
Energb(Ii) = X2(I+1)
Ii = Ii + 1
C
IF (X1(I).LT.0.1667D0) THEN
C *** Here we need to subdivide region i several times
200 CONTINUE
call reallocate_real_data(Energb, Ii, 100)
Energb(Ii) = Energb(Ii-1) + A
IF (Ii.GT.Maxkda) THEN
DO M=Ii-10,Ii
De = Energb(M) - Energb(M-1)
WRITE (6,40000) M, Energb(M), De
40000 FORMAT(' M, E, De=', I8, 1P2G20.10)
END DO
STOP '[#7 in Adjust_Auxil in dat/mdat5]'
END IF
Energb(Ii) = Energb(Ii-1) + A
Ii = Ii + 1
A = A*2.0D0
IF (Energb(Ii-1)+A .LT. X2(I+2)) GO TO 200
......@@ -363,9 +354,7 @@ C
END IF
C
220 CONTINUE
Kdatb = Ii - 1
IF (Kdatb.GT.Maxkda) WRITE (6,20100) Kdatb, Maxkda
20100 FORMAT (/' ########### WARNING ######### Kdatb, Maxkda =', 2I10)
Kdatb = Ii - 1
Size = B
C
C
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment