subroutine data_saveA1_R(pname,Y1) integer :: i,Np character(len=*) :: pname real(8),dimension(:) :: Y1 Np=size(Y1) open(free_unit(unit),file=reg(pname)) do i=1,Np write(unit,*)Y1(i) enddo close(unit) call file_gzip(reg(pname)) end subroutine data_saveA1_R subroutine data_saveA1_C(pname,Y1) integer :: i,Np character(len=*) :: pname complex(8),dimension(:) :: Y1 Np=size(Y1) open(free_unit(unit),file=reg(pname)) do i=1,Np write(unit,*)Y1(i) enddo close(unit) call file_gzip(reg(pname)) end subroutine data_saveA1_C !------------------------------------------------------------------! !------------------------------------------------------------------! !------------------------------------------------------------------! subroutine data_saveA2_R(pname,Y1,order,wspace) integer :: i,j,Ny1,Ny2 character(len=*) :: pname real(8),dimension(:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace Ny1=size(Y1,1) ; Ny2=size(Y1,2) open(free_unit(unit),file=reg(pname)) select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i=1,Ny1 do j=1,Ny2 write(unit,*)Y1(i,j) enddo if(wspace_)write(unit,*) enddo case ("C") do j=1,Ny2 do i=1,Ny1 write(unit,*)Y1(i,j) enddo if(wspace_)write(unit,*) enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA2_R subroutine data_saveA2_C(pname,Y1,order,wspace) integer :: i,j,Ny1,Ny2 character(len=*) :: pname complex(8),dimension(:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace Ny1=size(Y1,1) ; Ny2=size(Y1,2) open(free_unit(unit),file=reg(pname)) select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i=1,Ny1 do j=1,Ny2 write(unit,*)Y1(i,j) enddo if(wspace_)write(unit,*) enddo case ("C") do j=1,Ny2 do i=1,Ny1 write(unit,*)Y1(i,j) enddo if(wspace_)write(unit,*) enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA2_C !---------------------------- !---------------------------- !---------------------------- subroutine data_saveA3_R(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3 integer :: i1,i2,i3 character(len=*) :: pname real(8),dimension(:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace ! open(free_unit(unit),file=reg(pname)) ! Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) ! select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 write(unit,*)Y1(i1,i2,i3) enddo if(wspace_)write(unit,*) enddo enddo case ("C") do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3) enddo if(wspace_)write(unit,*) enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA3_R subroutine data_saveA3_C(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3 integer :: i1,i2,i3 character(len=*) :: pname complex(8),dimension(:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace ! open(free_unit(unit),file=reg(pname)) ! Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) ! select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 write(unit,*)Y1(i1,i2,i3) enddo if(wspace_)write(unit,*) enddo enddo case ("C") do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3) enddo if(wspace_)write(unit,*) enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA3_C !---------------------------- !---------------------------- !---------------------------- subroutine data_saveA4_R(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4 integer :: i1,i2,i3,i4 character(len=*) :: pname real(8),dimension(:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 write(unit,*)Y1(i1,i2,i3,i4) enddo if(wspace_)write(unit,*) enddo enddo enddo case ("C") do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4) enddo if(wspace_)write(unit,*) enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA4_R subroutine data_saveA4_C(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4 integer :: i1,i2,i3,i4 character(len=*) :: pname complex(8),dimension(:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 write(unit,*)Y1(i1,i2,i3,i4) enddo if(wspace_)write(unit,*) enddo enddo enddo case ("C") do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4) enddo if(wspace_)write(unit,*) enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA4_C !---------------------------- !---------------------------- !---------------------------- subroutine data_saveA5_R(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4,Ny5 integer :: i1,i2,i3,i4,i5 character(len=*) :: pname real(8),dimension(:,:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) Ny5=size(Y1,5) select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 do i5=1,Ny5 write(unit,*)Y1(i1,i2,i3,i4,i5) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo case ("C") do i5=1,Ny5 do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4,i5) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA5_R subroutine data_saveA5_C(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4,Ny5 integer :: i1,i2,i3,i4,i5 character(len=*) :: pname complex(8),dimension(:,:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) Ny5=size(Y1,5) select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 do i5=1,Ny5 write(unit,*)Y1(i1,i2,i3,i4,i5) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo case ("C") do i5=1,Ny5 do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4,i5) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA5_C !---------------------------- !---------------------------- !---------------------------- subroutine data_saveA6_R(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4,Ny5,Ny6 integer :: i1,i2,i3,i4,i5,i6 character(len=*) :: pname real(8),dimension(:,:,:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) Ny5=size(Y1,5) Ny6=size(Y1,6) ! select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 do i5=1,Ny5 do i6=1,Ny6 write(unit,*)Y1(i1,i2,i3,i4,i5,i6) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo case ("C") do i6=1,Ny6 do i5=1,Ny5 do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4,i5,i6) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA6_R subroutine data_saveA6_C(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4,Ny5,Ny6 integer :: i1,i2,i3,i4,i5,i6 character(len=*) :: pname complex(8),dimension(:,:,:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) Ny5=size(Y1,5) Ny6=size(Y1,6) ! select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 do i5=1,Ny5 do i6=1,Ny6 write(unit,*)Y1(i1,i2,i3,i4,i5,i6) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo case ("C") do i6=1,Ny6 do i5=1,Ny5 do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4,i5,i6) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA6_C !---------------------------- !---------------------------- !---------------------------- subroutine data_saveA7_R(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4,Ny5,Ny6,Ny7 integer :: i1,i2,i3,i4,i5,i6,i7 character(len=*) :: pname real(8),dimension(:,:,:,:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) Ny5=size(Y1,5) Ny6=size(Y1,6) Ny7=size(Y1,7) ! select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 do i5=1,Ny5 do i6=1,Ny6 do i7=1,Ny7 write(unit,*)Y1(i1,i2,i3,i4,i5,i6,i7) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo enddo case ("C") do i7=1,Ny7 do i6=1,Ny6 do i5=1,Ny5 do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4,i5,i6,i7) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA7_R subroutine data_saveA7_C(pname,Y1,order,wspace) integer :: Ny1,Ny2,Ny3,Ny4,Ny5,Ny6,Ny7 integer :: i1,i2,i3,i4,i5,i6,i7 character(len=*) :: pname complex(8),dimension(:,:,:,:,:,:,:) :: Y1 character(len=*),optional :: order logical,optional :: wspace character(len=1) :: order_ logical :: wspace_ order_ = "R" ; if(present(order))order_=trim(order(1:1)) wspace_= .true.; if(present(wspace))wspace_=wspace open(free_unit(unit),file=reg(pname)) Ny1=size(Y1,1) Ny2=size(Y1,2) Ny3=size(Y1,3) Ny4=size(Y1,4) Ny5=size(Y1,5) Ny6=size(Y1,6) Ny7=size(Y1,7) ! select case(order_) case default stop "save_array: order != Row-major, Col-major" case ("R") do i1=1,Ny1 do i2=1,Ny2 do i3=1,Ny3 do i4=1,Ny4 do i5=1,Ny5 do i6=1,Ny6 do i7=1,Ny7 write(unit,*)Y1(i1,i2,i3,i4,i5,i6,i7) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo enddo case ("C") do i7=1,Ny7 do i6=1,Ny6 do i5=1,Ny5 do i4=1,Ny4 do i3=1,Ny3 do i2=1,Ny2 do i1=1,Ny1 write(unit,*)Y1(i1,i2,i3,i4,i5,i6,i7) enddo if(wspace_)write(unit,*) enddo enddo enddo enddo enddo enddo end select close(unit) call file_gzip(reg(pname)) end subroutine data_saveA7_C