Changeset 7353

Timestamp:

02/05/11 08:53:26 (14 years ago)

Author:

Eric.Larour

Message:

Quick hack

File:

• issm/trunk/src/c/objects/Elements/Tria.cpp (modified) (3 diffs)

issm/trunk/src/c/objects/Elements/Tria.cpp ¶

@@ -4823 +4823 @@
         bool    elementonshelf = false;
 
+        double  shelf_melting=10;
+        double  yts;
+        int     swap;
+
+        /*recover parameters: */
+        parameters->FindParam(&yts,YtsEnum);
+
         /*Recover info at the vertices: */
         Input* surface_input =inputs->GetInput(SurfaceEnum); _assert_(surface_input);
@@ -4853 +4860 @@
 
         /*Now, update  shelf status of element. An element can only be on shelf if all its nodes are on shelf: */
+        swap=0;
         elementonshelf=false;
         for(i=0;i<3;i++){
@@ -4860 +4868 @@
                 }
         }
+        if(!this->IsOnShelf() && elementonshelf==true)swap=1;
     this->inputs->AddInput(new BoolInput(ElementOnIceShelfEnum,elementonshelf));
+
+        /*If this element just  became ungrounded, set its melting rate at 50 m/yr:*/
+        if(swap){
+                Input* melting_rate_input     =inputs->GetInput(MeltingRateEnum);     _assert_(melting_rate_input);
+                melting_rate_input->GetValuesPtr(&values,NULL); for(i=0;i<3;i++)values[i]=-shelf_melting/yts;
+        }
+
+
 }
 /*}}}*/