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

Bourtoulots laboratories in the place

Browse Source
This commit is contained in:
tTh
2026-05-31 23:44:58 +02:00
parent c3f27c9843
commit 8e14ec2c8d
6 changed files with 12 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
Fraktalism/voxelize.f90
View File

@@ -9,7 +9,7 @@ program voxelize
implicit none
integer, parameter :: DIMC = 320
integer, parameter :: DIMC = 451
integer, dimension(:,:,:), allocatable :: cube
type(t_point3d), dimension(:), allocatable :: points
integer :: errcode, foo, argc
@@ -19,7 +19,7 @@ program voxelize
double precision :: dmaxcube, delta
character(200) :: filename, string
write(0, *) "--- start of voxelize"
write(0, *) "--------- start of voxelize -------------------"
argc = IARGC()
if (2 .NE. argc) then
@@ -43,7 +43,7 @@ program voxelize
KA(1) = -1.3402 ; KA(2) = 1.5245
KA(3) = 1.0966 ; KA(4) = -2.3423
KB(1) = -1.2100 ; KB(2) = 1.3685
KB(1) = -1.2333 ; KB(2) = 1.3685
KB(3) = 1.3237 ; KB(4) = -2.3992
call interp4dp(KA, KB, KM, delta)
write(0, "(' --- coefs = ', 4F11.6)") KM
@@ -80,7 +80,7 @@ subroutine fcoor2icoor(in, out)
integer :: outvalue
invalue = (in + 2.0) / 2.0
outvalue = int(invalue * real(DIMC/2))
outvalue = int(invalue * real(DIMC/2.0))
! add molly-guard here
out = outvalue