From 0ca93a259796312d113dcda9c45f6c501ce7b318 Mon Sep 17 00:00:00 2001 From: ignis Date: Tue, 27 Oct 2020 08:06:01 +0900 Subject: [PATCH] temporarily fix viscosity ratio to 5 --- code/ysolve.f90 | 1 + 1 file changed, 1 insertion(+) diff --git a/code/ysolve.f90 b/code/ysolve.f90 index 3bbab83..4ed1c61 100644 --- a/code/ysolve.f90 +++ b/code/ysolve.f90 @@ -670,6 +670,7 @@ REAL FUNCTION diffusivity_sutherland(c) REAL, INTENT(IN) :: c REAL :: theta, As, Ts, T0, T1 + REAL, PARAMETER :: rvis=5.0d0 T0 = 1.0 T1 = (1.0 + bc)