Changeset 27873

Timestamp:

08/08/23 08:24:08 (20 months ago)

Author:

vverjans

Message:

BUG: fixing identification of basin for SMB lapse rates

Location:

issm/trunk-jpl

Files:

• src/c/classes/Elements/Element.cpp (modified) (1 diff)
• test/Archives/Archive257.arch (modified) ( previous)

issm/trunk-jpl/src/c/classes/Elements/Element.cpp ¶

@@ -2469 +2469 @@
            for(int v=0;v<numvertices;v++){
               /*Find elevation bin of Reference elevation and of Vertex*/
+                        bb1 = 0;
+                        bb2 = 0;
                         for(int ii=0;ii<(numelevbins-1);ii++){
-                                if(surfacelist[v]<=elevbins_b[ii]) bb1 = ii;    
-                                if(refelevation_b<=elevbins_b[ii]) bb2 = ii;
+                                if(surfacelist[v]>=elevbins_b[ii]) bb1 = ii+1;
+                                if(refelevation_b>=elevbins_b[ii]) bb2 = ii+1;
                         }
-                        /*Check for elevations above highest bin limit */
-                        if(surfacelist[v]>elevbins_b[numelevbins-1-1]) bb1 = numelevbins-1;
-                        if(refelevation_b>elevbins_b[numelevbins-1-1]) bb2 = numelevbins-1;
+
                         /*Vertex and Reference elevation in same elevation bin*/
                         if(bb1==bb2){