Changeset 21741

Timestamp:

05/23/17 10:12:54 (8 years ago)

Author:

Eric.Larour

Message:

CHG: ensure contiguous had to be enabled for the G_elastic and co arrays because they are quite big.

File:

• issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp

@@ -90 +90 @@
                 iomodel->FetchData(&degacc,"md.slr.degacc");
                 M=reCast<int,IssmDouble>(180./degacc+1.);
-                G_elastic=xNew<IssmDouble>(M);
-                U_elastic=xNew<IssmDouble>(M);
-                H_elastic=xNew<IssmDouble>(M);
+                G_elastic=xNew<IssmDouble>(M,"t");
+                U_elastic=xNew<IssmDouble>(M,"t");
+                H_elastic=xNew<IssmDouble>(M,"t");
                 
                 /*compute combined legendre + love number (elastic green function:*/
                 m=DetermineLocalSize(M,IssmComm::GetComm());
                 GetOwnershipBoundariesFromRange(&lower_row,&upper_row,m,IssmComm::GetComm());
-                G_elastic_local=xNew<IssmDouble>(m);
-                U_elastic_local=xNew<IssmDouble>(m);
-                H_elastic_local=xNew<IssmDouble>(m);
+                G_elastic_local=xNew<IssmDouble>(m,"t");
+                U_elastic_local=xNew<IssmDouble>(m,"t");
+                H_elastic_local=xNew<IssmDouble>(m,"t");
 
                 for(int i=lower_row;i<upper_row;i++){