streamlined non-Newtonian variables as defaults

Author: mrodrig6

src/forward_solver/m_imr_fd.m

@@ -11,6 +11,14 @@
         acos_opts, wave_opts, sigma_opts, thermal_opts, mass_opts] ...
         = f_call_params(varargin{:});
 
+    % defaults
+    
+    % non-Newtonian viscosity parameters
+    fnu             = 0;
+    intfnu          = 0;
+    dintfnu         = 0;
+    ddintfnu        = 0;
+    
     % equations settings
     radial          = eqns_opts(1);
     bubtherm        = eqns_opts(2);
@@ -97,12 +105,6 @@
     iDRe            = sigma_opts(13);
     iWe             = 1/We;
 
-    % viscosity parameters
-    fnu             = 0;
-    intfnu          = 0;
-    dintfnu         = 0;
-    ddintfnu        = 0;
-    
     % dimensionless thermal
     Foh             = thermal_opts(1);
     Br              = thermal_opts(2);

src/forward_solver/m_imr_spectral.m

@@ -11,6 +11,14 @@
         acos_opts, wave_opts, sigma_opts, thermal_opts, mass_opts] ...
         = f_call_params(varargin{:});
 
+    % defaults
+    
+    % viscosity parameters
+    fnu             = 0;
+    intfnu          = 0;
+    dintfnu         = 0;
+    ddintfnu        = 0;
+    
     % equations settings
     radial          = eqns_opts(1);
     bubtherm        = eqns_opts(2);
@@ -97,12 +105,6 @@
     iDRe            = sigma_opts(13);
     iWe             = 1/We;
 
-    % viscosity parameters
-    fnu             = 0;
-    intfnu          = 0;
-    dintfnu         = 0;
-    ddintfnu        = 0;
-    
     % dimensionless thermal
     Foh             = thermal_opts(1);
     Br              = thermal_opts(2);