Context Navigation

source: issm/trunk-jpl/src/c/analyses/SmbAnalysis.cpp@ 24253

Last change on this file since 24253 was 24253, checked in by Mathieu Morlighem, 5 years ago
CHG: not needed
File size: 20.0 KB

Rev	Line
[19528]	1	#include "./SmbAnalysis.h"
	2	#include "../toolkits/toolkits.h"
	3	#include "../classes/classes.h"
	4	#include "../shared/shared.h"
	5	#include "../modules/modules.h"
	6
[23317]	7	// FIX
	8	#include "./shared/io/Print/Print.h"
	9
[19528]	10	/Model processing/
	11	void SmbAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/{{{/
	12	/No constraints/
	13	}/}}}/
	14	void SmbAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/{{{/
	15	/No loads/
	16	}/}}}/
[23585]	17	void SmbAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel,bool isamr){/{{{/
[19528]	18	::CreateNodes(nodes,iomodel,SmbAnalysisEnum,P1Enum);
	19	}/}}}/
	20	int SmbAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/{{{/
	21	return 1;
	22	}/}}}/
	23	void SmbAnalysis::UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type){/{{{/
[23066]	24
[19528]	25	int smb_model;
[23317]	26	bool isdelta18o,ismungsm,isd18opd,issetpddfac,isprecipscaled,istemperaturescaled,isfirnwarming;
[23066]	27
[19528]	28	/Update elements: /
	29	int counter=0;
	30	for(int i=0;i<iomodel->numberofelements;i++){
	31	if(iomodel->my_elements[i]){
	32	Element* element=(Element*)elements->GetObjectByOffset(counter);
	33	element->Update(i,iomodel,analysis_counter,analysis_type,P1Enum);
	34	counter++;
	35	}
	36	}
[23066]	37
[19554]	38	/Figure out smb model: /
[20690]	39	iomodel->FindConstant(&smb_model,"md.smb.model");
[19528]	40	switch(smb_model){
	41	case SMBforcingEnum:
[20690]	42	iomodel->FetchDataToInput(elements,"md.smb.mass_balance",SmbMassBalanceEnum,0.);
[19528]	43	break;
[19554]	44	case SMBgembEnum:
[20690]	45	iomodel->FetchDataToInput(elements,"md.smb.Ta",SmbTaEnum);
	46	iomodel->FetchDataToInput(elements,"md.smb.V",SmbVEnum);
	47	iomodel->FetchDataToInput(elements,"md.smb.dswrf",SmbDswrfEnum);
	48	iomodel->FetchDataToInput(elements,"md.smb.dlwrf",SmbDlwrfEnum);
	49	iomodel->FetchDataToInput(elements,"md.smb.P",SmbPEnum);
	50	iomodel->FetchDataToInput(elements,"md.smb.eAir",SmbEAirEnum);
	51	iomodel->FetchDataToInput(elements,"md.smb.pAir",SmbPAirEnum);
	52	iomodel->FetchDataToInput(elements,"md.smb.zTop",SmbZTopEnum);
	53	iomodel->FetchDataToInput(elements,"md.smb.dzTop",SmbDzTopEnum);
	54	iomodel->FetchDataToInput(elements,"md.smb.dzMin",SmbDzMinEnum);
	55	iomodel->FetchDataToInput(elements,"md.smb.zY",SmbZYEnum);
	56	iomodel->FetchDataToInput(elements,"md.smb.zMax",SmbZMaxEnum);
	57	iomodel->FetchDataToInput(elements,"md.smb.zMin",SmbZMinEnum);
	58	iomodel->FetchDataToInput(elements,"md.smb.Tmean",SmbTmeanEnum);
[23468]	59	iomodel->FetchDataToInput(elements,"md.smb.Vmean",SmbVmeanEnum);
[20690]	60	iomodel->FetchDataToInput(elements,"md.smb.C",SmbCEnum);
	61	iomodel->FetchDataToInput(elements,"md.smb.Tz",SmbTzEnum);
	62	iomodel->FetchDataToInput(elements,"md.smb.Vz",SmbVzEnum);
[21232]	63	InputUpdateFromConstantx(elements,0.,SmbIsInitializedEnum);
	64	iomodel->FetchDataToInput(elements,"md.smb.Dzini",SmbDziniEnum);
	65	iomodel->FetchDataToInput(elements,"md.smb.Dini",SmbDiniEnum);
	66	iomodel->FetchDataToInput(elements,"md.smb.Reini",SmbReiniEnum);
	67	iomodel->FetchDataToInput(elements,"md.smb.Gdnini",SmbGdniniEnum);
	68	iomodel->FetchDataToInput(elements,"md.smb.Gspini",SmbGspiniEnum);
	69	iomodel->FetchDataToInput(elements,"md.smb.ECini",SmbECiniEnum);
	70	iomodel->FetchDataToInput(elements,"md.smb.Wini",SmbWiniEnum);
	71	iomodel->FetchDataToInput(elements,"md.smb.Aini",SmbAiniEnum);
	72	iomodel->FetchDataToInput(elements,"md.smb.Tini",SmbTiniEnum);
	73	iomodel->FetchDataToInput(elements,"md.smb.Sizeini",SmbSizeiniEnum);
[22450]	74	iomodel->FetchDataToInput(elements,"md.smb.aValue",SmbAValueEnum);
	75	iomodel->FetchDataToInput(elements,"md.smb.teValue",SmbTeValueEnum);
[19554]	76	break;
[19528]	77	case SMBpddEnum:
[20690]	78	iomodel->FindConstant(&isdelta18o,"md.smb.isdelta18o");
	79	iomodel->FindConstant(&ismungsm,"md.smb.ismungsm");
	80	iomodel->FetchDataToInput(elements,"md.thermal.spctemperature",ThermalSpctemperatureEnum);
	81	iomodel->FetchDataToInput(elements,"md.smb.s0p",SmbS0pEnum);
	82	iomodel->FetchDataToInput(elements,"md.smb.s0t",SmbS0tEnum);
[19528]	83	if(isdelta18o \|\| ismungsm){
[20690]	84	iomodel->FetchDataToInput(elements,"md.smb.temperatures_lgm",SmbTemperaturesLgmEnum);
	85	iomodel->FetchDataToInput(elements,"md.smb.temperatures_presentday",SmbTemperaturesPresentdayEnum);
	86	iomodel->FetchDataToInput(elements,"md.smb.precipitations_presentday",SmbPrecipitationsPresentdayEnum);
	87	iomodel->FetchDataToInput(elements,"md.smb.precipitations_lgm",SmbPrecipitationsLgmEnum);
[22852]	88	}else{
[20690]	89	iomodel->FetchDataToInput(elements,"md.smb.precipitation",SmbPrecipitationEnum);
	90	iomodel->FetchDataToInput(elements,"md.smb.monthlytemperatures",SmbMonthlytemperaturesEnum);
[19528]	91	}
	92	break;
[23317]	93	case SMBpddSicopolisEnum:
	94	iomodel->FetchDataToInput(elements,"md.smb.s0p",SmbS0pEnum);
	95	iomodel->FetchDataToInput(elements,"md.smb.s0t",SmbS0tEnum);
	96	iomodel->FindConstant(&isfirnwarming,"md.smb.isfirnwarming");
	97	iomodel->FetchDataToInput(elements,"md.smb.smb_corr",SmbSmbCorrEnum);
	98	iomodel->FetchDataToInput(elements,"md.smb.precipitation",SmbPrecipitationEnum);
	99	iomodel->FetchDataToInput(elements,"md.smb.monthlytemperatures",SmbMonthlytemperaturesEnum);
	100	iomodel->FetchDataToInput(elements,"md.smb.precipitation_anomaly",SmbPrecipitationsAnomalyEnum);
	101	iomodel->FetchDataToInput(elements,"md.smb.temperature_anomaly",SmbTemperaturesAnomalyEnum);
	102	break;
[19528]	103	case SMBd18opddEnum:
[22852]	104	iomodel->FindConstant(&istemperaturescaled,"md.smb.istemperaturescaled");
	105	iomodel->FindConstant(&isprecipscaled,"md.smb.isprecipscaled");
[20690]	106	iomodel->FindConstant(&ismungsm,"md.smb.ismungsm");
	107	iomodel->FindConstant(&isd18opd,"md.smb.isd18opd");
[22448]	108	iomodel->FindConstant(&issetpddfac,"md.smb.issetpddfac");
[20690]	109	iomodel->FetchDataToInput(elements,"md.thermal.spctemperature",ThermalSpctemperatureEnum);
	110	iomodel->FetchDataToInput(elements,"md.smb.s0p",SmbS0pEnum);
	111	iomodel->FetchDataToInput(elements,"md.smb.s0t",SmbS0tEnum);
	112	if(isd18opd){
	113	iomodel->FetchDataToInput(elements,"md.smb.temperatures_presentday",SmbTemperaturesPresentdayEnum);
	114	iomodel->FetchDataToInput(elements,"md.smb.precipitations_presentday",SmbPrecipitationsPresentdayEnum);
[22852]	115	if(!istemperaturescaled){
	116	iomodel->FetchDataToInput(elements,"md.smb.temperatures_reconstructed",SmbTemperaturesReconstructedEnum);
	117	}
	118	if(!isprecipscaled){
	119	iomodel->FetchDataToInput(elements,"md.smb.precipitations_reconstructed",SmbPrecipitationsReconstructedEnum);
	120	}
[19528]	121	}
[22448]	122	if(issetpddfac){
	123	iomodel->FetchDataToInput(elements,"md.smb.pddfac_snow",SmbPddfacSnowEnum,-1.);
	124	iomodel->FetchDataToInput(elements,"md.smb.pddfac_ice",SmbPddfacIceEnum,-1.);
	125	}
[19528]	126	break;
	127	case SMBgradientsEnum:
[20690]	128	iomodel->FetchDataToInput(elements,"md.smb.href",SmbHrefEnum);
	129	iomodel->FetchDataToInput(elements,"md.smb.smbref",SmbSmbrefEnum);
	130	iomodel->FetchDataToInput(elements,"md.smb.b_pos",SmbBPosEnum);
	131	iomodel->FetchDataToInput(elements,"md.smb.b_neg",SmbBNegEnum);
[19528]	132	break;
[21469]	133	case SMBgradientselaEnum:
	134	iomodel->FetchDataToInput(elements,"md.smb.ela",SmbElaEnum);
	135	iomodel->FetchDataToInput(elements,"md.smb.b_pos",SmbBPosEnum);
	136	iomodel->FetchDataToInput(elements,"md.smb.b_neg",SmbBNegEnum);
	137	iomodel->FetchDataToInput(elements,"md.smb.b_max",SmbBMaxEnum);
	138	iomodel->FetchDataToInput(elements,"md.smb.b_min",SmbBMinEnum);
	139	break;
[19528]	140	case SMBhenningEnum:
[20690]	141	iomodel->FetchDataToInput(elements,"md.smb.smbref",SmbSmbrefEnum,0.);
[19528]	142	break;
	143	case SMBcomponentsEnum:
[20690]	144	iomodel->FetchDataToInput(elements,"md.smb.accumulation",SmbAccumulationEnum,0.);
	145	iomodel->FetchDataToInput(elements,"md.smb.evaporation",SmbEvaporationEnum,0.);
	146	iomodel->FetchDataToInput(elements,"md.smb.runoff",SmbRunoffEnum,0.);
[19528]	147	break;
	148	case SMBmeltcomponentsEnum:
[20690]	149	iomodel->FetchDataToInput(elements,"md.smb.accumulation",SmbAccumulationEnum,0.);
	150	iomodel->FetchDataToInput(elements,"md.smb.evaporation",SmbEvaporationEnum,0.);
	151	iomodel->FetchDataToInput(elements,"md.smb.melt",SmbMeltEnum,0.);
	152	iomodel->FetchDataToInput(elements,"md.smb.refreeze",SmbRefreezeEnum,0.);
[19528]	153	break;
[23366]	154	case SMBgradientscomponentsEnum:
	155	iomodel->FetchDataToInput(elements,"md.smb.accualti",SmbAccualtiEnum);
	156	iomodel->FetchDataToInput(elements,"md.smb.accugrad",SmbAccugradEnum);
	157	iomodel->FetchDataToInput(elements,"md.smb.runoffalti",SmbRunoffaltiEnum);
	158	iomodel->FetchDataToInput(elements,"md.smb.runoffgrad",SmbRunoffgradEnum);
	159	break;
[23540]	160	case SMBsemicEnum:
	161	iomodel->FetchDataToInput(elements,"md.thermal.spctemperature",ThermalSpctemperatureEnum);
	162	iomodel->FetchDataToInput(elements,"md.smb.s0gcm",SmbS0gcmEnum);
	163	iomodel->FetchDataToInput(elements,"md.smb.dailysnowfall",SmbDailysnowfallEnum);
	164	iomodel->FetchDataToInput(elements,"md.smb.dailyrainfall",SmbDailyrainfallEnum);
	165	iomodel->FetchDataToInput(elements,"md.smb.dailydsradiation",SmbDailydsradiationEnum);
	166	iomodel->FetchDataToInput(elements,"md.smb.dailydlradiation",SmbDailydlradiationEnum);
	167	iomodel->FetchDataToInput(elements,"md.smb.dailywindspeed",SmbDailywindspeedEnum);
	168	iomodel->FetchDataToInput(elements,"md.smb.dailypressure",SmbDailypressureEnum);
	169	iomodel->FetchDataToInput(elements,"md.smb.dailyairdensity",SmbDailyairdensityEnum);
	170	iomodel->FetchDataToInput(elements,"md.smb.dailyairhumidity",SmbDailyairhumidityEnum);
	171	iomodel->FetchDataToInput(elements,"md.smb.dailytemperature",SmbDailytemperatureEnum);
	172	break;
[19528]	173	default:
	174	_error_("Surface mass balance model "<<EnumToStringx(smb_model)<<" not supported yet");
	175	}
	176
	177	}/}}}/
	178	void SmbAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/{{{/
	179
	180	int numoutputs;
	181	char** requestedoutputs = NULL;
[23317]	182	bool isdelta18o,ismungsm,isd18opd,issetpddfac,interp,isfirnwarming;
[24240]	183	int smb_model, smbslices;
[19528]	184	IssmDouble *temp = NULL;
	185	int N,M;
[23066]	186
[20690]	187	parameters->AddObject(iomodel->CopyConstantObject("md.smb.model",SmbEnum));
[23066]	188
[20690]	189	iomodel->FindConstant(&smb_model,"md.smb.model");
	190	iomodel->FindConstant(&interp,"md.timestepping.interp_forcings");
[23066]	191
[24240]	192	iomodel->FetchData(&smbslices,"md.smb.steps_per_step");
	193	parameters->AddObject(new IntParam(SmbStepsPerStepEnum,smbslices));
	194
[19528]	195	switch(smb_model){
	196	case SMBforcingEnum:
[23814]	197	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isclimatology",SmbIsclimatologyEnum));
[19528]	198	break;
[19554]	199	case SMBgembEnum:
[20690]	200	parameters->AddObject(iomodel->CopyConstantObject("md.smb.aIdx",SmbAIdxEnum));
	201	parameters->AddObject(iomodel->CopyConstantObject("md.smb.swIdx",SmbSwIdxEnum));
	202	parameters->AddObject(iomodel->CopyConstantObject("md.smb.denIdx",SmbDenIdxEnum));
[23468]	203	parameters->AddObject(iomodel->CopyConstantObject("md.smb.dsnowIdx",SmbDsnowIdxEnum));
[20690]	204	parameters->AddObject(iomodel->CopyConstantObject("md.smb.cldFrac",SmbCldFracEnum));
	205	parameters->AddObject(iomodel->CopyConstantObject("md.smb.t0wet",SmbT0wetEnum));
	206	parameters->AddObject(iomodel->CopyConstantObject("md.smb.t0dry",SmbT0dryEnum));
	207	parameters->AddObject(iomodel->CopyConstantObject("md.smb.K",SmbKEnum));
	208	parameters->AddObject(iomodel->CopyConstantObject("md.smb.aSnow",SmbASnowEnum));
	209	parameters->AddObject(iomodel->CopyConstantObject("md.smb.aIce",SmbAIceEnum));
	210	parameters->AddObject(iomodel->CopyConstantObject("md.smb.dt",SmbDtEnum));
	211	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isgraingrowth",SmbIsgraingrowthEnum));
	212	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isalbedo",SmbIsalbedoEnum));
	213	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isshortwave",SmbIsshortwaveEnum));
	214	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isthermal",SmbIsthermalEnum));
	215	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isaccumulation",SmbIsaccumulationEnum));
	216	parameters->AddObject(iomodel->CopyConstantObject("md.smb.ismelt",SmbIsmeltEnum));
	217	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isdensification",SmbIsdensificationEnum));
	218	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isturbulentflux",SmbIsturbulentfluxEnum));
[23808]	219	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isclimatology",SmbIsclimatologyEnum));
[20690]	220	parameters->AddObject(iomodel->CopyConstantObject("md.smb.InitDensityScaling",SmbInitDensityScalingEnum));
[22475]	221	parameters->AddObject(iomodel->CopyConstantObject("md.smb.ThermoDeltaTScaling",SmbThermoDeltaTScalingEnum));
[22482]	222	parameters->AddObject(iomodel->CopyConstantObject("md.smb.adThresh",SmbAdThreshEnum));
[19554]	223	break;
[19528]	224	case SMBpddEnum:
[20690]	225	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isdelta18o",SmbIsdelta18oEnum));
	226	parameters->AddObject(iomodel->CopyConstantObject("md.smb.ismungsm",SmbIsmungsmEnum));
[22448]	227	parameters->AddObject(iomodel->CopyConstantObject("md.smb.issetpddfac",SmbIssetpddfacEnum));
[20690]	228	iomodel->FindConstant(&isdelta18o,"md.smb.isdelta18o");
	229	iomodel->FindConstant(&ismungsm,"md.smb.ismungsm");
[19528]	230
	231	if(ismungsm){
[20690]	232	iomodel->FetchData(&temp,&N,&M,"md.smb.Pfac"); _assert_(N==2);
[19528]	233	parameters->AddObject(new TransientParam(SmbPfacEnum,&temp[0],&temp[M],interp,M));
[20690]	234	iomodel->DeleteData(temp,"md.smb.Pfac");
[23066]	235
[20690]	236	iomodel->FetchData(&temp,&N,&M,"md.smb.Tdiff"); _assert_(N==2);
[19528]	237	parameters->AddObject(new TransientParam(SmbTdiffEnum,&temp[0],&temp[M],interp,M));
[20690]	238	iomodel->DeleteData(temp,"md.smb.Tdiff");
[19528]	239
[20690]	240	iomodel->FetchData(&temp,&N,&M,"md.smb.sealev"); _assert_(N==2);
[19528]	241	parameters->AddObject(new TransientParam(SmbSealevEnum,&temp[0],&temp[M],interp,M));
[20690]	242	iomodel->DeleteData(temp,"md.smb.sealev");
[19528]	243	}
	244	if(isdelta18o){
[20690]	245	iomodel->FetchData(&temp,&N,&M,"md.smb.delta18o"); _assert_(N==2);
[19528]	246	parameters->AddObject(new TransientParam(SmbDelta18oEnum,&temp[0],&temp[M],interp,M));
[20690]	247	iomodel->DeleteData(temp,"md.smb.delta18o");
[19528]	248
[20690]	249	iomodel->FetchData(&temp,&N,&M,"md.smb.delta18o_surface"); _assert_(N==2);
[19528]	250	parameters->AddObject(new TransientParam(SmbDelta18oSurfaceEnum,&temp[0],&temp[M],interp,M));
[20690]	251	iomodel->DeleteData(temp,"md.smb.delta18o_surface");
[19528]	252	}
	253	break;
[23317]	254	case SMBpddSicopolisEnum:
[23328]	255	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isfirnwarming",SmbIsfirnwarmingEnum));
[23317]	256	break;
[19528]	257	case SMBd18opddEnum:
[20690]	258	parameters->AddObject(iomodel->CopyConstantObject("md.smb.ismungsm",SmbIsmungsmEnum));
	259	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isd18opd",SmbIsd18opdEnum));
[22448]	260	parameters->AddObject(iomodel->CopyConstantObject("md.smb.issetpddfac",SmbIssetpddfacEnum));
[20690]	261	iomodel->FindConstant(&ismungsm,"md.smb.ismungsm");
	262	iomodel->FindConstant(&isd18opd,"md.smb.isd18opd");
[22448]	263	iomodel->FindConstant(&issetpddfac,"md.smb.issetpddfac");
[19528]	264	if(isd18opd){
[22608]	265	parameters->AddObject(iomodel->CopyConstantObject("md.smb.f",SmbFEnum));
[22495]	266	parameters->AddObject(iomodel->CopyConstantObject("md.smb.istemperaturescaled",SmbIstemperaturescaledEnum));
[22852]	267	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isprecipscaled",SmbIsprecipscaledEnum));
[20690]	268	iomodel->FetchData(&temp,&N,&M,"md.smb.delta18o"); _assert_(N==2);
[19528]	269	parameters->AddObject(new TransientParam(SmbDelta18oEnum,&temp[0],&temp[M],interp,M));
[20690]	270	iomodel->DeleteData(temp,"md.smb.delta18o");
[19528]	271	}
	272	break;
	273	case SMBgradientsEnum:
	274	/Nothing to add to parameters/
	275	break;
[21469]	276	case SMBgradientselaEnum:
	277	/Nothing to add to parameters/
	278	break;
[19528]	279	case SMBhenningEnum:
	280	/Nothing to add to parameters/
	281	break;
	282	case SMBcomponentsEnum:
[23814]	283	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isclimatology",SmbIsclimatologyEnum));
[19528]	284	break;
	285	case SMBmeltcomponentsEnum:
[23814]	286	parameters->AddObject(iomodel->CopyConstantObject("md.smb.isclimatology",SmbIsclimatologyEnum));
[19528]	287	break;
[23366]	288	case SMBgradientscomponentsEnum:
	289	parameters->AddObject(iomodel->CopyConstantObject("md.smb.accualti",SmbAccualtiEnum));
	290	parameters->AddObject(iomodel->CopyConstantObject("md.smb.accugrad",SmbAccugradEnum));
	291	parameters->AddObject(iomodel->CopyConstantObject("md.smb.runoffalti",SmbRunoffaltiEnum));
	292	parameters->AddObject(iomodel->CopyConstantObject("md.smb.runoffgrad",SmbRunoffgradEnum));
	293
	294	iomodel->FetchData(&temp,&N,&M,"md.smb.accuref"); _assert_(N==2);
	295	parameters->AddObject(new TransientParam(SmbAccurefEnum,&temp[0],&temp[M],interp,M));
	296	iomodel->DeleteData(temp,"md.smb.accuref");
	297
	298	iomodel->FetchData(&temp,&N,&M,"md.smb.runoffref"); _assert_(N==2);
	299	parameters->AddObject(new TransientParam(SmbRunoffrefEnum,&temp[0],&temp[M],interp,M));
	300	iomodel->DeleteData(temp,"md.smb.runoffref");
	301	break;
[23540]	302	case SMBsemicEnum:
	303	/Nothing to add to parameters/
	304	break;
[19528]	305	default:
	306	_error_("Surface mass balance model "<<EnumToStringx(smb_model)<<" not supported yet");
	307	}
	308
[20690]	309	iomodel->FindConstant(&requestedoutputs,&numoutputs,"md.smb.requested_outputs");
[19528]	310	parameters->AddObject(new IntParam(SmbNumRequestedOutputsEnum,numoutputs));
	311	if(numoutputs)parameters->AddObject(new StringArrayParam(SmbRequestedOutputsEnum,requestedoutputs,numoutputs));
[20690]	312	iomodel->DeleteData(&requestedoutputs,numoutputs,"md.smb.requested_outputs");
[19528]	313
	314	}/}}}/
	315
	316	/Finite Element Analysis/
	317	void SmbAnalysis::Core(FemModel* femmodel){/{{{/
	318
	319	int smb_model;
	320
	321	/Figure out smb model: /
	322	femmodel->parameters->FindParam(&smb_model,SmbEnum);
[23066]	323
[19528]	324	/branch to correct module/
	325	switch(smb_model){
	326	case SMBforcingEnum:
[23814]	327	SmbForcingx(femmodel);
[19528]	328	break;
[19554]	329	case SMBgembEnum:
	330	Gembx(femmodel);
	331	break;
[19528]	332	case SMBpddEnum:
	333	bool isdelta18o,ismungsm;
	334	femmodel->parameters->FindParam(&isdelta18o,SmbIsdelta18oEnum);
	335	femmodel->parameters->FindParam(&ismungsm,SmbIsmungsmEnum);
	336	if(isdelta18o){
	337	if(VerboseSolution()) _printf0_(" call Delta18oParameterization module\n");
	338	Delta18oParameterizationx(femmodel);
[23366]	339	}
[19528]	340	if(ismungsm){
	341	if(VerboseSolution()) _printf0_(" call MungsmtpParameterization module\n");
	342	MungsmtpParameterizationx(femmodel);
[23366]	343	}
[19528]	344	if(VerboseSolution()) _printf0_(" call positive degree day module\n");
	345	PositiveDegreeDayx(femmodel);
	346	break;
[23317]	347	case SMBpddSicopolisEnum:
	348	if(VerboseSolution()) _printf0_(" call SICOPOLIS positive degree day module\n");
	349	PositiveDegreeDaySicopolisx(femmodel);
	350	break;
[19528]	351	case SMBd18opddEnum:
	352	bool isd18opd;
	353	femmodel->parameters->FindParam(&isd18opd,SmbIsd18opdEnum);
	354	if(isd18opd){
	355	if(VerboseSolution()) _printf0_(" call Delta18opdParameterization module\n");
	356	Delta18opdParameterizationx(femmodel);
	357	if(VerboseSolution()) _printf0_(" call positive degree day module\n");
	358	PositiveDegreeDayx(femmodel);
[23317]	359	}
[19528]	360	break;
	361	case SMBgradientsEnum:
	362	if(VerboseSolution())_printf0_(" call smb gradients module\n");
	363	SmbGradientsx(femmodel);
	364	break;
[21469]	365	case SMBgradientselaEnum:
	366	if(VerboseSolution())_printf0_(" call smb gradients ela module\n");
	367	SmbGradientsElax(femmodel);
	368	break;
[19528]	369	case SMBhenningEnum:
	370	if(VerboseSolution())_printf0_(" call smb Henning module\n");
	371	SmbHenningx(femmodel);
	372	break;
	373	case SMBcomponentsEnum:
	374	if(VerboseSolution())_printf0_(" call smb Components module\n");
	375	SmbComponentsx(femmodel);
	376	break;
	377	case SMBmeltcomponentsEnum:
	378	if(VerboseSolution())_printf0_(" call smb Melt Components module\n");
	379	SmbMeltComponentsx(femmodel);
	380	break;
	381	case SMBgcmEnum:
	382	/Nothing to be done/
	383	break;
[23366]	384	case SMBgradientscomponentsEnum:
	385	if(VerboseSolution())_printf0_(" call smb gradients components module\n");
	386	SmbGradientsComponentsx(femmodel);
	387	break;
[23540]	388	case SMBsemicEnum:
	389	#ifdef _HAVE_SEMIC_
	390	if(VerboseSolution())_printf0_(" call smb SEMIC module\n");
	391	SmbSemicx(femmodel);
	392	#else
	393	_error_("SEMIC not installed");
	394	#endif //_HAVE_SEMIC_
	395	break;
[19528]	396	default:
	397	_error_("Surface mass balance model "<<EnumToStringx(smb_model)<<" not supported yet");
	398	}
	399
	400	}/}}}/
	401	ElementVector* SmbAnalysis::CreateDVector(Element* element){/{{{/
	402	_error_("not implemented");
	403	}/}}}/
	404	ElementMatrix* SmbAnalysis::CreateJacobianMatrix(Element* element){/{{{/
	405	_error_("Not implemented");
	406	}/}}}/
	407	ElementMatrix* SmbAnalysis::CreateKMatrix(Element* element){/{{{/
	408	_error_("not implemented yet");
	409	}/}}}/
	410	ElementVector* SmbAnalysis::CreatePVector(Element* element){/{{{/
	411	_error_("not implemented yet");
	412	}/}}}/
	413	void SmbAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/{{{/
	414	_error_("not implemented yet");
	415	}/}}}/
	416	void SmbAnalysis::GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index){/{{{/
	417	_error_("Not implemented yet");
	418	}/}}}/
	419	void SmbAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/{{{/
	420	_error_("not implemented yet");
	421	}/}}}/
	422	void SmbAnalysis::UpdateConstraints(FemModel* femmodel){/{{{/
	423	/Default, do nothing/
	424	return;
	425	}/}}}/

Note: See TracBrowser for help on using the repository browser.

Download in other formats: