Changeset 927
- Timestamp:
- 06/11/09 16:47:05 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Pengrid.cpp
r926 r927 349 349 if (temperature<t_pmp){ //If T<Tpmp, there must be no melting. Therefore, melting should be constrained to 0 when T<Tpmp, instead of using spcs, use penalties 350 350 Ke[0][0]=kmax*pow(10,penalty_offset); 351 if (id==65){352 printf("Ke_terms[0] = %10.10lf\n",Ke[0]);353 printf("penalty_offset(K) = %g\n",penalty_offset);354 printf("kmax = %10.10lf\n",kmax);355 }356 351 } 357 352 … … 507 502 if (sub_analysis_type==SteadyAnalysisEnum()){ 508 503 P_terms[0]=melting_offset*pow(10,penalty_offset)*(temperature-t_pmp); 509 if (id==65){510 printf("penalty_offset(P) = %g\n",penalty_offset);511 printf("melting_offset(P) = %10.10lf\n",melting_offset);512 printf("P_terms[0] = %10.10lf\n",P_terms[0]);513 }514 504 } 515 505 else{
Note:
See TracChangeset
for help on using the changeset viewer.