[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 | }/*}}}*/
|
---|