Changeset 2330 for issm/trunk/src/c/parallel/ControlInitialization.cpp
- Timestamp:
- 09/25/09 17:00:28 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/ControlInitialization.cpp
r2316 r2330 34 34 35 35 /*flags: */ 36 int debug=0;36 int verbose=0; 37 37 int dim=-1; 38 38 int ishutter=0; … … 54 54 55 55 /*first recover parameters common to all solutions:*/ 56 model->FindParam(& debug,"debug");56 model->FindParam(&verbose,"verbose"); 57 57 model->FindParam(&dim,"dim"); 58 58 model->FindParam(&ishutter,"ishutter"); … … 83 83 84 84 //compute slopes 85 if( debug)_printf_("%s\n","computing bed slope (x and y derivatives)...");85 if(verbose)_printf_("%s\n","computing bed slope (x and y derivatives)..."); 86 86 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopeComputeAnalysisEnum(),BedXAnalysisEnum()); 87 87 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopeComputeAnalysisEnum(),BedYAnalysisEnum()); … … 97 97 98 98 //horizontal velocity 99 if( debug)_printf_("%s\n"," computing horizontal velocities...");99 if(verbose)_printf_("%s\n"," computing horizontal velocities..."); 100 100 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_dh,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum()); 101 if( debug)_printf_("%s\n"," extruding horizontal velocities...");101 if(verbose)_printf_("%s\n"," extruding horizontal velocities..."); 102 102 VecDuplicatePatch(&ug_horiz,ug); FieldExtrudex( ug_horiz,fem_dh->elements,fem_dh->nodes, fem_dh->loads,fem_dh-> materials,"velocity",1); 103 103 104 104 //vertical velocity 105 if( debug)_printf_("%s\n"," computing vertical velocities...");105 if(verbose)_printf_("%s\n"," computing vertical velocities..."); 106 106 inputs->Add("velocity",ug_horiz,numberofdofspernode_dh,numberofnodes); 107 107 diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum(),VertAnalysisEnum()); 108 108 109 109 //Create 3d u_g 110 if( debug)_printf_("%s\n"," combining horizontal and vertical velocities...");110 if(verbose)_printf_("%s\n"," combining horizontal and vertical velocities..."); 111 111 VecFree(&ug); ug=NewVec(numberofnodes*3); 112 112 xfree((void**)&dofset);dofset=dofsetgen(2,&dof01[0],3,numberofnodes*3); VecMerge(ug,ug_horiz,dofset,numberofnodes*2); … … 115 115 116 116 //Create 4d u_g 117 if( debug)_printf_("%s\n"," computing pressure according to Pattyn...");117 if(verbose)_printf_("%s\n"," computing pressure according to Pattyn..."); 118 118 ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->loads, fem_dh->materials, numberofnodes); 119 119 VecScale(pg,1.0/stokesreconditioning); … … 127 127 128 128 //update spcs 129 if( debug)_printf_("%s\n"," update boundary conditions for stokes using velocities previously computed...");129 if(verbose)_printf_("%s\n"," update boundary conditions for stokes using velocities previously computed..."); 130 130 xfree((void**)&dofset);dofset=dofsetgen(3,dof012,4,numberofnodes*4); VecMerge(fem_ds->yg->vector,ug,dofset,3*numberofnodes); 131 131 … … 134 134 135 135 //Compute Stokes velocities to speed up later runs 136 if( debug)_printf_("%s\n"," computing stokes velocities and pressure ...");136 if(verbose)_printf_("%s\n"," computing stokes velocities and pressure ..."); 137 137 VecFree(&ug); 138 138 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
Note:
See TracChangeset
for help on using the changeset viewer.