Changeset 25680
- Timestamp:
- 10/12/20 09:50:51 (4 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 25 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp
r25539 r25680 69 69 70 70 if(save_results){ 71 if(VerboseSolution()) _printf0_(" saving results\n");72 71 femmodel->RequestedOutputsx(&femmodel->results,&extvar_enum,1); 73 72 } -
issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
r25539 r25680 168 168 169 169 if(save_results){ 170 if(VerboseSolution()) _printf0_(" saving levelset results\n");171 170 int outputs[1] = {MaskIceLevelsetEnum}; 172 171 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],1); -
issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
r25611 r25680 1050 1050 1051 1051 if(isFS && !(isSSA || isHO || isL1L2 || isMLHO)){ 1052 if(VerboseSolution()) _printf0_(" computing velocities\n");1053 1054 1052 femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum); 1055 1053 … … 1075 1073 } 1076 1074 else if(!isFS && (isSSA || isHO || isL1L2 || isMLHO)){ 1077 if(VerboseSolution()) _printf0_(" computing velocities\n");1078 1079 1075 femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum); 1080 1076 if(newton>0) -
issm/trunk-jpl/src/c/cores/adjointbalancethickness2_core.cpp
r22178 r25680 33 33 /*Save results*/ 34 34 if(save_results || true){ 35 if(VerboseSolution()) _printf0_(" saving results\n");36 35 int outputs[1] = {AdjointEnum}; 37 36 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],1); -
issm/trunk-jpl/src/c/cores/adjointbalancethickness_core.cpp
r16518 r25680 33 33 /*Save results*/ 34 34 if(save_results){ 35 if(VerboseSolution()) _printf0_(" saving results\n");36 35 int outputs[1] = {AdjointEnum}; 37 36 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],1); -
issm/trunk-jpl/src/c/cores/balancethickness2_core.cpp
r22830 r25680 33 33 34 34 if(save_results){ 35 if(VerboseSolution()) _printf0_(" saving results\n");36 35 const int numoutputs = 1; 37 36 int outputs[numoutputs] = {ThicknessEnum}; -
issm/trunk-jpl/src/c/cores/balancethickness_core.cpp
r25239 r25680 36 36 37 37 if(save_results){ 38 if(VerboseSolution()) _printf0_(" saving results\n");39 38 int outputs = ThicknessEnum; 40 39 femmodel->RequestedOutputsx(&femmodel->results,&outputs,1); -
issm/trunk-jpl/src/c/cores/balancevelocity_core.cpp
r23380 r25680 32 32 33 33 if(save_results){ 34 if(VerboseSolution()) _printf0_(" saving results\n");35 34 int outputs[3] = {DrivingStressXEnum,DrivingStressYEnum,VelEnum}; 36 35 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3); -
issm/trunk-jpl/src/c/cores/bedslope_core.cpp
r24240 r25680 34 34 35 35 if(save_results){ 36 if(VerboseSolution()) _printf0_(" saving bedslopes results\n");37 36 if(domaintype!=Domain2DverticalEnum){ 38 37 int outputs[2] = {BedSlopeXEnum,BedSlopeYEnum}; -
issm/trunk-jpl/src/c/cores/damage_core.cpp
r24240 r25680 40 40 41 41 if(save_results){ 42 if(VerboseSolution()) _printf0_(" saving damage results\n");43 42 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); 44 43 } -
issm/trunk-jpl/src/c/cores/depthaverage_core.cpp
r21799 r25680 16 16 int inputenum,input_average_enum; 17 17 18 if(VerboseSolution()) _printf0_(" depth averaging solution...\n");19 20 18 /*Get parameters*/ 21 19 femmodel->parameters->FindParam(&domaintype,DomainTypeEnum); … … 30 28 } 31 29 30 if(VerboseSolution()) _printf0_(" depth averaging "<<EnumToStringx(inputenum)<<"\n"); 31 32 32 /*Special method for Penta, otherwise call solution sequence*/ 33 33 if(elementtype==PentaEnum){ -
issm/trunk-jpl/src/c/cores/esa_core.cpp
r25465 r25680 99 99 100 100 if(save_results){ 101 if(VerboseSolution()) _printf0_(" saving results\n");102 101 femmodel->parameters->FindParam(&requested_outputs,&numoutputs,EsaRequestedOutputsEnum); 103 102 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); -
issm/trunk-jpl/src/c/cores/extrudefrombase_core.cpp
r21799 r25680 15 15 int elementtype,domaintype; 16 16 17 if(VerboseSolution()) _printf0_(" extruding solution from base...\n");18 19 17 /*Get parameters*/ 20 18 femmodel->parameters->FindParam(&domaintype,DomainTypeEnum); … … 24 22 if(domaintype==Domain2DhorizontalEnum) return; 25 23 24 int inputenum; 25 femmodel->parameters->FindParam(&inputenum,InputToExtrudeEnum); 26 if(VerboseSolution()) _printf0_(" extruding "<<EnumToStringx(inputenum)<<" from base...\n"); 27 26 28 /*Special method for Penta, otherwise call solution sequence*/ 27 29 if(elementtype==PentaEnum){ 28 int inputenum; femmodel->parameters->FindParam(&inputenum,InputToExtrudeEnum);29 30 InputExtrudex(femmodel,inputenum,-1); 30 31 } -
issm/trunk-jpl/src/c/cores/gia_core.cpp
r25379 r25680 58 58 59 59 if(save_results){ 60 if(VerboseSolution()) _printf0_(" saving results\n");61 60 int outputs[2] = {UGiaEnum,UGiaRateEnum}; 62 61 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],2); -
issm/trunk-jpl/src/c/cores/groundingline_core.cpp
r25465 r25680 22 22 23 23 /*Move grounding line*/ 24 if(VerboseSolution()) _printf0_(" computing new grounding line position\n"); 24 25 GroundinglineMigrationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters); 25 26 … … 33 34 34 35 if(save_results){ 35 if(VerboseSolution()) _printf0_(" saving groundingline results\n");36 36 int outputs[3] = {SurfaceEnum,BaseEnum,MaskOceanLevelsetEnum}; 37 37 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3); -
issm/trunk-jpl/src/c/cores/hydrology_core.cpp
r25489 r25680 187 187 } 188 188 if(save_results){ 189 if(VerboseSolution()) _printf0_(" saving hydrology results \n");190 189 if(hydrology_model==HydrologydcEnum && ThawedNodes==0){ 191 190 if(VerboseSolution()) _printf0_(" No thawed node hydro is skiped \n");} -
issm/trunk-jpl/src/c/cores/masstransport_core.cpp
r24490 r25680 85 85 86 86 if(save_results){ 87 if(VerboseSolution()) _printf0_(" saving results\n");88 87 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); 89 88 } -
issm/trunk-jpl/src/c/cores/movingfront_core.cpp
r25241 r25680 48 48 } 49 49 50 50 51 /* start the work from here */ 52 if(VerboseSolution()) _printf0_(" computing calving and undercutting\n"); 51 53 Calvingx(femmodel); 52 54 FrontalForcingsx(femmodel); 53 if(VerboseSolution()) _printf0_(" computing level set transport\n");55 if(VerboseSolution()) _printf0_(" computing new ice front position\n"); 54 56 55 57 /* smoothen slope of lsf for computation of normal on ice domain*/ -
issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp
r25627 r25680 65 65 int numoutputs; 66 66 char **requested_outputs = NULL; 67 68 if(VerboseSolution()) _printf0_(" saving results\n");69 67 femmodel->parameters->FindParam(&requested_outputs,&numoutputs,SealevelriseRequestedOutputsEnum); 70 68 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); -
issm/trunk-jpl/src/c/cores/smb_core.cpp
r24790 r25680 103 103 104 104 if(save_results){ 105 if(VerboseSolution()) _printf0_(" saving smb results\n");106 105 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); 107 106 } -
issm/trunk-jpl/src/c/cores/steadystate_core.cpp
r22886 r25680 78 78 79 79 if(save_results){ 80 if(VerboseSolution()) _printf0_(" saving results\n");81 80 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); 82 81 } -
issm/trunk-jpl/src/c/cores/stressbalance_core.cpp
r25514 r25680 94 94 95 95 if(save_results){ 96 if(VerboseSolution()) _printf0_(" saving stressbalance results\n");97 96 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); 98 if(VerboseSolution()) _printf0_(" results saved\n");99 97 } 100 98 -
issm/trunk-jpl/src/c/cores/surfaceslope_core.cpp
r24240 r25680 38 38 39 39 if(save_results){ 40 if(VerboseSolution()) _printf0_("saving surface slopes results:\n");41 40 if(domaintype!=Domain2DverticalEnum){ 42 41 int outputs[2] = {SurfaceSlopeXEnum,SurfaceSlopeYEnum}; -
issm/trunk-jpl/src/c/cores/thermal_core.cpp
r24335 r25680 50 50 51 51 if(save_results){ 52 if(VerboseSolution()) _printf0_(" saving thermal results\n");53 52 femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs); 54 53 } -
issm/trunk-jpl/src/c/cores/transient_core.cpp
r25637 r25680 190 190 /* Using Hydrology dc coupled we need to compute smb in the hydrology inner time loop*/ 191 191 if(issmb) { 192 if(VerboseSolution()) _printf0_(" computing smb\n");193 192 smb_core(femmodel); 194 193 } 195 194 196 195 if(ishydrology){ 197 if(VerboseSolution()) _printf0_(" computing hydrology\n");198 196 int hydrology_model; 199 197 hydrology_core(femmodel); … … 203 201 204 202 if(isstressbalance && (step%sb_coupling_frequency==0 || step==1) ) { 205 if(VerboseSolution()) _printf0_(" computing stress balance\n");206 203 stressbalance_core(femmodel); 207 204 } 208 205 209 206 if(isdamageevolution) { 210 if(VerboseSolution()) _printf0_(" computing damage\n");211 207 damage_core(femmodel); 212 208 } 213 209 214 210 if(ismovingfront) { 215 if(VerboseSolution()) _printf0_(" computing moving front\n");216 211 movingfront_core(femmodel); 217 212 } … … 220 215 221 216 if(ismasstransport){ 222 if(VerboseSolution()) _printf0_(" computing mass transport\n");223 217 bmb_core(femmodel); 224 218 masstransport_core(femmodel); … … 227 221 228 222 if(isgroundingline){ 229 if(VerboseSolution()) _printf0_(" computing new grounding line position\n");230 223 groundingline_core(femmodel); 231 224 } 232 225 233 226 if(isgia){ 234 if(VerboseSolution()) _printf0_(" computing glacial isostatic adjustment\n");235 227 #ifdef _HAVE_GIA_ 236 228 gia_core(femmodel); … … 461 453 delete hdl_setadjoint; 462 454 455 if(VerboseSolution()) _printf0_(" computing gradient\n"); 463 456 tape_codi.evaluate(); 464 457
Note:
See TracChangeset
for help on using the changeset viewer.