aselimov
/
CACmb
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixes to parse_pos subroutine

Browse Source
master
Alex Selimov 5 years ago
parent 69ab46b6e3
commit 1db2b24499
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File

5
src/subroutines.f90
Unescape View File

@ -194,11 +194,12 @@ module subroutines
pos = box_bd(2*i) - pos pos = box_bd(2*i) - pos
else if ((index(pos_string,'+') > 0).and.(index(pos_string,'inf')>0)) then else if ((index(pos_string,'+') > 0).and.(index(pos_string,'inf')>0)) then
!Now extract the number we are reducing from infinity !Now extract the number we are reducing from infinity
if(index(pos_string,'inf') < index(pos_string,'-')) then if(index(pos_string,'inf') < index(pos_string,'+')) then
read(pos_string(index(pos_string,'+')+1:), *, iostat=iospara) pos read(pos_string(index(pos_string,'+')+1:), *, iostat=iospara) pos
else else
read(pos_string(1:index(pos_string,'+')-1), *, iostat=iospara) pos read(pos_string(1:index(pos_string,'+')-1), *, iostat=iospara) pos
end if end if
pos = box_bd(2*i-1) + pos
else if ((index(pos_string,'*') > 0).and.(index(pos_string,'inf')>0)) then else if ((index(pos_string,'*') > 0).and.(index(pos_string,'inf')>0)) then
!Now extract the number we are reducing from infinity !Now extract the number we are reducing from infinity
if(index(pos_string,'inf') < index(pos_string,'*')) then if(index(pos_string,'inf') < index(pos_string,'*')) then
@ -206,7 +207,7 @@ module subroutines
else else
read(pos_string(1:index(pos_string,'*')-1), *, iostat=iospara) pos read(pos_string(1:index(pos_string,'*')-1), *, iostat=iospara) pos
end if end if
pos = (box_bd(2*i)-box_bd(2*i-1))*pos pos = (box_bd(2*i)-box_bd(2*i-1))*pos + box_bd(2*i-1)
else else
read(pos_string, *, iostat=iospara) pos read(pos_string, *, iostat=iospara) pos
end if end if

Write Preview
Loading…
Cancel
Save