Changeset 23367
- Timestamp:
- 10/02/18 07:06:44 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
r23366 r23367 607 607 //thickness[i] = old_thickness[i]*(1.0+opening-closing)/(1.0-opening+closing); 608 608 609 if(element->nodes[i]->Sid()==2299){610 printf("for node %i \n",element->nodes[i]->Sid());611 printf(" old thickness is %g \n",old_thickness[i]);612 printf(" new thickness is %g \n",thickness[i]);613 printf(" closing is %g \n",((2.0*A*dt*pow(EPL_N,n))/(pow(n,n))));614 printf(" opening is %g \n",((rho_water*gravity*epl_conductivity*EPLgrad2*dt)/(rho_ice*latentheat)));615 }609 /* if(element->nodes[i]->Sid()==2299){ */ 610 /* printf("for node %i \n",element->nodes[i]->Sid()); */ 611 /* printf(" old thickness is %g \n",old_thickness[i]); */ 612 /* printf(" new thickness is %g \n",thickness[i]); */ 613 /* printf(" closing is %g \n",((2.0*A*dt*pow(EPL_N,n))/(pow(n,n)))); */ 614 /* printf(" opening is %g \n",((rho_water*gravity*epl_conductivity*EPLgrad2*dt)/(rho_ice*latentheat))); */ 615 /* } */ 616 616 /*Take care of otherthikening*/ 617 617 if(thickness[i]>max_thick){
Note:
See TracChangeset
for help on using the changeset viewer.