Changeset 3977

Timestamp:

05/27/10 15:45:16 (15 years ago)

Author:

seroussi

Message:

inputs melting and accumulation must be divided by md.yts

Location:

issm/trunk/src/c/objects/Elements

Files:

• Penta.cpp (modified) (1 diff)
• Tria.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -103 +103 @@
         }
         if (iomodel->melting_rate) {
-                for(i=0;i<6;i++)nodeinputs[i]=iomodel->melting_rate[penta_node_ids[i]-1];
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->melting_rate[penta_node_ids[i]-1]/iomodel->yts;
                 this->inputs->AddInput(new PentaVertexInput(MeltingRateEnum,nodeinputs));
         }
         if (iomodel->accumulation_rate) {
-                for(i=0;i<6;i++)nodeinputs[i]=iomodel->accumulation_rate[penta_node_ids[i]-1];
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->accumulation_rate[penta_node_ids[i]-1]/iomodel->yts;
                 this->inputs->AddInput(new PentaVertexInput(AccumulationRateEnum,nodeinputs));
         }

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -98 +98 @@
         }
         if (iomodel->melting_rate) {
-                for(i=0;i<3;i++)nodeinputs[i]=iomodel->melting_rate[tria_node_ids[i]-1];
+                for(i=0;i<3;i++)nodeinputs[i]=iomodel->melting_rate[tria_node_ids[i]-1]/iomodel->yts;
                 this->inputs->AddInput(new TriaVertexInput(MeltingRateEnum,nodeinputs));
         }
         if (iomodel->accumulation_rate) {
-                for(i=0;i<3;i++)nodeinputs[i]=iomodel->accumulation_rate[tria_node_ids[i]-1];
+                for(i=0;i<3;i++)nodeinputs[i]=iomodel->accumulation_rate[tria_node_ids[i]-1]/iomodel->yts;
                 this->inputs->AddInput(new TriaVertexInput(AccumulationRateEnum,nodeinputs));
         }