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

T^0.76 diffusivity

Browse source
This commit is contained in:
ignis 2019-03-24 13:56:57 +09:00
parent 805ca37b61
commit 5b5f958c0f
1 changed files with 3 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
code/ysolve.f90
View file

@ -134,7 +134,7 @@
IF (ux(1,i).gt.1.) ux(1,i)=1.
ENDDO
dm = (/ (diff * diffusivity_sigmoid(1.-y1(i,1,1)), i=1,nx) /)
dm = (/ (diff * diffusivity_T076(1.-y1(i,1,1)), i=1,nx) /)
nd=1
CALL dfnonp(nx,hx,ux(1,:),dux(1,:),nd,1)
@ -341,7 +341,7 @@
WRITE(otape,*) cdum,sigw
DO i = 0,100
WRITE(*,*) diffusivity_sigmoid(i / 100.)
WRITE(*,*) diffusivity_T076(i / 100.)
END DO
IF (.not.read_stdin) THEN
@ -432,7 +432,7 @@
ux(3,i)=r1(i,j,k) ! Y
ENDDO
dm = (/ (diff * diffusivity_sigmoid(1.-r1(i,1,1)), i=1,nx) /)
dm = (/ (diff * diffusivity_T076(1.-r1(i,1,1)), i=1,nx) /)
CALL dfnonp(nx,hx,ux(1:3,:),dux(1:3,:),3,1)
CALL dfnonp(nx,hx,dm,duy(1,:),1,1)