Changeset 27313

Timestamp:

10/17/22 07:25:14 (2 years ago)

Author:

jwilner

Message:

NEW: set calving rate to 0 where there is no ice

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (4 diffs)

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

@@ -569 +569 @@
         IssmDouble ds, db, da, dt, dw, r, R;
 
+        if(!IsIceInElement()){
+                for (int iv=0;iv<NUMVERTICES;iv++) calvingrate[iv]=0.;
+                this->AddInput(CalvingCalvingrateEnum,&calvingrate[0],P1DGEnum);
+                this->CalvingRateToVector();
+                return;
+        }
+
         /*Retrieve all inputs and parameters we will need*/
         IssmDouble rc        = FindParam(CalvingRcEnum);
@@ -618 +625 @@
 
          /*3. "Additional" crevasse opening*/
-         vel = sqrt(vx*vx + vy*vy)/(365.25*24*3600);
-         da = H* max(0., log(vel/1600.))/log(1.2);
+                        vel = sqrt(vx*vx + vy*vy);
+                        da = H* max(0., log(vel*365*24*3600/1600.))/log(1.2);
 
          /*4. deal with shallow ice*/
@@ -627 +634 @@
          dw = 0.;
          R = smb*365.25*24*3600; //convert from m/s to m/yr
-         if(R>1.5 && R<=3.){
+                        if(R>1.5 && R<=3.){
             dw = 4*1.5*(R - 1.5);
          }
@@ -636 +643 @@
          /*Total calving rate*/
          r = (ds+db+da+dt+dw)/H;
+                        //if(this->Id()==1){
+                        //      printf("rc = %g\n",rc);
+                        //      printf("ds = %g\n",ds);
+                        //}
                         calvingrate[iv]= mig_max * max(0., min(1., (r - rc)/(1 - rc))); //P&DC: mig_max = 3000 m/yr
                         _assert_(!xIsNan<IssmDouble>(calvingrate[iv]));