Changeset 27820

Timestamp:

07/03/23 12:54:20 (21 months ago)

Author:

Cheng Gong

Message:

BUG: convert time using yts in CalvingMIP EXP4

File:

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

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

@@ -1128 +1128 @@
         int         dim, domaintype;
         IssmDouble      vx, vy, vel, c, wrate;
-        IssmDouble  time, groundedice;
+        IssmDouble  time, groundedice, yts;
 
         /*Get problem dimension and whether there is moving front or not*/
         this->FindParam(&domaintype,DomainTypeEnum);
         this->FindParam(&time,TimeEnum);
+        this->FindParam(&yts,ConstantsYtsEnum);
 
         switch(domaintype){
@@ -1176 +1177 @@
                         case 4:
                                 /* 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)) {
+                                if((groundedice<0) && (time<=500.0*yts)) {
                                         wrate_input->GetInputValue(&wrate,&gauss);
                                 }