added input for variable diffusivity

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
 

code/test.sh

@@ -1,4 +1,4 @@
 cd sample/4pi-IC1
 ls
 ../../ex
-diff sfield.dat sfield.diff
+diff -y sfield.dat sfield.diff

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)