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

Add random position

Browse Source
This commit is contained in:
Alex Selimov 2020-01-28 11:35:14 -05:00
parent e91bcd5d1e
commit 282dcd3c58
2 changed files with 7 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/main.f90
View File

@ -23,6 +23,8 @@ program main
!Call initialization functions !Call initialization functions
call lattice_init call lattice_init
call box_init call box_init
call random_seed
end_mode_arg = 0 end_mode_arg = 0
! Command line parsing ! Command line parsing

5
src/subroutines.f90
View File

@ -178,6 +178,7 @@ module subroutines
real(kind=dp), intent(out) :: pos !The output parsed position value real(kind=dp), intent(out) :: pos !The output parsed position value
integer :: iospara integer :: iospara
real(kind=dp) :: rand
iospara = 0 iospara = 0
if(trim(adjustl(pos_string)) == 'inf') then if(trim(adjustl(pos_string)) == 'inf') then
@ -208,6 +209,10 @@ module subroutines
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 + box_bd(2*i-1) pos = (box_bd(2*i)-box_bd(2*i-1))*pos + box_bd(2*i-1)
else if (trim(adjustl(pos_string)) == 'rand') then
call random_number(rand)
pos = (box_bd(2*i)-box_bd(2*i-1))*rand + box_bd(2*i-1)
else else
read(pos_string, *, iostat=iospara) pos read(pos_string, *, iostat=iospara) pos
end if end if