readview.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	readview (ntr, adview, auview, fenv, nzsrad, ithermal, jobnamef)

Function/Subroutine Documentation

readview()

subroutine readview	(	integer	ntr,
		real*8, dimension(*)	adview,
		real*8, dimension(*)	auview,
		real*8, dimension(*)	fenv,
		integer	nzsrad,
		integer	ithermal,
		character*132, dimension(*)	jobnamef
	)

!     
!     reading the viewfactors from file
!
      implicit none
!     
      logical exi
!
      character*132 jobnamef(*),fnvw
!     
      integer ntr,nzsrad,ithermal,i,k
!
      real*8 adview(*),auview(*),fenv(*)
!     
      if(ithermal.eq.3) then
         write(*,*) '*WARNING in radmatrix: viewfactors are being'
         write(*,*) '         read from file for a thermomechani-'
         write(*,*) '         cal calculation: they will not be '
         write(*,*) '         recalculated in every iteration.'
      endif
!     
      write(*,*) 'Reading the viewfactors from file'
      write(*,*)
!     
      if(jobnamef(2)(1:1).eq.' ') then
         do i=1,132
            if(jobnamef(1)(i:i).eq.' ') exit
         enddo
         i=i-1
         fnvw=jobnamef(1)(1:i)//'.vwf'
      else
         fnvw=jobnamef(2)
      endif
      inquire(file=fnvw,exist=exi)
      if(exi) then
         open(10,file=fnvw,status='old',form='unformatted',
     &        access='sequential',err=10)
      else
         write(*,*) '*ERROR in radmatrix: viewfactor file ',fnvw
         write(*,*) 'does not exist'
         call exit(201)
      endif
!     
      read(10) (adview(k),k=1,ntr)
      read(10) (auview(k),k=1,2*nzsrad)
      read(10)(fenv(k),k=1,ntr)
!     
      close(10)
!     
      return
!
 10   write(*,*) '*ERROR in radmatrix: could not open file ',fnvw
      call exit(201)