CHG: fixed problem with Glads

ISSMteam · Jun 25, 2024 · 2063a05 · 2063a05
1 parent 1644d52
commit 2063a05
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/c/analyses/HydrologyGlaDSAnalysis.cpp b/src/c/analyses/HydrologyGlaDSAnalysis.cpp
@@ -455,8 +455,8 @@ ElementVector* HydrologyGlaDSAnalysis::CreatePVector(Element* element){/*{{{*/
 
 		/*Transient term if dt>0*/
 		if(dt>0.){
-			factor = gauss->weight*Jdet*e_v/(rho_water*g*dt)*phi_old;
 			phiold_input->GetInputValue(&phi_old,gauss);
+			factor = gauss->weight*Jdet*e_v/(rho_water*g*dt)*phi_old;
 			for(int i=0;i<numnodes;i++) pe->values[i] += factor*basis[i];
 		}
 	}