Mercurial > op > rk_sakamoti

comparison rkmain.f90 @ 5:0ca8347e4c47

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
add rk_Lorenz.R and rksub.f90
author "uncorrelated zombie" <uncorrelated@yahoo.co.jp>
date Tue, 23 Jul 2024 03:31:21 +0900
parents 2347cff808b0
children
comparison
equal deleted inserted replaced
4:8026ceb4de7a 5:0ca8347e4c47
1 program rkmain 1 program rkmain
2 use rkmod 2 use rkmod
3 implicit none 3 implicit none
4 real(kind=dpkind),dimension(3) :: x = [0.5, 0.5, 0.5] 4 real(kind=dpkind),dimension(3) :: x = [0.5, 0.5, 0.5]
5 integer :: max_i = 5000 5 integer :: max_i = 5000
6 real(kind=dpkind), allocatable, dimension(:, :) :: r 6 ! real(kind=dpkind), allocatable, dimension(:, :) :: r
7 real(kind=dpkind), dimension(5000, 4) :: r
7 8
8 allocate(r(max_i, 4)) 9 ! allocate(r(max_i, 4))
9 call calc(x, max_i, r) 10 call calc(x, max_i, r)
10 print *, r 11 open(17,file='rkmain.dat', form='unformatted', access = "stream", status='replace')
12 write(17) r
13 close(17)
11 end program 14 end program
12 15