tTh/Fortraneries
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

non working code inside

Browse Source
This commit is contained in:
tth 2022-03-18 22:08:08 +01:00
parent d987ff1685
commit 6d51ad3947
1 changed files with 14 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
BloubWorld/mergebloubs.f90
View File

@ -1,5 +1,9 @@
program mergebloubs program mergebloubs
!-------------------------------------------!
! THIS IS A BIG MESS OF BUGS !
!-------------------------------------------!
use bloubspace use bloubspace
implicit none implicit none
@ -21,7 +25,8 @@ program mergebloubs
call getarg(1, infile) call getarg(1, infile)
call getarg(2, outfile) call getarg(2, outfile)
write(0, '(2A20, I8)') trim(infile), trim(outfile), NB_MAX_BLOUBS write(0, '(A, 2A20, I8)') "*** mergebloubs ", &
trim(infile), trim(outfile), NB_MAX_BLOUBS
! --------------- allocate memory for the people ! --------------- allocate memory for the people
@ -61,6 +66,11 @@ program mergebloubs
STOP " : ERR READING FIRST BLOUB" STOP " : ERR READING FIRST BLOUB"
endif endif
call display_bloub (bloub, "first bloub") call display_bloub (bloub, "first bloub")
write(outu, iostat=errcode) bloub
if (0 .ne. errcode) then
STOP " : FIRST BLOUB, WRITE ERROR TO " // trim(outfile)
endif
les_bloubs(idx) = bloub les_bloubs(idx) = bloub
idx = idx + 1 idx = idx + 1
@ -69,7 +79,7 @@ program mergebloubs
do ! infinite loop do ! infinite loop
print *, "============ PASS ", idx ! print *, "============ PASS ", idx
if (idx .EQ. NB_MAX_BLOUBS) then if (idx .EQ. NB_MAX_BLOUBS) then
write(0, '(I8, A)') idx, " max number of bloubs reached" write(0, '(I8, A)') idx, " max number of bloubs reached"
exit exit
@ -112,7 +122,7 @@ program mergebloubs
les_bloubs(idx) = bloub les_bloubs(idx) = bloub
endif endif
write(outu, iostat=errcode) bloub ! no error control ? write(outu, iostat=errcode) bloub
if (0 .ne. errcode) then if (0 .ne. errcode) then
STOP " : WRITE ERROR TO " // trim(outfile) STOP " : WRITE ERROR TO " // trim(outfile)
endif endif
@ -149,7 +159,7 @@ contains
blr%vy = (bla%vy + blb%vy) / 2.0 blr%vy = (bla%vy + blb%vy) / 2.0
blr%vz = (bla%vz + blb%vz) / 2.0 blr%vz = (bla%vz + blb%vz) / 2.0
blr%radius = (bla%radius + blb%radius) / 1.414 blr%radius = (bla%radius + blb%radius) / 2.718
blr%age = min(bla%age, blb%age) blr%age = min(bla%age, blb%age)
! bring it to life ! ! bring it to life !