the problem is that How can i increase a (Fortran) allocated array size in a program?
in program like this:
... real(8), allocatable, dimension(:) :: x ... ! programs body N=100 allocate(x(N)) ... N=150 allocate(x(N)) ! but this is an error
there is two approach
First: using move_alloc
Second: using x=[x,(x(1),i=1,j-size(x))]