Index: ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (revision 17411) +++ ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (revision 17412) @@ -3021,10 +3021,10 @@ /*Retrieve all inputs and parameters*/ element->GetVerticesCoordinatesBase(&xyz_list_base); element->FindParam(&dt,TimesteppingTimeStepEnum); - if(dt==0) dt=pow(10,5); - IssmDouble rho_water =element->GetMaterialParameter(MaterialsRhoWaterEnum); - IssmDouble gravity =element->GetMaterialParameter(ConstantsGEnum); - Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input); + if(dt==0) dt=1.e+5; + IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoWaterEnum); + IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); + Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input); /* Start looping on the number of gaussian points: */ gauss=element->NewGaussBase(3);