Changeset 8392
- Timestamp:
- 05/23/11 11:52:10 (14 years ago)
- Location:
- issm/trunk
- Files:
-
- 114 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/examples/Bumps/Bump1_surface_bed/Square.par
r7810 r8392 21 21 md.eps_abs=10; %m/yr 22 22 md.penalty_offset=4; 23 md.penalty_ melting=10^7;23 md.penalty_basal_melting=10^7; 24 24 if md.numberofgrids<1000000, 25 25 md.sparsity=.001; … … 70 70 disp(' creating accumulation rates'); 71 71 md.accumulation_rate=ones(md.numberofgrids,1)/md.yts; %1m/a 72 md. melting_rate=0*ones(md.numberofgrids,1)/md.yts; %1m/a72 md.basal_melting_rate=0*ones(md.numberofgrids,1)/md.yts; %1m/a 73 73 74 74 %Deal with boundary conditions: -
issm/trunk/examples/Bumps/Bump1_surface_bed/runme.m
r8294 r8392 13 13 md.vy=PatchToVec(md.results.DiagnosticSolution.Vy); 14 14 md.accumulation_rate(:)=0; 15 md. melting_rate(:)=0;15 md.basal_melting_rate(:)=0; 16 16 md.thickness(:)=1; 17 17 md.bed=md.surface-md.thickness; -
issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h
r8376 r8392 226 226 MeltingOffsetEnum, 227 227 GroundingLineMeltingRateEnum, 228 MeltingRateEnum,228 BasalMeltingRateEnum, 229 229 MisfitEnum, 230 230 NodeOnBedEnum, … … 341 341 GravityEnum, 342 342 ThermalConductivityEnum, 343 MeltingRateCorrectionEnum,344 MeltingRateCorrectionApplyEnum,343 BasalMeltingRateCorrectionEnum, 344 BasalMeltingRateCorrectionApplyEnum, 345 345 MixedLayerCapacityEnum, 346 346 ThermalExchangeVelocityEnum, -
issm/trunk/src/c/EnumDefinitions/EnumToModelField.cpp
r8371 r8392 24 24 case GroundingLineMigrationEnum : return "gl_migration"; 25 25 case AccumulationRateEnum: return "accumulation_rate"; 26 case MeltingRateEnum: return "melting_rate";26 case BasalMeltingRateEnum: return "basal_melting_rate"; 27 27 default : _error_("No model field is associated to enum %s",EnumToStringx(en)); 28 28 } -
issm/trunk/src/c/modules/EnumToStringx/EnumToStringx.cpp
r8376 r8392 194 194 case MeltingOffsetEnum : return "MeltingOffset"; 195 195 case GroundingLineMeltingRateEnum : return "GroundingLineMeltingRate"; 196 case MeltingRateEnum : return "MeltingRate";196 case BasalMeltingRateEnum : return "BasalMeltingRate"; 197 197 case MisfitEnum : return "Misfit"; 198 198 case NodeOnBedEnum : return "NodeOnBed"; … … 295 295 case GravityEnum : return "Gravity"; 296 296 case ThermalConductivityEnum : return "ThermalConductivity"; 297 case MeltingRateCorrectionEnum : return "MeltingRateCorrection";298 case MeltingRateCorrectionApplyEnum : return "MeltingRateCorrectionApply";297 case BasalMeltingRateCorrectionEnum : return "BasalMeltingRateCorrection"; 298 case BasalMeltingRateCorrectionApplyEnum : return "BasalMeltingRateCorrectionApply"; 299 299 case MixedLayerCapacityEnum : return "MixedLayerCapacity"; 300 300 case ThermalExchangeVelocityEnum : return "ThermalExchangeVelocity"; -
issm/trunk/src/c/modules/ModelProcessorx/Balancethickness/UpdateElementsBalancethickness.cpp
r8330 r8392 29 29 IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx"); 30 30 IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy"); 31 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");31 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 32 32 IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate"); 33 33 IoModelFetchData(&iomodel->dhdt,NULL,NULL,iomodel_handle,"dhdt"); … … 58 58 xfree((void**)&iomodel->vx); 59 59 xfree((void**)&iomodel->vy); 60 xfree((void**)&iomodel-> melting_rate);60 xfree((void**)&iomodel->basal_melting_rate); 61 61 xfree((void**)&iomodel->accumulation_rate); 62 62 } -
issm/trunk/src/c/modules/ModelProcessorx/Balancevelocities/UpdateElementsBalancevelocities.cpp
r8330 r8392 29 29 IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx"); 30 30 IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy"); 31 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");31 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 32 32 IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate"); 33 33 … … 59 59 xfree((void**)&iomodel->vy); 60 60 xfree((void**)&iomodel->vz); 61 xfree((void**)&iomodel-> melting_rate);61 xfree((void**)&iomodel->basal_melting_rate); 62 62 xfree((void**)&iomodel->accumulation_rate); 63 63 } -
issm/trunk/src/c/modules/ModelProcessorx/CreateParameters.cpp
r8335 r8392 65 65 parameters->AddObject(new DoubleParam(HeatCapacityEnum,iomodel->heatcapacity)); 66 66 parameters->AddObject(new IntParam(ArtDiffEnum,iomodel->artdiff)); 67 parameters->AddObject(new IntParam( MeltingRateCorrectionApplyEnum,iomodel->melting_rate_correction_apply));67 parameters->AddObject(new IntParam(BasalMeltingRateCorrectionApplyEnum,iomodel->basal_melting_rate_correction_apply)); 68 68 parameters->AddObject(new DoubleParam(GroundingLineMeltingRateEnum,iomodel->gl_melting_rate)); 69 69 parameters->AddObject(new DoubleParam(PenaltyMeltingEnum,iomodel->penalty_melting)); -
issm/trunk/src/c/modules/ModelProcessorx/DiagnosticHoriz/UpdateElementsDiagnosticHoriz.cpp
r8330 r8392 46 46 IoModelFetchData(&iomodel->vz,NULL,NULL,iomodel_handle,"vz"); 47 47 IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate"); 48 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");48 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 49 49 IoModelFetchData(&iomodel->nodeonstokes,NULL,NULL,iomodel_handle,"nodeonstokes"); 50 50 } … … 82 82 xfree((void**)&iomodel->elementonwater); 83 83 xfree((void**)&iomodel->accumulation_rate); 84 xfree((void**)&iomodel-> melting_rate);84 xfree((void**)&iomodel->basal_melting_rate); 85 85 xfree((void**)&iomodel->nodeonstokes); 86 86 xfree((void**)&iomodel->vx); -
issm/trunk/src/c/modules/ModelProcessorx/DiagnosticVert/UpdateElementsDiagnosticVert.cpp
r8330 r8392 33 33 IoModelFetchData(&iomodel->elementonsurface,NULL,NULL,iomodel_handle,"elementonsurface"); 34 34 IoModelFetchData(&iomodel->elementonwater,NULL,NULL,iomodel_handle,"elementonwater"); 35 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");35 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 36 36 IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate"); 37 37 IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx"); … … 59 59 xfree((void**)&iomodel->elementonsurface); 60 60 xfree((void**)&iomodel->elementonwater); 61 xfree((void**)&iomodel-> melting_rate);61 xfree((void**)&iomodel->basal_melting_rate); 62 62 xfree((void**)&iomodel->accumulation_rate); 63 63 xfree((void**)&iomodel->vx); -
issm/trunk/src/c/modules/ModelProcessorx/Hydrology/UpdateElementsHydrology.cpp
r8330 r8392 37 37 IoModelFetchData(&iomodel->pressure,NULL,NULL,iomodel_handle,"pressure"); 38 38 IoModelFetchData(&iomodel->temperature,NULL,NULL,iomodel_handle,"temperature"); 39 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");39 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 40 40 IoModelFetchData(&iomodel->watercolumn,NULL,NULL,iomodel_handle,"watercolumn"); 41 41 IoModelFetchData(&iomodel->geothermalflux,NULL,NULL,iomodel_handle,"geothermalflux"); … … 73 73 xfree((void**)&iomodel->pressure); 74 74 xfree((void**)&iomodel->temperature); 75 xfree((void**)&iomodel-> melting_rate);75 xfree((void**)&iomodel->basal_melting_rate); 76 76 xfree((void**)&iomodel->watercolumn); 77 77 xfree((void**)&iomodel->geothermalflux); -
issm/trunk/src/c/modules/ModelProcessorx/Melting/UpdateElementsMelting.cpp
r8330 r8392 39 39 IoModelFetchData(&iomodel->rheology_n,NULL,NULL,iomodel_handle,"rheology_n"); 40 40 IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate"); 41 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");41 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 42 42 IoModelFetchData(&iomodel->pressure,NULL,NULL,iomodel_handle,"pressure"); 43 43 … … 69 69 xfree((void**)&iomodel->rheology_n); 70 70 xfree((void**)&iomodel->accumulation_rate); 71 xfree((void**)&iomodel-> melting_rate);71 xfree((void**)&iomodel->basal_melting_rate); 72 72 xfree((void**)&iomodel->pressure); 73 73 } -
issm/trunk/src/c/modules/ModelProcessorx/Prognostic/UpdateElementsPrognostic.cpp
r8330 r8392 28 28 IoModelFetchData(&iomodel->elementonwater,NULL,NULL,iomodel_handle,"elementonwater"); 29 29 IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate"); 30 IoModelFetchData(&iomodel-> melting_rate,NULL,NULL,iomodel_handle,"melting_rate");31 if(iomodel-> melting_rate_correction_apply)IoModelFetchData(&iomodel->melting_rate_correction,NULL,NULL,iomodel_handle,"melting_rate_correction");30 IoModelFetchData(&iomodel->basal_melting_rate,NULL,NULL,iomodel_handle,"basal_melting_rate"); 31 if(iomodel->basal_melting_rate_correction_apply)IoModelFetchData(&iomodel->basal_melting_rate_correction,NULL,NULL,iomodel_handle,"basal_melting_rate_correction"); 32 32 IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx"); 33 33 IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy"); … … 60 60 xfree((void**)&iomodel->elementonwater); 61 61 xfree((void**)&iomodel->accumulation_rate); 62 xfree((void**)&iomodel-> melting_rate);63 xfree((void**)&iomodel-> melting_rate_correction);62 xfree((void**)&iomodel->basal_melting_rate); 63 xfree((void**)&iomodel->basal_melting_rate_correction); 64 64 xfree((void**)&iomodel->vx); 65 65 xfree((void**)&iomodel->vy); -
issm/trunk/src/c/modules/StringToEnumx/StringToEnumx.cpp
r8376 r8392 192 192 else if (strcmp(name,"MeltingOffset")==0) return MeltingOffsetEnum; 193 193 else if (strcmp(name,"GroundingLineMeltingRate")==0) return GroundingLineMeltingRateEnum; 194 else if (strcmp(name," MeltingRate")==0) returnMeltingRateEnum;194 else if (strcmp(name,"BasalMeltingRate")==0) return BasalMeltingRateEnum; 195 195 else if (strcmp(name,"Misfit")==0) return MisfitEnum; 196 196 else if (strcmp(name,"NodeOnBed")==0) return NodeOnBedEnum; … … 293 293 else if (strcmp(name,"Gravity")==0) return GravityEnum; 294 294 else if (strcmp(name,"ThermalConductivity")==0) return ThermalConductivityEnum; 295 else if (strcmp(name," MeltingRateCorrection")==0) returnMeltingRateCorrectionEnum;296 else if (strcmp(name," MeltingRateCorrectionApply")==0) returnMeltingRateCorrectionApplyEnum;295 else if (strcmp(name,"BasalMeltingRateCorrection")==0) return BasalMeltingRateCorrectionEnum; 296 else if (strcmp(name,"BasalMeltingRateCorrectionApply")==0) return BasalMeltingRateCorrectionApplyEnum; 297 297 else if (strcmp(name,"MixedLayerCapacity")==0) return MixedLayerCapacityEnum; 298 298 else if (strcmp(name,"ThermalExchangeVelocity")==0) return ThermalExchangeVelocityEnum; -
issm/trunk/src/c/objects/Elements/Penta.cpp
r8386 r8392 4003 4003 4004 4004 } 4005 if (iomodel-> melting_rate) {4006 for(i=0;i<6;i++)nodeinputs[i]=iomodel-> melting_rate[penta_vertex_ids[i]-1]/iomodel->yts;4007 this->inputs->AddInput(new PentaVertexInput( MeltingRateEnum,nodeinputs));4005 if (iomodel->basal_melting_rate) { 4006 for(i=0;i<6;i++)nodeinputs[i]=iomodel->basal_melting_rate[penta_vertex_ids[i]-1]/iomodel->yts; 4007 this->inputs->AddInput(new PentaVertexInput(BasalMeltingRateEnum,nodeinputs)); 4008 4008 } 4009 4009 if (iomodel->accumulation_rate) { … … 4197 4197 } 4198 4198 else if (analysis_type==MeltingAnalysisEnum){ 4199 InputUpdateFromSolutionOneDof(solution, MeltingRateEnum);4199 InputUpdateFromSolutionOneDof(solution,BasalMeltingRateEnum); 4200 4200 } 4201 4201 else{ … … 5306 5306 name==SurfaceSlopeXEnum || 5307 5307 name==SurfaceSlopeYEnum || 5308 name== MeltingRateEnum ||5308 name==BasalMeltingRateEnum || 5309 5309 name==AccumulationRateEnum || 5310 5310 name==GeothermalFluxEnum || -
issm/trunk/src/c/objects/Elements/Tria.cpp
r8386 r8392 1671 1671 GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES); 1672 1672 Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); _assert_(accumulation_input); 1673 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);1673 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 1674 1674 Input* dhdt_input=inputs->GetInput(DhDtEnum); _assert_(dhdt_input); 1675 1675 … … 1714 1714 GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES); 1715 1715 Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); _assert_(accumulation_input); 1716 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);1716 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 1717 1717 Input* dhdt_input=inputs->GetInput(DhDtEnum); _assert_(dhdt_input); 1718 1718 … … 1757 1757 GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES); 1758 1758 Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); _assert_(accumulation_input); 1759 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);1759 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 1760 1760 1761 1761 /* Start looping on the number of gaussian points: */ … … 1802 1802 inputs->GetParameterValue(&approximation,ApproximationEnum); 1803 1803 Input* bed_input=inputs->GetInput(BedEnum); _assert_(bed_input); 1804 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);1804 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 1805 1805 Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input); 1806 1806 Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input); … … 2401 2401 GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES); 2402 2402 this->parameters->FindParam(&dt,DtEnum); 2403 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);2403 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 2404 2404 Input* old_watercolumn_input=inputs->GetInput(WaterColumnOldEnum);_assert_(old_watercolumn_input); 2405 2405 … … 2447 2447 this->parameters->FindParam(&dt,DtEnum); 2448 2448 Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); _assert_(accumulation_input); 2449 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);2450 Input* melting_correction_input=inputs->GetInput( MeltingRateCorrectionEnum);2449 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 2450 Input* melting_correction_input=inputs->GetInput(BasalMeltingRateCorrectionEnum); 2451 2451 Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input); 2452 2452 … … 2495 2495 GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES); 2496 2496 Input* accumulation_input=inputs->GetInput(AccumulationRateEnum); _assert_(accumulation_input); 2497 Input* melting_input=inputs->GetInput( MeltingRateEnum); _assert_(melting_input);2497 Input* melting_input=inputs->GetInput(BasalMeltingRateEnum); _assert_(melting_input); 2498 2498 Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input); 2499 2499 … … 3872 3872 this->inputs->AddInput(new TriaVertexInput(ThicknessObsEnum,nodeinputs)); 3873 3873 } 3874 if (iomodel-> melting_rate) {3875 for(i=0;i<3;i++)nodeinputs[i]=iomodel-> melting_rate[tria_vertex_ids[i]-1]/iomodel->yts;3876 this->inputs->AddInput(new TriaVertexInput( MeltingRateEnum,nodeinputs));3877 } 3878 if (iomodel-> melting_rate_correction) {3879 for(i=0;i<3;i++)nodeinputs[i]=iomodel-> melting_rate_correction[tria_vertex_ids[i]-1]/iomodel->yts;3880 this->inputs->AddInput(new TriaVertexInput( MeltingRateCorrectionEnum,nodeinputs));3874 if (iomodel->basal_melting_rate) { 3875 for(i=0;i<3;i++)nodeinputs[i]=iomodel->basal_melting_rate[tria_vertex_ids[i]-1]/iomodel->yts; 3876 this->inputs->AddInput(new TriaVertexInput(BasalMeltingRateEnum,nodeinputs)); 3877 } 3878 if (iomodel->basal_melting_rate_correction) { 3879 for(i=0;i<3;i++)nodeinputs[i]=iomodel->basal_melting_rate_correction[tria_vertex_ids[i]-1]/iomodel->yts; 3880 this->inputs->AddInput(new TriaVertexInput(BasalMeltingRateCorrectionEnum,nodeinputs)); 3881 3881 } 3882 3882 if (iomodel->watercolumn){ … … 4528 4528 name==SurfaceSlopeXEnum || 4529 4529 name==SurfaceSlopeYEnum || 4530 name== MeltingRateEnum ||4530 name==BasalMeltingRateEnum || 4531 4531 name==WaterColumnEnum || 4532 4532 name==AccumulationRateEnum || … … 5148 5148 /*If this element just became ungrounded, set its melting rate at 50 m/yr:*/ 5149 5149 if(swap){ 5150 Input* melting_rate_input =inputs->GetInput(MeltingRateEnum); _assert_(melting_rate_input);5151 melting_rate_input->GetValuesPtr(&values,NULL); for(i=0;i<3;i++)values[i]=gl_melting_rate/yts;5150 Input* basal_melting_rate_input =inputs->GetInput(BasalMeltingRateEnum); _assert_(basal_melting_rate_input); 5151 basal_melting_rate_input->GetValuesPtr(&values,NULL); for(i=0;i<3;i++)values[i]=gl_melting_rate/yts; 5152 5152 } 5153 5153 -
issm/trunk/src/c/objects/IoModel.cpp
r8365 r8392 85 85 xfree((void**)&this->edges); 86 86 xfree((void**)&this->geothermalflux); 87 xfree((void**)&this-> melting_rate);87 xfree((void**)&this->basal_melting_rate); 88 88 xfree((void**)&this->watercolumn); 89 xfree((void**)&this-> melting_rate_correction);89 xfree((void**)&this->basal_melting_rate_correction); 90 90 xfree((void**)&this->accumulation_rate); 91 91 xfree((void**)&this->forcingtypes); … … 219 219 IoModelFetchData(&this->mixed_layer_capacity,iomodel_handle,"mixed_layer_capacity"); 220 220 IoModelFetchData(&this->thermal_exchange_velocity,iomodel_handle,"thermal_exchange_velocity"); 221 IoModelFetchData(&this-> melting_rate_correction_apply,iomodel_handle,"melting_rate_correction_apply");221 IoModelFetchData(&this->basal_melting_rate_correction_apply,iomodel_handle,"basal_melting_rate_correction_apply"); 222 222 IoModelFetchData(&this->gl_melting_rate,iomodel_handle,"gl_melting_rate"); 223 223 IoModelFetchData(&this->rheology_law,iomodel_handle,"rheology_law"); … … 285 285 this->temperature=NULL; 286 286 this->gl_melting_rate=0; 287 this-> melting_rate=NULL;287 this->basal_melting_rate=NULL; 288 288 this->watercolumn=NULL; 289 this-> melting_rate_correction=NULL;290 this-> melting_rate_correction_apply=0;289 this->basal_melting_rate_correction=NULL; 290 this->basal_melting_rate_correction_apply=0; 291 291 this->geothermalflux=NULL; 292 292 this->elementonbed=NULL; -
issm/trunk/src/c/objects/IoModel.h
r8365 r8392 189 189 190 190 /*basal: */ 191 double* melting_rate;191 double* basal_melting_rate; 192 192 double* watercolumn; 193 193 double gl_melting_rate; 194 double* melting_rate_correction;195 int melting_rate_correction_apply;194 double* basal_melting_rate_correction; 195 int basal_melting_rate_correction_apply; 196 196 double* accumulation_rate; 197 197 double* dhdt; -
issm/trunk/src/c/shared/Numerics/UnitConversion.cpp
r8365 r8392 68 68 case VelObsEnum: scale=yts;break; //m/yr 69 69 case DhDtEnum: scale=yts;break; //m/yr 70 case MeltingRateEnum: scale=yts;break; //m/yr70 case BasalMeltingRateEnum: scale=yts;break; //m/yr 71 71 case AccumulationRateEnum: scale=yts;break; //m/yr 72 72 case MisfitEnum: scale=pow(yts,2);break; //(m/yr)^2 -
issm/trunk/src/c/solutions/steadystate_core.cpp
r8073 r8392 62 62 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,PressureEnum); 63 63 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,TemperatureEnum); 64 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters, MeltingRateEnum);64 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BasalMeltingRateEnum); 65 65 } 66 66 } -
issm/trunk/src/c/solutions/thermal_core.cpp
r6412 r8392 50 50 _printf_(VerboseSolution()," saving results\n"); 51 51 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,TemperatureEnum,i+1,time); 52 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters, MeltingRateEnum,i+1,time);52 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BasalMeltingRateEnum,i+1,time); 53 53 } 54 54 -
issm/trunk/src/c/solutions/transient2d_core.cpp
r8365 r8392 87 87 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BedEnum,step,time); 88 88 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AccumulationRateEnum,step,time); 89 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters, MeltingRateEnum,step,time);89 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BasalMeltingRateEnum,step,time); 90 90 if(gl_migration!=NoneEnum)InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,ElementOnIceShelfEnum,step,time); 91 91 -
issm/trunk/src/c/solutions/transient3d_core.cpp
r8363 r8392 85 85 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BedEnum,step,time); 86 86 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,TemperatureEnum,step,time); 87 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters, MeltingRateEnum,step,time);87 InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BasalMeltingRateEnum,step,time); 88 88 89 89 /*unload results*/ -
issm/trunk/src/m/classes/model.m
r8372 r8392 266 266 temperature=NaN; %temperature solution vector 267 267 gl_melting_rate=NaN; 268 melting_rate=NaN;269 melting_rate_correction=NaN;270 melting_rate_correction_apply=0;268 basal_melting_rate=NaN; 269 basal_melting_rate_correction=NaN; 270 basal_melting_rate_correction_apply=0; 271 271 pressure=NaN; 272 272 … … 425 425 if isfield(structmd,'B'), md.rheology_B=structmd.B; end 426 426 if isfield(structmd,'n'), md.rheology_n=structmd.n; end 427 if isfield(structmd,'melting'), md. melting_rate=structmd.melting; end427 if isfield(structmd,'melting'), md.basal_melting_rate=structmd.melting; end 428 428 if isfield(structmd,'accumulation'), md.accumulation_rate=structmd.accumulation; end 429 429 if isfield(structmd,'numberofgrids'), md.numberofnodes=structmd.numberofgrids; end … … 617 617 md.latentheat=3.34*10^5; 618 618 619 % melting rate correction:619 %basal melting rate correction: 620 620 md.gl_melting_rate=0; 621 md. melting_rate_correction_apply=0; %do not apply621 md.basal_melting_rate_correction_apply=0; %do not apply 622 622 623 623 %ice thermal conductivity lamda (W/m/K) -
issm/trunk/src/m/enum/EnumToModelField.m
r8383 r8392 22 22 case GroundingLineMigrationEnum(), string='gl_migration'; return 23 23 case AccumulationRateEnum(), string='accumulation_rate'; return 24 case MeltingRateEnum(), string='melting_rate'; return24 case BasalMeltingRateEnum(), string='basal_melting_rate'; return 25 25 otherwise, error(['Enum ' num2str(enum) ' not found associated to any model field']); 26 26 -
issm/trunk/src/m/model/collapse.m
r8298 r8392 76 76 77 77 %special for thermal modeling: 78 md. melting_rate=project2d(md,md.melting_rate,1);78 md.basal_melting_rate=project2d(md,md.basal_melting_rate,1); 79 79 md.observed_temperature=DepthAverage(md,md.observed_temperature); 80 80 md.geothermalflux=project2d(md,md.geothermalflux,1); %bedrock only gets geothermal flux -
issm/trunk/src/m/model/display/displaybc.m
r8298 r8392 27 27 disp(sprintf('\n thermal:')); 28 28 fielddisplay(md,'spctemperature','constraints flag list (first column) and values (second column)'); 29 fielddisplay(md,' melting_rate','melting rate [m/a]');29 fielddisplay(md,'basal_melting_rate','basal melting rate [m/a]'); -
issm/trunk/src/m/model/display/displaythermal.m
r5054 r8392 17 17 disp(sprintf('\n boundary conditions:')); 18 18 fielddisplay(md,'spctemperature','constraints flag list (first column) and values (second column)'); 19 fielddisplay(md,' melting_rate','melting rate [m/a]');19 fielddisplay(md,'basal_melting_rate','basal melting rate [m/a]'); 20 20 21 21 disp(sprintf('\n transient:')); -
issm/trunk/src/m/model/extrude.m
r8369 r8392 243 243 244 244 %special for thermal modeling: 245 md. melting_rate=project3d(md,md.melting_rate,'node',1);245 md.basal_melting_rate=project3d(md,md.basal_melting_rate,'node',1); 246 246 md.observed_temperature=project3d(md,md.observed_temperature,'node'); 247 247 if ~isnan(md.geothermalflux) -
issm/trunk/src/m/model/ismodelselfconsistent.m
r8365 r8392 92 92 %}}} 93 93 %SIZE NUMBEROFNODES {{{1 94 fields={'x','y','z','rheology_B','drag_coefficient',' melting_rate','accumulation_rate','surface','thickness','bed','nodeonbed','nodeonsurface'};94 fields={'x','y','z','rheology_B','drag_coefficient','basal_melting_rate','accumulation_rate','surface','thickness','bed','nodeonbed','nodeonsurface'}; 95 95 checksize(md,fields,[md.numberofnodes 1]); 96 96 %}}} … … 505 505 506 506 %INITIAL TEMPERATURE, MELTING AND ACCUMULATION 507 fields={'temperature','accumulation_rate',' melting_rate'};507 fields={'temperature','accumulation_rate','basal_melting_rate'}; 508 508 checksize(md,fields,[md.numberofnodes 1]); 509 509 checknan(md,fields); … … 521 521 % {{{2 522 522 %VELOCITIES MELTING AND ACCUMULATION 523 fields={'vx','vy','accumulation_rate',' melting_rate','dhdt'};523 fields={'vx','vy','accumulation_rate','basal_melting_rate','dhdt'}; 524 524 checksize(md,fields,[md.numberofnodes 1]); 525 525 checknan(md,fields); … … 540 540 % {{{2 541 541 %VELOCITIES MELTING AND ACCUMULATION 542 fields={'vx','vy','accumulation_rate',' melting_rate'};542 fields={'vx','vy','accumulation_rate','basal_melting_rate'}; 543 543 checksize(md,fields,[md.numberofnodes 1]); 544 544 checknan(md,fields); -
issm/trunk/src/m/model/marshall.m
r8372 r8392 91 91 WriteData(fid,md.accumulation_rate,'Mat','accumulation_rate'); 92 92 WriteData(fid,md.gl_melting_rate,'Scalar','gl_melting_rate'); 93 WriteData(fid,md. melting_rate,'Mat','melting_rate');94 WriteData(fid,md. melting_rate_correction_apply,'Integer','melting_rate_correction_apply');93 WriteData(fid,md.basal_melting_rate,'Mat','basal_melting_rate'); 94 WriteData(fid,md.basal_melting_rate_correction_apply,'Integer','basal_melting_rate_correction_apply'); 95 95 96 96 %deal with forcings … … 103 103 case 'accumulation_rate' 104 104 forcingtypes(i)=AccumulationRateEnum; 105 case ' melting_rate'106 forcingtypes(i)= MeltingRateEnum;105 case 'basal_melting_rate' 106 forcingtypes(i)=BasalMeltingRateEnum; 107 107 otherwise 108 108 error(['forcing ' forcingnames{i} ' not supported yet']); … … 126 126 end 127 127 128 if md. melting_rate_correction_apply,129 WriteData(fid,md. melting_rate_correction,'Mat','melting_rate_correction');128 if md.basal_melting_rate_correction_apply, 129 WriteData(fid,md.basal_melting_rate_correction,'Mat','basal_melting_rate_correction'); 130 130 end 131 131 WriteData(fid,md.dhdt,'Mat','dhdt'); -
issm/trunk/src/m/model/modeldefault/defaultparams.m
r8298 r8392 94 94 95 95 disp(' thermal model'); 96 md. melting_rate=zeros(md.numberofnodes,1);96 md.basal_melting_rate=zeros(md.numberofnodes,1); 97 97 md.observed_temperature=md.temperature; 98 98 -
issm/trunk/src/m/model/structtomodel.m
r8323 r8392 25 25 if isfield(structmd,'B'), md.rheology_B=structmd.B; end 26 26 if isfield(structmd,'n'), md.rheology_n=structmd.n; end 27 if isfield(structmd,'melting'), md. melting_rate=structmd.melting; end27 if isfield(structmd,'melting'), md.basal_melting_rate=structmd.melting; end 28 28 if isfield(structmd,'accumulation'), md.accumulation_rate=structmd.accumulation; end 29 29 if isfield(structmd,'numberofgrids'), md.numberofnodes=structmd.numberofgrids; end -
issm/trunk/src/m/model/tres.m
r8298 r8392 84 84 md.pressure=PatchToVec(md.results.SteadystateSolution.Pressure); 85 85 md.temperature=PatchToVec(md.results.SteadystateSolution.Temperature); 86 md. melting_rate=PatchToVec(md.results.SteadystateSolution.MeltingRate);86 md.basal_melting_rate=PatchToVec(md.results.SteadystateSolution.BasalMeltingRate); 87 87 88 88 if md.control_analysis==1, … … 94 94 elseif strcmpi(string,'thermal'), 95 95 md.temperature=PatchToVec(md.results.ThermalSolution.Temperature); 96 md. melting_rate=PatchToVec(md.results.ThermalSolution.MeltingRate);96 md.basal_melting_rate=PatchToVec(md.results.ThermalSolution.BasalMeltingRate); 97 97 elseif strcmpi(string,'hydrology'), 98 98 md.watercolumn=PatchToVec(md.results.HydrologySolution.WaterColumn); -
issm/trunk/src/m/solutions/steadystate_core.m
r8075 r8392 46 46 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum); 47 47 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum); 48 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters, MeltingRateEnum);48 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BasalMeltingRateEnum); 49 49 end 50 50 -
issm/trunk/src/m/solutions/thermal_core.m
r6323 r8392 31 31 issmprintf(VerboseSolution,'\n%s',[' saving results']); 32 32 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,i,time); 33 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters, MeltingRateEnum,i,time);33 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BasalMeltingRateEnum,i,time); 34 34 end 35 35 end -
issm/trunk/src/m/solutions/transient3d_core.m
r8363 r8392 65 65 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time); 66 66 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,step,time); 67 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters, MeltingRateEnum,step,time);67 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BasalMeltingRateEnum,step,time); 68 68 end 69 69 end -
issm/trunk/src/m/utils/BC/SetIceSheetBC.m
r8367 r8392 28 28 end 29 29 30 %Create zeros melting_rate and accumulation_rate if not specified30 %Create zeros basal_melting_rate and accumulation_rate if not specified 31 31 if isnan(md.accumulation_rate), 32 32 md.accumulation_rate=zeros(md.numberofnodes,1); … … 34 34 disp(' no accumulation_rate specified: values set as zero'); 35 35 end 36 if isnan(md. melting_rate),37 md. melting_rate=zeros(md.numberofnodes,1);38 disp(' no melting_rate specified: values set as zero');36 if isnan(md.basal_melting_rate), 37 md.basal_melting_rate=zeros(md.numberofnodes,1); 38 disp(' no basal_melting_rate specified: values set as zero'); 39 39 end 40 40 if isnan(md.dhdt), -
issm/trunk/src/m/utils/BC/SetIceShelfBC.m
r8367 r8392 59 59 md.pressureload=pressureload; 60 60 61 %Create zeros melting_rate and accumulation_rate if not specified61 %Create zeros basal_melting_rate and accumulation_rate if not specified 62 62 if isnan(md.accumulation_rate), 63 63 md.accumulation_rate=zeros(md.numberofnodes,1); … … 65 65 disp(' no accumulation_rate specified: values set as zero'); 66 66 end 67 if isnan(md. melting_rate),68 md. melting_rate=zeros(md.numberofnodes,1);69 disp(' no melting_rate specified: values set as zero');67 if isnan(md.basal_melting_rate), 68 md.basal_melting_rate=zeros(md.numberofnodes,1); 69 disp(' no basal_melting_rate specified: values set as zero'); 70 70 end 71 71 if isnan(md.dhdt), -
issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
r8367 r8392 69 69 md.pressureload=pressureload; 70 70 71 %Create zeros melting_rate and accumulation_rate if not specified71 %Create zeros basal_melting_rate and accumulation_rate if not specified 72 72 if isnan(md.accumulation_rate), 73 73 md.accumulation_rate=zeros(md.numberofnodes,1); … … 75 75 disp(' no accumulation_rate specified: values set as zero'); 76 76 end 77 if isnan(md. melting_rate),78 md. melting_rate=zeros(md.numberofnodes,1);79 disp(' no melting_rate specified: values set as zero');77 if isnan(md.basal_melting_rate), 78 md.basal_melting_rate=zeros(md.numberofnodes,1); 79 disp(' no basal_melting_rate specified: values set as zero'); 80 80 end 81 81 if isnan(md.dhdt), -
issm/trunk/test/NightlyRun/test115.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test116.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test117.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature',' MeltingRate'};10 field_names ={'Temperature','BasalMeltingRate'}; 11 11 field_tolerances={1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution.Temperature),... 14 PatchToVec(md.results.ThermalSolution. MeltingRate),...14 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 15 15 }; -
issm/trunk/test/NightlyRun/test118.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature',' MeltingRate'};10 field_names ={'Temperature','BasalMeltingRate'}; 11 11 field_tolerances={1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution.Temperature),... 14 PatchToVec(md.results.ThermalSolution. MeltingRate),...14 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 15 15 }; -
issm/trunk/test/NightlyRun/test1208.m
r8384 r8392 26 26 27 27 %Fields and tolerances to track changes 28 field_names ={'Vx','Vy','Vz','Vel','Pressure','Thickness','Bed','Surface','Temperature',' MeltingRate'};28 field_names ={'Vx','Vy','Vz','Vel','Pressure','Thickness','Bed','Surface','Temperature','BasalMeltingRate'}; 29 29 field_tolerances={1e-08,1e-08,1e-07,1e-08,1e-08,1e-08,1e-08,1e-08,1e-07,1e-07}; 30 30 field_values={... … … 38 38 PatchToVec(md.results.Transient3DSolution(end).Surface),... 39 39 PatchToVec(md.results.Transient3DSolution(end).Temperature),... 40 PatchToVec(md.results.Transient3DSolution(end). MeltingRate),...40 PatchToVec(md.results.Transient3DSolution(end).BasalMeltingRate),... 41 41 }; -
issm/trunk/test/NightlyRun/test121.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test122.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test1301.m
r8304 r8392 29 29 %plot results 30 30 comp_melting=zeros(md.numberofnodes,1); 31 comp_melting(md.results.ThermalSolution. MeltingRate.index)=md.results.ThermalSolution.MeltingRate.value;31 comp_melting(md.results.ThermalSolution.BasalMeltingRate.index)=md.results.ThermalSolution.BasalMeltingRate.value; 32 32 relative=abs((comp_melting-melting)./melting)*100; 33 33 relative(find(comp_melting==melting))=0; … … 43 43 44 44 %Fields and tolerances to track changes 45 field_names ={' Melting'};45 field_names ={'BasalMelting'}; 46 46 field_tolerances={1e-13}; 47 47 field_values ={comp_melting}; -
issm/trunk/test/NightlyRun/test211.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test212.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test213.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature1','MeltingRate1','Temperature2',' MeltingRate2','Temperature3','MeltingRate3'};10 field_names ={'Temperature1','MeltingRate1','Temperature2','BasalMeltingRate2','Temperature3','BasalMeltingRate3'}; 11 11 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution(1).Temperature),... 14 PatchToVec(md.results.ThermalSolution(1). MeltingRate),...14 PatchToVec(md.results.ThermalSolution(1).BasalMeltingRate),... 15 15 PatchToVec(md.results.ThermalSolution(2).Temperature),... 16 PatchToVec(md.results.ThermalSolution(2). MeltingRate),...16 PatchToVec(md.results.ThermalSolution(2).BasalMeltingRate),... 17 17 PatchToVec(md.results.ThermalSolution(3).Temperature),... 18 PatchToVec(md.results.ThermalSolution(3). MeltingRate),...18 PatchToVec(md.results.ThermalSolution(3).BasalMeltingRate),... 19 19 }; -
issm/trunk/test/NightlyRun/test214.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature1',' MeltingRate1','Temperature2','MeltingRate2','Temperature3','MeltingRate3'};10 field_names ={'Temperature1','BasalMeltingRate1','Temperature2','BasalMeltingRate2','Temperature3','BasalMeltingRate3'}; 11 11 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution(1).Temperature),... 14 PatchToVec(md.results.ThermalSolution(1). MeltingRate),...14 PatchToVec(md.results.ThermalSolution(1).BasalMeltingRate),... 15 15 PatchToVec(md.results.ThermalSolution(2).Temperature),... 16 PatchToVec(md.results.ThermalSolution(2). MeltingRate),...16 PatchToVec(md.results.ThermalSolution(2).BasalMeltingRate),... 17 17 PatchToVec(md.results.ThermalSolution(3).Temperature),... 18 PatchToVec(md.results.ThermalSolution(3). MeltingRate),...18 PatchToVec(md.results.ThermalSolution(3).BasalMeltingRate),... 19 19 }; -
issm/trunk/test/NightlyRun/test217.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test218.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test219.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test220.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={... 14 14 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,... … … 25 25 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 26 26 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 27 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...27 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vx),... 29 29 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 35 35 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 36 36 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 37 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...37 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 38 38 PatchToVec(md.results.Transient3DSolution(3).Vx),... 39 39 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 45 45 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 46 46 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 47 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...47 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 48 48 }; -
issm/trunk/test/NightlyRun/test221.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 14 14 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08}; 15 15 field_values={... … … 23 23 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 24 24 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 25 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...25 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vx),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 33 33 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 34 34 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 35 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...35 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vx),... 37 37 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 43 43 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 44 44 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 45 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...45 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 46 46 }; -
issm/trunk/test/NightlyRun/test222.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 14 14 field_tolerances={... 15 15 1e-08,1e-08,1e-06,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,... … … 26 26 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 27 27 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 28 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...28 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 29 29 PatchToVec(md.results.Transient3DSolution(2).Vx),... 30 30 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 36 36 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 37 37 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 38 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...38 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 39 39 PatchToVec(md.results.Transient3DSolution(3).Vx),... 40 40 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 46 46 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 47 47 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 48 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...48 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 49 49 }; -
issm/trunk/test/NightlyRun/test247.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 14 14 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06}; 15 15 field_values={... … … 23 23 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 24 24 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 25 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...25 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vx),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 33 33 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 34 34 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 35 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...35 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vx),... 37 37 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 43 43 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 44 44 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 45 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...45 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 46 46 }; -
issm/trunk/test/NightlyRun/test248.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 13 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 14 14 field_tolerances={... 15 15 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,... … … 26 26 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 27 27 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 28 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...28 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 29 29 PatchToVec(md.results.Transient3DSolution(2).Vx),... 30 30 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 36 36 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 37 37 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 38 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...38 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 39 39 PatchToVec(md.results.Transient3DSolution(3).Vx),... 40 40 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 46 46 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 47 47 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 48 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...48 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 49 49 }; -
issm/trunk/test/NightlyRun/test249.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 15 15 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06}; 16 16 field_values={... … … 24 24 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 25 25 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 26 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...26 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vx),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 34 34 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 35 35 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 36 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...36 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 37 37 PatchToVec(md.results.Transient3DSolution(3).Vx),... 38 38 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 44 44 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 45 45 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 46 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...46 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 47 47 }; -
issm/trunk/test/NightlyRun/test250.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 15 15 field_tolerances={... 16 16 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,... … … 27 27 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 28 28 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 29 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...29 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 30 30 PatchToVec(md.results.Transient3DSolution(2).Vx),... 31 31 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 37 37 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 38 38 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 39 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...39 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 40 40 PatchToVec(md.results.Transient3DSolution(3).Vx),... 41 41 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 47 47 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 48 48 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 49 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...49 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 50 50 }; -
issm/trunk/test/NightlyRun/test253.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 15 15 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06}; 16 16 field_values={... … … 24 24 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 25 25 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 26 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...26 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vx),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 34 34 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 35 35 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 36 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...36 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 37 37 PatchToVec(md.results.Transient3DSolution(3).Vx),... 38 38 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 44 44 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 45 45 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 46 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...46 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 47 47 }; -
issm/trunk/test/NightlyRun/test254.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 15 15 field_tolerances={... 16 16 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,... … … 27 27 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 28 28 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 29 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...29 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 30 30 PatchToVec(md.results.Transient3DSolution(2).Vx),... 31 31 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 37 37 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 38 38 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 39 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...39 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 40 40 PatchToVec(md.results.Transient3DSolution(3).Vx),... 41 41 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 47 47 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 48 48 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 49 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...49 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 50 50 }; -
issm/trunk/test/NightlyRun/test323.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test324.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test325.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature',' MeltingRate'};10 field_names ={'Temperature','BasalMeltingRate'}; 11 11 field_tolerances={1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution.Temperature),... 14 PatchToVec(md.results.ThermalSolution. MeltingRate),...14 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 15 15 }; -
issm/trunk/test/NightlyRun/test326.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature',' MeltingRate'};10 field_names ={'Temperature','BasalMeltingRate'}; 11 11 field_tolerances={1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution.Temperature),... 14 PatchToVec(md.results.ThermalSolution. MeltingRate),...14 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 15 15 }; -
issm/trunk/test/NightlyRun/test329.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test330.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-13, 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,..., 14 14 5*1e-11,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... … … 24 24 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 25 25 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 26 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...26 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vx),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 34 34 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 35 35 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 36 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...36 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 37 37 PatchToVec(md.results.Transient3DSolution(3).Vx),... 38 38 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 44 44 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 45 45 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 46 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...46 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 47 47 }; -
issm/trunk/test/NightlyRun/test333.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test334.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... 14 14 1e-09,1e-09,1e-10,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... … … 24 24 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 25 25 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 26 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...26 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vx),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 34 34 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 35 35 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 36 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...36 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 37 37 PatchToVec(md.results.Transient3DSolution(3).Vx),... 38 38 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 44 44 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 45 45 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 46 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...46 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 47 47 }; -
issm/trunk/test/NightlyRun/test335.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'}; 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 13 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 13 14 }; … … 19 20 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 21 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...22 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 23 }; -
issm/trunk/test/NightlyRun/test336.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test347.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};14 %'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 14 %'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 15 15 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 16 16 field_values={... … … 24 24 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 25 25 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 26 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...26 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 27 27 PatchToVec(md.results.Transient3DSolution(2).Vx),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 34 34 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 35 35 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 36 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...36 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 37 37 % PatchToVec(md.results.Transient3DSolution(3).Vx),... 38 38 % PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 44 44 % PatchToVec(md.results.Transient3DSolution(3).Thickness),... 45 45 % PatchToVec(md.results.Transient3DSolution(3).Temperature),... 46 % PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...46 % PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 47 47 }; -
issm/trunk/test/NightlyRun/test348.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};14 %'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};12 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 14 %'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 15 15 field_tolerances={1e-13, 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,..., 16 16 5*1e-11,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... … … 26 26 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 27 27 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 28 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...28 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 29 29 PatchToVec(md.results.Transient3DSolution(2).Vx),... 30 30 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 36 36 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 37 37 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 38 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...38 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 39 39 % PatchToVec(md.results.Transient3DSolution(3).Vx),... 40 40 % PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 46 46 % PatchToVec(md.results.Transient3DSolution(3).Thickness),... 47 47 % PatchToVec(md.results.Transient3DSolution(3).Temperature),... 48 % PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...48 % PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 49 49 }; -
issm/trunk/test/NightlyRun/test411.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test412.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature',' MeltingRate'};11 field_names ={'Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution.Temperature),... 15 PatchToVec(md.results.ThermalSolution. MeltingRate),...15 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 16 16 }; -
issm/trunk/test/NightlyRun/test413.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature',' MeltingRate'};10 field_names ={'Temperature','BasalMeltingRate'}; 11 11 field_tolerances={1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution.Temperature),... 14 PatchToVec(md.results.ThermalSolution. MeltingRate),...14 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 15 15 }; -
issm/trunk/test/NightlyRun/test414.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature',' MeltingRate'};10 field_names ={'Temperature','BasalMeltingRate'}; 11 11 field_tolerances={1e-13,1e-13}; 12 12 field_values={... 13 13 PatchToVec(md.results.ThermalSolution.Temperature),... 14 PatchToVec(md.results.ThermalSolution. MeltingRate),...14 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 15 15 }; -
issm/trunk/test/NightlyRun/test417.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 PatchToVec(md.results.Transient3DSolution(3).Vx),... 36 36 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 42 42 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 43 43 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 44 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...44 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 45 45 }; -
issm/trunk/test/NightlyRun/test418.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2', ...12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3',' MeltingRate3'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2', ... 12 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalMeltingRate3'}; 13 13 field_tolerances={... 14 14 1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,1e-05,... … … 25 25 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 26 26 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 27 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...27 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 28 28 PatchToVec(md.results.Transient3DSolution(2).Vx),... 29 29 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 35 35 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 36 36 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 37 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...37 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 38 38 PatchToVec(md.results.Transient3DSolution(3).Vx),... 39 39 PatchToVec(md.results.Transient3DSolution(3).Vy),... … … 45 45 PatchToVec(md.results.Transient3DSolution(3).Thickness),... 46 46 PatchToVec(md.results.Transient3DSolution(3).Temperature),... 47 PatchToVec(md.results.Transient3DSolution(3). MeltingRate),...47 PatchToVec(md.results.Transient3DSolution(3).BasalMeltingRate),... 48 48 }; -
issm/trunk/test/NightlyRun/test419.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test420.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test421.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06}; 13 13 field_values={... … … 18 18 PatchToVec(md.results.SteadystateSolution.Pressure),... 19 19 PatchToVec(md.results.SteadystateSolution.Temperature),... 20 PatchToVec(md.results.SteadystateSolution. MeltingRate),...20 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 21 21 }; -
issm/trunk/test/NightlyRun/test422.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-05 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test429.m
r8304 r8392 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' ' MeltingRate'};26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalMeltingRate'}; 27 27 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 28 28 field_values={... … … 36 36 PatchToVec(md.results.SteadystateSolution.Vz),... 37 37 PatchToVec(md.results.SteadystateSolution.Temperature),... 38 PatchToVec(md.results.SteadystateSolution. MeltingRate)38 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate) 39 39 }; -
issm/trunk/test/NightlyRun/test430.m
r8304 r8392 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' ' MeltingRate'};26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalMeltingRate'}; 27 27 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 28 28 field_values={... … … 36 36 PatchToVec(md.results.SteadystateSolution.Vz),... 37 37 PatchToVec(md.results.SteadystateSolution.Temperature),... 38 PatchToVec(md.results.SteadystateSolution. MeltingRate)38 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate) 39 39 }; -
issm/trunk/test/NightlyRun/test431.m
r8304 r8392 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' ' MeltingRate'};26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalMeltingRate'}; 27 27 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06}; 28 28 field_values={... … … 36 36 PatchToVec(md.results.SteadystateSolution.Vz),... 37 37 PatchToVec(md.results.SteadystateSolution.Temperature),... 38 PatchToVec(md.results.SteadystateSolution. MeltingRate)38 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate) 39 39 }; -
issm/trunk/test/NightlyRun/test432.m
r8304 r8392 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' ' MeltingRate'};26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalMeltingRate'}; 27 27 field_tolerances={1e-08,1e-07,1e-08,1e-08,1e-08,1e-08,1e-08,1e-07,1e-08,1e-05}; 28 28 field_values={... … … 36 36 PatchToVec(md.results.SteadystateSolution.Vz),... 37 37 PatchToVec(md.results.SteadystateSolution.Temperature),... 38 PatchToVec(md.results.SteadystateSolution. MeltingRate)38 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate) 39 39 }; -
issm/trunk/test/NightlyRun/test509.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 12 12 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12}; 13 13 field_values={... … … 21 21 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 22 22 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 23 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...23 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 24 24 PatchToVec(md.results.Transient3DSolution(2).Vx),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 31 31 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 32 32 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 33 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...33 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 34 34 }; -
issm/trunk/test/NightlyRun/test510.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 12 12 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-12,1e-12,1e-12,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 13 13 field_values={... … … 21 21 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 22 22 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 23 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...23 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 24 24 PatchToVec(md.results.Transient3DSolution(2).Vx),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 31 31 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 32 32 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 33 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...33 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 34 34 }; -
issm/trunk/test/NightlyRun/test511.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 12 12 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 13 13 field_values={... … … 21 21 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 22 22 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 23 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...23 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 24 24 PatchToVec(md.results.Transient3DSolution(2).Vx),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 31 31 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 32 32 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 33 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...33 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 34 34 }; -
issm/trunk/test/NightlyRun/test512.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};10 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 11 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 12 12 field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-12,1e-11,1e-12,1e-11,1e-12,1e-12,... 13 13 1e-11,1e-11,1e-09,1e-11,1e-11,1e-10,1e-11,1e-10,1e-11,1e-08}; … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 }; -
issm/trunk/test/NightlyRun/test513.m
r8295 r8392 9 9 %Fields and tolerances to track changes 10 10 %The thickness wants to be lower than 1 so we contrain it to 1 11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 13 13 field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 }; -
issm/trunk/test/NightlyRun/test514.m
r8295 r8392 9 9 %Fields and tolerances to track changes 10 10 %The thickness wants to be lower than 1 so we contrain it to 1 11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1',' MeltingRate1', ...12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2',' MeltingRate2'};11 field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalMeltingRate1', ... 12 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalMeltingRate2'}; 13 13 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06,1e-06,1e-06,1e-06,1e-06,1e-06,1e-06,1e-06,1e-06,1e-06}; 14 14 field_values={... … … 22 22 PatchToVec(md.results.Transient3DSolution(1).Thickness),... 23 23 PatchToVec(md.results.Transient3DSolution(1).Temperature),... 24 PatchToVec(md.results.Transient3DSolution(1). MeltingRate),...24 PatchToVec(md.results.Transient3DSolution(1).BasalMeltingRate),... 25 25 PatchToVec(md.results.Transient3DSolution(2).Vx),... 26 26 PatchToVec(md.results.Transient3DSolution(2).Vy),... … … 32 32 PatchToVec(md.results.Transient3DSolution(2).Thickness),... 33 33 PatchToVec(md.results.Transient3DSolution(2).Temperature),... 34 PatchToVec(md.results.Transient3DSolution(2). MeltingRate),...34 PatchToVec(md.results.Transient3DSolution(2).BasalMeltingRate),... 35 35 }; -
issm/trunk/test/NightlyRun/test515.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test516.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-08 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test517.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test518.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test519.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test520.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature',' MeltingRate'};11 field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalMeltingRate'}; 12 12 field_tolerances={1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-06 13 13 }; … … 19 19 PatchToVec(md.results.SteadystateSolution.Pressure),... 20 20 PatchToVec(md.results.SteadystateSolution.Temperature),... 21 PatchToVec(md.results.SteadystateSolution. MeltingRate),...21 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate),... 22 22 }; -
issm/trunk/test/NightlyRun/test525.m
r8304 r8392 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' ' MeltingRate'};26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalMeltingRate'}; 27 27 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 28 28 field_values={... … … 36 36 PatchToVec(md.results.SteadystateSolution.Vz),... 37 37 PatchToVec(md.results.SteadystateSolution.Temperature),... 38 PatchToVec(md.results.SteadystateSolution. MeltingRate)38 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate) 39 39 }; -
issm/trunk/test/NightlyRun/test526.m
r8304 r8392 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' ' MeltingRate'};26 field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalMeltingRate'}; 27 27 field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-6,1e-6,1e-6,1e-6,1e-6,1e-06}; 28 28 field_values={... … … 36 36 PatchToVec(md.results.SteadystateSolution.Vz),... 37 37 PatchToVec(md.results.SteadystateSolution.Temperature),... 38 PatchToVec(md.results.SteadystateSolution. MeltingRate)38 PatchToVec(md.results.SteadystateSolution.BasalMeltingRate) 39 39 }; -
issm/trunk/test/NightlyRun/test529.m
r8295 r8392 8 8 9 9 %Fields and tolerances to track changes 10 field_names ={'Temperature1',' MeltingRate1', ...11 'Temperature2',' MeltingRate2'};10 field_names ={'Temperature1','BasalMeltingRate1', ... 11 'Temperature2','BasalMeltingRate2'}; 12 12 field_tolerances={1e-13,1e-9,1e-13,1e-9}; 13 13 field_values={... 14 14 PatchToVec(md.results.ThermalSolution(1).Temperature),... 15 PatchToVec(md.results.ThermalSolution(1). MeltingRate),...15 PatchToVec(md.results.ThermalSolution(1).BasalMeltingRate),... 16 16 PatchToVec(md.results.ThermalSolution(2).Temperature),... 17 PatchToVec(md.results.ThermalSolution(2). MeltingRate),...17 PatchToVec(md.results.ThermalSolution(2).BasalMeltingRate),... 18 18 }; -
issm/trunk/test/NightlyRun/test530.m
r8295 r8392 9 9 10 10 %Fields and tolerances to track changes 11 field_names ={'Temperature1',' MeltingRate1', ...12 'Temperature2',' MeltingRate2'};11 field_names ={'Temperature1','BasalMeltingRate1', ... 12 'Temperature2','BasalMeltingRate2'}; 13 13 field_tolerances={1e-13,1e-8,1e-13,1e-8}; 14 14 field_values={... 15 15 PatchToVec(md.results.ThermalSolution(1).Temperature),... 16 PatchToVec(md.results.ThermalSolution(1). MeltingRate),...16 PatchToVec(md.results.ThermalSolution(1).BasalMeltingRate),... 17 17 PatchToVec(md.results.ThermalSolution(2).Temperature),... 18 PatchToVec(md.results.ThermalSolution(2). MeltingRate),...18 PatchToVec(md.results.ThermalSolution(2).BasalMeltingRate),... 19 19 }; -
issm/trunk/test/NightlyRun/test531.m
r8295 r8392 10 10 11 11 %Fields and tolerances to track changes 12 field_names ={'Temperature',' MeltingRate'};12 field_names ={'Temperature','BasalMeltingRate'}; 13 13 field_tolerances={1e-13,1e-12}; 14 14 field_values={... 15 15 PatchToVec(md.results.ThermalSolution.Temperature),... 16 PatchToVec(md.results.ThermalSolution. MeltingRate),...16 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 17 17 }; -
issm/trunk/test/NightlyRun/test532.m
r8295 r8392 11 11 12 12 %Fields and tolerances to track changes 13 field_names ={'Temperature',' MeltingRate'};13 field_names ={'Temperature','BasalMeltingRate'}; 14 14 field_tolerances={1e-11,1e-11}; 15 15 field_values={... 16 16 PatchToVec(md.results.ThermalSolution.Temperature),... 17 PatchToVec(md.results.ThermalSolution. MeltingRate),...17 PatchToVec(md.results.ThermalSolution.BasalMeltingRate),... 18 18 }; -
issm/trunk/test/Par/79North.par
r8368 r8392 24 24 25 25 %Ice shelf melting and accumulation 26 md. melting_rate=zeros(md.numberofnodes,1);26 md.basal_melting_rate=zeros(md.numberofnodes,1); 27 27 pos=zeros(md.numberofnodes,1); 28 28 pos(md.elements(find(md.elementoniceshelf),:))=1; 29 md. melting_rate(find(pos))=10;29 md.basal_melting_rate(find(pos))=10; 30 30 md.accumulation_rate=15*ones(md.numberofnodes,1); 31 31 md.forcings.accumulation_rate=[15*ones(md.numberofnodes,1);1]; -
issm/trunk/test/Par/RoundSheetShelf.par
r8368 r8392 48 48 md.accumulation_rate=-10*ones(md.numberofnodes,1); 49 49 md.forcings.accumulation_rate=[-10*ones(md.numberofnodes,1);1]; 50 md. melting_rate=zeros(md.numberofnodes,1);51 pos=find(md.nodeoniceshelf);md. melting_rate(pos)=10;50 md.basal_melting_rate=zeros(md.numberofnodes,1); 51 pos=find(md.nodeoniceshelf);md.basal_melting_rate(pos)=10; 52 52 53 53 %Friction -
issm/trunk/test/Par/SquareEISMINT.par
r8368 r8392 29 29 md.accumulation_rate=0.2*ones(md.numberofnodes,1); %0m/a 30 30 md.forcings.accumulation_rate=[0.2*ones(md.numberofnodes,1);1]; %0m/a 31 md. melting_rate=0*ones(md.numberofnodes,1); %0m/a31 md.basal_melting_rate=0*ones(md.numberofnodes,1); %0m/a 32 32 33 33 disp(' boundary conditions '); -
issm/trunk/test/Par/SquareSheetShelf.par
r8368 r8392 30 30 md.accumulation_rate=10*ones(md.numberofnodes,1); 31 31 md.forcings.accumulation_rate=[10*ones(md.numberofnodes,1);1]; 32 md. melting_rate=5*ones(md.numberofnodes,1);32 md.basal_melting_rate=5*ones(md.numberofnodes,1); 33 33 34 34 %Friction -
issm/trunk/test/Par/SquareShelfConstrained.par
r8368 r8392 27 27 md.accumulation_rate=10*ones(md.numberofnodes,1); 28 28 md.forcings.accumulation_rate=[10*ones(md.numberofnodes,1);1]; 29 md. melting_rate=5*ones(md.numberofnodes,1);29 md.basal_melting_rate=5*ones(md.numberofnodes,1); 30 30 31 31 %Friction -
issm/trunk/test/Par/SquareThermal.par
r8368 r8392 36 36 md.accumulation_rate=ones(md.numberofnodes,1)/md.yts; %1m/a 37 37 md.forcings.accumulation_rate=[ones(md.numberofnodes,1)/md.yts;1]; %1m/a 38 md. melting_rate=0*ones(md.numberofnodes,1)/md.yts; %1m/a38 md.basal_melting_rate=0*ones(md.numberofnodes,1)/md.yts; %1m/a 39 39 40 40 %Deal with boundary conditions:
Note:
See TracChangeset
for help on using the changeset viewer.