Changeset 16905 for issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

Timestamp:

11/23/13 13:44:24 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: in 2d zmin is actually ymin (same for max)

File:

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

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

@@ -2586 +2586 @@
 
         /*Initialize gauss points*/
-        IssmDouble zmax=xyz_list[0*3+2]; for(int i=1;i<numvertices;i++) if(xyz_list[i*3+2]>zmax) zmax=xyz_list[i*3+2];
-        IssmDouble zmin=xyz_list[0*3+2]; for(int i=1;i<numvertices;i++) if(xyz_list[i*3+2]<zmin) zmin=xyz_list[i*3+2];
-        if(zmax>0 && zmin<0) gauss=element->NewGauss(xyz_list,xyz_list_front,3,30);//refined in vertical because of the sea level discontinuity
-        else                 gauss=element->NewGauss(xyz_list,xyz_list_front,3,3);
+        IssmDouble zmax=xyz_list[0*3+(dim-1)]; for(int i=1;i<numvertices;i++) if(xyz_list[i*3+(dim-1)]>zmax) zmax=xyz_list[i*3+(dim-1)];
+        IssmDouble zmin=xyz_list[0*3+(dim-1)]; for(int i=1;i<numvertices;i++) if(xyz_list[i*3+(dim-1)]<zmin) zmin=xyz_list[i*3+(dim-1)];
+        if(zmax>0. && zmin<0.) gauss=element->NewGauss(xyz_list,xyz_list_front,3,30);//refined in vertical because of the sea level discontinuity
+        else                   gauss=element->NewGauss(xyz_list,xyz_list_front,3,3);
 
         /* Start  looping on the number of gaussian points: */