[17802] | 1 | Index: ../trunk-jpl/src/c/classes/Elements/Tria.cpp
|
---|
| 2 | ===================================================================
|
---|
| 3 | --- ../trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 17031)
|
---|
| 4 | +++ ../trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 17032)
|
---|
| 5 | @@ -4731,7 +4731,7 @@
|
---|
| 6 | // wh_trans=old_transfer[i]+0.66*(wh_trans-old_transfer[i]);
|
---|
| 7 |
|
---|
| 8 | /* if(this->Id()==27){ */
|
---|
| 9 | - /* 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]); */
|
---|
| 10 | + /* 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]); */
|
---|
| 11 | /* } */
|
---|
| 12 |
|
---|
| 13 | /*Assign output pointer*/
|
---|
| 14 | @@ -4896,8 +4896,8 @@
|
---|
| 15 | /*And proceed to the real thing*/
|
---|
| 16 | 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)));
|
---|
| 17 | // thickness[i] = old_thickness[i]+0.66*(thickness[i]-old_thickness[i]);
|
---|
| 18 | - if(this->Id()==27){
|
---|
| 19 | - 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]);
|
---|
| 20 | + /* if(this->Id()==27){ */
|
---|
| 21 | + /* 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]); */
|
---|
| 22 | }
|
---|
| 23 | }
|
---|
| 24 | }
|
---|