[cig-commits] commit: Make sinker.input work for viscosity differences of up to 10^4 on 4 levels. Does not work for 5 levels. :(

Mercurial hg at geodynamics.org
Wed Apr 20 16:56:05 PDT 2011


changeset:   181:6f94da5f348e
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Apr 20 10:34:35 2011 -0700
files:       input/sinker.input
description:
Make sinker.input work for viscosity differences of up to 10^4 on 4 levels.  Does not work for 5 levels. :(


diff -r 949397d18c76 -r 6f94da5f348e input/sinker.input
--- a/input/sinker.input	Wed Apr 20 10:33:08 2011 -0700
+++ b/input/sinker.input	Wed Apr 20 10:34:35 2011 -0700
@@ -34,16 +34,16 @@ FACStokes {
         enable_logging = TRUE   // Bool flag to switch logging on/off
         max_cycles = 100         // Max number of FAC cycles to use
         residual_tol = 1e-5     // Residual tolerance to solve for
-        num_pre_sweeps = 0      // Number of presmoothing sweeps to use
-        num_post_sweeps = 0     // Number of postsmoothing sweeps to use
+        num_pre_sweeps = 100      // Number of presmoothing sweeps to use
+        num_post_sweeps = 100     // Number of postsmoothing sweeps to use
         smoothing_choice = "Tackley"
         coarse_solver_choice = "Tackley"
         // smoothing_choice = "Gerya"
         // coarse_solver_choice = "Gerya"
-        coarse_solver_max_iterations = 250
+        coarse_solver_max_iterations = 2500
         coarse_solver_tolerance = 1e-8
-        // p_prolongation_method = "P_MDPI_REFINE"
-        p_prolongation_method = "P_REFINE"
+        p_prolongation_method = "P_MDPI_REFINE"
+        // p_prolongation_method = "P_REFINE"
         v_prolongation_method = "V_REFINE"
     }
     bc_coefs {



More information about the CIG-COMMITS mailing list