source: issm/oecreview/Archive/16554-17801/ISSM-17031-17032.diff

../trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -4731 +4731 @@
                                         //                                      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 +4896 @@
                                         /*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]); */
                                         }
                         }
                 }