pdf_read_2d Subroutine

subroutine pdf_read_2d(self, pfile)

Arguments

Type IntentOptional Attributes Name
type(pdf_kernel_2d) :: self
character(len=*) :: pfile

Calls

proc~~pdf_read_2d~~CallsGraph proc~pdf_read_2d pdf_read_2d free_unit free_unit proc~pdf_read_2d->free_unit pdf_allocate pdf_allocate proc~pdf_read_2d->pdf_allocate

Source Code

subroutine pdf_read_2d(self,pfile)
  type(pdf_kernel_2d) :: self    
  character(len=*) :: pfile
  integer          :: i,N1,N2,unit
  if(.not.self%status)then
     print*,"PDF_READ: PDF not allocated"
  endif
  open(free_unit(unit),file=trim(pfile))
  read(unit,*)N1,N2
  call pdf_allocate(self,[N1,N2])
  read(unit,*)self%xmin
  read(unit,*)self%xmax
  read(unit,*)self%dx
  read(unit,*)self%Ndata   
  read(unit,*)self%x
  read(unit,*)self%y
  read(unit,*)self%pdf
  read(unit,*)self%sigma
  read(unit,*)self%status
  read(unit,*)self%variance
  read(unit,*)self%rescale
  close(unit)
end subroutine pdf_read_2d