Index: LMDZ6/trunk/libf/phylmd/orografi_strato.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/orografi_strato.F90	(revision 3586)
+++ LMDZ6/trunk/libf/phylmd/orografi_strato.F90	(revision 3587)
@@ -1868,5 +1868,5 @@
   IF (dist_min/=dist_min_glo) cell=-1
 !ym in future find the point at equator close to (0,0) coordinates.
-  !PRINT *, 'SUGWD distmin dist_min_glo cell=', dist_min,dist_min_glo,cell
+  PRINT *, 'SUGWD distmin dist_min_glo cell=', dist_min,dist_min_glo,cell
 
   nktopg_tmp=nktopg
@@ -1875,6 +1875,9 @@
   IF (cell/=-1) THEN
 
+    !print*,'SUGWD shape ',shape(pplay),cell+1
+
     DO jk = 1, nlev
-      zpm1r = pplay(cell+1, jk)/paprs(cell+1, 1)
+      !zpm1r = pplay(cell+1, jk)/paprs(cell+1, 1)
+      zpm1r = pplay(cell, jk)/paprs(cell, 1)
       IF (zpm1r>=zsigt) THEN
         nktopg_tmp = jk