Mercurial > op > rk_sakamoti

diff 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
line wrap: on
line diff
--- a/rkmain.f90	Tue Jul 23 01:48:31 2024 +0900
+++ b/rkmain.f90	Tue Jul 23 03:31:21 2024 +0900
@@ -3,10 +3,13 @@
   implicit none
   real(kind=dpkind),dimension(3) :: x = [0.5, 0.5, 0.5]
   integer :: max_i = 5000
-  real(kind=dpkind), allocatable, dimension(:, :) :: r
+!  real(kind=dpkind), allocatable, dimension(:, :) :: r
+  real(kind=dpkind), dimension(5000, 4) :: r
 
-  allocate(r(max_i, 4))
+!  allocate(r(max_i, 4))
   call calc(x, max_i, r)
-  print *, r
+  open(17,file='rkmain.dat', form='unformatted', access = "stream", status='replace')
+  write(17) r
+  close(17)
 end program