Changeset 27828

Timestamp:

07/14/23 11:28:42 (20 months ago)

Author:

Cheng Gong

Message:

CHG: set the calving rate for CalvingMIP EXP4 as a function of time

File:

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

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

@@ -1178 +1178 @@
                                 /* Exp 4: set c=v-wrate(given), for the first 500 years, then c=0 for the second 500 years*/
                                 if((groundedice<0) && (time<=500.0*yts)) {
-                                        wrate_input->GetInputValue(&wrate,&gauss);
+                                //      wrate_input->GetInputValue(&wrate,&gauss);
+                                        wrate = -750*sin(2.0*M_PI*time/yts/1000)/yts;  // m/a -> m/s
                                 }
                                 else {