data_saveA2_C Subroutine

subroutine data_saveA2_C(pname, Y1, order, wspace)

Arguments

Type IntentOptional Attributes Name
character(len=*) :: pname
complex(kind=8), dimension(:,:) :: Y1
character(len=*), optional :: order
logical, optional :: wspace

Calls

proc~~data_savea2_c~2~~CallsGraph proc~data_savea2_c~2 data_saveA2_C file_gzip file_gzip proc~data_savea2_c~2->file_gzip free_unit free_unit proc~data_savea2_c~2->free_unit reg reg proc~data_savea2_c~2->reg

Source Code

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