writeinletfile Subroutine

public subroutine writeinletfile()

Uses

  • proc~~writeinletfile~~UsesGraph proc~writeinletfile modinlet::writeinletfile module~modglobal modglobal proc~writeinletfile->module~modglobal module~modmpi modmpi proc~writeinletfile->module~modmpi mpi mpi module~modmpi->mpi

Arguments

None

Called by

proc~~writeinletfile~~CalledByGraph proc~writeinletfile modinlet::writeinletfile proc~inletgen modinlet::inletgen proc~inletgen->proc~writeinletfile proc~inletgennotemp modinlet::inletgennotemp proc~inletgennotemp->proc~writeinletfile proc~boundary modboundary::boundary proc~boundary->proc~inletgen proc~boundary->proc~inletgennotemp proc~readinitfiles modstartup::readinitfiles proc~readinitfiles->proc~boundary program~dalesurban DALESURBAN program~dalesurban->proc~boundary proc~startup modstartup::startup program~dalesurban->proc~startup proc~startup->proc~readinitfiles

Contents

Source Code

writeinletfile

Source Code

  subroutine writeinletfile
    use modglobal, only : jb,je,kb,ke,cexpnr,ifoutput,nstore,ltempeq
    use modmpi,    only : cmyid,myid
!    use modinletdata, only : storeu0inletbc,storev0inletbc,storew0inletbc,nfile

    implicit none
    integer fileid
    integer j,k,n
    character(24) name

      name = 'inlet/inlet_    k   .'
      write (name(13:16)  ,'(i4.4)') nfile
      name(18:20)= cmyid
      name(22:24)= cexpnr

      write(6,*) 'Writing Inlet velocity: ', name
      open  (ifoutput,file=name,form='unformatted',position='append')

      write(ifoutput)  (((storeu0inletbc (j,k,n),j=jb,je),k=kb,ke),  n=1,nstore)
      write(ifoutput)  (((storev0inletbc (j,k,n),j=jb,je),k=kb,ke),  n=1,nstore)
      write(ifoutput)  (((storew0inletbc (j,k,n),j=jb,je),k=kb,ke+1),n=1,nstore)
      close (ifoutput)

      if (ltempeq ) then
        name = 'inlet/itemp_    k   .'
        write (name(13:16)  ,'(i4.4)') nfile
        name(18:20)= cmyid
        name(22:24)= cexpnr

        write(6,*) 'Writing Inlet temperature: ', name
        open  (ifoutput,file=name,form='unformatted',position='append')

        write(ifoutput)  (((storet0inletbc (j,k,n),j=jb,je),k=kb,ke),  n=1,nstore)
        close (ifoutput)
      end if

  end subroutine writeinletfile