Changeset 20093

Timestamp:

02/08/16 15:15:15 (9 years ago)

Author:

jbondzio

Message:

CHG: compute in steadystate_core 1.velocity 2.temperature, since temperature is highly sensitive to velocity. This bugfix also makes steadystate basal BCs for enthalpy work.

Location:

issm/trunk-jpl

Files:

• src/c/cores/steadystate_core.cpp (modified) (1 diff)
• test/Archives/Archive318.nc (modified) ( previous)
• test/Archives/Archive410.nc (modified) ( previous)
• test/Archives/Archive411.nc (modified) ( previous)
• test/Archives/Archive415.nc (modified) ( previous)
• test/Archives/Archive416.nc (modified) ( previous)
• test/Archives/Archive431.nc (modified) ( previous)
• test/Archives/Archive432.nc (modified) ( previous)
• test/Archives/Archive436.nc (modified) ( previous)
• test/Archives/Archive508.nc (modified) ( previous)
• test/Archives/Archive509.nc (modified) ( previous)
• test/Archives/Archive510.nc (modified) ( previous)
• test/Archives/Archive513.nc (modified) ( previous)

issm/trunk-jpl/src/c/cores/steadystate_core.cpp

@@ -48 +48 @@
         for(;;){
 
-                if(VerboseSolution()) _printf0_("   computing temperature and velocity for step: " << step << "\n");
-                thermal_core(femmodel);
-                if(!isenthalpy)femmodel->SetCurrentConfiguration(ThermalAnalysisEnum);/*Could be MeltingAnalysis...*/
-                GetSolutionFromInputsx(&tg,femmodel);
+                /* Compute first velocity, then temperature due to high sensitivity of temperature to velocity. */
+                if(VerboseSolution()) _printf0_("   computing velocity and temperature for step: " << step << "\n");
 
                 if(VerboseSolution()) _printf0_("   computing new velocity\n");
                 stressbalance_core(femmodel);
                 GetSolutionFromInputsx(&ug,femmodel);
+
+                thermal_core(femmodel);
+                if(!isenthalpy)femmodel->SetCurrentConfiguration(ThermalAnalysisEnum);/*Could be MeltingAnalysis...*/
+                GetSolutionFromInputsx(&tg,femmodel);
 
                 if(step>1){