Changeset 21224 for issm/trunk-jpl/src/c/classes/Elements/Element.cpp
- Timestamp:
- 09/22/16 16:36:34 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/classes/Elements/Element.cpp ¶
r21218 r21224 2246 2246 2247 2247 /*Intermediary variables: {{{*/ 2248 IssmDoubleisinitialized;2249 IssmDoubleisrestart;2248 bool isinitialized; 2249 bool isrestart; 2250 2250 IssmDouble zTop,dzTop,zMax,zMin,zY,dzMin; 2251 2251 IssmDouble Tmean; … … 2331 2331 parameters->FindParam(&isturbulentflux,SmbIsturbulentfluxEnum); 2332 2332 parameters->FindParam(&init_scaling,SmbInitDensityScalingEnum); 2333 2333 2334 /*}}}*/ 2334 2335 /*Retrieve inputs: {{{*/ … … 2366 2367 Tz_input->GetInputValue(&Tz,gauss); 2367 2368 Vz_input->GetInputValue(&Vz,gauss); 2368 isinitialized_input->GetInputValue(&isinitialized ,gauss);2369 isrestart_input->GetInputValue(&isrestart ,gauss);2369 isinitialized_input->GetInputValue(&isinitialized); 2370 isrestart_input->GetInputValue(&isrestart); 2370 2371 /*}}}*/ 2371 2372 … … 2391 2392 2392 2393 /*Flag the initialization:*/ 2393 this->AddInput(new DoubleInput(SmbIsInitializedEnum,false));2394 this->AddInput(new BoolInput(SmbIsInitializedEnum,false)); 2394 2395 } 2395 2396 else if(isrestart){ //Retrieve the snow properties from previous run … … 2435 2436 2436 2437 /*Flag the initialization:*/ 2437 this->AddInput(new DoubleInput(SmbIsrestartEnum,false));//*CL* add2438 this->AddInput(new BoolInput(SmbIsrestartEnum,false)); 2438 2439 } 2439 2440 else{ 2440 2441 2441 /*Recover inputs: */ 2442 2442 DoubleArrayInput* dz_input= dynamic_cast<DoubleArrayInput*>(this->GetInput(SmbDzEnum)); _assert_(dz_input); … … 2449 2449 DoubleArrayInput* a_input= dynamic_cast<DoubleArrayInput*>(this->GetInput(SmbAEnum));_assert_(a_input); 2450 2450 DoubleArrayInput* T_input= dynamic_cast<DoubleArrayInput*>(this->GetInput(SmbTEnum));_assert_(T_input); 2451 2452 2451 2452 /*Recover arrays: */ 2453 2453 dz_input->GetValues(&dz,&m); 2454 2454 d_input->GetValues(&d,&m); 2455 2455 re_input->GetValues(&re,&m); 2456 2456 gdn_input->GetValues(&gdn,&m); … … 2460 2460 a_input->GetValues(&a,&m); 2461 2461 T_input->GetValues(&T,&m); 2462 2463 2462 2463 //fixed lower temperatuer bounday condition - T is fixed 2464 2464 T_bottom=T[m-1]; 2465 2465 2466 2466 } /*}}}*/ 2467 2467 2468 2468 // determine initial mass [kg] … … 2578 2578 count++; 2579 2579 } //for (t=time;t<time+dt;t=t+smb_dt) 2580 2581 2580 2581 /*Save generated inputs: */ 2582 2582 this->AddInput(new DoubleArrayInput(SmbDzEnum,dz,m)); 2583 2583 this->AddInput(new DoubleArrayInput(SmbDEnum,d,m)); … … 2595 2595 this->AddInput(new DoubleInput(SmbM_addEnum,sumMassAdd/yts)); 2596 2596 2597 2597 /*Free allocations:{{{*/ 2598 2598 xDelete<IssmDouble>(dz); 2599 2599 xDelete<IssmDouble>(d);
Note:
See TracChangeset
for help on using the changeset viewer.