source: issm/oecreview/Archive/16554-17801/ISSM-17706-17707.diff

../trunk-jpl/src/c/modules/GroundinglineMigrationx/GroundinglineMigrationx.cpp

@@ -94 +94 @@
                         phi[i]=-1.;
                 }
                 else{
-                        _error_("not supported (vertexfloating="<<serial_vertexfloating[i]<<" vertexgrounded="<<serial_vertexgrounded[i]<<")");
+                        phi[i]=9999.;
                 }
         }
 

../trunk-jpl/src/c/cores/transient_core.cpp

@@ -166 +166 @@
                 if(isgroundingline){
                         if(VerboseSolution()) _printf0_("   computing new grounding line position\n");
                         GroundinglineMigrationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
+
+                        if(groundingline_migration==ContactEnum){
+                                femmodel->parameters->SetParam(MaskGroundediceLevelsetEnum,InputToExtrudeEnum);
+                                extrudefrombase_core(femmodel);
+                        }
                         if(save_results){
                                 int outputs[3] = {SurfaceEnum,BaseEnum,MaskGroundediceLevelsetEnum};
                                 femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3);

../trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -4367 +4367 @@
         density     = rho_ice/rho_water;
         oldfloating = this->IsFloating(); //figure out if element is floating before we start to change everything
 
-        if(migration_style==ContactEnum){
-                for(i=0;i<NUMVERTICES;i++) phi[i]=phi_ungrounding[vertices[i]->Pid()];
+        if(migration_style == ContactEnum){
+                for(i = 0;i < NUMVERTICES;i++) phi[i] = phi_ungrounding[vertices[i]->Pid()];
                 this->inputs->AddInput(new TriaInput(MaskGroundediceLevelsetEnum,&phi[0],P1Enum));
+
+                /*go through vertices, and update inputs, considering them to be TriaVertex type: */
+                for(i = 0;i < NUMVERTICES;i++){
+                        /*Ice shelf: if bed below bathymetry, impose it at the bathymetry and update surface, elso do nothing */
+                        if(phi[i] >= 0.){
+                                if(b[i] < r[i]){
+                                        b[i]  = r[i];
+                                }
+                        }
+                }
+
+                /*Add basal melting rate if element just ungrounded*/
+                for(i=0;i<NUMVERTICES;i++){
+                        if(phi[i]<0.){
+                                floatingelement=true;
+                                break;
+                        }
+                }
+                if(!oldfloating && floatingelement==true){
+                        for(i=0;i<NUMVERTICES;i++) melting[i]=gl_melting_rate/yts;
+                        this->inputs->AddInput(new TriaInput(BasalforcingsMeltingRateEnum,&melting[0],P1Enum));
+                } 
+
+                /*Update inputs*/
+                this->inputs->AddInput(new TriaInput(BaseEnum,&b[0],P1Enum));
                 return;
         }
 
@@ -4437 +4462 @@
                 }
         }
 
-   /*Add basal melting rate if element just ungrounded*/
+        /*Add basal melting rate if element just ungrounded*/
         if(!oldfloating && floatingelement==true){
                 for(i=0;i<NUMVERTICES;i++)melting[i]=gl_melting_rate/yts;
                 this->inputs->AddInput(new TriaInput(BasalforcingsMeltingRateEnum,&melting[0],P1Enum));