Index: ../trunk-jpl/src/c/classes/Elements/Tria.cpp =================================================================== --- ../trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 17031) +++ ../trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 17032) @@ -4731,7 +4731,7 @@ // wh_trans=old_transfer[i]+0.66*(wh_trans-old_transfer[i]); /* if(this->Id()==27){ */ - /* printf("old %e, transfer is %e ,eplhead %e, sedhead %e, factor %e \n",old_transfer[i],wh_trans,epl_head[i],sed_head[i],epl_specificstoring[i]*epl_thickness[i]*sed_trans[i]); */ + /* printf("old %e, transfer is %e ,eplhead %e, sedhead %e, factor %e \n",old_transfer[i],wh_trans,epl_head[i],sed_head[i],epl_specificstoring[i]*epl_thickness[i]*sed_trans[i]); */ /* } */ /*Assign output pointer*/ @@ -4896,8 +4896,8 @@ /*And proceed to the real thing*/ thickness[i] = old_thickness[i]*(1+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n))); // thickness[i] = old_thickness[i]+0.66*(thickness[i]-old_thickness[i]); - if(this->Id()==27){ - printf("old %e, thickness is %e ,eplhead %e, slopeX %e, slopeY %e\n",old_thickness[i],thickness[i],eplhead[i],epl_slopeX[i],epl_slopeY[i]); + /* if(this->Id()==27){ */ + /* printf("old %e, thickness is %e ,eplhead %e, slopeX %e, slopeY %e\n",old_thickness[i],thickness[i],eplhead[i],epl_slopeX[i],epl_slopeY[i]); */ } } }