From d6261282a12370a610c8fbf78794f8be7e0ddcf9 Mon Sep 17 00:00:00 2001 From: Yeongdo Park Date: Fri, 15 Mar 2019 15:17:32 -0400 Subject: [PATCH] added input for variable diffusivity --- code/sample/4pi-IC1/itape | 2 ++ code/test.sh | 2 +- code/ysolve.f90 | 5 +++++ 3 files changed, 8 insertions(+), 1 deletion(-) diff --git a/code/sample/4pi-IC1/itape b/code/sample/4pi-IC1/itape index 3e343cc..993a20d 100644 --- a/code/sample/4pi-IC1/itape +++ b/code/sample/4pi-IC1/itape @@ -16,4 +16,6 @@ c_ref 0.01 !0.003 ! min_wr 0. ! 5.0e-14 prof_wr 1.0 + rvis 1.0 + cstar 1.0 diff --git a/code/test.sh b/code/test.sh index 0744735..f05c4f4 100644 --- a/code/test.sh +++ b/code/test.sh @@ -1,4 +1,4 @@ cd sample/4pi-IC1 ls ../../ex -diff sfield.dat sfield.diff +diff -y sfield.dat sfield.diff diff --git a/code/ysolve.f90 b/code/ysolve.f90 index 284bddb..6a9886a 100644 --- a/code/ysolve.f90 +++ b/code/ysolve.f90 @@ -4,6 +4,7 @@ PRIVATE REAL, PARAMETER :: pi=3.14159265358979323846 REAL :: hx,dt,vis,sc,diff,pre,ac,bc,tf,t_now,t_uf,dt_uf + REAL :: rvis,cstar REAL :: c_cut,c_ref,refwr,minf,tar_lo,u0,ctmp,l_0,lo_flm=0. REAL :: er_lof=0., erdot=0.,min_wr,prof_wr REAL :: pflame,pflold,oldu @@ -291,6 +292,10 @@ WRITE(otape,*) cdum,min_wr READ(itape,*) cdum,prof_wr WRITE(otape,*) cdum,prof_wr + READ(itape,*) cdum,rvis + WRITE(otape,*) cdum,rvis + READ(itape,*) cdum,cstar + WRITE(otape,*) cdum,cstar CLOSE(itape)