main Program

solve double mach reflection problem x方向格子点数 y方向格子点数 時間刻み幅 比熱比と関連定数 定数(計算速度性能に寄与するかはよくわからない) パラメータ等設定 includeを用いるのは本来好ましくない 最大反復数 終了時間


Calls

program~~main~~CallsGraph program~main main proc~init init program~main->proc~init df df program~main->df proc~outf outf program~main->proc~outf proc~stbc stbc program~main->proc~stbc proc~intp intp program~main->proc~intp proc~grid grid program~main->proc~grid proc~c2p c2p program~main->proc~c2p proc~iflx iflx program~main->proc~iflx proc~intg intg program~main->proc~intg cons cons program~main->cons proc~init->cons y y proc~init->y pris pris proc~init->pris dsign dsign proc~init->dsign x x proc~init->x proc~outf->cons proc~stbc->cons proc~stbc->dsign proc~stbc->x proc~intp->pris pyl pyl proc~intp->pyl proc~intp->dsign pxr pxr proc~intp->pxr pxl pxl proc~intp->pxl pyr pyr proc~intp->pyr proc~grid->y proc~grid->x proc~c2p->cons proc~c2p->pris proc~iflx->pyl proc~iflx->pxr proc~iflx->pxl proc~iflx->pyr dnflx dnflx proc~iflx->dnflx proc~slau slau proc~iflx->proc~slau conl conl proc~iflx->conl conr conr proc~iflx->conr xflx xflx proc~iflx->xflx yflx yflx proc~iflx->yflx proc~intg->df proc~intg->xflx proc~intg->yflx proc~slau->dsign proc~slau->dnflx proc~slau->conl proc~slau->conr

Contents


Common Blocks

格子点座標 0とmax+1は便宜上 計算領域は0-max

TypeAttributesNameInitial
real :: x(0:jmax+1,0:kmax+1)
real :: y(0:jmax+1,0:kmax+1)

TypeAttributesNameInitial
real :: dltx
real :: dlty

保存量 

TypeAttributesNameInitial
real :: cons(4,0:jmax+1,0:kmax+1)

基本量 

TypeAttributesNameInitial
real :: pris(4,0:jmax+1,0:kmax+1)

t:現在時間, tt:サブステップの時間

TypeAttributesNameInitial
real :: tn
real :: tt

反復数

TypeAttributesNameInitial
integer :: im

before shock state

TypeAttributesNameInitial
real :: r0
real :: p0
real :: u0
real :: v0
real :: ru0
real :: rv0
real :: e0

after shock state

TypeAttributesNameInitial
real :: r1
real :: p1
real :: u1
real :: v1
real :: ru1
real :: rv1
real :: e1

非粘性流束

TypeAttributesNameInitial
real :: xflx(4,jmax,kmax)
real :: yflx(4,jmax,kmax)

流束評価面に補間された値 l:left, r:right

TypeAttributesNameInitial
real :: pxl(4,0:jmax+1,0:kmax+1)
real :: pxr(4,0:jmax+1,0:kmax+1)
real :: pyl(4,0:jmax+1,0:kmax+1)
real :: pyr(4,0:jmax+1,0:kmax+1)

時間増分

TypeAttributesNameInitial
real :: df(4,jmax,kmax)

非粘性流束評価で用いる値

TypeAttributesNameInitial
real :: conL(4)
real :: conR(4)
real :: dnflx(4)

流束評価面の法線ベクトル

TypeAttributesNameInitial
real :: dnx
real :: dny

Variables

TypeAttributesNameInitial
double precision :: conn(4,0:jmax+1,0:kmax+1)

n stepの保存量の一時的置き場