research/incomp-flame-post
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

removed version info from m_terms.f90

Browse source
This commit is contained in:
ignis 2019-05-23 11:47:42 +09:00
parent eec864f2c7
commit 05a6fb92d5
3 changed files with 0 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
code/code_gen/code_gen.py
View file

@ -92,8 +92,6 @@ def test():
with open("terms.input") as inputfile:
terms_raw = ((inputfile.read()))
vinfo = VersionInfo(terms_raw)
tree = Lark(calc_grammar, parser='lalr' ).parse(terms_raw)
@ -107,8 +105,6 @@ def test():
cdict = ir4.print_program()
cdict["module_data"] = vinfo.fparams() + "\n" + cdict["module_data"]
print mod_form.format( cdict )
# print mod_form.format( ir4.print_program() )

3
code/m_terms.f90
View file

@ -7,9 +7,6 @@ use m_calculate
implicit none
character (len=*), parameter :: build_date="Fri May 10 04:36:01 2019"
character (len=*), parameter :: build_base="b4db0e2f8a3669de8450525bdabbfeeb7c215005"
character (len=*), parameter :: term_info="[u, v, w, y]"//char(10)//""//char(10)//"c = 1.0 - y"//char(10)//"wrate = rxn_rate(c)"//char(10)//"fsd = sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))"//char(10)//""//char(10)//"sd = ((d2dx(c) + d2dy(c) + d2dz(c)) * $rod + wrate) / fsd"//char(10)//""//char(10)//"nx = - ddx(c) / fsd"//char(10)//"ny = - ddy(c) / fsd"//char(10)//"nz = - ddz(c) / fsd"//char(10)//""//char(10)//"vn = u * nx + v * ny + w * nz"//char(10)//""//char(10)//"divn = ddx(nx) + ddy(ny) + ddz(nz)"//char(10)//""//char(10)//"absk = abs(divn)"//char(10)//""//char(10)//"lapc = d2dx(c) + d2dy(c) + d2dz(c)"//char(10)//""//char(10)//"t1 = (vn + sd) * nx"//char(10)//""//char(10)//"t2 = t1' * nx'"//char(10)//""//char(10)//"t3 = t1' * divn'"//char(10)//""//char(10)//"t4 = t1' * absk'"//char(10)//""//char(10)//"t5 = sqr(t1')"//char(10)//""//char(10)//"t6 = sqr(nx')"//char(10)//""//char(10)//"t7 = sqr(divn')"//char(10)//""//char(10)//"t8 = sqr(absk')"//char(10)//""//char(10)//"avg { c, u, v, w, fsd, lapc, ddx_c, ddy_c, ddz_c, d2dx_c, d2dy_c, d2dz_c }"//char(10)//""//char(10)//"avg fsd { u, v, w, nx, ny, nz, t1, t2, t3, t4, t5, t6, t7, t8, absk, divn, sd }"//char(10)//""
character (len = *), parameter :: output_header="x avg_c avg_d2dx_c avg_d2dy_c avg_d2dz_c avg_ddx_c avg_ddy_c avg_ddz_c avg_fsd avg_lapc avg_u avg_v avg_w fsd_avg_absk fsd_avg_divn fsd_avg_nx fsd_avg_ny fsd_avg_nz fsd_avg_sd fsd_avg_t1 fsd_avg_u fsd_avg_v fsd_avg_w fsd_avg_t2 fsd_avg_t3 fsd_avg_t4 fsd_avg_t5 fsd_avg_t6 fsd_avg_t7 fsd_avg_t8"
real*8, allocatable, dimension(:) :: fsd_avg_divn
real*8, allocatable, dimension(:) :: fsd_avg_ny

5
code/post.f90
View file

@ -100,11 +100,6 @@
!========================================================================================
SUBROUTINE PRINT_BANNER
WRITE(*,*) ' build date: ', build_date
WRITE(*,*) ' build base: ', build_base
WRITE(*,*) ' terms to process '
WRITE(*,*) term_info
WRITE(*,*) ' This program, x-edge-cold-bc-5-hybrid, is written by D. Kim, 2018'
WRITE(*,*) ' It is to study the statistics of the flame parameters at the leading edge'
WRITE(*,*) ' in turbulent premixed flames.'