Changeset 27797

Timestamp:

06/21/23 14:23:42 (22 months ago)

Author:

Cheng Gong

Message:

ADD: CalvingMIP Exp 2 and 4, second half with no calving, use md.calvingcalvingmip.experiment=0

Location:

issm/trunk-jpl/src

Files:

• c/classes/Elements/Tria.cpp (modified) (1 diff)
• m/classes/calvingcalvingmip.m (modified) (1 diff)

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

@@ -1170 +1170 @@
                                 /*Get wrate*/
                                 wrate_input->GetInputValue(&wrate,gauss);
+                                break;
+                        case 0:
+                                /* This is for second half of Exp 2 and 4: to set c=0*/
+                                wrate = vel;
                                 break;
                         default:

issm/trunk-jpl/src/m/classes/calvingcalvingmip.m

@@ -41 +41 @@
 
                         md = checkfield(md,'fieldname','calving.min_thickness','>=',0,'NaN',1,'Inf',1,'numel',1);
-                        md = checkfield(md,'fieldname','calving.experiment','values',[1, 2, 3, 4, 5]);
+                        md = checkfield(md,'fieldname','calving.experiment','values',[0, 1, 2, 3, 4, 5]);
                 end % }}}
                 function disp(self) % {{{