TURBO Q VSTAGE Initial File Format




Turbo Q VSTAGE Initial File Format (QU*0)


Pseudo Code:

ni1 = ni+1
nj1 = nj+1
nk1 = nk+1

read iter
read (((p(i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((q(1,i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((q(2,i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((q(3,i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((q(4,i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((q(5,i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((dmut(i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((qt1(i,j,k),i=1,ni1),j=1,nj1),k=1,nk1),
(((qt2(i,j,k),i=1,ni1),j=1,nj1),k=1,nk1)
read num_k_tsbc
if (num_k_tsbc != 0)
   read nlast
   read (trec(mt),mt=1,nlast)
   do n=1,num_k_tsbc
     read is,ie,js,je
     do mt=1,nlast
       read ((((qpk[n](i,j,k,mt,l),i=is,ie),j=js,je),k=1,2),l=1,7)
     enddo
   enddo
   if (IminBC eq slide_ts_i (402))
     read js,je,ks,ke
     do mt=1,nlast
       read ((((qpi(i,j,k,mt,l),i=1,2),j=js,je),k=ks,ke),l=1,7)
     enddo
   endif
   if (ImaxBC eq slide_ts_i (402))
     read js,je,ks,ke
     do mt=1,nlast
       read ((((qpi(i,j,k,mt,l),i=1,2),j=js,je),k=ks,ke),l=1,7)
     enddo
   endif
endif
if (IminBC eq wbsteady_in (203))
   read ((( pus(i,j,k),i=1,2),j=2,nj),k=2,nk),
((( rus(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((ruus(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((rvus(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((rwus(i,j,k),i=1,2),j=2,nj),k=2,nk),
((( eus(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((akus(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((epus(i,j,k),i=1,2),j=2,nj),k=2,nk)
else if (IminBC eq wb_unsteady_in (204))
   read mkus
   read ((( pus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
((( rus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
(((ruus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
(((rvus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
(((rwus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
((( eus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
((( zus(i,j,k),i=1,3),j=1,nj),k=1,mkus),
(((akus(i,j,k),i=1,2),j=2,nj),k=2,mkus),
(((epus(i,j,k),i=1,2),j=2,nj),k=2,mkus)
endif
if (ImaxBC eq wb_steady_exit (302))
   read ((( pbas(i,j,k),i=1,2),j=2,nj),k=2,nk),
((( rbas(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((rubas(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((rvbas(i,j,k),i=1,2),j=2,nj),k=2,nk),
(((rwbas(i,j,k),i=1,2),j=2,nj),k=2,nk),
((( ebas(i,j,k),i=1,2),j=2,nj),k=2,nk)
else if (ImaxBC eq wb_unsteady_exit (302))
   read mkds
   read ((( pbas(i,j,k),i=1,2),j=2,nj),k=2,mkds),
((( rbas(i,j,k),i=1,2),j=2,nj),k=2,mkds),
(((rubas(i,j,k),i=1,2),j=2,nj),k=2,mkds),
(((rvbas(i,j,k),i=1,2),j=2,nj),k=2,mkds),
(((rwbas(i,j,k),i=1,2),j=2,nj),k=2,mkds),
((( ebas(i,j,k),i=1,2),j=2,nj),k=2,mkds),
((( zbas(i,j,k),i=1,3),j=1,nj),k=1,mkds)
endif
read num_source
do n=1,num_source
   read isr(n),jsr(n),ksr(n),sfr(n),sfru(n),sfrv(n),sfrw(n),sfe(n),sfk(n),sfep(n)
enddo

Description:

ni,nj,nk [integer]: block dimensions
iter [integer]: Set to 0
p,q,dmut,qt1,qt2 [double]: cell centered flow variables
num_k_tsbc [integer]: number of time shift k patches (>0)
or
-(number of time shift i patches) (<0)
or
zero (= 0)
nlast [integer]: number of time steps stored
trec [float]: time record
is,ie,js,je,ks,ke [integer]: starting and ending indices of time shift patches
qpk,qpi [float]: time shift data on patch
mkus [integer]: number of tangential cells (us)
pus,rus,ruus,rvus,rwus,eus,zus,akus,epus [double]: inlet profile data
mkds [integer]: number of tangential cells (ds)
pbas,rbas,rubas,rvbas,rwbas,ebas,zbas [double]: exit profile data
num_source [integer]: number of source points
isr,jsr,ksr [integer]: cell center location
sfr,sfru,sfrv,sfrw,sfe,sfk,sfep [double]: source values

Notes:

Written per block. Dimensions and boundary conditions are known.