From 5b5f958c0f0fb92ab6bb273a2ba30d1427b1147d Mon Sep 17 00:00:00 2001 From: ignis Date: Sun, 24 Mar 2019 13:56:57 +0900 Subject: [PATCH] T^0.76 diffusivity --- code/ysolve.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/ysolve.f90 b/code/ysolve.f90 index 18c277b..8d835ff 100644 --- a/code/ysolve.f90 +++ b/code/ysolve.f90 @@ -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)