Fortraneries/GravityField/essai.f90

83 lines
2.2 KiB
Fortran
Raw Normal View History

!-----------------------------------------------------------------------
2022-11-30 03:24:24 +01:00
!-
! test program for the project "gravity field"
!-
!-----------------------------------------------------------------------
2022-11-24 00:52:13 +01:00
program essai
use realfield
use spitpgm ! XXX
use pixrgb
2022-11-24 00:52:13 +01:00
implicit none
2022-11-30 22:46:44 +01:00
call init_random()
2022-11-24 00:52:13 +01:00
2023-02-11 17:00:58 +01:00
call essai_near_planet(2048, 2048)
2022-11-28 13:47:44 +01:00
2022-11-30 22:46:44 +01:00
STOP 'BECAUSE YOLO'
2022-11-24 00:52:13 +01:00
contains
!-----------------------------------------------------------------------
!-
2022-12-12 00:57:00 +01:00
! computation of thr nearest planet
2022-11-30 22:46:44 +01:00
!-
2022-12-12 00:57:00 +01:00
subroutine essai_near_planet(nbplanets, szfield)
integer, intent(in) :: nbplanets, szfield
type(t_pixrgb), dimension(:,:), allocatable :: cmap
integer :: ix, iy
real :: fx, fy, dx, dy
integer :: near, ipl, errcode
real :: curdist, smalldist
type(massbody) :: planets(nbplanets)
2022-12-12 00:57:00 +01:00
print *, "near planets test", nbplanets, szfield
2022-11-30 22:46:44 +01:00
allocate(cmap(szfield, szfield), stat=errcode)
! map = -1
2022-12-12 00:57:00 +01:00
! create some random bodies
do ipl=1, nbplanets
planets(ipl)%posx = rand() * szfield
planets(ipl)%posy = rand() * szfield
planets(ipl)%serial = ipl
end do
! call save_bodies_to_txt_file(planets, "planets.txt")
! loop over all the location of the field
do ix=1, szfield
fx = real(ix)
do iy=1, szfield
fy = real(iy)
near = -1
smalldist = 1e37
! loop over all the planet's bodies
do ipl=1, nbplanets
! compute the "fake" distance
dx = fx - planets(ipl)%posx
dy = fy - planets(ipl)%posy
curdist = (dx*dx) + (dy*dy)
if (curdist .LT. smalldist) then
near = ipl
smalldist = curdist
endif
end do ! loop on ipl
2022-12-16 21:16:02 +01:00
cmap(ix, iy)%r = mod(near*4, 255)
cmap(ix, iy)%g = mod(near*7, 255)
cmap(ix, iy)%b = mod(near*11, 255)
2022-12-12 00:57:00 +01:00
enddo
write(0, *) "row", ix, " on", szfield
2022-11-30 22:46:44 +01:00
enddo
2022-11-24 00:52:13 +01:00
2023-02-11 17:00:58 +01:00
call rgbpix_spit_as_pnm_8(cmap, "rgb.pnm")
2022-11-24 00:52:13 +01:00
2022-11-28 13:47:44 +01:00
end subroutine
!-----------------------------------------------------------------------
2022-11-24 00:52:13 +01:00
end program