Save Constants in SSA material models, else Budd law can't find gravity

Author: RupertGladstone

elmerice/Utils/SSAMaterialModels.F90

@@ -83,7 +83,7 @@ FUNCTION SSAEffectiveFriction(Element,nn,Basis,ub,SEP,PartlyGrounded,h,rho,rhow,
    TYPE(Solver_t), POINTER :: pSolver => NULL()
 
    SAVE FirstTime
-   SAVE :: GLnIP, GMSol, BedRockSol, pSolver, MaskName
+   SAVE :: GLnIP, GMSol, BedRockSol, pSolver, MaskName, Constants
 
    Material => GetMaterial(Element)