2022-02-07 09:45:08 +11:00
|
|
|
program movebloubs
|
|
|
|
|
|
|
|
|
|
use bloubspace
|
|
|
|
|
use povstuff
|
|
|
|
|
|
|
|
|
|
implicit none
|
|
|
|
|
|
|
|
|
|
character(200) :: infile, outfile
|
|
|
|
|
integer :: inu, outu, errcode, i
|
|
|
|
|
type(t_bloubs) :: bloub
|
|
|
|
|
double precision :: bx, by, bz
|
2022-02-08 12:53:49 +11:00
|
|
|
logical :: add_new_bloub = .FALSE.
|
2022-02-07 09:45:08 +11:00
|
|
|
|
|
|
|
|
i = IARGC()
|
|
|
|
|
if (i .ne. 2) then
|
|
|
|
|
STOP ": BAD ARGS ON COMMAND LINE"
|
|
|
|
|
endif
|
|
|
|
|
call getarg(1, infile)
|
|
|
|
|
call getarg(2, outfile)
|
|
|
|
|
|
|
|
|
|
inu = 42 ; outu = 51
|
|
|
|
|
|
|
|
|
|
write (0, '(A)') &
|
|
|
|
|
"*** moving bloubs from "//trim(infile)//" to "//trim(outfile)
|
|
|
|
|
|
|
|
|
|
open(unit=inu, file=trim(infile), form='unformatted', &
|
|
|
|
|
iostat=errcode, &
|
|
|
|
|
action='read', status='old')
|
|
|
|
|
if (0 .ne. errcode) then
|
|
|
|
|
STOP " : CAN'T OPEN FILE " // trim(infile)
|
|
|
|
|
endif
|
|
|
|
|
|
|
|
|
|
open(unit=outu, file=trim(outfile), form='unformatted', &
|
|
|
|
|
iostat=errcode, &
|
|
|
|
|
action='write', status='replace')
|
|
|
|
|
if (0 .ne. errcode) then
|
2022-02-08 12:53:49 +11:00
|
|
|
STOP " : CAN'T OPEN " // trim(outfile) // "FOR WRITE"
|
2022-02-07 09:45:08 +11:00
|
|
|
endif
|
|
|
|
|
|
|
|
|
|
bx = 0.0; by = 0.0; bz = 0.0
|
|
|
|
|
|
|
|
|
|
do
|
|
|
|
|
read (unit=inu, iostat=errcode) bloub
|
|
|
|
|
if (0 .ne. errcode) then
|
|
|
|
|
exit
|
|
|
|
|
endif
|
|
|
|
|
|
2022-02-08 12:53:49 +11:00
|
|
|
! moving and boundingboxing
|
|
|
|
|
call move_bloub (bloub, 0.11)
|
|
|
|
|
call bound_a_blob(bloub)
|
2022-02-07 09:45:08 +11:00
|
|
|
|
2022-02-08 12:53:49 +11:00
|
|
|
! calcul du barycentre
|
2022-02-07 09:45:08 +11:00
|
|
|
bx = bx + bloub%px
|
|
|
|
|
by = by + bloub%py
|
|
|
|
|
bz = bz + bloub%pz
|
|
|
|
|
|
2022-02-08 12:53:49 +11:00
|
|
|
write(outu, iostat=errcode) bloub ! no error control ?
|
|
|
|
|
if (0 .ne. errcode) then
|
|
|
|
|
STOP " : WRITE ERROR TO " // trim(outfile)
|
2022-02-07 12:08:17 +11:00
|
|
|
endif
|
2022-02-07 09:45:08 +11:00
|
|
|
|
|
|
|
|
enddo
|
|
|
|
|
|
2022-02-08 12:53:49 +11:00
|
|
|
if (add_new_bloub) then
|
|
|
|
|
! and now, we inject a new bloub in the universe
|
|
|
|
|
bloub%nick = 'newbie '
|
|
|
|
|
bloub%alive = .TRUE.
|
|
|
|
|
call random_pv(bloub)
|
|
|
|
|
bloub%radius = 0.042
|
|
|
|
|
bloub%seq = 0
|
|
|
|
|
!
|
|
|
|
|
! where is the bug ?
|
|
|
|
|
!
|
|
|
|
|
call display_bloub (bloub, "new bloub")
|
|
|
|
|
write(outu) bloub ! no error control ?
|
|
|
|
|
endif
|
2022-02-07 09:45:08 +11:00
|
|
|
|
2022-02-08 12:53:49 +11:00
|
|
|
! ther was some strange bias in this data
|
|
|
|
|
! may be a random not symetric around 0.0 ?
|
|
|
|
|
write (0, '(A,3(F11.3,3X))') "barycentre : ", bx, by, bz
|
2022-02-07 09:45:08 +11:00
|
|
|
|
|
|
|
|
close(inu) ; close(outu)
|
|
|
|
|
|
|
|
|
|
end program
|
